Transcriptome analysis identifies genes involved in the somatic embryogenesis of Eucalyptus

Genome-Wide Identification and Expression Analysis of Nitrate Transporter (NRT) Gene Family in Eucalyptus grandis

Eucalyptus grandis is an important planted hardwood tree worldwide with fast growth and good wood performance. The nitrate transporter (NRT) gene family is a major core involved in nitrogen (N) absorption and utilization in plants, but the comprehensive characterization of NRT genes in E. grandis remains mostly elusive. In this study, a total of 75 EgNRT genes were identified from the genome of E. grandis that were distributed unevenly across ten chromosomes, except Chr9. A phylogenetic analysis showed that the EgNRT proteins could be divided into three classes, namely NRT1, NRT2 and NRT3, which contained 69, 4 and 2 members, respectively. The cis-regulatory elements in the promoter regions of EgNRT genes were mainly involved in phytohormone and stress response. The transcriptome analysis indicated that the differentially expressed genes of leaf and root in E. grandis under different N supply conditions were mainly involved in the metabolic process and plant hormone signal transduction. In addition, the transcriptome-based and RT-qPCR analysis revealed that the expression of 13 EgNRT genes, especially EgNRT1.3, EgNRT1.38, EgNRT1.39 and EgNRT1.52, was significantly upregulated in the root under low-N-supply treatment, suggesting that those genes might play a critical role in root response to nitrate deficiency. Taken together, these results would provide valuable information for characterizing the roles of EgNRTs and facilitate the clarification of the molecular mechanism underlying EgNRT-mediated N absorption and distribution in E. grandis.

Reactive oxygen species act as the key signaling molecules mediating light-induced anthocyanin biosynthesis in Eucalyptus

Light plays a pivotal role in regulating anthocyanin biosynthesis in plants, and the early light-responsive signals that initiate anthocyanin biosynthesis remain to be elucidated. In this study, we showed that the anthocyanin biosynthesis in Eucalyptus is hypersensitive to increased light intensity. The combined transcriptomic and metabolomic analyses were conducted on Eucalyptus leaves after moderate (ML; 100 μmol m-2 s-1) and high (HL; 300 μmol m-2 s-1) light intensity treatments. The results identified 1940, 1096, 1173, and 2756 differentially expressed genes at 6, 12, 24, and 36 h after HL treatment, respectively. The metabolomic results revealed the primary anthocyanin types, and other differentially accumulated flavonoids and phenylpropane intermediates that were produced in response to HL, which well aligned with the transcriptome results. Moreover, biochemical analysis showed that HL inhibited peroxidase activity and increased the ROS level in Eucalyptus leaves. ROS depletion through co-application of the antioxidants rutin, uric acid, and melatonin significantly reduced, and even abolished, anthocyanin biosynthesis induced by HL treatment. Additionally, exogenous application of hydrogen peroxide efficiently induced anthocyanin biosynthesis within 24 h, even under ML conditions, suggesting that ROS played a major role in activating anthocyanin biosynthesis. A HL-responsive MYB transcription factor EgrMYB113 was identified to play an important role in regulating anthocyanin biosynthesis by targeting multiple anthocyanin biosynthesis genes. Additionally, the results demonstrated that gibberellic acid and sugar signaling contributed to HL-induced anthocyanin biosynthesis. Conclusively, these results suggested that HL triggers multiple signaling pathways to induce anthocyanin biosynthesis, with ROS acting as indispensable mediators in Eucalyptus.

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.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Genomic survey and expression analysis of LcARFs reveal multiple functions to somatic embryogenesis in Liriodendron

BACKGROUND

Auxin response factors (ARFs) are critical transcription factors that mediate the auxin signaling pathway and are essential for regulating plant growth. However, there is a lack of understanding regarding the ARF gene family in Liriodendron chinense, a vital species in landscaping and economics. Thus, further research is needed to explore the roles of ARFs in L. chinense and their potential applications in plant development.

RESULT

In this study, we have identified 20 LcARF genes that belong to three subfamilies in the genome of L. chinense. The analysis of their conserved domains, gene structure, and phylogeny suggests that LcARFs may be evolutionarily conserved and functionally similar to other plant ARFs. The expression of LcARFs varies in different tissues. Additionally, they are also involved in different developmental stages of somatic embryogenesis. Overexpression of LcARF1, LcARF2a, and LcARF5 led to increased activity within callus. Additionally, our promoter-GFP fusion study indicated that LcARF1 may play a role in embryogenesis. Overall, this study provides insights into the functions of LcARFs in plant development and embryogenesis, which could facilitate the improvement of somatic embryogenesis in L. chinense.

CONCLUSION

The research findings presented in this study shed light on the regulatory roles of LcARFs in somatic embryogenesis in L. chinense and may aid in accelerating the breeding process of this tree species. By identifying the specific LcARFs involved in different stages of somatic embryogenesis, this study provides a basis for developing targeted breeding strategies aimed at optimizing somatic embryogenesis in L. chinense, which holds great potential for improving the growth and productivity of this economically important species.

Cytokinin promotes anthocyanin biosynthesis via regulating sugar accumulation and MYB113 expression in Eucalyptus

Anthocyanins are flavonoid-like substances that play important roles in plants' adaptation to various environmental stresses. In this research, we discovered that cytokinin (CK) alone could effectively induce the anthocyanin biosynthesis in Eucalyptus and many other perennial woody plant species, but not in tobacco and Arabidopsis, suggesting a diverse role of CK in regulating anthocyanin biosynthesis in different species. Transcriptomic and metabolomic strategies were used to further clarify the specific role of CK in regulating anthocyanin biosynthesis in Eucalyptus. The results showed that 801 and 2241 genes were differentially regulated at 6 and 24 h, respectively, after CK treatment. Pathway analysis showed that most of the differentially expressed genes were categorized into pathways related to cellular metabolism or transport of metabolites, including amino acids and sugars. The metabolomic results well supported the transcriptome data, which showed that most of the differentially regulated metabolites were related to the metabolism of sugar, amino acids and flavonoids. Moreover, CK treatment significantly induced the accumulation of sucrose in the CK-treated leaves, while sugar starvation mimicked by either defoliation or shading treatment of the basal leaves significantly reduced the sugar increase of the CK-treated leaves and thus inhibited CK-induced anthocyanin biosynthesis. The results of in vitro experiment also suggested that CK-induced anthocyanin in Eucalyptus was sugar-dependent. Furthermore, we identified an early CK-responsive transcription factor MYB113 in Eucalyptus, the expression of which was significantly upregulated by CK treatment in Eucalyptus, but was inhibited in Arabidopsis. Importantly, the overexpression of EgrMYB113 in the Eucalyptus hairy roots was associated with significant anthocyanin accumulation and upregulation of most of the anthocyanin biosynthetic genes. In conclusion, our study demonstrates a key role of CK in the regulation of anthocyanin biosynthesis in Eucalyptus, providing a molecular basis for further understanding the regulatory mechanism and diversity of hormone-regulated anthocyanin biosynthesis in different plant species.

Integrated Metabolomic and Transcriptomic Analyses Reveal the Potential Regulation of Flavonoids in the Production of Embryogenic Cultures during Early Somatic Embryogenesis of Longan (Dimocarpus longan Lour.)

Embryogenic cultures of longan (Dimocarpus longan Lour.) contain various metabolites with pharmacological properties that may function in the regulation of somatic embryogenesis (SE). In this study, based on widely targeted metabolomics, 501 metabolites were obtained from the embryogenic calli, incomplete compact proembryogenic cultures, and globular embryos during early SE of longan, among which 41 flavonoids were differentially accumulated during the SE. Using RNA sequencing, 36 flavonoid-biosynthesis-related genes and 43 MYB and 52 bHLH transcription factors were identified as differentially expressed genes. Furthermore, Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the flavonoid metabolism-related pathways were significantly enriched during the early SE. These results suggested that the changes in flavonoid levels in the embryogenic cultures of longan were mediated by MYBs and bHLHs via regulating flavonoid-biosynthesis-related genes, thus potentially regulating early SE. The identified metabolites in the embryogenic cultures of longan can be used to develop pharmaceutical ingredients.

Molecular mechanism of somatic embryogenesis in paeonia ostii 'Fengdan' based on transcriptome analysis combined histomorphological observation and metabolite determination

BACKGROUND

Tree peony (Paeonia sect. Moutan DC.) is a famous flower native to China with high ornamental, medicinal, and oil value. However, the low regeneration rate of callus is one of the main constraints for the establishment of a genetic transformation system in tree peony. By histomorphological observation, transcriptomic analysis and metabolite determination, we investigated the molecular mechanism of somatic embryogenesis after the establishment of a culture system and the induction of somatic embryo(SE) formation.

RESULTS

We found that SE formation was successfully induced when cotyledons were used as explants. A total of 3185 differentially expressed genes were screened by comparative transcriptomic analysis of embryogenic callus (EC), SE, and non-embryogenic callus (NEC). Compared to NEC, the auxin synthesis-related genes GH3.6 and PCO2 were up-regulated, whereas cytokinin dehydrogenase (CKX6) and CYP450 family genes were down-regulated in somatic embryogenesis. In SE, the auxin content was significantly higher than the cytokinin content. The methyltransferase-related gene S-adenosylmethionine synthase (SAMS) and the flavonoid biosynthesis-related gene (ANS and F3'5'H) were down-regulated in somatic embryogenesis. The determination of flavonoids showed that rhoifolin and hyperoside had the highest content in SE. The results of transcriptome analysis were consistent with the relative expression of 8 candidate genes by quantitative polymerase chain reaction analysis.

CONCLUSION

The results revealed that auxin and cytokinin may play a key role in 'Fengdan' somatic embryogenesis. The genes related to somatic embryogenesis were revealed, which has partly elucidated the molecular mechanism of somatic embryogenesis in 'Fengdan'.

Changes and transcriptome regulation of endogenous hormones during somatic embryogenesis in Ormosia henryi Prain

Introduction

Ormosia henryi is a rare and endangered plant growing in southern China. Somatic embryo culture is an effective measure for the rapid propagation of O. henryi. It has not been reported how regulatory genes induce somatic embryogenesis by regulating endogenous hormone changes during the process of somatic embryogenesis in O. henryi.

Methods

In this study, we analysed the endogenous hormone levels and transcriptome data of nonembryogenic callus (NEC), embryogenic callus (EC), globular embryo (GE) and cotyledon embryo (CE) in O. henryi.

Results

The results showed that the indole-3-acetic acid (IAA) content was higher and the cytokinins (CKs) content was lower in EC than in NEC, and the gibberellins (GAs) and abscisic acid (ABA) contents were significantly higher in NEC than in EC. The contents of IAA, CKs, GAs and ABA increased significantly with EC development. The expression patterns of differentially expressed genes (DEGs) involved in the biosynthesis and signal transduction of auxin (AUX) (YUCCA and SAUR), CKs (B-ARR), GAs (GA3ox, GA20ox, GID1 and DELLA) and ABA (ZEP, ABA2, AAO3, CYP97A3, PYL and ABF) were consistent with the levels of endogenous hormones during somatic embryogenesis (SE). In this study, 316 different transcription factors (TFs) regulating phytohormones were detected during SE. AUX/IAA were downregulated in the process of EC formation and GE differentiation into CE, but other TFs were upregulated and downregulated.

Conclusion

Therefore, we believe that relatively high IAA content and low CKs, GAs and ABA contents contribute to EC formation. The differential expression of AUX, CKs, GAs and ABA biosynthesis and signal transduction genes affected the endogenous hormone levels at different stages of SE in O. henryi. The downregulated expression of AUX/IAA inhibited NEC induction, promoted EC formation and GE differentiation into CE.

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.

Comparative transcriptomics and bioinformatics analysis of genes related to photosynthesis in Eucalyptus camaldulensis

The timber species Eucalyptus camaldulensis is one of the most important in southern China. Therefore, it is essential to understand the photosynthetic pattern in eucalyptus leaves. In the present study, eighteen photosynthesis-related genes were analyzed using bioinformatics methods. The results indicated that there were ten differentially expressed ribose-5-phosphate isomerase genes (RPI), and six of them were up-regulated in the mature leaves compared to the young leaves, while others were down-regulated. The differential expression of four rubisco methyltransferase genes (RBCMT) were observed. Two of them were up-regulated, while two were down-regulated in mature leaves compared to young leaves. Furthermore, two ribulose-phosphate-3-epimerase genes (RPE) were up-regulated in the mature leaves compared to the young leaves. In contrast, two genes involved in triosephosphate isomerase (TIM) were down-regulated in mature leaves compared with young leaves. The current study provides basic information about the transcriptome of E. camaldulensis and lays a foundation for further research in developing and utilizing important photosynthetic genes.

BREVIPEDICELLUS and ERECTA control the expression of AtPRX17 to prevent Arabidopsis callus browning

Efficient in vitro callus generation is required for tissue culture propagation, a process that allows for plant regeneration and transgenic breeding for desired phenotypes. Identifying genes and regulatory elements that prevent impaired callus growth and callus browning is essential for the development of in vitro callus systems. Here, we show that the BREVIPEDICELLUS and ERECTA pathways in Arabidopsis calli converge to prevent callus browning, and positively regulate the expression of the isoperoxidase gene AtPRX17 in rapidly growing calli. Loss-of-function mutations in both BREVIPEDICELLUS and ERECTA resulted in markedly increased callus browning. Transgenic lines expressing 35S pro::AtPRX17 in the bp-5 er105 double mutant background fully rescued this phenotypic abnormality. Using in vivo (chromatin immunoprecipitation-PCR and transient expression) and in vitro (electrophoretic mobility shift assays) assays, we observed that the BREVIPEDICELLUS protein binds directly to the upstream sequence of AtPRX17 to promote its transcription during callus growth. ERECTA is a ubiquitous factor required for cell proliferation and growth. We show that ERECTA positively regulates the expression of the transcription factor WRKY6, which directly binds to a separate site on the AtPRX17 promoter, further increasing its expression. Our data reveal an important molecular mechanism involved in the regulation of peroxidase isozyme expression to reduce Arabidopsis callus browning.

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.

Global Transcriptome and Coexpression Network Analyses Reveal New Insights Into Somatic Embryogenesis in Hybrid Sweetgum (Liquidambar styraciflua × Liquidambar formosana)

Somatic embryogenesis (SE) is a process of somatic cells that dedifferentiate to totipotent embryonic stem cells and generate embryos in vitro. Despite recent scientific headway in deciphering the difficulties of somatic embryogenesis, the overall picture of key genes, pathways, and co-expression networks regulating SE is still fragmented. Therefore, deciphering the molecular basis of somatic embryogenesis of hybrid sweetgum remains pertinent. In the present study, we analyzed the transcriptome profiles and gene expression regulation changes via RNA sequencing from three distinct developmental stages of hybrid sweetgum: non-embryogenic callus (NEC), embryogenic callus (EC), and redifferentiation. Comparative transcriptome analysis showed that 19,957 genes were differentially expressed in ten pairwise comparisons of SE. Among these, plant hormone signaling-related genes, especially the auxin and cytokinin signaling components, were significantly enriched in NEC and EC early. The K-means method was used to identify multiple transcription factors, including HB-WOX, B3-ARF, AP2/ERF, and GRFs (growth regulating factors). These transcription factors showed distinct stage- or tissue-specific expression patterns mirroring each of the 12 superclusters to which they belonged. For example, the WOX transcription factor family was expressed only at NEC and EC stages, ARF transcription factor was expressed in EC early, and GRFs was expressed in late SE. It was noteworthy that the AP2/ERF transcription factor family was expressed during the whole SE process, but almost not in roots, stems and leaves. A weighted gene co-expression network analysis (WGCNA) was used in conjunction with the gene expression profiles to recognize the genes and modules that may associate with specific tissues and stages. We constructed co-expression networks and revealed 22 gene modules. Four of these modules with properties relating to embryonic potential, early somatic embryogenesis, and somatic embryo development, as well as some hub genes, were identified for further functional studied. Through a combination analysis of WGCNA and K-means, SE-related genes including AUX22, ABI3, ARF3, ARF5, AIL1, AIL5, AGL15, WOX11, WOX9, IAA29, BBM1, MYB36, LEA6, SMR4 and others were obtained, indicating that these genes play an important role in the processes underlying the progression from EC to somatic embryos (SEs) morphogenesis. The transcriptome information provided here will form the foundation for future research on genetic transformation and epigenetic control of plant embryogenesis at a molecular level. In follow-up studies, these data could be used to construct a regulatory network for SE; Key genes obtained from coexpression network analysis at each critical stage of somatic embryo can be considered as potential candidate genes to verify these networks.

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.