Integration analysis of transcriptome and proteome profiles brings new insights of somatic embryogenesis of two eucalyptus species.

BACKGROUND: Somatic embryogenesis (SE) is recognized as a promising technology for plant vegetative propagation. Although previous studies have identified some key regulators involved in the SE process in plant, our knowledge about the molecular changes in the SE process and key regulators associated with high embryogenic potential is still poor, especially in the important fiber and energy source tree - eucalyptus. RESULTS: In this study, we analyzed the transcriptome and proteome profiles of E. camaldulensis (with high embryogenic potential) and E. grandis x urophylla (with low embryogenic potential) in SE process: callus induction and development. A total of 12,121 differentially expressed genes (DEGs) and 3,922 differentially expressed proteins (DEPs) were identified in the SE of the two eucalyptus species. Integration analysis identified 1,353 (131 to 546) DEGs/DEPs shared by the two eucalyptus species in the SE process, including 142, 13 and 186 DEGs/DEPs commonly upregulated in the callus induction, maturation and development, respectively. Further, we found that the trihelix transcription factor ASR3 isoform X2 was commonly upregulated in the callus induction of the two eucalyptus species. The SOX30 and WRKY40 TFs were specifically upregulated in the callus induction of E. camaldulensis. Three TFs (bHLH62, bHLH35 isoform X2, RAP2-1) were specifically downregulated in the callus induction of E. grandis x urophylla. WGCNA identified 125 and 26 genes/proteins with high correlation (Pearson correlation > 0.8 or < -0.8) with ASR3 TF in the SE of E. camaldulensis and E. grandis x urophylla, respectively. The potential target gene expression patterns of ASR3 TF were then validated using qRT-PCR in the material. CONCLUSIONS: This is the first time to integrate multiple omics technologies to study the SE of eucalyptus. The findings will enhance our understanding of molecular regulation mechanisms of SE in eucalyptus. The output will also benefit the eucalyptus breeding program.

Identification of MYB gene family in medicinal tea tree Melaleuca alternifolia (Maiden and Betche) cheel and analysis of members regulating terpene biosynthesis.

BACKGROUND: The tea tree (Melaleuca alternifolia) is renowned for its production of tea tree oil, an essential oil primarily composed of terpenes extracted from its shoot. MYB transcription factors, which are one of the largest TF families, play a crucial role in regulating primary and secondary metabolite synthesis. However, knowledge of the MYB gene family in M. alternifolia is limited. METHODS AND RESULTS: Here, we conducted a comprehensive genome-wide analysis of MYB genes in M. alternifolia, referred to as MaMYBs, including phylogenetic relationships, structures, promoter regions, and GO annotations. Our findings classified 219 MaMYBs into four subfamilies: one 5R-MYB, four 3R-MYBs, sixty-one MYB-related, and the remaining 153 are all 2R-MYBs. Seven genes (MYB189, MYB146, MYB44, MYB29, MYB175, MYB162, and MYB160) were linked to terpenoid synthesis based on GO annotation. Phylogenetic analysis with Arabidopsis homologous MYB genes suggested that MYB193 and MYB163 may also be involved in terpenoid synthesis. Additionally, through correlation analysis of gene expression and metabolite content, we identified 42 MYB genes associated with metabolite content. CONCLUSION: The results provide valuable insights into the importance of MYB transcription factors in essential oil production in M. alternifolia. These findings lay the groundwork for a better understanding of the MYB regulatory network and the development of novel strategies to enhance essential oil synthesis in M. alternifolia.

Proteome profiles during early stage of somatic embryogenesis of two Eucalyptus species.

BACKGROUND: Somatic embryogenesis (SE) was recognized as an important tool for plants to propagate. However, our knowledge about the proteins involved in early SE including the callus dedifferentiation is still limited, especially in the economic woody tree - Eucalyptus. RESULTS: We used the data-independent acquisition mass-spectrometry to study the different proteome profiles of early SE of two Eucalyptus species-E. camaldulensis (high regeneratively potential) and E. grandis x urophylla (low regenerative potential). Initially, 35,207 peptides and 7,077 proteins were identified in the stem and tissue-culture induced callus of the two Eucalyptus species. MSstat identified 2,078 and 2,807 differentially expressed proteins (DEPs) in early SE of E. camaldulensis and E. grandis x urophylla, respectively. They shared 760 upregulated and 420 downregulated proteins, including 4 transcription factors, 31 ribosomal proteins, 1 histone, 3 zinc finger proteins (ZFPs), 16 glutathione transferases, 10 glucosyltransferases, ARF19, WOX8 and PIN1. These proteins might be involved in the early SE of Eucalyptus. By combining the miRNA and RNA-Seq results, some miRNA ~ gene/protein regulatory networks were identified in early SE of Eucalyptus, such as miR160 ~ TPP2, miR164 ~ UXS2, miR169 ~ COX11 and miR535 ~ Eucgr.E01067. Further, we found SERK, WRKY, ZFP and ABC transporter might be related with high SE potential. CONCLUSIONS: Overall, our study identified proteins involved in the early SE and related to the high regeneration potential of Eucalyptus. It greatly enhanced our understanding of the early SE and the SE capacity of Eucalyptus.

Gene expression programs during callus development in tissue culture of two Eucalyptus species.

BACKGROUND: Eucalyptus is a highly diverse genus of the Myrtaceae family and widely planted in the world for timber and pulp production. Tissue culture induced callus has become a common tool for Eucalyptus breeding, however, our knowledge about the genes related to the callus maturation and shoot regeneration is still poor. RESULTS: We set up an experiment to monitor the callus induction and callus development of two Eucalyptus species - E. camaldulensis (high embryogenic potential) and E. grandis x urophylla (low embryogenic potential). Then, we performed transcriptome sequencing for primary callus, mature callus, shoot regeneration stage callus and senescence callus. We identified 707 upregulated and 694 downregulated genes during the maturation process of the two Eucalyptus species and most of them were involved in the signaling pathways like plant hormone and MAPK. Next, we identified 135 and 142 genes that might play important roles during the callus development of E. camaldulensis and E. grandis x urophylla, respectively. Further, we found 15 DEGs shared by these two Eucalyptus species during the callus development, including Eucgr.D00640 (stem-specific protein TSJT1), Eucgr.B00171 (BTB/POZ and TAZ domain-containing protein 1), Eucgr.C00948 (zinc finger CCCH domain-containing protein 20), Eucgr.K01667 (stomatal closure-related actinbinding protein 3), Eucgr.C00663 (glutaredoxin-C10) and Eucgr.C00419 (UPF0481 protein At3g47200). Interestingly, the expression patterns of these genes displayed "N" shape in the samples. Further, we found 51 genes that were dysregulated during the callus development of E. camaldulensis but without changes in E. grandis x urophylla, such as Eucgr.B02127 (GRF1-interacting factor 1), Eucgr.C00947 (transcription factor MYB36), Eucgr.B02752 (laccase-7), Eucgr.B03985 (transcription factor MYB108), Eucgr.D00536 (GDSL esterase/lipase At5g45920) and Eucgr.B02347 (scarecrow-like protein 34). These 51 genes might be associated with the high propagation ability of Eucalyptus and 22 might be induced after the dedifferentiation. Last, we performed WGCNA to identify the co-expressed genes during the callus development of Eucalyptus and qRT-PCR experiment to validate the gene expression patterns. CONCLUSIONS: This is the first time to globally study the gene profiles during the callus development of Eucalyptus. The results will improve our understanding of gene regulation and molecular mechanisms in the callus maturation and shoot regeneration.

The complete chloroplast genome of Castanopsis hystrix Hook. f. & Thomson ex A. DC. 1863 (Fagaceae).

Castanopsis hystrix Hook. f. & Thomson ex A. DC. 1863 (Fagaceae) is an evergreen broad-leaved tree with high economic and ecological value. In this study, the complete chloroplast genome of C. hystrix was sequenced, assembled and annotated. The plastome (plastid genome) of C. hystrix was 160,624 bp in size, consisting of a pair of inverted repeats (IRs, 25,699 bp), a large-single-copy (LSC, 90,276 bp) region, and a small-single-copy (SSC, 18,950 bp). The overall GC content of C. hystrix was 36.8%. A total of 133 genes were annotated, including 88 protein-coding genes (PCG), 37 transfer RNA genes (tRNA), and eight ribosomal RNA genes (rRNA). A maximum likelihood analysis showed that the Castanopsis species form a monophyletic clade. C. hystrix is most closely related to C. tibetana with 100% bootstrap support value. The result enriches the genomic data for the genus Castanopsis, which will contribute to future studies in phylogenetics and evolution.

The complete chloroplast genome sequence of chaya (Cnidoscolus aconitifolius) and phylogenetic analysis.

Cnidoscolus aconitifolius is a leafy green heathy vegetable and medicinal plant belongs to the family Euphorbiaceae. In the present study, we sequenced the complete chloroplast genome of C. aconitifolius, which is 158,658 bp in length and consisted of two copies of inverted repeat (IR) of 26,982 bp separated by a large single copy (LSC) of 87,022 bp and a small single copy (SSC) of 17,672 bp. The GC content of C. aconitifolius was 36.3%. A total of 130 genes were predicted, including 86 protein-coding genes, 36 tRNAs and 8 rRNAs. The plastid phylogenomic analysis support C. aconitifolius is closely related to Manihot esculenta.

The complete chloroplast genome sequence and phylogenetic analysis of Nanhaia speciosa (Fabaceae).

Nanhaia speciosa (Fabaceae) is a woody perennial vine used as an important traditionally Chinese medicine. In this study, the complete chloroplast genome of Nanhaia speciosa was sequenced and assembled. The chloroplast genome of N. speciosa was 132,551bp in length including only one copy of the inverted repeat (IR). It encoded a total of 110 genes, containing 76 protein-coding genes, 30 tRNA and 4 rRNA. The overall GC content was 34.1%. Phylogenetic analysis using a matrix of 69 protein-coding genes illustrated that N. speciosa is most closely related to Wisteriopsis reticulata of tribe Wisterieae.

Characterization of the complete chloroplast genome of the Taiwan alder Alnus formosana (Betulaceae) based on next-generation sequencing technology.

Alnus formosana (Betulaceae) is an important ecological and economic deciduous tree species widely distributed throughout subtropical regions of Taiwan province, China. At the present study, the complete chloroplast genome of A. formosana was assumbled using next-generation sequencing technology. The complete chloroplast sequence is 161,029 bp in length, which consisted of a large single copy (LSC, 89,720 bp) and a small single copy (SSC; 19,205 bp) separated a pair of inverted repeats (IRs; 26,052 bp). The overall guanine-cytosine (GC) content was 36.4%. A total of 131 genes were annotated, including 85 protein-coding genes, 37 tRNAs, eight rRNAs and one pseudogene (ψycf1). The phylogenetic analysis fully resolved A. formosana in a clade with A. japonica. The plastome of A. formosana will provide informative genomic resources for further phylogenetic application and genetic improvement.

Genome-wide identification of microRNAs involved in the somatic embryogenesis of Eucalyptus.

MicroRNAs (miRNAs) are small noncoding RNAs (18-24 nt) and function in many biological processes in plants. Although Eucalyptus trees are widely planted across the world, our understanding of the miRNA regulation in the somatic embryogenesis (SE) of Eucalyptus is still poor. Here we reported, for the first time, the miRNA profiles of differentiated and dedifferentiated tissues of two Eucalyptus species and identified miRNAs involved in SE of Eucalyptus. Stem and tissue culture-induced callus were obtained from the subculture seedlings of E. camaldulensis and E. grandis x urophylla and were used as differentiated and dedifferentiated samples, respectively. Small RNA sequencing generated 304.2 million clean reads for the Eucalyptus samples (n = 3) and identified 888 miRNA precursors (197 known and 691 novel) for Eucalyptus. These miRNAs were mainly distributed in chromosomes Chr03, Chr05, and Chr08 and can produce 46 miRNA clusters. Then, we identified 327 and 343 differentially expressed miRNAs (DEmiRs) in the dedifferentiation process of E. camaldulensis and E. grandis x urophylla, respectively. DEmiRs shared by the two Eucalyptus species might be involved in the development of embryonic callus, such as MIR156, MIR159, MIR160, MIR164, MIR166, MIR169, MIR171, MIR399, and MIR482. Notably, we identified 81 upregulated and 67 downregulated miRNAs specific to E. camaldulensis, which might be associated with the high embryogenic potential. Target prediction and functional analysis showed that they might be involved in longevity regulating and plant hormone signal transduction pathways. Further, using the gene expression profiles, we observed the negative regulation of miRNA-target pairs, such as MIR160~ARF18, MIR396~GRF6, MIR166~ATHB15/HD-ZIP, and MIR156/MIR157~SPL1. Interestingly, transcription factors such as WRKY, MYB, GAMYB, TCP4, and PIL1 were found to be regulated by the DEmiRs. The genes encoding PIL1 and RPS21C, regulated by upregulated miRNAs (e.g., egd-N-miR63-5p, egd-N-miR63-5p, and MIR169,) were downregulated exclusively in the dedifferentiation of E. camaldulensis. This is the first time to study the miRNA regulation in the dedifferentiation process of Eucalyptus and it will provide a valuable resource for future studies. More importantly, it will improve our understanding of miRNA regulation during the somatic embryogenesis of Eucalyptus and benefit the Eucalyptus breeding program.