Transcriptome analysis of embryonic domains in Norway spruce reveals potential regulators of suspensor cell death

The terminal differentiation and elimination of the embryo-suspensor is the earliest manifestation of programmed cell death (PCD) during plant ontogenesis. Molecular regulation of suspensor PCD remains poorly understood. Norway spruce (Picea abies) embryos provide a powerful model for studying embryo development because of their large size, sequenced genome, and the possibility to obtain a large number of embryos at a specific developmental stage through somatic embryogenesis. Here, we have carried out global gene expression analysis of the Norway spruce embryo-suspensor versus embryonal mass (a gymnosperm analogue of embryo proper) using RNA sequencing. We have identified that suspensors have enhanced expression of the NAC domain-containing transcription factors, XND1 and ANAC075, previously shown to be involved in the initiation of developmental PCD in Arabidiopsis. The analysis has also revealed enhanced expression of Norway spruce homologues of the known executioners of both developmental and stress-induced cell deaths, such as metacaspase 9 (MC9), cysteine endopeptidase-1 (CEP1) and ribonuclease 3 (RNS3). Interestingly, a spruce homologue of bax inhibitor-1 (PaBI-1, for Picea abies BI-1), an evolutionarily conserved cell death suppressor, was likewise up-regulated in the embryo-suspensor. Since Arabidopsis BI-1 so far has been implicated only in the endoplasmic reticulum (ER)-stress induced cell death, we investigated its role in embryogenesis and suspensor PCD using RNA interference (RNAi). We have found that PaBI-1-deficient lines formed a large number of abnormal embryos with suppressed suspensor elongation and disturbed polarity. Cytochemical staining of suspensor cells has revealed that PaBI-1 deficiency suppresses vacuolar cell death and induces necrotic type of cell death previously shown to compromise embryo development. This study demonstrates that a large number of cell-death components are conserved between angiosperms and gymnosperms and establishes a new role for BI-1 in the progression of vacuolar cell death.

Introduction Plant embryogenesis starts with the asymmetric division of the zygote in the plane perpendicular to the future apical-basal axis of the embryo. This division generates a small apical cell and a large basal cell, the progenitors of two structurally and functionally distinct domains: embryo proper (in angiosperms) or embryonal mass (EM, in gymnosperms) and suspensor, respectively [1]. The apical domain gives rise to the plant, whereas the suspensor functions as a conduit of growth factors and nutrients to the growing apical domain and is gradually eliminated through programmed cell death (PCD).
The terminal differentiation and elimination of the embryo-suspensor is the earliest manifestation of PCD in plant life. In Norway spruce (Picea abies L. Karst.), the suspensor contains several files of elongated cells, derived through a series of asymmetric cell divisions in the EM. Once produced, these cells undergo terminal differentiation and embark on the PCD pathway. Generation of new layers of suspensor cells thus results in a gradient of PCD stages along apical-basal axis of an embryo. While suspensor cells adjacent to the EM are at the commitment stage of PCD, the cells at the lower layers of the suspensor are characterized by increased degree of dismantling. Therefore, position of the cell within the suspensor of spruce embryos can be used as a marker of PCD stage [2,3,4].
Most examples of plant developmental PCD, including the death of the embryo-suspensor, belong to the class of vacuolar cell death [5]. During vacuolar cell death, the cell contents are removed completely by a combination of autophagy-like engulfment of the cytoplasm and organelles and vacuolar collapse. Necrosis is another major class of plant PCD characterized by mitochondrial dysfunction and early rupture of plasma membrane, resulting in incomplete removal of cell contents [5]. It has been shown that genetic suppression of vacuolar PCD in the terminally-differentiated cells can trigger necrosis [6].
Our understanding of the molecular machinery regulating developmental PCD in plants is advancing, yet remains limited compared to animal-specific apoptosis. During terminal differentiation, the plant cell achieves the competency for death through expression of transcription factors (TFs) that regulate expression of genes controlling PCD triggers and executioner [7,8]. Ethylene, reactive oxygen species (ROS), calcium influx and a decrease in pH have all been implicated as potential PCD triggers [7,9].
Autophagy and activity of hydrolytic enzymes, such as cysteine, serine and aspartic proteases and nucleases execute PCD and are directly responsible for cell dismantling and morphology of cell corpse. In Arabidopsis, a type II metacaspase AtMC9, XYLEM CYSTEINE PEPTIDASES 1 and 2 (XCP1 and XCP2) are involved in post mortem clearance of root xylem cell contents [10]. In Norway spruce, execution of vacuolar PCD in the embryo-suspensor requires activity of a type II metacaspase mcII-Pa and autophagy [4,6,11]. A S1-P1 type nuclease BIFUNCTIONAL NUCLEASE 1 (BFN1) participates in DNA degradation during lateral root cap cell death [12].
The large size of Norway spruce embryos, which contain several millimetre-long suspensors (as compared to the Arabidopsis suspensor being composed of a single file of 6-9 small cells), the use of somatic embryogenesis to provide an unlimited number of genetically identical embryos at a specific developmental stage and the sequenced genome make somatic embryos of Norway spruce a powerful model system for studying molecular mechanisms of developmental PCD. Here, we took advantage of this system to compare transcriptomes of the living (EM) and dying (embryo-suspensor) domains of plant embryos using high-throughput RNA sequencing (RNA-Seq). Our analysis revealed a subset of genes highly expressed in the suspensor and therefore representing potential PCD initiators and executioners. Among these genes, we have found a spruce homologue of BAX INHIBITOR 1 (BI-1), which has previously been shown to act as a suppressor of ER stress-mediated cell death in Arabidopsis [13]. Silencing of Norway spruce BI-1 (PaBI-1, for P. abies BI-1) by RNA interference (RNAi) induced a switch from vacuolar to necrotic cell death in the suspensor leading to abnormal embryo development. Our findings thus not only define that an anti-necrotic role of BI-1 is conserved between angiosperms and gymnosperms but also connect this role to the regulation of developmental PCD.

Norway spruce somatic embryo culture
Two embryogenic cell lines (88:22 and 11:18) originated from zygotic embryos of two independent fertilization events were cultured as previously described [14]. Briefly, the lines were maintained by weekly subculture on half-strength solidified LP medium containing growth regulators 9.0 μM 2,4-diclorophenoxyacetic acid (2,4-D) and 4.4 μM 6-benzylaminopurine (6-BA) (hereafter referred as proliferation medium). To stimulate embryo development, the cell lines were transferred onto solidified half-strength LP medium without growth regulators for one week followed by an additional week on solidified embryo maturation medium DKM supplemented with 30 μM abscisic acid (ABA). The cultures were grown in the dark at 22 o C.
Individual somatic embryos during the transition from early embryogeny to late embryogeny were split into EMs and suspensors under a Zeiss STEMI SV8 stereoscope (Germany). Briefly, the EM was hold in a petri plate with fine forceps and dissected by cutting through the first layer of the suspensor using a Feather surgical blade (Japan). The separated embryonic domains (S1 Fig) were treated with RNA stabilization reagent RNAlater (Qiagen) in a screwcap tube and then snap-frozen in liquid nitrogen before storage at -80 o .

RNA extraction, cDNA synthesis and RNA-seq
Total RNA was extracted from the EMs and suspensors using RNAqueous Micro kit (Ambion). RNA quality was assessed in terms of RNA integrity number (RIN) by Bioanalyzer (Agilent 2100 expert). To obtain sufficient amount of RNA required for RNA-seq, the RNA extracted from the embryos of cell line 11:18 was subsequently amplified using MessageAm-pIIaRNA Kit (Ambion).
cDNA library preparation and subsequent sequencing were performed at the SciLifeLab (Stockholm, Sweden). Strand-specific cDNA libraries were prepared with TruSeq Stranded mRNA Sample prep kit of 96 dual indexes (Illumina, CA, USA) according to the manufacturer's instructions except for the following changes. The protocols were automated in Agilent NGS workstation (Agilent, CA, USA) using purification steps as previously described [15,16]. Clonal clusters were generated using cBot (Illumina) and sequenced on HiSeq2500 (Illumina) according to manufacturer's instructions. Bcl to Fastq conversion was performed with bcl2Fastq v1.8.3 from the CASAVA software suite. The quality scale was Sanger / phred33 / Illumina 1.9. The obtained data were deposited to the European Nucleotide Archive (ENA) and is accessible under the accession number PRJEB22154.
Statistical analysis of single-gene differential expression between EM and suspensor was performed in R (v3.3.2; R Core Team 2015) using the Bioconductor (v3.4; [22]) DESeq2 package (v1.14.1; [23]). FDR adjusted p-values were used to assess significance; a common FDR threshold of 1% was used throughout. For the data quality assessment (QA) and visualization, the read counts were normalized using a variance stabilizing transformation (VST) as implemented in DESeq2. The biological relevance of the data e.g. biological replicates similarity was assessed by Principal Component Analysis (PCA) and other visualizations (e.g. heatmaps), using custom R scripts. An overview of the data, including raw and post-QC read counts and alignment rates is given in S1 Table.

Annotation of differentially expressed genes (DEGs) and gene ontology (GO) enrichment analysis
Fasta sequences for the DEGs were extracted from ConGenIE.org [19] and imported to Blas-t2GO for Gene Ontology (GO) enrichment analyses. For gene annotation, Basic Local Alignment Search Tool (http://www.ncbi.nlm.nih.gov/BLAST) of Blast2GO Version 4 [24] was used to find sequences similar to the DEGs. GO mapping option of Blast2GO was used to retrieve GO terms associated to the hits obtained by the BLAST search. For functional enrichment analysis, the numbers of annotated sequences in each GO term were counted by Blast2GO. Categorization of TFs was done by comparing the DEGs with the P. abies TFs listed in the Plant Transcription Factor Database (PlantTFDB) version 4.0 [25].

Quantitative real-time PCR (qRT-PCR)
cDNA was synthesized from 500 μg of RNA isolated from embryogenic cell line 11:18 using Maxima First Strand cDNA synthesis kit (Thermo Scientific). A twentieth part concentration of each cDNA sample was utilized for the analysis using Dynamo Flash SYBR Green kit (Thermo Scientific) in a CFX PCR thermal cycler (sequences of all primers used in this study are listed in S2 Table). ΔΔC T method was used to measure the fold expression of five genes of interest normalized to the expression of two reference genes: cell division control 2 (CDC2) and phosphoglucomutase which were selected on the basis of their stability tested by Biogazelle-qbase+ software.

Cloning of PaBI-1 hairpin construct
For creating the hairpin of Norway spruce PaBI-1 the corresponding parts of PaBI-1 cDNA were amplified using attB1_AS_BI1_F/AsBI1_R_HindIII and attB2_S_BI1_R/S_BI1_F_Hin-dIII using Phusion DNA polymerase (Thermo Scientific) and cut with HindIII (Thermo scientific), which produced 360 bp and 466 bp fragments of PaBI-1 coding DNA sequences, respectively. The fragments were ligated with T4 DNA Ligase (Thermo Scientific) and recombined into pDONR/Zeocin (Zeo) vector by BP Clonase (Invitrogen) followed by recombination into modified pMDC32 vector containing suspensor-specific NIP (nodulin-like intrinsic protein) promoter [26] by LR Clonase (Invitrogen). The resulting recombined vector was checked by digestion with restriction enzymes followed by sequencing.

Transformation
To suppress PaBI-1 expression with the hairpin construct, Agrobacterium-mediated transformation of embryogenic cell lines 11:18 was performed as described previously [27]. pMDC32:: PaBI-1 was transformed into Agrobacterium tumefaciens strain GV3101. Overnight culture of transformed Agrobacterium was centrifuged and resuspended in infiltration buffer (10 mM MES, 10 mM MgCl 2 , pH 5.5) to an OD 600 of 10. Resuspended bacteria were mixed with 5-dayold spruce proliferating culture in 1:9 ratio. Acetosyringone was added to the mixed culture to a final concentration of 150 μM and the culture was left for 5 h with regular shaking. Thereafter the culture was transferred onto a filter paper placed on solidified proliferation medium and left for 2 days at 22˚C in the dark after which, the filter paper was transferred to a fresh proliferation medium supplemented with 400 μg/ml Timentin and 250 μg/ml cefotaxime and cultured for additional 5 days under the same conditions. Next, the cells on filter paper were transferred to proliferation medium with 5 μg/ml hygromycin, 400 μg/ml Timentin and 250 μg/ml cefotaxime and incubated for 7 days under the same condition. At the next and each following subculture with 7-day intervals, the concentration of hygromycin was increased to 7.5 μg/ml, while other components remained the same. Once transformants emerged, they were transferred to a fresh proliferation medium with increasing concentration of hygromycin (10 μg/ml and 15 μg/ml) every week and finally maintained on the proliferation medium with 15 μg/ml hygromycin.

Embryo staining and confocal microscopy
To detect the presence of dead cells, 10 ml of 5-day-old spruce culture resuspended in liquid proliferation medium devoid of growth regulators was stained with 0.0025% (w/v) Evans blue solution for 10 min at room temperature as described previously [6]. An aliquot of stained embryos was transferred onto a glass slide and imaged with a bright field microscopes (Zeiss Axioplan) with 2.5X Plan-Neofluar objective, NA 0.075. Images were taken using a DFC295 camera and LAS AF v3.2 software. The length of the suspensor was measured with ImageJ software.
For detection of necrotic cells, 5-day-old embryogenic suspension culture grown in liquid medium without growth regulators was stained first with 2 μg/ml FDA for 15 min and then counterstained with 1 μM FM4-64 for maximum 10 min as described previously [6]. The cells were imaged within 10 min after staining with FM4-64 using the sequential mode of a Zeiss LSM 780 confocal microscope, with a 20X objective, NA 0.8, excitation at 488 nm and 561 nm and emission at 490-587 nm and 582-754 nm for detection of FDA and FM4-64 staining, respectively.

Overview of Norway spruce somatic embryo transcriptome
To identify genes differentially expressed in the EM versus suspensor, RNA was isolated from the corresponding domains of somatic embryos originated from two unrelated genotypes (viz. cell lines). Thus, four RNA-Seq libraries were sequenced using Illumina technology which generated a total of~410 million raw reads, the expression proxy of 47,552 genes (Fig 1A and 1B). Expression of 32.8% percent (23,184) of the 70,736 Norway spruce genes was not detected in the embryos. Accounting for the difference in library size (i.e. sequencing depth) and the technical and biological variability of the samples, a total of 451 genes were found to be differentially expressed between the two embryonic domains at a 1% adjusted p-value cutoff (False Discovery Rate). Of these 451 genes, 53 and 398 were up-regulated in the EM and suspensor, respectively (S1 and S2 Datasets; Fig 1A and 1B).

Genes up-regulated in the EM
About 85% (45 out of 53) of the genes up-regulated in the EM could be mapped to GO ( Fig  1A; S1 Dataset). Using Blast2GO PRO analysis we have identified that a significant fraction of these genes (13 out of 45) encode enzymes possessing oxidoreductase, transferase, hydrolase or isomerase activity ( Fig 1C, Table 1). Of these enzymes, seven fall into flavonoid biosynthesis pathway [28] (Fig 2), suggesting enhanced ROS scavenging activity maintained in the meristematic cells of the EM [29]. Using the Plant Transcription Factor Database [25] we identified six MYB and MYB-like TFs up-regulated in the EM ( Table 2). The function of most of these TFs in Arabidopsis was linked to the regulation of flavonoid biosynthesis [37,38]. In particular, MYB12 was shown to activate promoters of several flavonoid biosynthesis enzymes, such as chalcone synthase (TT4), flavanone 3-hydroxylase (TT6), flavonol synthase (FLS) and chalcone flavanone isomerase (CHI) [39]. In agreement, TT4, TT6 and CHI were transcriptionally upregulted in the EM (Fig 2; S2 Fig). Polar auxin transport is essential for the correct apical-basal patterning of conifer embryos [41]. Flavonoids modulate the activity of auxin-transporting P-glycoproteins and may also be involved in modulating the activity of regulatory proteins e.g. phosphatases and kinases [42]. Thus, up-regulation of flavonoid biosynthesis genes and related TFs may promote spruce embryo development by regulating polar auxin transport and functions of housekeeping enzymes, in addition to more general role of flavonoids in ROS scavenging and abiotic stress resistance.

Genes up-regulated in the embryo-suspensor
Altogether, about 80% (318 out of 398) of the genes up-regulated in the suspensor were assigned to a GO category (Fig 1A; S2 Dataset), functional enrichment analysis of which revealed that a significant fraction of these genes (74 out of 318) encode enzymes possessing oxidoreductase, transferase, hydrolase or lyase activity, with hydrolases being clearly predominating ( Fig 1C; Table 1).
Genes related to cell elongation and cell wall modification. The suspensor cells in Norway spruce achieve highly-elongated morphology accompanied by growth of lytic vacuoles prior to tonoplast collapse and complete clearance of cellular contents [1]. This developmental process requires abundance of structural materials (e.g. carbohydrates and lipids), as well as activity of cell wall modifying enzymes. We have found transcriptional up-regulation of aquaporins, which accelerate water uptake to facilitate cell expansion [43] and choline kinase (CK) (as detected by qRT-PCR, S2 Fig), which is involved in the biosynthesis of phosphatidylcholine [44], the major component of plasma membrane and tonoplast [45]. Enzymes participating in cell wall loosening and reorganization, such as xyloglucan endotransglucosylase/hydrolase, galactosidases and pectinesterase were likewise up-regulated in the embryo-suspensor (Table 3).
Transcription factors. In animals, besides a major role of post-translational regulation, developmental apoptosis is also controlled at the transcriptional level [50]. Our analysis of the transcriptome of the Norway spruce embryo-suspensor suggests that transcriptional control   (Table 4).

Cell-death triggers and stress-responsive genes.
Reactive oxygen species (ROS) act as signalling molecules to control plant developmental PCD [7]. We observed enhanced expression of three homologs of L-ascorbate oxidase in the Norway spruce embryo-suspensor (Table 5). It has been shown that enhanced expression of L-ascorbate oxidase causes accumulation of H 2 O 2 and changes redox homeostasis [59]. A gene encoding another H 2 O 2 producing enzyme, germin, which is known to act at early stages of stress-induced cell death [60], was likewise up-regulated in the suspensor (Table 5). Germin has previously been shown to be transcriptionally up-regulated in the suspensor of Larix marschlinsii somatic embryos, where it has been suggested to also participate in cell wall remodelling [61].
Apart from the potential cell-death triggers, we have also found increased expression of a number of genes associated with stress response (Table 5). These include a cytochrome p450 involved in plant response to osmotic stress [62], hypersensitive response (HR) to bacteria and senescence [63]; alcohol oxidase involved in anoxia tolerance through alcoholic fermentation [64] and two genes encoding heat shock proteins (HSPs). Interestingly, we also found enhanced expression of a homologue of evolutionary conserved ER-stress induced cell-death suppressor BAX INHIBITOR 1 (BI-1), PaBI-1 (Table 5; S2 Fig) [13] and Bcl2-associated athanogene 1 (BAG1) implicated in proteasomal degradation of misfolded protein [65].
Proteases and nucleases. Among hydrolytic enzymes transcriptionally up-regulated in the embryo-suspensor, several cysteine peptidases have previously been reported to be involved in PCD and development of different plant tissues and organs [10,11,67] (Table 5). For example, the Norway spruce gene MA_103463g0010 is a homologue of Arabidopsis papain-like cysteine protease CEP1, which mediates tapetal PCD and pollen development [67], whereas gene MA_616703g0010 is a homologue of Arabidopsis METACASPASE 9 (AtMc9) participating in post-mortem xylem vessel clearance [10]. The Arabidopsis homologues of two up-regulated genes for cathepsin B-like cysteine protease were previously implicated in senescence, HR, UV, oxidative stress-and ER stress-mediated PCD [68][69][70]. Noteworthy, RNA-Seq analysis did not reveal transcriptional up-regulation of another type-II metacaspase, mcII-Pa, which is known to play a major role in PCD of Norway spruce embryo-suspensor [4,11]. Changes of mcII-Pa gene are too subtle for us to observe them, or the changes in that gene are not common to all individual cells within a pool and hence the effect is diluted or neutralised-i.e. unobservable at that resolution. Post-translational modification of mcII-Pa, i.e., autoprocessing and calcium binding, required to convert mcII-Pa zymogen to active enzyme [11] could be another explanation of this finding.

PaBI-1 is required for embryo patterning and suppression of necrotic cell death
BI-1 is an evolutionary conserved cell-death suppressor localized to the ER membrane and regulating calcium and lipid dynamics under ER stress [13,76]. Although BI-1 has been implicated to suppress chemically induced ER stress-mediated cell death, as well as cell death induced by necrotrophic fungi and heat stress [13,73], the function of BI-1 in plant developmental PCD remains unknown.
To investigate the role of PaBI-1 in embryo development and associated PCD, we suppressed PaBI-1 expression using RNAi (Fig 3A and 3B). The gene silencing impaired apical-basal Table 5. Examples of potential anti-and pro-cell-death genes up-regulated in the Norway spruce embryo-suspensor and functions of their Arabidopsis homologues. patterning ( Fig 3A) through suppression of anisotropic expansion of the suspensor cells, as revealed by measuring the length of Evan's blue positive cells (Fig 3C). This accounted for the increased frequency of aberrant early embryos lacking properly formed suspensors in the PaBI-1 RNAi lines (Fig 3D) and the decreased number of cotyledonary embryos (Fig 3E and 3F). Based on the morphological criteria, most examples of plant cell death can be divided in two main classes: vacuolar cell death and necrosis [5]. Silencing of PaBI-1 switched the mode of cell death in the suspensor from vacuolar to necrotic, as indicated by the presence of shrunken and largely undigested protoplast (Fig 4A). Quantification of this necrotic hallmark using double staining with FDA and FM4-64 has shown ten-fold increase in the frequency of necrotic cells in PaBI-1 RNAi lines, as compared with the control line ( Fig 4B).

ConGenIE ID
In animals, BI-1 functions as a suppressor of apoptosis induced by cell-death effector Bax. Plants lack direct homologues of Bax and other members of Bcl-2 family, but surprisingly, a homologue of BI-1 is present in plants [77]. The mechanism by which PaBI-1 suppresses necrosis and sustains progression of vacuolar cell death in spruce embryo-suspensor remains to be identified. Yet, one can propose that this mechanism is reliant on either maintenance of ER homeostasis or interaction with autophagy pathway or a combination of both for the following two reasons.
Firstly, ER stress is induced by accumulation of misfolded or unfolded proteins. In Arabidopsis, BI-1 has been shown to be transcriptionally up-regulated upon ER stress-induced unfolded protein response (UPR). AtBI-1 keeps the cell alive until the ER homeostasis is reestablished by the activity of ER chaperons such as Bip2 [13]. The enhanced expression of PaBI-1 in the suspensor suggests that this developmental PCD may likewise implicate ER stress. Vacuolar cell death is believed to be a slow process [5] and the cell remains metabolically active until vacuolar collapse. We speculate that down-regulation of PaBI-1 compromises ER homeostasis in the terminally-differentiated suspensor cells, thus triggering necrosis.
Secondly, it has been recently reported that in Nicotiana benthamiana, tobacco BI-1 interacts with autophagy-related protein ATG6 and silencing of BI-1 reduced autophagic flux and enhanced N-mediated hypersensitive response (HR) cell death [78]. In spruce, autophagy sustains vacuolar cell death and prevents necrosis of the suspensor cells [6]. Therefore, PaBI-1 deficiency-induced necrotic death of the suspensor cells might be due to compromised autophagy, the possibility that calls for further investigations.

Conclusions
In this study, we performed global transcriptome analysis of the two embryonic domains of Norway spruce to identify potential regulators of suspensor PCD. Since the embryo-suspensor in Norway spruce is composed of several cell layers representing a gradient of successive stages of PCD [1][2][3], we attempted to ascribe identified suspensor-specific DEGs to the sequence of events underlying terminal cell differentiation, initiation and execution of PCD (Fig 5). Our observations suggest that suspensor cells adjacent to the EM and undergoing terminal Thereafter cells initiate PCD by expressing transcription factors and stress-responsive genes. The latter might act to suppress rapid cell death and to allow gradual cell dismantling characteristic for vacuolar PCD. At the basal end of the suspensor, hydrolytic enzymes (proteases and nucleases) execute PCD by processing protein and nucleic acid substrates. Note that some gene names correspond to Arabidopsis homologues. https://doi.org/10.1371/journal.pone.0192945.g005 Transcriptome analysis of spruce embryo suspensor cell death differentiation express genes controlling loosening and reorganization of the cell wall, as well as aquaporin and CK, thus facilitating cell expansion and vacuole enlargement. Next, the elongated suspensor cells initiate PCD by expressing a subset of dedicated TFs, which in turn confer simultaneous expression of PCD triggers (e.g. L-ascorbate oxidase and germins) and stressresponse genes preventing rapid cell collapse (e.g. PaBI-1) and sustaining slow progression of vacuolar cell death. It is conceivable that all or some of these TFs, in particular XND1 and ANAC075, are already highly expressed in the first suspensor cell layer composed of cells undergoing terminal cell differentiation to orchestrate expression of genes regulating cell wall loosening and reorganization. Finally, execution of PCD at the basal end of the suspensor requires enhanced expression of hydrolytic enzymes acting to degrade proteins and nucleic acids (Fig 5). Clearly, a higher resolution, single cell layer-specific mRNA isolation procedure assisted by laser capture microdissection is required to verify the proposed sequence of transcriptional changes occurring throughout Norway spruce embryo-suspensor.
Olvera-Carrillo and colleagues [74] have recently performed a comparative bioinformatics analysis of transcriptomes of developmentally-regulated cell deaths in Arabidopsis and concluded that various examples of developmental PCD share a common set of regulators, including RNS3, BFN1, PASPA3, AtMc9 and SCPL48. Availability of conifer genomes enables evolutionary insight into molecular regulation of plant development in gymnosperm versus angiosperm lineages, and our transcriptomics data suggest significant conservation of developmental cell death components between the lineages. This conservation encompasses both TFs (e.g. ANAC072, ANAC075 and XND1) and downstream effectors, such as RNS3 and metacaspase-9. Identification of anti-cell death protein BI-1 as a potential regulator of Norway spruce suspensor PCD indicates that there might yet be some interesting differences between angiosperms and gymnosperms. Further functional studies are required to link PaBI-1 and other candidate components of suspensor PCD into biochemical pathways.