The Larix kaempferi genome reveals new insights into wood properties.

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Estimating the distribution and productivity characters of Larix kaempferi in response to climate change.

Understanding the distribution, net primary productivity (NPP) and environmental constraints of Larix kaempferi is crucial to predict how global climate change will affect its growth and future dynamics. We simulated future changes in the globally suitable distribution patterns and the NPP dynamics under different representative concentration pathways (RCPs) using MaxEnt and Physiological Principles in Predicting Growth (3-PG) models. The results showed that suitable distribution areas for Larix kaempferi were concentrated in Europe and Asia, followed by North America, under current climate conditions. Globally, about 33.75% of the suitable area was in China. Suitable areas decreased and shifted northward in Asia, Europe and China in the RCP scenarios. Larix kaempferi could adapt or move to higher latitudes/altitudes to mitigate the negative impacts of climate change. The NPP of Larix kaempferi in China was 241.85-863.57 g m-2 a-1 simulated by the 3-PG model after local parameterization, which was consistent with the measured NPP. Changes in NPP were predicted in future climates. When the correlations between climate factors and NPP were examined, under the more optimistic scenarios, NPP would increase significantly. The key parameters of the 3-PG model were the optimal temperature for growth, forest age, and the number of days of lost productivity in each frost period. Therefore, climate change has a quantitative and significant impact on the distribution and productivity of L. kaempferi, which was estimated successfully with the two modeling approaches. Our results will contribute to the improved cultivation, environment and management of L. kaempferi and potentially of other deciduous gymnosperms.

Multi-omics sequencing provides insight into floral transition in Catalpa bungei. C.A. Mey.

BACKGROUND: Floral transition plays an important role in development, and proper time is necessary to improve the value of valuable ornamental trees. The molecular mechanisms of floral transition remain unknown in perennial woody plants. "Bairihua" is a type of C. bungei that can undergo floral transition in the first planting year. RESULTS: Here, we combined short-read next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to provide a more complete view of transcriptome regulation during floral transition in C. bungei. The circadian rhythm-plant pathway may be the critical pathway during floral transition in early flowering (EF) C. bungei, according to horizontal and vertical analysis in EF and normal flowering (NF) C. bungei. SBP and MIKC-MADS-box were seemingly involved in EF during floral transition. A total of 61 hub genes were associated with floral transition in the MEturquoise model with Weighted Gene Co-expression Network Analysis (WGCNA). The results reveal that ten hub genes had a close connection with the GASA homologue gene (Cbu.gene.18280), and the ten co-expressed genes belong to five flowering-related pathways. Furthermore, our study provides new insights into the complexity and regulation of alternative splicing (AS). The ratio or number of isoforms of some floral transition-related genes is different in different periods or in different sub-genomes. CONCLUSIONS: Our results will be a useful reference for the study of floral transition in other perennial woody plants. Further molecular investigations are needed to verify our sequencing data.

Potential function of CbuSPL and gene encoding its interacting protein during flowering in Catalpa bungei.

BACKGROUND: "Bairihua", a variety of the Catalpa bungei, has a large amount of flowers and a long flowering period which make it an excellent material for flowering researches in trees. SPL is one of the hub genes that regulate both flowering transition and development. RESULTS: SPL homologues CbuSPL9 was cloned using degenerate primers with RACE. Expression studies during flowering transition in "Bairihua" and ectopic expression in Arabidopsis showed that CbuSPL9 was functional similarly with its Arabidopsis homologues. In the next step, we used Y2H to identify the proteins that could interact with CbuSPL9. HMGA, an architectural transcriptional factor, was identified and cloned for further research. BiFC and BLI showed that CbuSPL9 could form a heterodimer with CbuHMGA in the nucleus. The expression analysis showed that CbuHMGA had a similar expression trend to that of CbuSPL9 during flowering in "Bairihua". Intriguingly, ectopic expression of CbuHMGA in Arabidopsis would lead to aberrant flowers, but did not effect flowering time. CONCLUSIONS: Our results implied a novel pathway that CbuSPL9 regulated flowering development, but not flowering transition, with the participation of CbuHMGA. Further investments need to be done to verify the details of this pathway.

Quantitative Phosphoproteomic and Physiological Analyses Provide Insights into the Formation of the Variegated Leaf in Catalpa fargesii.

Variegated plants are valuable materials for investigating leaf color regulated mechanisms. To unveil the role of posttranslational modification in the variegated phenotype, we conducted global quantitative phosphoproteomic analysis on different leaf color sectors of Maiyuanjinqiu and the corresponding of Catalpa fargesii using Ti4+-IMAC phosphopeptide enrichment. A total of 3778 phosphorylated sites assigned to 1646 phosphoproteins were identified, and 3221 in 1434 proteins were quantified. Differential phosphoproteins (above 1.5 or below 1/1.5) in various leaf color sectors were selected for functional enrichment analyses. Gene ontology (GO) enrichment revealed that processes of photosynthesis, regulation of the generation of precursor metabolites, response to stress, homeostasis, amino acid metabolism, transport-related processes, and most of the energy metabolisms might contribute to leaf color. KEGG pathway enrichment analysis was performed based on differential phosphoproteins (DPs) in different organelles. The result showed that most enriched pathways were located in the chloroplasts and cytosol. The phosphorylation levels of glycometabolism enzymes might greatly affect leaf variegation. Measurements of fluorescence parameters and enzyme activities confirmed that protein phosphorylation could affect plant physiology by regulating enzyme activity. These results provide new clues for further study the formation mechanisms of naturally variegated phenotype.

Genome-wide analysis of long non-coding RNAs in Catalpa bungei and their potential function in floral transition using high-throughput sequencing.

BACKGROUND: Long non-coding RNAs (lncRNAs) have crucial roles in various biological regulatory processes. However, the study of lncRNAs is limited in woody plants. Catalpa bungei is a valuable ornamental tree with a long cultivation history in China, and a deeper understanding of the floral transition mechanism in C. bungei would be interesting from both economic and scientific perspectives. RESULTS: In this study, we categorized C. bungei buds from early flowering (EF) and normal flowering (NF) varieties into three consecutive developmental stages. These buds were used to systematically study lncRNAs during floral transition using high-throughput sequencing to identify molecular regulatory networks. Quantitative real-time PCR was performed to study RNA expression changes in different stages. In total, 12,532 lncRNAs and 26,936 messenger RNAs (mRNAs) were detected. Moreover, 680 differentially expressed genes and 817 differentially expressed lncRNAs were detected during the initiation of floral transition. The results highlight the mRNAs and lncRNAs that may be involved in floral transition, as well as the many lncRNAs serving as microRNA precursors. We predicted the functions of lncRNAs by analysing the relationships between lncRNAs and mRNAs. Seven lncRNA-mRNA interaction pairs may participate in floral transition. CONCLUSIONS: This study is the first to identify lncRNAs and their potential functions in floral transition, providing a starting point for detailed determination of the functions of lncRNAs in C. bungei.

Proteomic analysis of stress-related proteins and metabolic pathways in Picea asperata somatic embryos during partial desiccation.

Partial desiccation treatment (PDT) stimulates germination and enhances the conversion of conifer somatic embryos. To better understand the mechanisms underlying the responses of somatic embryos to PDT, we used proteomic and physiological analyses to investigate these responses during PDT in Picea asperata. Comparative proteomic analysis revealed that, during PDT, stress-related proteins were mainly involved in osmosis, endogenous hormones, antioxidative proteins, molecular chaperones and defence-related proteins. Compared with those in cotyledonary embryos before PDT, these stress-related proteins remained at high levels on days 7 (D7) and 14 (D14) of PDT. The proteins that differentially accumulated in the somatic embryos on D7 were mapped to stress and/or stimuli. They may also be involved in the glyoxylate cycle and the chitin metabolic process. The most significant difference in the differentially accumulated proteins occurred in the metabolic pathways of photosynthesis on D14. Furthermore, in accordance with the changes in stress-related proteins, analyses of changes in water content, abscisic acid, indoleacetic acid and H2 O2 levels in the embryos indicated that PDT is involved in water-deficit tolerance and affects endogenous hormones. Our results provide insight into the mechanisms responsible for the transition from morphologically mature to physiologically mature somatic embryos during the PDT process in P. asperata.

A computational framework for mapping the timing of vegetative phase change.

Phase change plays a prominent role in determining the form of growth and development. Although considerable attention has been focused on identifying the regulatory control mechanisms of phase change, a detailed understanding of the genetic architecture of this phenomenon is still lacking. We address this issue by deriving a computational model. The model is founded on the framework of functional mapping aimed at characterizing the interplay between quantitative trait loci (QTLs) and development through biologically meaningful mathematical equations. A multiphasic growth equation was implemented into functional mapping, which, via a series of hypothesis tests, allows the quantification of how QTLs regulate the timing and pattern of vegetative phase transition between independently regulated, temporally coordinated processes. The model was applied to analyze stem radial growth data of an interspecific hybrid family derived from two Populus species during the first 24 yr of ontogeny. Several key QTLs related to phase change have been characterized, most of which were observed to be in the adjacent regions of candidate genes. The identification of phase transition QTLs, whose expression is regulated by endogenous and environmental signals, may enhance our understanding of the evolution of development in changing environments.

A quantitative model of transcriptional differentiation driving host-pathogen interactions.

Despite our expanding knowledge about the biochemistry of gene regulation involved in host-pathogen interactions, a quantitative understanding of this process at a transcriptional level is still limited. We devise and assess a computational framework that can address this question. This framework is founded on a mixture model-based likelihood, equipped with functionality to cluster genes per dynamic and functional changes of gene expression within an interconnected system composed of the host and pathogen. If genes from the host and pathogen are clustered in the same group due to a similar pattern of dynamic profiles, they are likely to be reciprocally co-evolving. If genes from the two organisms are clustered in different groups, this means that they experience strong host-pathogen interactions. The framework can test the rates of change for individual gene clusters during pathogenic infection and quantify their impacts on host-pathogen interactions. The framework was validated by a pathological study of poplar leaves infected by fungal Marssonina brunnea in which co-evolving and interactive genes that determine poplar-fungus interactions are identified. The new framework should find its wide application to studying host-pathogen interactions for any other interconnected systems.

Transcriptomic and proteomic analyses of embryogenic tissues in Picea balfouriana treated with 6-benzylaminopurine.

The cytokinin 6-benzylaminopurine (6-BAP) influences the embryogenic capacity of the tissues of Picea balfouriana during long subculture (after 3 months). Tissues that proliferate in 3.6 and 5 µM 6-BAP exhibit the highest and lowest embryogenic capacity, respectively, generating 113 ± 6 and 23 ± 3 mature embryos per 100 mg of tissue. In this study, a comparative transcriptomic and proteomic approach was applied to characterize the genes and proteins that are differentially expressed among tissues under the influence of different levels of 6-BAP. A total of 51 375 unigenes and 2617 proteins were obtained after quality filtering. There were 2770 transcripts for proteins found among these unigenes. Gene ontology (GO) analysis of the differentially expressed unigenes and proteins showed that they were involved in cell and binding activity and were enriched in ribosome and glutathione metabolism pathways. Ribosomal proteins, glutathione S-transferase proteins, germin-like proteins and calmodulin-independent protein kinases were up-regulated in the embryogenic tissues with the highest embryogenic ability (treated with 3.6 µM 6-BAP), which was validated via quantitative real-time polymerase chain reaction (qRT-PCR) analysis, and these proteins might serve as molecular markers of embryogenic ability. Data are available via Sequence Read Archive (SRA) and ProteomeXchange with identifier SRP042246 and PXD001022, respectively.

Transcriptome analysis of callus from Picea balfouriana.

BACKGROUND: Picea likiangensis var. balfouriana (Rehd. et Wils.) Hillier ex Slavin (also known as Picea balfouriana) is an ecologically and economically important conifer that grows rapidly under optimum conditions and produces high-quality wood. It has a wide geographic distribution and is prevalent in southwest and eastern regions of China. Under suboptimal conditions, P. balfouriana grows slowly, which restricts its cultivation. Somatic embryogenesis has been used in the mass propagation of commercial species. However, low initiation rates are a common problem and the mechanisms involved in the induction of somatic embryogenesis are not fully understood. To understand the molecular mechanisms regulating somatic embryogenesis in P. balfouriana, high-throughput RNA-seq technology was used to investigate the transcriptomes of embryogenic and non-embryogenic tissues from three P. balfouriana genotypes. We compared the genes expressed in these tissues to identify molecular markers with embryogenic potential. RESULTS: A total of 55,078,846 nucleotide sequence reads were obtained for the embryogenic and non-embryogenic tissues of P. balfouriana, and 49.56% of them uniquely matched 22,295 (84.3%) of the 26,437 genes in the Picea abies genome database (Nature 497: 579-584, 2013). Differential gene expression analysis identified 1,418 differentially expressed genes (false discovery rate <0.0001; fold change ≥2) in the embryogenic tissues relative to the non-embryogenic tissues, including 431 significantly upregulated and 987 significantly downregulated genes. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis revealed that the most significantly altered genes were involved in plant hormone signal transduction, metabolic pathways (starch and sucrose metabolism), and phenylalanine metabolism. CONCLUSIONS: We found that the initiation of embryogenic tissues affected gene expression in many KEGG pathways, but predominantly in plant hormone signal transduction, plant-pathogen interaction, and starch and sucrose metabolism. The changes in multiple pathways related to induction in the P. balfouriana embryogenic tissues described here, will contribute to a more comprehensive understanding of the mechanisms involved in the initiation of somatic embryogenesis. Additionally, we found that somatic embryogenesis receptor kinase (SERK), arabinogalactan proteins, and members of the WUS-related homeobox protein family may play important roles and could act as molecular markers in the early stage of somatic embryogenesis, as reported previously.

Transcriptome and proteome profiling of adventitious root development in hybrid larch (Larix kaempferi × Larix olgensis).

BACKGROUND: Hybrids of larch (Larix kaempferi × Larix olgensis) are important afforestation species in northeastern China. They are routinely propagated via rooted stem cuttings. Despite the importance of rooting, little is known about the regulation of adventitious root development in larch hybrids. 454 GS FLX Titanium technology represents a new method for characterizing the transcriptomes of non-model species. This method can be used to identify differentially expressed genes, and then two-dimensional difference gel electrophoresis (2D-DIGE) and matrix-assisted laser desorption-ionization time-of-flight mass spectrometry (MALDI-TOF/TOF MS) analyses can be used to analyze their corresponding proteins. In this study, we analyzed semi-lignified cuttings of two clones of L. kaempferi × L. olgensis with different rooting capacities to study the molecular basis of adventitious root development. RESULTS: We analyzed two clones; clone 25-5, with strong rooting capacity, and clone 23-12, with weak rooting capacity. We constructed four cDNA libraries from 25-5 and 23-12 at two development stages. Sequencing was conducted using the 454 pyrosequencing platform. A total of 957832 raw reads was produced; 95.07% were high-quality reads, and were assembled into 45137 contigs and 61647 singletons. The functions of the unigenes, as indicated by their Gene Ontology annotation, included diverse roles in the molecular functions, biological processes, and cellular component categories. We analyzed 75 protein spots (-fold change ≥ 2, P ≤ 0.05) by 2D-DIGE, and identified the differentially expressed proteins using MALDI-TOF/TOF MS. A joint analysis of transcriptome and proteome showed genes related to two pathways, polyamine synthesis and stress response, might play an important role on adventitious root development. CONCLUSIONS: These results provide fundamental and important information for research on the molecular mechanism of adventitious root development. We also demonstrated for the first time the combined use of two important technologies as a powerful approach to advance research on non-model, but otherwise important, larch species.

Age-related trends in genetic parameters for Larix kaempferi and their implications for early selection.

BACKGROUND: Japanese larch (Larix kaempferi) has been introduced in China at the end of the 19th century, and as one successful exotic species, is becoming the preferred coniferous in northern China and sub-tropical alpine region. The rotation age is about 25-28 years for L. kaempferi as pulpwood in Henan province. Waiting for even one-half rotation age for final evaluation will be inefficient due to accumulated testing costs and delayed return on investment, which suggests that selection at an early age is highly desirable for L. kaempferi improvement programs in Henan province. In this study, we determined age trends of genetic parameters and evaluated early selection efficiency for L. kaempferi in Henan province to find out the appropriate trait for early selection and its selection age. RESULTS: Growth traits of 78 clones were measured periodically from age 2 to age 15 in a clonal trial of Larix kaempferi establishted at Son town, Henan Province. The genetic variation among clones, age-age correlations, and age trends in genetic parameters for growth traits were analyzed. Variant analysis revealed that tree height (HGT) and diameter at breast (DBH) were significant (1% level) among clones at every ages. The clonal repeatability of growth traits varied year-by-year, reaching the highest levels at different ages for different traits (0.77 at age 2 for HGT, 0.70 at age 5 for DBH and 0.66 from age 8 to age 10 for volume, respectively). The age-age genetic correlations ranged from 0.904 to 1.000 for HGT, and from 0943 to 1.000 for DBH. DBH at different ages was more genetically correlated to volume-15 than HGT. At the phenotypic level, HGT was always less correlated to volume-15 than DBH. With the estimates of efficiencies of early selection, the recommendation from present study was that the optimum age of early selection was age 2 for HGT and age 5 for DBH. CONCLUSIONS: Our study showed that there were significant (1% level) on growth traits among clones at every ages. The genetic parameters for growth traits varied from age to age. We found dual trait selection was more efficient than single trait selection for early selection.

A genome-wide survey of microRNA truncation and 3' nucleotide addition events in larch (Larix leptolepis).

MicroRNAs (miRNAs) play essential roles in numerous developmental and metabolic processes in animals and plants. Although the framework of miRNA biogenesis and function is established, the mechanism of miRNA degradation or modification remains to be investigated in plants. Mature miRNAs may be truncated or added nucleotides to generate variants. A detailed analysis of small RNA deep sequencing data sets resulted in the cloning of a large number of variants derived from larch miRNAs. Many 5'- and/or 3'-end truncated versions of miRNAs suggested that larch miRNAs might be degraded through either 5'-3' or 3'-5'. The relative abundance of variants truncated from 3'-end was higher than that of 5'-end for most miRNAs. The addition of adenine, uridine, and cytidine to the 3'-end of miRNAs was globally present, and the subtle variability in isomiR abundance might be regulated and biologically meaningful. It is the first report for cytidine addition in plant, and our examination of published small RNA deep sequencing data sets of Arabidopsis, rice, and moss suggests that cytidine addition to miRNA 3'-end exists broadly in plants. In addition, the nucleotide addition might be associated with 3'-5' miRNA degradation. Our results provide valuable information for a genome-wide survey of miRNA truncation and modification in larch or plants.

Dynamic expression of small RNA populations in larch (Larix leptolepis).

Small RNAs (sRNAs) are emerging as essential regulators of biological processes. However, several studies have reported that gymnosperms do not express appreciable amounts of 24-nt sRNAs, and conifers in particular may have a unique sRNA-silencing signature. Here, we compared the sRNA transcriptomes of Japanese larch somatic embryos (SE) and seedlings. SE sRNAs exhibited a length bias toward 24 nt, while seedlings showed a bias toward a 21-nt length. We also confirmed that larch is capable of producing 24-nt sRNAs based on a polyacrylamide gel analysis. The sRNA expression patterns varied according to developmental stage, which might be associated with Dicer-like 3 and RNA-dependent RNA polymerase2 (RDR2) levels. Our data suggest that many MIR loci that produce canonical microRNAs (miRNAs, 20-22 nt) and long sRNAs (23-26 nt) have dual functions; the latter were preferentially produced in SE compared to seedlings. However, the ratio of miRNAs to total sRNAs in seedlings was higher than in SE, and most miRNAs were upregulated in seedlings. Trans-acting small interfering RNAs (ta-siRNAs) generated from TAS3 triggered by miR390 were identified, and levels of the three detected ta-siRNAs peaked in mature embryos, which was not consistent with the lowest RDR6 level. These findings indicate that larch, and possibly other gymnosperms, shares a common sRNA pathway with other land plants, and that the sRNA distribution pattern varies according to developmental stage, which may be attributable to the expression of sRNA pathway genes.

[Mechanism on differential gene expression and heterosis formation].

Despite the rediscovery of heterosis about a century ago and the suggestion of various genetic models to explain this phenomenon, little consensus has yet been reached about the genetic basis of heterosis. Following the genome organization variation and gene effects, an understanding of gene differential expression in hybrids and its parents provides a new opportunity to speculate on mechanisms that might lead to heterosis. Investigation on allele-specific gene expression in hybrid and gene differential expression between hybrids and its parents might contribute to improve our understanding of the molecular basis of heterosis and eventually guide breeding practices. In this review, we discussed the recent researches on allelic-specific expression in hybrid which was frequently observed in recent studies and analyzed its regulatory mechanism. All possible modes of gene action, including additivity, high- and low-parent dominance, underdominance, and over-dominance, were observed when investigating gene differential expression between hybrids and its parents. Data from transcriptomic studies screened several heterosis-associated genes and highlighted the importance of certain key biochemical pathways that may prove to be quintessential for the manifestation of heterosis. So far, no uniform global expression pat-terns were observed in these gene expression studies. Most heterosis-associated gene expression analyses have not revealed a predominant functional category to which differentially expressed genes belong. However, these gene expression profiling studies represent a first step towards the definition of the complex gene expression networks that might be relevant in the context of heterosis. New technique on gene expression profile and advancements in bioinformatics will facilitate our understanding of the genetic basis of heterosis at the gene-expression level.

Sequencing the genome of Marssonina brunnea reveals fungus-poplar co-evolution.

BACKGROUND: The fungus Marssonina brunnea is a causal pathogen of Marssonina leaf spot that devastates poplar plantations by defoliating susceptible trees before normal fall leaf drop. RESULTS: We sequence the genome of M. brunnea with a size of 52 Mb assembled into 89 scaffolds, representing the first sequenced Dermateaceae genome. By inoculating this fungus onto a poplar hybrid clone, we investigate how M. brunnea interacts and co-evolves with its host to colonize poplar leaves. While a handful of virulence genes in M. brunnea, mostly from the LysM family, are detected to up-regulate during infection, the poplar down-regulates its resistance genes, such as nucleotide binding site domains and leucine rich repeats, in response to infection. From 10,027 predicted proteins of M. brunnea in a comparison with those from poplar, we identify four poplar transferases that stimulate the host to resist M. brunnea. These transferas-encoding genes may have driven the co-evolution of M. brunnea and Populus during the process of infection and anti-infection. CONCLUSIONS: Our results from the draft sequence of the M. brunnea genome provide evidence for genome-genome interactions that play an important role in poplar-pathogen co-evolution. This knowledge could help to design effective strategies for controlling Marssonina leaf spot in poplar.

Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis.

MicroRNAs (miRNAs) are emerging as essential regulators of biological processes. Somatic embryogenesis is one of the most important techniques for gymnosperm-breeding programs, but there is little understanding of its underlying mechanism. To investigate the roles of miRNAs during somatic embryogenesis in larch, we constructed a small RNA library from somatic embryos. High-throughput sequencing of the library identified 83 conserved miRNAs from 35 families, 16 novel miRNAs, and 14 plausible miRNA candidates, with a high proportion specific to larch or gymnosperms. qRT-PCR analysis demonstrated that both the conserved and novel or candidate miRNAs were expressed in larch. Several miRNA precursor sequences were obtained via RACE. We predicted 110 target genes using bioinformatics, and validated 9 of them by 5' RACE. 11 conserved miRNA families including 17 miRNAs with critical functions in plant development and six target mRNAs were detected by qRT-PCR in the larch SE. Stage-specific expression of miRNAs and their targets indicate their possible modulation on SE of larch: miR171a/b might exert function on PEMs, while miR171c acts in the induction process of larch SE; miR397 and miR398 mainly involved in modulation of PEM propagation and transition to single embryo; miR162 and miR168 exert their regulatory function during total SE process, especially during stages 5-8; miR156, miR159, miR160, miR166, miR167, and miR390 might play regulatory roles during cotyledonary embryo development. These findings indicate that larch and possibly other gymnosperms have complex mechanisms of gene regulation involving specific and common miRNAs operating post-transcriptionally during embryogenesis.

Transcriptome profiling and in silico analysis of somatic embryos in Japanese larch (Larix leptolepis).

UNLABELLED: Japanese larch (Larix leptolepis) is an ecologically and economically important species mainly grown in northeastern China, Japan and Europe. However, erratic flowering and poor germplasm resources caused by high embryo abortion rates have hampered breeding of Larix species. Somatic embryogenesis (SE) is an effective tool for the production of L. leptolepis with desirable characteristics, such as expression of totipotency, preparation of synthetic seeds, and genetic transformation. However, public genomic resources for this species are limited. We sequenced 591,759 raw expressed sequence tags (ESTs) from a 454 sequencing cDNA library of L. leptolepis somatic embryos, resulting in 572,403 high-quality reads. These reads were assembled into 70,927 unique sequences (UniGenes), including 32,321 contigs and 38,606 singletons. After removal of low-quality sequences, 65,115 UniGenes were annotated using the UniProtKB program. Based on their sequence similarity with known proteins, the matched 30,372 sequences from 664 species were estimated to represent approximately 19,000 unique genes. Gene ontology analysis revealed 21,324 UniGenes assigned to 51 categories. By Kyoto Encyclopedia of Genes and Genomes mapping, 25,773 transcripts were associated with 160 biochemical pathways. Further analysis screened four signal transduction pathways represented by 337 enzymes and 17 secondary metabolites. In silico analysis reveals that 207 UniESTs in Larix are homologous to MAPKs genes identified from other model plants, which may be involved in regulating SE development. This study provides an initial insight into the Larix transcriptomes of the pro-embryogenic mass and is a sound basis for future studies. KEY MESSAGE: We constructed a large, full-length 454 sequencing cDNA library of Larix leptolepis during somatic embryogenesis. More than 590,000 sequences were obtained and a deep-coverage EST database was constructed.

Integrated transcriptomic and metabolic analyses provide insights into the maintenance of embryogenic potential and the biosynthesis of phenolic acids and flavonoids involving transcription factors in Larix kaempferi (Lamb.) Carr.

Somatic embryogenesis (SE) techniques have been established for micropropagation or basic research related to plant development in many conifer species. The frequent occurrence of non-embryogenic callus (NEC) during SE has impose constraints on the application of somatic embryogenesis SE in Larix kaempferi (Lamb.) Carr, but the potential regulatory mechanisms are poorly understood. In this study, integrated transcriptomic and metabolomic analyses were performed in embryogenic callus (EC) and NEC originating from a single immature zygotic embryo to better decipher the key molecular and metabolic mechanisms required for embryogenic potential maintenance. The results showed that a total of 13,842 differentially expressed genes (DEGs) were found in EC and NEC, among which many were enriched in plant hormone signal transduction, starch and sucrose metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis, and the biosynthesis of amino acids pathways. Metabolite profiling showed that 441 differentially accumulated metabolites (DAMs) were identified in EC and NEC. Both EC and NEC had vigorous primary metabolic activities, while most secondary metabolites were upregulated in NEC. Many totipotency-related transcription factor (TF) genes such as BBMs, WUSs, and LEC1 showed higher expression levels in EC compared with NEC, which may result in the higher accumulation of indole 3-acetic acid (IAA) in EC. NEC was characterized by upregulation of genes and metabolites associated with stress responses, such as DEGs involved in jasmonic acid (JA) and ethylene (ETH) biosynthesis and signal transduction pathways, and DEGs and DAMs related to phenylpropanoid and flavonoid biosynthesis. We predicted and analyzed TFs that could target several key co-expressed structural DEGs including two C4H genes, two CcoAOMT genes and three HCT genes involved in phenylpropanoid and flavonoid biosynthesis. Based on the targeted relationship and the co-expression network, two ERFs (Lk23436 and Lk458687), one MYB (Lk34626) and one C2C2-dof (Lk37167) may play an important role in regulating phenolic acid and flavonoid biosynthesis by transcriptionally regulating the expression of these structural genes. This study shows an approach involving integrated transcriptomic and metabolic analyses to obtain insights into molecular events underlying embryogenic potential maintenance and the biosynthesis mechanisms of key metabolites involving TF regulation, which provides valuable information for the improvement of SE efficiency in L. kaempferi.