Genome-wide identification and expression analysis of GRAS gene family in Eucalyptus grandis.

BACKGROUND: The GRAS gene family is a class of plant-specific transcription factors with important roles in many biological processes, such as signal transduction, disease resistance and stress tolerance, plant growth and development. So far, no information available describes the functions of the GRAS genes in Eucalyptus grandis. RESULTS: A total of 82 GRAS genes were identified with amino acid lengths ranging from 267 to 817 aa, and most EgrGRAS genes had one exon. Members of the GRAS gene family of Eucalyptus grandis are divided into 9 subfamilies with different protein structures, while members of the same subfamily have similar gene structures and conserved motifs. Moreover, these EgrGRAS genes expanded primarily due to segmental duplication. In addition, cis-acting element analysis showed that this family of genes was involved involved in the signal transduction of various plant hormones, growth and development, and stress response. The qRT-PCR data indicated that 18 EgrGRAS genes significantly responded to hormonal and abiotic stresses. Among them, the expression of EgrGRAS13, EgrGRAS68 and EgrGRAS55 genes was significantly up-regulated during the treatment period, and it was hypothesised that members of the EgrGRAS family play an important role in stress tolerance. CONCLUSIONS: In this study, the phylogenetic relationship, conserved domains, cis-elements and expression patterns of GRAS gene family of Eucalyptus grandis were analyzed, which filled the gap in the identification of GRAS gene family of Eucalyptus grandis and laid the foundation for analyzing the function of EgrGRAS gene in hormone and stress response.

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.

Xylem cell size regulation is a key adaptive response to water deficit in Eucalyptus grandis.

Future climatic scenarios forecast increasingly frequent droughts that will pose substantial consequences on tree mortality. In light of this, drought-tolerant eucalypts have been propagated; however, the severity of these conditions will invoke adaptive responses, impacting the commercially valuable wood properties. To determine what mechanisms govern the wood anatomical adaptive response, highly controlled drought experiments were conducted in Eucalyptus grandis W. Hill ex Maiden, with the tree physiology and transcriptome closely monitored. In response to water deficit, E. grandis displays an isohydric stomatal response to conserve water and enable stem growth to continue, albeit at a reduced rate. Maintaining gaseous exchange is likely a critical short-term response that drives the formation of hydraulically safer xylem. For instance, the development of significantly smaller fibers and vessels was found to increase cellular density, thereby promoting drought tolerance through improved functional redundancy, as well as implosion and cavitation resistance. The transcriptome was explored to identify the molecular mechanisms responsible for controlling xylem cell size during prolonged water deficit. Downregulation of genes associated with cell wall remodeling and the biosynthesis of cellulose, hemicellulose and pectin appeared to coincide with a reduction in cellular enlargement during drought. Furthermore, transcript levels of NAC and MYB transcription factors, vital for cell wall component biosynthesis, were reduced, while those linked to lignification increased. The upregulation of EgCAD and various peroxidases under water deficit did not correlate with an increased lignin composition. However, with the elevated cellular density, a higher lignin content per xylem cross-sectional area was observed, potentially enhancing hydraulic safety. These results support the requirement for higher density, drought-adapted wood as a long-term adaptive response in E. grandis, which is largely influenced by the isohydric stomatal response coupled with cellular expansion-related molecular processes.

Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus.

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Insights into Leptocybe invasa resistance in Eucalyptus: phenotyping, genotyping and in silico approaches.

The gall wasp, Leptocybe invasa, poses a significant global threat to Eucalyptus cultivation, by causing substantial economic losses. The objective of this study was to differentiate between resistant and susceptible genotypes by morphological characteristics using image analysis based on the damage caused by the gall wasp. In addition, consensus sequences derived from transposable elements (TEs) and the genome of Eucalyptus spp. Were identified by in silico analysis. Furthermore, another objective was to discriminate Eucalyptus genotypes in response to Leptocybe invasa by conducting molecular analyses involving transposable elements and inter simple sequence markers. For image analysis, the GroundEye ® system was used to collect images of 60 leaves from six genotypes, three of which were resistant and three susceptible. Eucalyptus spp. sequences were obtained from the GenBank database by in silico analysis and pairwise alignments with TE sequences were conducted using BLASTN. Multiple sequence alignment was performed with Clustal Omega, followed by the identification of conserved regions in Jalview. A motif signature was generated using Weblogo. For molecular characterization using ISSR markers and TEs, samples of young leaves were obtained from a total of 80 Eucalyptus seedlings, of which 50 were classified as resistant and 30 as susceptible to L. invasa. It was possible to distinguish gall wasp susceptible and resistant genotypes by image analysis. In silico analysis enabled the identification of conserved regions in the Eucalyptus spp. genome, which were associated with proteins involved in secondary metabolite production, e.g., terpenes, which play a role in the response to L. invasa. The discrimination capacity of TEs and ISSR primers was demonstrated and bands were generated that could be used to identify resistant genotypes. However, increasing the number of markers required to discriminate genotypes in both cases is suggested.

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.

Evolution of rarity and phylogeny determine above- and belowground biomass in plant-plant interactions.

Rare species are often considered inferior competitors due to occupancy of small ranges, specific habitats, and small local populations. However, the phylogenetic relatedness and rarity level (level 1-7 and common) of interacting species in plant-plant interactions are not often considered when predicting the response of rare plants in a biotic context. We used a common garden of 25 species of Tasmanian Eucalyptus, to differentiate non-additive patterns in the biomass of rare versus common species when grown in mixtures varying in phylogenetic relatedness and rarity. We demonstrate that rare species maintain progressively positive non-additive responses in biomass when interacting with phylogenetically intermediate, less rare and common species. This trend is not reflected in common species that out-performed in monocultures compared to mixtures. These results offer predictability as to how rare species' productivity will respond within various plant-plant interactions. However, species-specific interactions, such as those involving E. globulus, yielded a 97% increase in biomass compared to other species-specific interaction outcomes. These results are important because they suggest that plant rarity may also be shaped by biotic interactions, in addition to the known environmental and population factors normally used to describe rarity. Rare species may utilize potentially facilitative interactions with phylogenetically intermediate and common species to escape the effects of limiting similarity. Biotically mediated increases in rare plant biomass may have subsequent effects on the competitive ability and geographic occurrence of rare species, allowing rare species to persist at low abundance across plant communities. Through the consideration of species rarity and evolutionary history, we can more accurately predict plant-plant interaction dynamics to preserve unique ecosystem functions and fundamentally challenge what it means to be "rare".

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence.

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis.

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

High-quality genome resource of Lasiodiplodia pseudotheobromae associated with die-back on Eucalyptus trees.

OBJECTIVES: Lasiodiplodia pseudotheobromae is an important fungal pathogen associated with die-back, canker and shoot blight in many plant hosts with a wide geographic distribution. The aim of our study was to provide high-quality genome assemblies and sequence annotation resources of L. pseudotheobromae, to facilitate future studies on the systematics, population genetics and genomics of the fungal pathogen L. pseudotheobromae. DATA DESCRIPTION: High-quality genomes of five L. pseudotheobromae isolates were sequenced based on Oxford Nanopore technology (ONT) and Illumina HiSeq sequencing platform. The total size of each assembly ranged from 43 Mb to 43.86 Mb and over 11,000 protein-coding genes were predicted from each genome. The proteins of predicted genes were annotated using multiple public databases, among the annotated protein-coding genes, more than 4,300 genes were predicted as potential virulence genes by the Pathogen Host Interactions (PHI) database. Moreover, the genome comparative analysis among L. pseudotheobromae and other closely related species revealed that 7,408 gene clusters were shared among them and 152 gene clusters unique to L. pseudotheobromae. This genome and associated datasets provided here will serve as a useful resource for further analyses of this fungal pathogen species.

Reference genes for Eucalyptus spp. under Beauveria bassiana inoculation and subsequently infestation by the galling wasp Leptocybe invasa.

Relative gene expression analysis through RT-qPCR is an important molecular technique that helps understanding different molecular mechanisms, such as the plant defense response to insect pests. However, the use of RT-qPCR for gene expression analysis can be affected by factors that directly affect the reliability of the results. Among these factors, the appropriate choice of reference genes is crucial and can strongly impact RT-qPCR relative gene expression analyses, highlighting the importance in correctly choosing the most suitable genes for the success of the analysis. Thus, this study aimed to select and validate reference genes for relative gene expression studies through RT-qPCR in hybrids of Eucalyptus tereticornis × Eucalyptus camaldulensis (drought tolerant and susceptible to Leptocybe invasa) under conditions of inoculation by the Beauveria bassiana fungus and subsequent infestation by L. invasa. The expression level and stability of eleven candidate genes were evaluated. Stability was analyzed using the RefFinder tool, which integrates the geNorm, NormFinder, BestKeeper, and Delta-Ct algorithms. The selected reference genes were validated through the expression analysis of the transcriptional factor EcDREB2 (dehydration-responsive element-binding protein 2). For all treatments evaluated, EcPTB, EcPP2A-1, and EcEUC12 were the best reference genes. The triplets EcPTB/EcEUC12/EcUBP6, EcPP2A-1/EcEUC12/EcPTB, EcIDH/EcSAND/Ecα-TUB, EcPP2A-1/Ecα-TUB/EcPTB, and EcPP2A-1/EcUPL7/EcSAND were the best reference genes for the control plants, mother plants, plants inoculated with B. bassiana, plants infested with L. invasa, and plants inoculated with B. bassiana and subsequently infested with L. invasa, respectively. The best determined reference genes were used to normalize the RT-qPCR expression data for each experimental condition evaluated. The results emphasize the importance of this type of study to ensure the reliability of relative gene expression analyses. Furthermore, the findings of this study can be used as a basis for future research, comprising gene expression analysis of different eucalyptus metabolic pathways.

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.

Identification of WUSCHEL-related homeobox gene and truncated small peptides in transformation efficiency improvement in Eucalyptus.

BACKGROUND: The WUSCHEL-related Homeobox (WOX) genes, which encode plant-specific homeobox (HB) transcription factors, play crucial roles in regulating plant growth and development. However, the functions of WOX genes are little known in Eucalyptus, one of the fastest-growing tree resources with considerable widespread cultivation worldwide. RESULTS: A total of nine WOX genes named EgWOX1-EgWOX9 were retrieved and designated from Eucalyptus grandis. From the three divided clades marked as Modern/WUS, Intermediate and Ancient, the largest group Modern/WUS (6 EgWOXs) contains a specific domain with 8 amino acids: TLQLFPLR. The collinearity, cis-regulatory elements, protein-protein interaction network and gene expression analysis reveal that the WUS proteins in E. grandis involve in regulating meristems development and regeneration. Furthermore, by externally adding of truncated peptides isolated from WUS specific domain, the transformation efficiency in E. urophylla × E. grandis DH32-29 was significant enhanced. The transcriptomics data further reveals that the use of small peptides activates metabolism pathways such as starch and sucrose metabolism, phenylpropanoid biosynthesis and flavonoid biosynthesis. CONCLUSIONS: Peptides isolated from WUS protein can be utilized to enhance the transformation efficiency in Eucalyptus, thereby contributing to the high-efficiency breeding of Eucalyptus.

Genes expression profiles in vascular cambium of Eucalyptus urophylla × Eucalyptus grandis at different ages.

BACKGROUND: Wood is a secondary xylem generated by vascular cambium. Vascular cambium activities mainly include cambium proliferation and vascular tissue formation through secondary growth, thereby producing new secondary phloem inward and secondary xylem outward and leading to continuous tree thickening and wood formation. Wood formation is a complex biological process, which is strictly regulated by multiple genes. Therefore, molecular level research on the vascular cambium of different tree ages can lead to the identification of both key and related genes involved in wood formation and further explain the molecular regulation mechanism of wood formation. RESULTS: In the present study, RNA-Seq and Pac-Bio Iso-Seq were used for profiling gene expression changes in Eucalyptus urophylla × Eucalyptus grandis (E. urograndis) vascular cambium at four different ages. A total of 59,770 non-redundant transcripts and 1892 differentially expressed genes (DEGs) were identified. The expression trends of the DEGs related to cell division and differentiation, cell wall biosynthesis, phytohormone, and transcription factors were analyzed. The DEGs encoding expansin, kinesin, cycline, PAL, GRP9, KNOX, C2C2-dof, REV, etc., were highly expressed in E. urograndis at three years old, leading to positive effects on growth and development. Moreover, some gene family members, such as NAC, MYB, HD-ZIP III, RPK, and RAP, play different regulatory roles in wood formation because of their sophisticated transcriptional network and function redundantly. CONCLUSIONS: These candidate genes are a potential resource to further study wood formation, especially in fast-growing and adaptable eucalyptus. The results may also serve as a basis for further research to unravel the molecular mechanism underlying wood formation.

Safety Assessment of the CP4 EPSPS and NPTII Proteins in Eucalyptus.

Glyphosate herbicide treatment is essential to sustainable Eucalyptus plantation management in Brazil. Eucalyptus is highly sensitive to glyphosate, and Suzano/FuturaGene has genetically modified eucalyptus to tolerate glyphosate, with the aim of both protecting eucalyptus trees from glyphosate application damage and improving weed management. This study presents the biosafety results of the glyphosate-tolerant eucalyptus event 751K032, which expresses the selection marker neomycin phosphotransferase II (NPTII) enzyme and CP4-EPSPS, a glyphosate-tolerant variant of plant 5-enolpyruvyl-shikimate-3-phosphate synthase enzyme. The transgenic genetically modified (GM) event 751K032 behaved in the plantations like conventional non-transgenic eucalyptus clone, FGN-K, and had no effects on arthropods and soil microorganisms. The engineered NPTII and CP4 EPSPS proteins were heat-labile, readily digestible, and according to the bioinformatics analyses, unlikely to cause an allergenic or toxic reaction in humans or animals. This assessment of the biosafety of the glyphosate-tolerant eucalyptus event 751K032 concludes that it is safe to be used for wood production.

Functional responses to multiple sequential abiotic stress (waterlogging-drought) in three woody taxa with different root systems and stress tolerance.

There is generally a trade-off in the resistance to drought and to waterlogging. However, several species are sequentially subjected to both stressors in many environments. We evaluated the ecophysiological strategies to cope with multiple sequential stress of waterlogging and drought (W + D) of three taxa differing in stress resistance and root morphology: the phreatophic Eucalyptus camaldulensis (Ec) and two shallow-rooted willow clones: Salix matsudana x Salix alba (SmxSa) and Salix nigra (Sn4). Individuals of the three taxa were grown in pots and assigned to either of four treatments: Control (well-watered plants), well-watered followed by drought (C + D); waterlogged for 15 days followed by drought (W15d + D) and waterlogged for 30 days followed by drought (W30d + D). Biomass allocation, growth (diameter, height, length of leaves, and roots), specific leaf area, stomatal conductance, water potential, hydraulic conductivity of roots and branches, leaf C13 and root cortical aerenchyma formation were determined at different stages of the experiment. Ec growth was not affected by W + D, developing tolerance strategies at leaf and whole plant levels. Differential effects of W + D were observed in both Salix clones depending on the time of waterlogging. In Sn4 and SmxSa, the root biomass was affected in W15d + D treatment, but a root tolerance response (aerenchyma and adventitious root formation) was observed in W30d + D. In the three taxa, and contrary to expectations, the previous exposure to a waterlogging period did not increase the susceptibility of the plants to a subsequent drought event. On the contrary, we found tolerance, which depended on the time of waterlogging exposure.

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses.

Metal tolerance proteins (MTPs) as Me2+/H+(K+) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd2+ and Cu2+ exposure might be responsible for metal translocation from roots to leaves.

Survivorship and food consumption of immatures and adults of Apis mellifera and Scaptotrigona bipunctata exposed to genetically modified eucalyptus pollen.

Eucalyptus comprises the largest planted area of cultivated production forest in Brazil. Genetic modification (GM) of eucalyptus can provide additional characteristics for increasing productivity and protecting wood yield, as well as potentially altering fiber for a diversity of industrial uses. However, prior to releasing a new GM plant, risk assessments studies with non-target organisms must be undertaken. Bees are prominent biological models since they play an important role in varied ecosystems, including for Eucalyptus pollination. The main goal of this study was to evaluate whether a novel event (Eucalyptus 751K032), which carries the cp4-epsps gene that encodes the protein CP4-EPSPS and nptII gene that encodes the protein NPTII, might adversely affect honey bees (Apis mellifera) and stingless bees (Scaptotrigona bipunctata). The experiments were performed in southern Brazil, as follows: (i) larvae and adults were separately investigated, (ii) three or four different pollen diets were offered to bees, depending on larval or adult status, and (iii) two biological attributes, i.e., survivorship of larvae and adults and food intake by adults were evaluated. The diets were prepared with pollen from GM Eucalyptus 751K032; pollen from conventional Eucalyptus clone FGN-K, multifloral pollen or pure larval food. The insecticide dimethoate was used to evaluate the sensitivity of bees to toxic substances. Datasets were analyzed with Chi-square test, survival curves and repeated measures ANOVA. Results indicated no evidence of adverse effects of Eucalyptus pollen 751K032 on either honey bees or stingless bees assessed here. Therefore, the main findings suggest that the novel event may be considered harmless to these organisms since neither survivorship nor food consumption by bees were affected by it.

E. urophylla × E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus.

BACKGROUND: Eucalyptus urophylla × Eucalyptus grandis, an economically important forest tree, provides important raw material for energy and reduces damage to native forests. However, the absence of a high-quality E. urophylla × E. grandis reference genome has significantly hindered its evolution and genetic analysis. RESULTS: We successfully presented a high-quality reference genome of E. urophylla × E. grandis (545.75 Mb; scaffold N50, 51.62 Mb) using a combination of the Illumina, PacBio HiFi, and Hi-C sequencing platforms. A total of 34,502 genes and 58.56% of the repetitive sequences in this genome were annotated. Using genome evolution analyses, we identified a recent whole-genome duplication (WGD) event in E. urophylla × E. grandis. We further found that gene families associated with starch and sucrose metabolism, flavonoid biosynthesis, and plant-pathogen interaction were significantly expanded in E. urophylla × E. grandis. Moreover, comparative genomic and evolutionary analyses showed large structural variations among the different chromosomes of the 34 Eucalyptus accessions, which were divided into six clades. CONCLUSIONS: Overall, our findings provide a valuable resource for expanding our understanding of the E. urophylla × E. grandis genome evolution, genetic improvement, and its comparative biology.

Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations.

Gibberellins (GAs) are the key regulators controlling plant growth, wood production and the stress responses in perennial woody plants. The role of GA in regulating the above-mentioned processes in Eucalyptus remain largely unclear. There is still a lack of systematic identification and functional characterization of GA-related genes in Eucalyptus. In this study, a total of 59,948 expressed genes were identified from the major vegetative tissues of the E. grandis × E. urophylla using transcriptome sequencing. Then, the key gene families in each step of GA biosynthesis, degradation and signaling were investigated and compared with those of Arabidopsis, rice, and Populus. The expression profile generated using Real-time quantitative PCR showed that most of these genes exhibited diverse expression patterns in different vegetative organs and in response to abiotic stresses. Furthermore, we selectively overexpressed EguGA20ox1, EguGA20ox2 and EguGA2ox1 in both Arabidopsis and Eucalyptus via Agrobacterium tumefaciens or A. rhizogenes-mediated transformation. Though both Arabidopsis EguGA20ox1- and EguGA20ox2-overexpressing (OE) lines exhibited better vegetative growth performance, they were more sensitive to abiotic stress, unlike EguGA2ox1-OE plants, which exhibited enhanced stress resistance. Moreover, overexpression of EguGA20ox in Eucalyptus roots caused significantly accelerated hairy root initiation and elongation and improved root xylem differentiation. Our study provided a comprehensive and systematic study of the genes of the GA metabolism and signaling and identified the role of GA20ox and GA2ox in regulating plant growth, stress tolerance, and xylem development in Eucalyptus; this could benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.