Medium term water deficit elicits distinct transcriptome responses in Eucalyptus species of contrasting environmental origin

Transcriptional reprogramming during recovery from drought stress in Eucalyptus grandis

The importance of drought as a constraint to agriculture and forestry is increasing with climate change. Genetic improvement of plants' resilience is one of the mitigation strategies to curb this threat. Although recovery from drought stress is important to long-term drought adaptation and has been considered as an indicator of dehydration tolerance in annual crops, this has not been well explored in forest trees. Thus, we aimed to investigate the physiological and transcriptional changes during drought stress and rewatering in Eucalyptus grandis W. Hill ex Maiden. We set up a greenhouse experiment where we imposed drought stress on 2-year-old seedlings and rewatered the recovery group after 17 days of drought. Our measurement of leaf stomatal conductance (gs) showed that, while gs was reduced by drought stress, it fully recovered after 5 days of rewatering. The RNA-seq analysis from stem samples revealed that genes related to known stress responses such as phytohormone and reactive oxygen species signaling were upregulated, while genes involved in metabolism and growth were downregulated due to drought stress. We observed reprogramming of signal transduction pathways and metabolic processes at 1 day of rewatering, indicating a quick response to rewatering. Our results suggest that recovery from drought stress may entail alterations in the jasmonic acid, salicylic acid, ethylene and brassinosteroid signaling pathways. Using co-expression network analysis, we identified hub genes, including the putative orthologs of ABI1, ABF2, ABF3, HAI2, BAM1, GolS2 and SIP1 during drought and CAT2, G6PD1, ADG1 and FD-1 during recovery. Taken together, by highlighting the molecular processes and identifying key genes, this study gives an overview of the mechanisms underlying the response of E. grandis to drought stress and recovery that trees may face repeatedly throughout their long life cycle. This provides a useful reference to the identification and further investigation of signaling pathways and target genes for future tree improvement.

The roles of divergent and parallel molecular evolution contributing to thermal adaptive strategies in trees

Local adaptation is a driver of biological diversity, and species may develop analogous (parallel evolution) or alternative (divergent evolution) solutions to similar ecological challenges. We expect these adaptive solutions would culminate in both phenotypic and genotypic signals. Using two Eucalyptus species (Eucalyptus grandis and Eucalyptus tereticornis) with overlapping distributions grown under contrasting 'local' temperature conditions to investigate the independent contribution of adaptation and plasticity at molecular, physiological and morphological levels. The link between gene expression and traits markedly differed between species. Divergent evolution was the dominant pattern driving adaptation (91% of all significant genes); but overlapping gene (homologous) responses were dependent on the determining factor (plastic, adaptive or genotype by environment interaction). Ninety-eight percent of the plastic homologs were similarly regulated, while 50% of the adaptive homologs and 100% of the interaction homologs were antagonistical. Parallel evolution for the adaptive effect in homologous genes was greater than expected but not in favour of divergent evolution. Heat shock proteins for E. grandis were almost entirely driven by adaptation, and plasticity in E. tereticornis. These results suggest divergent molecular evolutionary solutions dominated the adaptive mechanisms among species, even in similar ecological circumstances. Suggesting that tree species with overlapping distributions are unlikely to equally persist in the future.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Insight of transcriptional regulators reveals the tolerance mechanism of carpet-grass (Axonopus compressus) against drought

BACKGROUND

Carpet grass [Axonopus compressus (L.)] is an important warm-season perennial grass around the world and is known for its adaptability to varied environmental conditions. However, Carpet grass lacks enough data in public data banks, which confined our comprehension of the mechanisms of environmental adaptations, gene discovery, and development of molecular markers. In current study, the DEGs (differentially expressed genes) in Axonopus compressus under drought stress (DS) were identified and compared with CK (control) by RNA-Seq.

RESULTS

A total of 263,835 unigenes were identified in Axonopus compressus, and 201,303 (also added to the numbers of the remaining 2 databases) a sequence of unigenes significantly matched in at least one of the seven databases. A total of 153,697 (58.25%) unigenes classified to 144 KEGG pathways, and 7444 unigenes were expressed differentially between DS and CK, of which 4249 were up-regulated and 3195 were down-regulated unigenes. Of the 50 significantly enriched GO terms, 18, 6, and 14 items were related to BP, CC, and MF respectively. Analysis of KEGG enrichment revealed 2569 DEGs involved in 143 different pathways, under drought stress. 2747 DEGs were up-regulated and 2502 DEGs were down-regulated. Moreover, we identified 352 transcription factors (TFs) in Axonopus compressus, of which 270 were differentially expressed between CK and DS. The qRT-PCR validation experiment also supports the transcriptional response of Axonopus compressus against drought. Accuracy of transcriptome unigenes of Axonopus compressus was assessed with BLAST, which showed 3300 sequences of Axonopus compressus in the NCBI.

CONCLUSION

The 7444 unigenes were found to be between DS and CK treatments, which indicate the existence of a strong mechanism of drought tolerance in Axonopus compressus. The current findings provide the first framework for further investigations for the particular roles of these unigenes in Axonopus compressus in response to drought.

Haplotype- and SNP-Based GWAS for Growth and Wood Quality Traits in Eucalyptus cladocalyx Trees under Arid Conditions

The agricultural and forestry productivity of Mediterranean ecosystems is strongly threatened by the adverse effects of climate change, including an increase in severe droughts and changes in rainfall distribution. In the present study, we performed a genome-wide association study (GWAS) to identify single-nucleotide polymorphisms (SNPs) and haplotype blocks associated with the growth and wood quality of Eucalyptus cladocalyx, a tree species suitable for low-rainfall sites. The study was conducted in a progeny-provenance trial established in an arid site with Mediterranean patterns located in the southern Atacama Desert, Chile. A total of 87 SNPs and 3 haplotype blocks were significantly associated with the 6 traits under study (tree height, diameter at breast height, slenderness coefficient, first bifurcation height, stem straightness, and pilodyn penetration). In addition, 11 loci were identified as pleiotropic through Bayesian multivariate regression and were mainly associated with wood hardness, height, and diameter. In general, the GWAS revealed associations with genes related to primary metabolism and biosynthesis of cell wall components. Additionally, associations coinciding with stress response genes, such as GEM-related 5 and prohibitin-3, were detected. The findings of this study provide valuable information regarding genetic control of morphological traits related to adaptation to arid environments.

Molecular characterization and expression analysis of small heat shock protein 17.3 gene from Sorbus pohuashanensis (Hance) Hedl. in response to abiotic stress

Sorbus pohuashanensis, a native tree species in China that is distributed at high altitudes. However, the problem of adaptability when introducing S. pohuashanensis to low altitude areas has not been solved. sHSPs can respond and play an essential role when exposing to abiotic stresses for plants. In this study, we aimed to investigate the expression patterns underlying the abiotic stress response of the small heat shock protein 17.3 gene from S. pohuashanensis (SpHSP17.3) at growing low altitude. 1 to 4 years old seedlings of S. pohuashanensis were used as materials for the gene cloning, the tissue-specific expression and the expression analysis underlying the response to abiotic stress using the transgentic methods and qPCR. We identified the open reading frame (ORF) sequence of SpHSP17.3 of 471 bp, which encodes a 17.3 kD protein of 156 amino acids that is located in cytoplasmic. We found that SpHSP17.3 had the highest expression in the stem, followed sequentially by fruit, root, and flower. The expression level of SpHSP17.3 in the leaves was significantly induced by the high temperature (42 °C), NaCl salt and drought stress of S. pohuashanensis. Notably, the same SpHSP17.3 expression trend was detected in the SpHSP17.3-overexpressing homozygous transgenic Arabidopsis underlying the high temperature, NaCl salt and drought stress, and the SpHSP17.3-overexpressing homozygous transgenic Arabidopsis also showed higher seed germination rates under the NaCl salt stress conditions. Our results suggested that SpHSP17.3 is involved in the response to high temperature, Nacl salt, and drought stress which would play a certain effect in the adaptability of introduction and domestication of S. pohuashanensis.

Hardwood Tree Genomics: Unlocking Woody Plant Biology

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

Reconfigured Cyanogenic Glucoside Biosynthesis in Eucalyptus cladocalyx Involves a Cytochrome P450 CYP706C55

Cyanogenic glucosides are a class of specialized metabolites widespread in the plant kingdom. Cyanogenic glucosides are α-hydroxynitriles, and their hydrolysis releases toxic hydrogen cyanide, providing an effective chemical defense against herbivores. Eucalyptus cladocalyx is a cyanogenic tree, allocating up to 20% of leaf nitrogen to the biosynthesis of the cyanogenic monoglucoside, prunasin. Here, mass spectrometry analyses of E. cladocalyx tissues revealed spatial and ontogenetic variations in prunasin content, as well as the presence of the cyanogenic diglucoside amygdalin in flower buds and flowers. The identification and biochemical characterization of the prunasin biosynthetic enzymes revealed a unique enzyme configuration for prunasin production in E. cladocalyx This result indicates that a multifunctional cytochrome P450 (CYP), CYP79A125, catalyzes the initial conversion of l-phenylalanine into its corresponding aldoxime, phenylacetaldoxime; a function consistent with other members of the CYP79 family. In contrast to the single multifunctional CYP known from other plant species, the conversion of phenylacetaldoxime to the α-hydroxynitrile, mandelonitrile, is catalyzed by two distinct CYPs. CYP706C55 catalyzes the dehydration of phenylacetaldoxime, an unusual CYP reaction. The resulting phenylacetonitrile is subsequently hydroxylatedby CYP71B103 to form mandelonitrile. The final glucosylation step to yield prunasin is catalyzed by a UDP-glucosyltransferase, UGT85A59. Members of the CYP706 family have not been reported previously to participate in the biosynthesis of cyanogenic glucosides, and the pathway structure in E. cladocalyx represents an example of convergent evolution in the biosynthesis of cyanogenic glucosides in plants.