In this study, 28,691 genome sequences and 16,566 expressed sequence tags (ESTs) of Eucalyptus were derived from the GenBank database. A total of 2292 SSR loci were sought out from 1785 effective sequences. Through analyses of SSR loci information, the SSR motif length was negatively correlated with the abundance of the SSRs. In the EST sequences of Eucalyptus, triplet repeat motifs were the most abundant, and dinucleotide repeats motifs had the highest frequencies. Subsequently, 395 pairs of primers were designed based on the SSR loci. Using optimized SSR-PCR conditions, 340 pairs of primers were successfully screened, with a success rate of 86.1%. By construction of a maximum likelihood phylogenetic tree of six eucalypt species, represented by five species of the genus Eucalyptus and one of the genus Corymbia, the genetic relationships of Eucalyptus urophylla and E. camaldulensis suggested by this tree was found to differ from that suggested by traditional morphological taxonomy. The results provide insights for evaluating genetic diversity of Eucalyptus and analysis of Eucalyptus phylogenetics using SSR markers.

Acanthus ilicifolius is an important medicinal plant in mangrove forests, which is rich in secondary metabolites with various biological activities. In this study, we used transcriptomic analysis to obtain differentially expressed genes in the flavonoid metabolic pathway and metabolomic methods to detect changes in the types and content in the flavonoid metabolic synthesis pathway. The results showed that DEGs were identified in the mature roots vs. leaves comparison (9001 up-regulated and 8910 down-regulated), mature roots vs. stems comparison (5861 up-regulated and 7374 down-regulated), and mature stems vs. leaves comparison (10,837 up-regulated and 11,903 down-regulated). Furthermore, two AiCHS genes and four AiCHI genes were up-regulated in the mature roots vs. stems of mature A. ilicifolius, and were down-regulated in mature stems vs. leaves, which were highly expressed in the A. ilicifolius stems. A total of 215 differential metabolites were found in the roots vs. leaves of mature A. ilicifolius, 173 differential metabolites in the roots vs. stems, and 228 differential metabolites in the stems vs. leaves. The metabolomic results showed that some flavonoids in A. ilicifolius stems were higher than in the roots. A total of 18 flavonoid differential metabolites were detected in the roots, stems, and leaves of mature A. ilicifolius. In mature leaves, quercetin-3-O-glucoside-7-O-rhamnoside, gossypitrin, isoquercitrin, quercetin 3,7-bis-O-β-D-glucoside, and isorhamnetin 3-O-β-(2″-O-acetyl-β-D-glucuronide) were found in a high content, while in mature roots, di-O-methylquercetin and isorhamnetin were the major compounds. The combined analysis of the metabolome and transcriptome revealed that DEGs and differential metabolites were related to flavonoid biosynthesis. This study provides a theoretical basis for analyzing the molecular mechanism of flavonoid synthesis in A. ilicifolius and provides a reference for further research and exploitation of its medicinal value.

Cellulose is the world’s most abundant renewable energy resource, and a variety of cellulose synthase genes are involved in the biosynthesis of cellulose. In the process of cellulose synthesis, all cellulose synthases are interrelated and act synergistically. In this study, we analyzed the contents of cellulose, hemicellulose, and lignin in the different parts and tissues of E. grandis. The results showed that the cellulose content had greater differences among three different heights. On this basis, we carried out the transcriptome-wide profiling of gene expression patterns using RNA sequencing. A total of 2066 differentially expressed genes were identified for three pairwise comparisons between three different heights, most of which were related to the programmed photosynthetic membrane and photosystem. A total of 100 transcripts of CSs (58 CesA and 42 Csl) were obtained from transcriptome libraries. The expression pattern of these genes indicated that different CS genes had a wide range of expression profiles. A phylogenetic analysis of 135 reference CS genes showed that the CSs of E. grandis were clustered into six major groups (CesA1-9, CslA, CslB/H, CslD, CslE, and CslG). Based on the weighted gene co-expression network analysis, a dual-directional regulation mechanism between Csl and CesA proteins in the cellulose biosynthesis was identified. The gene expression profile analysis, using qRT-PCR in different tissues of E. grandis, demonstrated that the CSs were highly expressed in xylem, and CesAs had a higher relative expression than Csls. The analysis of sequence similarity combined with the expression pattern indicated that the CesA1, 3, and 6 transcripts were associated with the biosynthesis of the secondary cell wall, and CesA4, 5, and 7 transcripts were more likely to associate with the biosynthesis of the primary cell wall. Finally, the qRT-PCR analysis confirmed the expression of 11 selected CSs in three different parts of E. grandis.

Cellulose, an essential structural component in the plant cell wall and a renewable biomass resource, plays a significant role in nature. Eucalyptus’s excellent timber tree species (including Eucalyptus grandis Hill) provide many raw materials for the paper and wood industries. The synthesis of cellulose is a very complex process involving multiple genes and regulated by various biological networks. However, research on regulating associated genes and non-coding RNAs during cellulose synthesis in E. grandis remains lacking. In this study, the wood anatomical characteristics and chemical indexes of E. grandis were analyzed by taking three different parts (diameter at breast height (DBH), middle and upper part of the trunk) from the main stem of E. grandis as raw materials. The role of non-coding RNAs (Long non-coding RNA, lncRNA; Micro RNA, miRNA; Circle RNA, circRNA) on regulating candidate genes was presented, and the network map of ceRNA (Competing endogenous RNA) regulation during wood cellulose biosynthesis of E. grandis was constructed. The transcriptome sequencing of nine samples obtained from the trunk of the immature xylem in E. grandis at DBH, middle and upper parts had a 95.81 G clean reading, 57,480 transcripts, 7365 lncRNAs, and 5180 circRNAs. Each sample had 172–306 known miRNAs and 1644–3508 new miRNAs. A total of 190 DE-lncRNAs (Differentially expressed long non-coding RNAs), 174 DE-miRNAs (Differentially expressed micro RNAs), and 270 DE-circRNAs (Differentially expressed circle RNAs) were obtained by comparing transcript expression levels. Four lncRNAs and nine miRNAs were screened out, and the ceRNA regulatory network was constructed. LncRNA1 and lncRNA4 regulated the genes responsible for cellulose synthesis in E. grandis, which were overexpressed in 84K (Populus Alba × Populus glandulosa) poplar. The cellulose and lignin content in lncRNA4-oe were significantly higher than wild type 84K poplar and lncRNA1-oe. The average plant height, middle and basal part of the stem diameter in lncRNA4-oe were significantly higher than the wild type. However, there was no significant difference between the growth of lncRNA1-oe and the wild type. Further studies are warranted to explore the molecular regulatory mechanism of cellulose biosynthesis in Eucalyptus species.

The WUSCHEL-related homeobox (WOX) gene family are key players in the rooting process. Eucalyptus is an important plant species of artificial forests in China. It is mainly grown through tissue culture of many excellent clonal materials, in which rooting is a key step. In the present study, by using the genome data of Eucalyptus grandis, Corymbia citriodora, E. pellita, and E. urophylla × E. grandis, the members of the eucalyptus WOX gene family were identified and analyzed by bioinformatics techniques. The eucalyptus WOX gene family members are unstable proteins, with 7 acidic proteins and 24 basic proteins, and no signal peptide region was detected. Subcellular localization prediction indicated that all these proteins are localized in the nucleus. Motif analysis showed that eucalyptus WOX genes share the same motifs. Phylogenetic tree and gene expression analyses revealed that the eucalyptus WOX genes are highly conserved during the evolution process. Moreover, the WOX protein sequences are also highly conserved within the species, with higher similarity between woody plants. The EupWOX gene showed tissue-specific expression, with EupWOX1 and EupWOX11 specifically expressed in the roots of E. urophylla × E. pellita clonal tissue culture during the late-stage rooting. This finding suggests that EupWOX1 may be a key regulatory gene induced by the root primordium and is critically related to the rooting rate. EupWOX1, EupWOX5, and EupWOX13 could be the key regulatory genes for adventitious root formation. EupWOX1, EupWOX5, and EupWOX13 could be the key regulatory genes for the elongation of adventitious roots and the growth of adventitious lateral roots. EupWOX5 and EupWOX13 could play a critical role, not only in the formation of adventitious roots and adventitious lateral roots of E. urophylla clonal tissue culture but also in the root elongation process. These results will help us understand the complexity of rooting in different lines and provide valuable information for future functional characterization of specific genes in eucalyptus clones.

Eucalyptus camaldulensis Dehn. is one of the most morphologically and genetically variable Eucalyptus species. Growth, Leptocybe invasa Fisher & La Salle susceptibility, pilodyn penetration and other traits up to age 36 months were assessed in a seed source/family trial in China comprising 112 seedlots representing five natural stand and six exotic seed sources. Genetic diversity and population structure of this trial population were also analyzed using 48 simple sequence repeat (SSR) markers. The key objective was to examine whether the genomic data could provide value over information obtained from just quantitative trait data. Significant genetic variation was found among seed sources and among families within seed sources for most quantitative traits. The ratio of variance among seed sources to variance among families within seed sources, based on variances estimated from quantitative trait data, varied from 0.1% (height at 9 months) up to 75.2% (bark thickness). Equivalent ratios estimated from the AMOVA on SSR loci data were similar for height (ages 24 and 36 months) and also pilodyn penetration at 36 months, but not for 9-month height or 36-month bark thickness. From 48 SSR loci examined, the genetic differentiation coefficient (among seed sources) was 0.086, indicating low genetic differentiation among seed sources. While overall genetic diversity in the trial population examined was high, the levels within the different seed sources varied markedly. Prior to this study, genetic distances among families from the three exotic seed sources (from domesticated Indian populations) in the trial, along with their genetic distances from, and relatedness to, families from five natural stand seed sources (Australian) in the trial were unknown. The SSR loci data removed uncertainties and revealed that the exotic sources increased the breadth of genetic origins represented in the trial population—information that could not have been obtained from just the quantitative trait data.