WUSCHEL-RELATED HOMEOBOX 2 is important for protoderm and suspensor development in the gymnosperm Norway spruce
© Zhu et al. 2016
Received: 19 August 2015
Accepted: 6 January 2016
Published: 19 January 2016
Distinct expression domains of WUSCHEL-RELATED HOMEOBOX (WOX) gene family members are involved in patterning and morphogenesis of the early embryo in Arabidopsis. However, the role of WOX genes in other taxa, including gymnosperms, remains elusive. Here, we use somatic embryos and reverse genetics for studying expression and function of PaWOX2, the corresponding homolog of AtWOX2 in the gymnosperm Picea abies (Pa; Norway spruce).
The mRNA level of PaWOX2 was transiently up-regulated during early and late embryogeny. PaWOX2 mRNA in early and early late embryos was detected both in the embryonal mass and in the upper part of the suspensor. Down-regulation of PaWOX2 during development of early embryos resulted in aberrant early embryos, which failed to form a proper protoderm. Cells on the surface layer of the embryonal mass became vacuolated, and new embryogenic tissue differentiated from the embryonal mass. In addition, the aberrant early embryos lacked a distinct border between the embryonal mass, and the suspensor and the length of the suspensor cells was reduced. Down-regulation of PaWOX2 in the beginning of embryo development, before late embryos were formed, caused a significant decrease in the yield of mature embryos. On the contrary, down-regulation of PaWOX2 after late embryos were formed had no effect on further embryo development and maturation.
Our data suggest an evolutionarily conserved function of WOX2 in protoderm formation early during embryo development among seed plants. In addition, PaWOX2 might exert a unique function in suspensor expansion in gymnosperms.
The basic plant body pattern is set up during embryogenesis. In seed plants, this body plan has been described as the superimposition of two patterns: an apical-basal and a radial pattern . The molecular processes that establish this primary body plan have mainly been studied in the angiosperm model species Arabidopsis (Arabidopsis thaliana). In contrast, knowledge about the molecular regulation of embryo development in gymnosperms is limited, partly owing to the lack of identified zygotic embryo defective mutants and genetic tractability. However, by using somatic embryos and reverse genetics it has been possible to study the regulation of embryo development in conifers . We are studying the early stages of embryo development in Norway spruce and especially the role of members of the WUSCHEL-RELATED HOMEOBOX (WOX) gene family, which encode transcription factors that play important roles in the determination of cell fate during embryogenesis in angiosperms [3, 4].
Angiosperms and gymnosperms separated approximately 300 million years ago and expectedly their patterning during embryogenesis differs significantly. For convenience, embryogenesis in Arabidopsis can be divided into three general phases, described as proembryogeny, early embryogeny (globular-stage to heart-stage transition) and late embryogeny . Proembryogeny begins after fertilization. An asymmetric cell division of the zygote generates a smaller apical cell and a larger basal cell. The apical cell is the founder of the embryo proper, while the basal cell develops into the suspensor. The apical cell undergoes several rounds of stereotyped asymmetric divisions, giving rise to cells with positionally determined cell fate. At the beginning of early embryogeny, in the 8-cell embryo proper, a single round of tangential divisions separate the outer layer of eight cells from the eight inner cells . The inner cells are founder cells of the ground tissue and vascular elements. The outer cells form the protoderm which will become the epidermis. Plant epidermis is characterized by the secretion of lipids and waxes to its outer cell wall . The continuous hydrophobic layer forms the cuticle. Owing to the periclinal cell division pattern in the protoderm, the protodermal cells remain essentially separated from the inner cells throughout embryogenesis. Root and shoot meristems are established during the transition from globular to heart-stage. After the heart-stage the suspensor is degraded by programmed cell death . During late embryogeny, there is a switch from pattern formation to storage product accumulation.
In conifers, the sequence of embryo development can also be described by three phases: proembryogeny, early embryogeny, and late embryogeny . Proembryogeny begins when the zygote undergoes several rounds of nuclear duplications. The produced free nuclei are first arranged in a tier before cellularization. After a cell division, two tiers are formed. The cells in the upper tier elongate to form a functional suspensor and the cells in the lower tier divide creating the embryonal mass (analogous to the embryo proper in angiosperms). Early embryogeny begins with the elongation of the embryonal suspensor. Cells in the outer layer of the embryonal mass divide mainly anticlinally, but also periclinally giving rise to additional internal layers , unlike angiosperms where only anticlinal cell divisions take place. Nevertheless, the outer cell layer in the embryonal mass in conifer embryos defines a functional protoderm [9, 10]. Late embryogeny is a period of histogenesis and organogenesis. Early during this phase, the suspensor cells are dismantled by programmed cell death , and the root and shoot apical meristems are delineated.
Members of the WOX gene family play important roles in determining cell fate during plant development. Phylogenetic analyses have identified three major clades in the WOX gene family: the modern (WUS and WOX1-7), the intermediate (WOX8, 9, 11, and 12) and the ancient clade (WOX10, 13, and 14) . In the gymnosperm Norway spruce (Picea abies), 11 WOX genes that belong to all major clades, have been identified .
All WOX genes examined show very specific expression patterns, both spatially and temporally, which are important for their molecular functions . WOX2, WOX8 and WOX9 have been implicated in patterning and morphogenesis of the early embryo in Arabidopsis [3, 15, 16]. AtWOX2 is expressed in the egg cell and the zygote . After the first cell division, AtWOX2 marks the apical cell, while, AtWOX8 and AtWOX9 mark the basal cell. AtWOX8 and AtWOX9 share redundant functions . The Arabidopsis wox8wox9 double mutant shows aberrant orientation of cell division planes in the embryo . Previously, we showed that the Norway spruce gene PaWOX8/9 is the orthologue of AtWOX8 and AtWOX9 . Similarly to Arabidopsis, down-regulation of PaWOX8/9 causes disturbed orientation of the cell division plane and cell fate determination during early embryo pattern formation suggesting that PaWOX8/9 exerts an evolutionarily conserved function.
It has previously been shown that PaWOX2, which is highly similar in sequence to AtWOX2, is expressed during embryo development in gymnosperms [13, 18–20]. Herein we report that PaWOX2 is transiently expressed in the embryonal mass and in the upper part of the suspensor in early and early late embryos. Furthermore, by using reverse genetics, we functionally characterize PaWOX2. Down-regulation of PaWOX2 during development of early embryos results in aberrant early embryos, which lack distinct embryonal mass and suspensor domains, a proper protoderm, and have shorter suspensor cells. This suggests an evolutionarily conserved function of WOX2 in protoderm formation early during embryo development among seed plants. In addition, PaWOX2 might exert a unique function in proper suspensor cell expansion in gymnosperms.
Materials and methods
The embryogenic cell line 61:21 of Norway spruce (Picea abies L. Karst) has been used in this study. The cell line was established as described by Högberg et al. . The cell line was stored in liquid nitrogen and thawed a couple of months before the start of experiments. After thawing the cell cultures were maintained as described previously . Briefly, proembryogenic masses (PEMs) were maintained on solidified proliferation medium containing the plant growth regulators (PGRs) auxin and cytokinin. To stimulate development of somatic embryos the cultures were first transferred to pre-maturation medium lacking PGRs for one week and then to maturation medium containing 30 μM abscisic acid (ABA). Early embryos (EEs) differentiated after one week on maturation medium; early late embryos (LE1s) and LE2s developed after two and three weeks on maturation medium, respectively; maturing embryos (ME1s), characterized by the initiation of cotyledons, developed after five weeks on maturation medium. ME2s (almost fully matured embryos) and ME3s (fully matured embryos) developed after about eight weeks on maturation medium.
RNA extraction, cDNA synthesis and quantitative real-time PCR
Samples for analyzing the mRNA level of PaWOX2 (accession number: AM286747) during embryo development, were collected from nine sequential developmental stages: PEM1 (after seven days on proliferation medium), PEM2 and PEM3 (after three and seven days on pre-maturation medium respectively), EE, LE1, LE2, ME1, ME2 and ME3 (after one to eight weeks on maturation medium). The sampling was performed at mid-day and samples were frozen in liquid nitrogen and stored at −80 °C after collection.
Total RNA was isolated using the Spectrum Plant Total RNA kit (Sigma-Aldrich, USA) according to the manufacturer’s instructions. For each sample, 1 μg of total RNA was reverse transcribed with RevertAid H Minus First Strand cDNA Synthesis Kit (Fermentas, Thermo Scientific, Sweden) using an equimolar ratio of random and oligo-dT primers according to the manufacturer’s instructions.
Quantitative real-time PCR (qRT-PCR) was performed as described previously . Three reference genes, CELL DIVISION CONTROL2 (PaCDC2), ELONGATION FACTOR 1 (PaEF1) and PHOSPHOGLUCOMUTASE (PaPHOS) were used . Two to three biological replicates, each with three technical replicates were performed for each test. The primer sequences are presented in Additional file 1: Table S1. Statistical analysis was done by t-test.
RNA in situ hybridization
For RNA in situ hybridization (ISH) the following materials were used: ovules from cones collected in the end of June and somatic embryos (EEs, LE1s and ME1s). The ovules were fixed and embedded as described by Karlgren, et al. . The somatic embryos were fixed in 3.7 % formaldehyde, 5.0 % acetic acid and 50 % ethanol overnight and embedded in Technovit 8100. A gene-specific fragment was used as a probe. The probes were prepared with DIG RNA Labeling Kit (see primer sequences in Additional file 1: Table S1) (Sigma-Aldrich, USA). In situ hybridization was performed essentially as described by Karlgren et al. . Sections of 10 μm were hybridized to digoxigenin-labeled RNA probes. The pictures were processed using Adobe Photoshop CS6 13.0 software.
RNA interference vector construction
The coding sequence (CDS) of PaWOX2 was amplified from a cDNA library of early somatic embryos of Norway spruce . The full-length CDS was subcloned into the pJET1.2/blunt cloning vector using the CloneJET™ PCR Cloning Kit (Fermentas, Thermo Scientific, Sweden). To obtain RNA interference (RNAi) constructs, two overlapping fragments of PaWOX2 were amplified and fused to form a hairpin structure for PaWOX2 (Additional file 2: Figure S1). To fuse these fragments, EcoRI and BamHI digestion sites were added on forward primers as linkers. The hairpin was confirmed by sequencing. Primers are presented in Additional file 3: Table S2.
Hairpin structures were introduced into pENTR™/D-TOPO® (Invitrogen, Carlsbad, CA, USA) and then inserted by att site LR recombination into the destination vector pMDC7 [LexA-VP16-ER (XVE) β-estradiol inducible promoter, which is derived from the pER8 vector and contains the estrogen receptor-based transactivator XVE] [24, 25] or pMDC32 (35S constitutive promoter) . Hairpin structures were confirmed by sequencing. Vectors were introduced by electroporation into Agrobacterium tumefaciens strain GV3101.
Transgenic cell lines
Embryogenic cultures were transformed by co-cultivation with A. tumefaciens as described previously . Stably transformed lines were selected after four weeks. Genomic DNA was isolated from PEMs from selected lines by using the DNeasy plant mini kit (Qiagen, Germany), according to the manufacturer’s instructions. Transformed lines were confirmed by PCR.
The mRNA level of PaWOX2 in PaWOX2 RNAi lines was analyzed by qRT-PCR. In the case of the inducible XVE-WOX2i lines, cultures were induced with β-estradiol (10 μM) for 48 h before the analysis. The lines 35S:WOX2i.2, 35S:WOX2i.3, 35S:WOX2i.4 and XVE-WOX2i.12 were selected for further studies (Additional file 4: Figure S2). The untransformed 61:21 line was used as a control. In addition, in the time-laps tracking experiments we also included a transformed control line (T-control), expressing the reporter GUS (β-glucuronidase) under the 35S promoter.
To study if PaWOX2 regulates cell division, the mRNA level of nine cell-cycle-regulating genes [PaRETINOBLASTOMA-RELATED PROTEIN-LIKE (PaRBRL), PaEXTRA SPINDLE POLES (PaESP), two E2F family genes (PaE2FABL) and five CYCLIN-LIKE (PaCYCLs) genes] were analysed in EEs from the control and line 35S:WOX2i.4 by qRT-PCR as previously described .
Samples of EEs for morphological analysis were collected from the control and lines 35S:WOX2i.2, 35S:WOX2i.3, 35S:WOX2i.4 and XVE-WOX2i.12 after one week on maturation medium. The samples were embedded by mixing with 2 ml of 1.2 % (w/v) Seaplaque agarose (FMC BioProducts, USA) in 60 mm Petri dishes. The length and width of suspensor cells of about 40 EEs from the control and 35S:WOX2i lines were measured using the ImageJ software (ver. 1.48 g) .
To further study embryo morphology, 26 LE1s from the control and line 35S:WOX2i.4 were scanned with a Zeiss 780 confocal microscope (Carl Zeiss AG), using the 488 nm Argon laser line, and the 20x objective (NA = 0.80).
For histological analysis, EEs from the control and line 35S:WOX2i.4 were fixed in 3.7 % formaldehyde, 5.0 % acetic acid and 50 % ethanol overnight. Subsequently, samples were dehydrated in 50, 75, 90 and 100 % ethanol series. Finally, the samples were embedded in Technovit 8100 (Kulzer, Wehrheim, Germany). The embryos were processed for serial sectioning (10 μM) on a Zeiss HM 355 microtome.
The cuticle of untreated LE1s from the control and line 35S:WOX2i.4 was stained in freshly prepared Oil Red, 0.2 % (w/v) in water, for 5 min and then washed in water . The stained embryos were hand-sectioned and examined under a Zeiss Axioplan microscope in dark field with a 5x objective (NA = 0.12).
Time-lapse tracking analysis was performed to examine in great detail the developmental pattern from EE to ME. EEs from the controls (both untransformed and T-control) and lines 35S:WOX2i.2, 35S:WOX2i.3 and 35S:WOX2i.4 (50 embryos per line) were sampled after one week on maturation medium and transferred to fresh maturation medium. Embryo morphology was examined every second day for 15 days.
To study the effect of PaWOX2 on the maturation process, after one week on pre-maturation medium cultures from the control and lines 35S:WOX2i.2, 35S:WOX2i.3, 35S:WOX2i.4 and XVE-WOX2i.12 were re-suspended in liquid pre-maturation medium and plated out as a thin layer on filter paper placed on maturation medium. For the XVE-WOX2i line, β-estradiol (10 μM) was added to the maturation medium, either from the start or after two weeks on maturation medium, when LE2s had already developed. The development of embryos was recorded for 14 days on maturation medium. The time points examined were 1, 3, 6, 10 and 14 days. The number of ME3s developed per initial gram of tissue was estimated after seven weeks on maturation medium.
Expression of PaWOX2 during embryo development
PaWOX2 is required for proper protoderm formation
In order to examine the function of PaWOX2, we constructed stably transformed PaWOX2 RNAi lines, using constitutive (35S; referred as 35S:WOX2i) or inducible (XVE; referred as XVE-WOX2i) promoters to drive expression of a hairpin RNA used to promote PaWOX2 specific RNAi. The down-regulation of PaWOX2 was confirmed by qRT-PCR (Additional file 4: Figure S2). The transcript level of PaWOX2 in non-induced tissue was lower than in the control. In accordance, it was previously shown that the XVE promoter is partially active in Norway spruce even in the absence of β-estradiol .
PaWOX2 is not regulating cell division at the transcriptional level
Previously, we showed that PaWOX8/9, which specifies the cell division plane orientation in the basal part of the embryonal mass, regulates cell division at the transcriptional level . In order to elucidate whether PaWOX2 exerts a similar function, we compared the mRNA abundance of nine cell-cycle-regulating genes in EEs from the control and from line 35S:WOX2i.4 by qRT-PCR. No significant differences in mRNA level of the tested cell-cycle-regulating genes could be observed between control and the RNAi line (Additional file 8: Figure S4). The previous suggests that PaWOX2 and PaWOX8/9 regulate the expression of different targets during embryo development.
PaWOX2 is important for the development of mature embryos
For elucidating when during embryo development PaWOX2 is important, we analyzed development of embryos of the inducible XVE-WOX2i lines. Initially 12 XVE-WOX2i lines treated with β-estradiol were screened for their potential to develop MEs. In 8 out of these 12 lines, the cultures continued to proliferate on maturation medium and the yield of MEs was low. Cell line XVE-WOX2i.12 was choosen for further studies. The frequency of EEs with normal morphology was significantly lower in line XVE-WOX2i.12 induced with β-estradiol than in the control (Fig. 3b). For examining the role of PaWOX2 in the formation of MEs, line XVE-WOX2i.12 was induced with β-estradiol either from the first or from the third week on maturation medium. When down-regulation of PaWOX2 was induced from the first week on maturation medium, the yield of MEs was significantly decreased (Fig. 7c, Additional file 11: Table S7). However, the yield of MEs was not affected when down-regulation of PaWOX2 was induced at the third week on maturation medium, when LE2s had formed. These results suggest that the transient high expression of PaWOX2 during the formation of EEs and LE1s is important for further development of the embryos.
In this study we show that the mRNA level of PaWOX2 in somatic embryos is transiently up-regulated during early and late embryogeny. Results obtained by qRT-PCR analyses, were in good agreement with results obtained by in situ mRNA hybridization. The expression of PaWOX2 in EEs and LE1s was detected in the embryonal mass and in suspensor cells that are in proximity to the embryonal mass. In situ mRNA analyses were also conducted on zygotic embryos. Owing to the difficulty to obtain high quality sections of the zygotic embryos we could not localize the signal to specific cell-types. However, we could detect signal in late embryos but not in maturing embryos. A comparable expression pattern for PaWOX2 in somatic and zygotic embryos of Norway spruce has been reported by Palovaara, et al. . Taken together, these results indicate that PaWOX2 is expressed in a similar way in somatic and zygotic embryos during embryo development.
Hybridization signals could be detected in the region of the mega gametophyte that resides in front of the growing embryo. In Scots pine (Pinus sylvestris) it has been shown that a cone-shaped zone is formed in the mega gametophyte in front of the dominant embryo . The cells in this zone are degenerating to make room for the actively growing embryo. We assume that the hybridization signals detected in our analyses correspond to the region in the mega gametophyte where the cells are committed to undergo cell death. However, at present we do not know why these cells are PaWOX2-positive.
To investigate the function of PaWOX2, embryogenic RNAi lines were established to down-regulate mRNA level of PaWOX2. EEs and LEs in35S:WOX2i lines showed aberrant morphology characterized by a severe disturbance of development in the border between the embryonal mass and the suspensor, decrease in the length of the suspensor cells and lack of the protoderm. Vacuolated cells formed in the outer cell layer of the embryonal mass and new embryogenic tissue was initiated from the embryonal mass. The previous result suggests that PaWOX2 has a role in the regulation of early embryo development, both in the embryonal mass and the suspensor.
Normal EEs and LEs have a distinct border between the compact globular embryonal mass in the apical part and elongated suspensor cells in the basal part. In contrast, the border between the embryonal mass and the suspensor in embryos from 35S:WOX2i lines is severely disturbed. In Arabidopsis, AtWOX2 shares redundant functions with AtWOX8 in regulating embryonic apical patterning [15, 31]. Previously, we showed that PaWOX8/9 is crucial for embryo patterning by regulating the orientation of the cell division plane at the basal part of the embryonal mass during early and late embryo development in Norway spruce. Expression level of cell-cycle-regulating genes, such as E2F is reduced in PaWOX8/9 RNAi lines . However, we could not detect any changes in mRNA levels of the tested cell-cycle-regulating genes in EEs from 35S:WOX2i lines. The previous suggests that PaWOX2 may not control cell division at the transcriptional level as PaWOX8/9 does, and that PaWOX2 and PaWOX8/9 might be involved in early embryo patterning through distinct pathways.
In Norway spruce, the most basally situated cells in the embryonal mass undergo asymmetric divisions to form the suspensor cells. Thereby, several files of suspensor cells are continuously formed. The suspensor cells do not divide but they elongate and become highly vacuolated . In contrast, many suspensor cells in embryos from 35:WOX2i lines fail to elongate. The defect in the suspensor cells is consistent with the expression pattern of PaWOX2. PaWOX2 is expressed in the embryonal mass and in suspensor cells that are in proximity to the embryonal mass. In Arabidopsis, expression of AtWOX2 is limited to the upper tier  and no aberrant phenotype has been reported for the basal lineage of the embryos from the wox2 mutant. Taken together, PaWOX2 is functioning both in the embryonal mass and the suspensor while AtWOX2 is only functioning in the apical part of the embryo proper.
In addition to the apical-basal polarity along the shoot-root axis, the basic body plan of plant embryos also shows a radial organization of primary tissue players: epidermis (derived from the protoderm), cortex and endodermis (derived from the ground tissue), pericycle and vascular tissues (derived from the procambium). Protoderm differentiation, which is the earliest event of radial pattern formation in plant embryogenesis [5, 33], is essential for normal patterning during early embryo development in both gymnosperms and angiosperms [9, 10, 34, 35]. In Norway spruce, normal EEs and LE1s have a compact globular embryonal mass covered by the protoderm. In contrast, the embryonal mass in EEs and LE1s from 35S:WOX2i lines with reduced PaWOX2 mRNA level, lack a distinct protoderm. In a similar way, embryos from the Arabidopsis wox2 mutant fail to correctly separate the protodermal layer by periclinal divisions in the upper part of the embryo . We have previously shown that differentiation of the outer cell layer in the embryonal mass is also regulated by PICEA ABIES HOMEOBOX1 (PaHB1), which encodes a protein highly similar to those from the HD-GL2 angiosperm counterparts . The GL2-type homeobox genes, e.g. ARABIDOPSIS THALIANA MERISTEM LAYER 1 (ATML1) from Arabidopsis and OVULE39 (O39) from Phalaenopsis, are considered to be a specialized subgroup of plant homeobox genes that exert a specific role in protoderm development . This implies that the determination of the protoderm is regulated in a similar way in seed plants.
The epidermal cells secrete lipids and waxes to their outer cell wall. The previous results in the formation of a cuticle layer early during embryo development, after differentiation of the protoderm [7, 37]. The protoderm in normal LE1s is stained with Oil Red suggesting that the protoderm is covered by a cuticularized layer . In contrast, LE1s from 35S:WOX2i lines either lack a cuticularized layer or the layer is thin and patchy, indicating that the function of some cells in the outer cell layer is disturbed.
Approximately one-third of the wox2 Arabidopsis mutant embryos fail to correctly separate the protoderm layer by periclinal divisions in one or more cells of the prospective shoot region . Intriguingly, at a later stage, most of these mutants (92 %) recover and give rise to fertile plants [3, 16]. In our study, about one-third of the EEs from line 35S:WOX2i.4 showed an aberrant phenotype that is unique for 35S:WOX2i lines. Furthermore, the yield of MEs in the 35S:WOX2i lines is significantly reduced compared to the control. By using an inducible promoter we could show that down regulation of PaWOX2 when LEs had already developed has no effect on further development of the embryos. Therefore, the transient high expression of PaWOX2 during development of EEs and LEs is crucial for further development of LEs. The low recovering rate implies a more essential role for PaWOX2 during early embryo development in Norway spruce compared to that of AtWOX2 in Arabidopsis. However, an alternative scenario could be based on the previously suggested functional redundancy of AtWOX1, AtWOX3 and AtWOX5 with AtWOX2 . It has been shown that the redundant functions of AtWOX1, AtWOX3 and AtWOX5 become essential for normal development of the apical lineage in the absence of AtWOX2 activity . WOX1 homologs have not been found in Norway spruce . More research is required before we know if WOX genes in conifers exert redundant functions.
PaWOX2 is transiently expressed in the embryonal mass and in the upper part of the suspensor in early and early late embryos. Our results suggest that PaWOX2 is required for the establishment of embryonal mass and suspensor domains, expansion of suspensor cells and proper protoderm formation during early embryo development. EEs that fail to form a protoderm cannot develop further and degenerate. Taken together, our results suggest that WOX2 exerts a conserved role in protoderm development in both gymnosperms and angiosperm, whereas its function in the expansion of suspensor cells may be restricted to the gymnosperm lineage.
Prof. Ueli Grossniklaus and Dr. Nam-hai Chua are acknowledged for providing the vector pMDC7. Prof. Ueli Grossniklaus is acknowledged for providing the vector pMDC32. This work was financially supported by The Swedish governmental research program Trees and Crops for the Future (TC4F) received by SvA (salary for TZ), The Excellence Research program at The Swedish University of Agricultural Sciences received by SvA (salary for JA and JS), Knut and Alice Wallenberg Foundation and The Swedish Foundation for Strategic Research received by P Bozhkov (salary for PM).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Mayer U, Ruiz RAT, Berleth T, Misera S, Jurgens G. Mutations affecting body organization in the Arabidopsis embryo. Nature. 1991;353(6343):402–7.View ArticleGoogle Scholar
- von Arnold S, Clapham D. Spruce embryogenesis. Methods Mol Biol. 2008;427:31–47.View ArticleGoogle Scholar
- Haecker A, Gross-Hardt R, Geiges B, Sarkar A, Breuninger H, Herrmann M, et al. Expression dynamics of WOX genes mark cell fate decisions during early embryonic patterning in Arabidopsis thaliana. Development. 2004;131(3):657–68.PubMedView ArticleGoogle Scholar
- Nardmann J, Ji JB, Werr W, Scanlon MJ. The maize duplicate genes narrow sheath1 and narrow sheath2 encode a conserved homeobox gene function in a lateral domain of shoot apical meristems. Development. 2004;131(12):2827–39.PubMedView ArticleGoogle Scholar
- Goldberg RB, de Paiva G, Yadegari R. Plant embryogenesis: zygote to seed. Science. 1994;266(5185):605–14.PubMedView ArticleGoogle Scholar
- De Smet I, Beeckman T. Asymmetric cell division in land plants and algae: the driving force for differentiation. Nat Rev Mol Cell Bio. 2011;12(3):177–88.View ArticleGoogle Scholar
- Javelle M, Vernoud V, Rogowsky PM, Ingram GC. Epidermis: the formation and functions of a fundamental plant tissue. New Phytol. 2011;189(1):17–39.PubMedView ArticleGoogle Scholar
- Hardev S. Embryology of gymnosperms. Berlin: Gerbrüder Borntraeger; 1978.Google Scholar
- Ingouff M, Farbos I, Lagercrantz U, von Arnold S. PaHB1 is an evolutionary conserved HD-GL2 homeobox gene expressed in the protoderm during Norway spruce embryo development. Genesis. 2001;30(4):220–30.PubMedView ArticleGoogle Scholar
- Sabala I, Elfstrand M, Farbos I, Clapham D, von Arnold S. Tissue-specific expression of Pa18, a putative lipid transfer protein gene, during embryo development in Norway spruce (Picea abies). Plant Mol Biol. 2000;42(3):461–78.PubMedView ArticleGoogle Scholar
- Filonova LH, Bozhkov PV, Brukhin VB, Daniel G, Zhivotovsky B, von Arnold S. Two waves of programmed cell death occur during formation and development of somatic embryos in the gymnosperm, Norway spruce. J Cell Sci. 2000;113(Pt 24):4399–411.PubMedGoogle Scholar
- van der Graaff E, Laux T, Rensing SA. The WUS homeobox-containing (WOX) protein family. Genome Biol. 2009;10(12):248.PubMedPubMed CentralView ArticleGoogle Scholar
- Hedman H, Zhu T, von Arnold S, Sohlberg JJ. Analysis of the WUSCHEL-RELATED HOMEOBOX gene family in the conifer Picea abies reveals extensive conservation as well as dynamic patterns. BMC Plant Biol. 2013;13:89.PubMedPubMed CentralView ArticleGoogle Scholar
- Ueda M, Zhang Z, Laux T. Transcriptional activation of Arabidopsis axis patterning genes WOX8/9 links zygote polarity to embryo development. Dev Cell. 2011;20(2):264–70.PubMedView ArticleGoogle Scholar
- Wu X, Chory J, Weigel D. Combinations of WOX activities regulate tissue proliferation during Arabidopsis embryonic development. Dev Biol. 2007;309(2):306–16.PubMedPubMed CentralView ArticleGoogle Scholar
- Breuninger H, Rikirsch E, Hermann M, Ueda M, Laux T. Differential expression of WOX genes mediates apical-basal axis formation in the Arabidopsis embryo. Dev Cell. 2008;14(6):867–76.PubMedView ArticleGoogle Scholar
- Zhu T, Moschou PN, Alvarez JM, Sohlberg JJ, von Arnold S. WUSCHEL-RELATED HOMEOBOX 8/9 is important for proper embryo patterning in the gymnosperm Norway spruce. J Exp Bot. 2014;65(22):6543–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Nardmann J, Werr W. Symplesiomorphies in the WUSCHEL clade suggest that the last common ancestor of seed plants contained at least four independent stem cell niches. New Phytol. 2013;199(4):1081–92.PubMedView ArticleGoogle Scholar
- Palovaara J, Hakman I. WOX2 and polar auxin transport during spruce embryo pattern formation. Plant Signal Behav. 2009;4(2):153–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Palovaara J, Hallberg H, Stasolla C, Hakman I. Comparative expression pattern analysis of WUSCHEL-RELATED HOMEOBOX2 (WOX2) and WOX8/9 in developing seeds and somatic embryos of the gymnosperm Picea abies. New Phytol. 2010;188(1):122–35.PubMedView ArticleGoogle Scholar
- Högberg K-A, Ekberg I, Norell L, von Arnold S. Integration of somatic embryogenesis in a tree breeding programme – a case study with Picea abies. Can J For Res. 1998;28:1536–45.View ArticleGoogle Scholar
- Vestman D, Larsson E, Uddenberg D, Cairney J, Clapham D, Sundberg E, et al. Important processes during differentiation and early development of somatic embryos of Norway spruce as revealed by changes in global gene expression. Tree Genet Genomes. 2011;7(2):347–62.View ArticleGoogle Scholar
- Karlgren A, Carlsson J, Gyllenstrand N, Lagercrantz U, Sundstrom JF. Non-radioactive in situ hybridization protocol applicable for Norway spruce and a range of plant species. Journal of visualized experiments: JoVE 2009(26):1-13.Google Scholar
- Brand L, Horler M, Nuesch E, Vassalli S, Barrell P, Yang W, et al. A versatile and reliable two-component system for tissue-specific gene induction in Arabidopsis. Plant Physiol. 2006;141(4):1194–204.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuo J, Niu QW, Chua NH. Technical advance: an estrogen receptor-based transactivator XVE mediates highly inducible gene expression in transgenic plants. Plant J. 2000;24(2):265–73.PubMedView ArticleGoogle Scholar
- Curtis MD, Grossniklaus U. A gateway cloning vector set for high-throughput functional analysis of genes in planta. Plant Physiol. 2003;133(2):462–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.PubMedView ArticleGoogle Scholar
- Akin DE, Rigsby LL, Morrison WH. Oil Red as a histochemical stain for natural fibers and plant cuticle. Ind Crop Prod. 2004;19(2):119–24.View ArticleGoogle Scholar
- Alvarez JM, Sohlberg J, Engstrom P, Zhu T, Englund M, Moschou PN, et al. The WUSCHEL-RELATED HOMEOBOX 3 gene PaWOX3 regulates lateral organ formation in Norway spruce. New Phytol. 2015;208(4):1078–88.PubMedView ArticleGoogle Scholar
- Filonova LH, von Arnold S, Daniel G, Bozhkov PV. Programmed cell death eliminates all but one embryo in a polyembryonic plant seed. Cell Death Differ. 2002;9(10):1057–62.PubMedView ArticleGoogle Scholar
- Lie C, Kelsom C, Wu X. WOX2 and STIMPY-LIKE/WOX8 promote cotyledon boundary formation in Arabidopsis. Plant J. 2012;72(4):674–82.PubMedView ArticleGoogle Scholar
- Smertenko AP, Bozhkov PV, Filonova LH, von Arnold S, Hussey PJ. Re-organisation of the cytoskeleton during developmental programmed cell death in Picea abies embryos. Plant J. 2003;33(5):813–24.PubMedView ArticleGoogle Scholar
- Laux T, Jurgens G. Embryogenesis: a new start in life. Plant Cell. 1997;9(7):989–1000.PubMedPubMed CentralView ArticleGoogle Scholar
- Dodeman VL, Ducreux G, Kreis M. Zygotic embryogenesis versus somatic embryogenesis. J Exp Bot. 1997;48(313):1493–509.View ArticleGoogle Scholar
- Watanabe M, Tanaka H, Watanabe D, Machida C, Machida Y. The ACR4 receptor-like kinase is required for surface formation of epidermis-related tissues in Arabidopsis thaliana. Plant J. 2004;39(3):298–308.PubMedView ArticleGoogle Scholar
- Ito M, Sentoku N, Nishimura A, Hong SK, Sato Y, Matsuoka M. Position dependent expression of GL2-type homeobox gene, Roc1: significance for protoderm differentiation and radial pattern formation in early rice embryogenesis. Plant J. 2002;29(4):497–507.PubMedView ArticleGoogle Scholar
- Yeats TH, Rose JK. The formation and function of plant cuticles. Plant Physiol. 2013;163(1):5–20.PubMedPubMed CentralView ArticleGoogle Scholar