Transcript profiling for early stages during embryo development in Scots pine

Arabidopsis ribosomal RNA processing meerling mutants exhibit suspensor-derived polyembryony due to direct reprogramming of the suspensor

Embryo development in Arabidopsis (Arabidopsis thaliana) starts off with an asymmetric division of the zygote to generate the precursors of the embryo proper and the supporting extraembryonic suspensor. The suspensor degenerates as the development of the embryo proper proceeds beyond the heart stage. Until the globular stage, the suspensor maintains embryonic potential and can form embryos in the absence of the developing embryo proper. We report a mutant called meerling-1 (mrl-1), which shows a high penetrance of suspensor-derived polyembryony due to delayed development of the embryo proper. Eventually, embryos from both apical and suspensor lineages successfully develop into normal plants and complete their life cycle. We identified the causal mutation as a genomic rearrangement altering the promoter of the Arabidopsis U3 SMALL NUCLEOLAR RNA-ASSOCIATED PROTEIN 18 (UTP18) homolog that encodes a nucleolar-localized WD40-repeat protein involved in processing 18S preribosomal RNA. Accordingly, root-specific knockout of UTP18 caused growth arrest and accumulation of unprocessed 18S pre-rRNA. We generated the mrl-2 loss-of-function mutant and observed asynchronous megagametophyte development causing embryo sac abortion. Together, our results indicate that promoter rearrangement decreased UTP18 protein abundance during early stage embryo proper development, triggering suspensor-derived embryogenesis. Our data support the existence of noncell autonomous signaling from the embryo proper to prevent direct reprogramming of the suspensor toward embryonic fate.

Identification of Metabolic Pathways Differentially Regulated in Somatic and Zygotic Embryos of Maritime Pine

Embryogenesis is a complex phase of conifer development involving hundreds of genes, and a proper understanding of this process is critical not only to produce embryos with different applied purposes but also for comparative studies with angiosperms. A global view of transcriptome dynamics during pine somatic and zygotic embryogenesis is currently missing. Here, we present a genome-wide transcriptome analysis of somatic and zygotic embryos at three developmental stages to identify conserved biological processes and gene functions during late embryogenesis. Most of the differences became more significant as the developmental process progressed from early to cotyledonary stages, and a higher number of genes were differentially expressed in somatic than in zygotic embryos. Metabolic pathways substantially affected included those involved in amino acid biosynthesis and utilization, and this difference was already observable at early developmental stages. Overall, this effect was found to be independent of the line (genotype) used to produce the somatic embryos. Additionally, transcription factors differentially expressed in somatic versus zygotic embryos were analyzed. Some potential hub regulatory genes were identified that can provide clues as to what transcription factors are controlling the process and to how the observed differences between somatic and zygotic embryogenesis in conifers could be regulated.

Transcriptome analysis reveals key developmental and metabolic regulatory aspects of oil palm (Elaeis guineensis Jacq.) during zygotic embryo development

BACKGROUND

Oil palm is the most efficient oil-producing crop in the world, and the yield of palm oil is associated with embryonic development. However, a comprehensive understanding of zygotic embryo development at the molecular level remains elusive. In order to address this issue, we report the transcriptomic analysis of zygotic embryo development in oil palm, specifically focusing on regulatory genes involved in important biological pathways.

RESULTS

In this study, three cDNA libraries were prepared from embryos at S1 (early-stage), S2 (middle-stage), and S3 (late-stage). There were 16,367, 16,500, and 18,012 genes characterized at the S1, S2, and S3 stages of embryonic development, respectively. A total of 1522, 2698, and 142 genes were differentially expressed in S1 vs S2, S1 vs S3, and S2 vs S3, respectively. Using Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to identify key genes and pathways. In the hormone signaling pathway, genes related to auxin antagonize the output of cytokinin which regulates the development of embryo meristem. The genes related to abscisic acid negatively regulating the synthesis of gibberellin were strongly up-regulated in the mid-late stage of embryonic development. The results were reported the early synthesis and mid-late degradation of sucrose, as well as the activation of the continuous degradation pathway of temporary starch, providing the nutrients needed for differentiation of the embryonic cell. Moreover, the transcripts of genes involved in fatty acid synthesis were also abundantly accumulated in the zygotic embryos.

CONCLUSION

Taken together, our research provides a new perspective on the developmental and metabolic regulation of zygotic embryo development at the transcriptional level in oil palm.

Transcriptomic Time-Series Analyses of Gene Expression Profile During Zygotic Embryo Development in Picea mongolica

Zygotic embryogenesis is a critical process during seed development in gymnosperms. However, knowledge on the genome-wide transcriptional activation that guides this process in conifers is limited, especially in Picea mongolica. This tree species is endemic to semiarid habitats of Inner Mongolia in China. To extend what is known about the molecular events underpinning its zygotic embryogenesis, comparative transcriptomic analyses of gene expression in zygotic embryos were performed by RNA sequencing in P. mongolica. Our results showed that most changes in transcript levels occurred in the early embryonic pattering determination and formation of mature embryos. Transcripts related to embryogenic competence, cell division pattern, hormones, and stress response genes were identified during embryogenesis. Auxin is essential for early embryo patterning and pre-cotyledon embryonic formation. However, ABA is a major regulator of embryo maturation. Moreover, we found that methylation-related gene expression is associated with activation of early-stage embryos, late embryogenesis abundant proteins, and storage/energy-related genes with late and mature embryos. Furthermore, network analysis revealed stage-specific and multistage gene expression clusters during embryogenesis. WOX, MYB, AP2, and HLH transcription factors seem more relevant to embryogenesis in different stages. Our results provide large-scale and comprehensive transcriptome data for embryo development in P. mongolica. These data will lay a foundation for the protection and utilization of P. mongolica resources.

Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

Despite their ecological and economical importance, conifers genomic resources are limited, mainly due to the large size and complexity of their genomes. Additionally, the available genomic resources lack complete structural and functional annotation. Transcriptomic resources have been commonly used to compensate for these deficiencies, though for most conifer species they are limited to a small number of tissues, or capture only a fraction of the genes present in the genome. Here we provide an atlas of gene expression patterns for conifer Pinus sylvestris across five tissues: embryo, megagametophyte, needle, phloem and vegetative bud. We used a wide range of tissues and focused our analyses on the expression profiles of genes at tissue level. We provide comprehensive information of the per-tissue normalized expression level, indication of tissue preferential upregulation and tissue-specificity of expression. We identified a total of 48,001 tissue preferentially upregulated and tissue specifically expressed genes, of which 28% have annotation in the Swiss-Prot database. Even though most of the putative genes identified do not have functional information in current biological databases, the tissue-specific patterns discovered provide valuable information about their potential functions for further studies, as for example in the areas of plant physiology, population genetics and genomics in general. As we provide information on tissue specificity at both diploid and haploid life stages, our data will also contribute to the understanding of evolutionary rates of different tissue types and ploidy levels.

In silico characterization of putative gene homologues involved in somatic embryogenesis suggests that some conifer species may lack LEC2, one of the key regulators of initiation of the process

BACKGROUND

Somatic embryogenesis (SE) is the process in which somatic embryos develop from somatic tissue in vitro on medium in most cases supplemented with growth regulators. Knowledge of genes involved in regulation of initiation and of development of somatic embryos is crucial for application of SE as an efficient tool to enable genetic improvement across genotypes by clonal propagation.

RESULTS

Current work presents in silico identification of putative homologues of central regulators of SE initiation and development in conifers focusing mainly on key transcription factors (TFs) e.g. BBM, LEC1, LEC1-LIKE, LEC2 and FUSCA3, based on sequence similarity using BLASTP. Protein sequences of well-characterised candidates genes from Arabidopsis thaliana were used to query the databases (Gymno PLAZA, Congenie, GenBank) including whole-genome sequence data from two representative species from the genus Picea (Picea abies) and Pinus (Pinus taeda), for finding putative conifer homologues, using BLASTP. Identification of corresponding conifer proteins was further confirmed by domain search (Conserved Domain Database), alignment (MUSCLE) with respective sequences of Arabidopsis thaliana proteins and phylogenetic analysis (Phylogeny.fr).

CONCLUSIONS

This in silico analysis suggests absence of LEC2 in Picea abies and Pinus taeda, the conifer species whose genomes have been sequenced. Based on available sequence data to date, LEC2 was also not detected in the other conifer species included in the study. LEC2 is one of the key TFs associated with initiation and regulation of the process of SE in angiosperms. Potential alternative mechanisms that might be functional in conifers to compensate the lack of LEC2 are discussed.

Regulation of Somatic Embryo Development in Norway Spruce

Somatic embryogenesis in Norway spruce combined with reverse genetics can be used as a model to study the regulation of embryo development in conifers. The somatic embryo system includes a sequence of developmental stages, which are similar in morphology to their zygotic counterparts. The system can be sufficiently synchronized to enable the collection and study of a large number of somatic embryos at each developmental stage.Here we describe a protocol for establishing transgenic cell lines in which genes of interest are upregulated or downregulated. Furthermore, we present methods for comparing embryo morphology and development in transgenic and control cell lines, including phenotyping the embryos, histological analysis, and tracking embryo development. The expression pattern of different genes is determined by GUS reporter assays.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Comprehensive assembly and analysis of the transcriptome of maritime pine developing embryos

BACKGROUND

There are clear differences in embryo development between angiosperm and gymnosperm species. Most of the current knowledge on gene expression and regulation during plant embryo development has derived from studies on angiosperms species, in particular from the model plant Arabidopsis thaliana. The few published studies on transcript profiling of conifer embryogenesis show the existence of many putative embryo-specific transcripts without an assigned function. In order to extend the knowledge on the transcriptomic expression during conifer embryogenesis, we sequenced the transcriptome of zygotic embryos for several developmental stages that cover most of Pinus pinaster (maritime pine) embryogenesis.

RESULTS

Total RNA samples collected from five zygotic embryo developmental stages were sequenced with Illumina technology. A de novo transcriptome was assembled as no genome sequence is yet published for Pinus pinaster. The transcriptome of reference for the period of zygotic embryogenesis in maritime pine contains 67,429 transcripts, which likely encode 58,527 proteins. The annotation shows a significant percentage, 31%, of predicted proteins exclusively present in pine embryogenesis. Functional categories and enrichment analysis of the differentially expressed transcripts evidenced carbohydrate transport and metabolism over-representation in early embryo stages, as highlighted by the identification of many putative glycoside hydrolases, possibly associated with cell wall modification, and carbohydrate transport transcripts. Moreover, the predominance of chromatin remodelling events was detected in early to middle embryogenesis, associated with an active synthesis of histones and their post-translational modifiers related to increased transcription, as well as silencing of transposons.

CONCLUSIONS

Our results extend the understanding of gene expression and regulation during zygotic embryogenesis in conifers and are a valuable resource to support further improvements in somatic embryogenesis for vegetative propagation of conifer species. Specific transcripts associated with carbohydrate metabolism, monosaccharide transport and epigenetic regulation seem to play an important role in pine early embryogenesis and may be a source of reliable molecular markers for early embryogenesis.

Metabolome and transcriptome profiling reveal new insights into somatic embryo germination in Norway spruce (Picea abies)

Transcriptome, metabolome and histological profiling were performed on normal and aberrant somatic embryo germinants of Norway spruce (Picea abies L. Karst) providing a simplistic systems biology description of conifer germination. Aberrant germinants (AGs) formed periderm-like tissue at the apical pole and lacked shoot growth above the cotyledons. Transcriptome profiling (RNA-Sequencing) revealed a total of 370 differentially expressed genes at ≥1 or ≤-1 log2-fold change, where 92% were down-regulated in AGs compared with normal germinants (NGs). Genes associated with shoot apical meristem formation were down-regulated in AGs, or not differentially expressed between AGs and NGs. Genes involved in hormone signaling and transport were also down-regulated. Metabolite profiling by gas chromatography-mass spectrometry (MS) and liquid chromatography-MS revealed biochemical difference between AGs and NGs, notably increased levels of sugars including glucose in AGs. Genes involved in glucose signaling were down-regulated and genes involved in starch biosynthesis were up-regulated, suggesting involvement of sugar signaling during late embryo development and germination. The overall results provide new data enabling further studies to confirm potential markers for a normal germination process in conifers.