Characterization of the Gene Expression Profile Response to Drought Stress in Populus ussuriensis Using PacBio SMRT and Illumina Sequencing

Allelic variations of WAK106-E2Fa-DPb1-UGT74E2 module regulate fibre properties in Populus tomentosa

Wood formation, intricately linked to the carbohydrate metabolism pathway, underpins the capacity of trees to produce renewable resources and offer vital ecosystem services. Despite their importance, the genetic regulatory mechanisms governing wood fibre properties in woody plants remain enigmatic. In this study, we identified a pivotal module comprising 158 high-priority core genes implicated in wood formation, drawing upon tissue-specific gene expression profiles from 22 Populus samples. Initially, we conducted a module-based association study in a natural population of 435 Populus tomentosa, pinpointing PtoDPb1 as the key gene contributing to wood formation through the carbohydrate metabolic pathway. Overexpressing PtoDPb1 led to a 52.91% surge in cellulose content, a reduction of 14.34% in fibre length, and an increment of 38.21% in fibre width in transgenic poplar. Moreover, by integrating co-expression patterns, RNA-sequencing analysis, and expression quantitative trait nucleotide (eQTN) mapping, we identified a PtoDPb1-mediated genetic module of PtoWAK106-PtoDPb1-PtoE2Fa-PtoUGT74E2 responsible for fibre properties in Populus. Additionally, we discovered the two PtoDPb1 haplotypes that influenced protein interaction efficiency between PtoE2Fa-PtoDPb1 and PtoDPb1-PtoWAK106, respectively. The transcriptional activation activity of the PtoE2Fa-PtoDPb1 haplotype-1 complex on the promoter of PtoUGT74E2 surpassed that of the PtoE2Fa-PtoDPb1 haplotype-2 complex. Taken together, our findings provide novel insights into the regulatory mechanisms of fibre properties in Populus, orchestrated by PtoDPb1, and offer a practical module for expediting genetic breeding in woody plants via molecular design.

Integration of genome-wide association studies, metabolomics, and transcriptomics reveals phenolic acid- and flavonoid-associated genes and their regulatory elements under drought stress in rapeseed flowers

Introduction

Biochemical and metabolic processes help plants tolerate the adverse effects of drought. In plants accumulating bioactive compounds, understanding the genetic control of the biosynthesis of biochemical pathways helps the discovery of candidate gene (CG)-metabolite relationships.

Methods

The metabolic profile of flowers in 119 rapeseed (Brassica napus) accessions was assessed over two irrigation treatments, one a well-watered (WW) condition and the other a drought stress (DS) regime. We integrated information gained from 52,157 single-nucleotide polymorphism (SNP) markers, metabolites, and transcriptomes to identify linked SNPs and CGs responsible for the genetic control of flower phenolic compounds and regulatory elements.

Results

In a genome-wide association study (GWAS), of the SNPs tested, 29,310 SNPs were qualified to assess the population structure and linkage disequilibrium (LD), of which several SNPs for radical scavenging activity (RSA) and total flavanol content (TFLC) were common between the two irrigation conditions and pleiotropic SNPs were found for chlorogenic and coumaric acids content. The principal component analysis (PCA) and stepwise regression showed that chlorogenic acid and epicatechin in WW and myricetin in DS conditions were the most important components for RSA. The hierarchical cluster analysis (HCA) showed that vanillic acid, myricetin, gallic acid, and catechin were closely associated in both irrigation conditions. Analysis of GWAS showed that 60 CGs were identified, of which 18 were involved in stress-induced pathways, phenylpropanoid pathway, and flavonoid modifications. Of the CGs, PAL1, CHI, UGT89B1, FLS3, CCR1, and CYP75B137 contributed to flavonoid biosynthetic pathways. The results of RNA sequencing (RNA-seq) revealed that the transcript levels of PAL, CHI, and CYP75B137 known as early flavonoid biosynthesis-related genes and FLS3, CCR1, and UGT89B1 related to the later stages were increased during drought conditions. The transcription factors (TFs) NAC035 and ERF119 related to flavonoids and phenolic acids were upregulated under drought conditions.

Discussion

These findings expand our knowledge on the response mechanisms to DS, particularly regarding the regulation of key phenolic biosynthetic genes in rapeseed. Our data also provided specific linked SNPs for marker-assisted selection (MAS) programs and CGs as resources toward realizing metabolomics-associated breeding of rapeseed.

Effects of drought and rehydration on root gene expression in seedlings of Pinus massoniana Lamb

The mechanisms underlying plant response to drought involve the expression of numerous functional and regulatory genes. Transcriptome sequencing based on the second- and/or third-generation high-throughput sequencing platforms has proven to be powerful for investigating the transcriptional landscape under drought stress. However, the full-length transcriptomes related to drought responses in the important conifer genus Pinus L. remained to be delineated using the third-generation sequencing technology. With the objectives of identifying the candidate genes responsible for drought and/or rehydration and clarifying the expression profile of key genes involved in drought regulation, we combined the third- and second-generation sequencing techniques to perform transcriptome analysis on seedling roots under drought stress and rewatering in the drought-tolerant conifer Pinus massoniana Lamb. A sum of 294,114 unique full-length transcripts were produced with a mean length of 3217 bp and N50 estimate of 5075 bp, including 279,560 and 124,438 unique full-length transcripts being functionally annotated and Gene Ontology enriched, respectively. A total of 4076, 6295 and 18,093 differentially expressed genes (DEGs) were identified in three pair-wise comparisons of drought-treatment versus control transcriptomes, including 2703, 3576 and 8273 upregulated and 1373, 2719 and 9820 downregulated DEGs, respectively. Moreover, 157, 196 and 691 DEGs were identified as transcription factors in the three transcriptome comparisons and grouped into 26, 34 and 44 transcription factor families, respectively. Gene Ontology enrichment analysis revealed that a remarkable number of DEGs were enriched in soluble sugar-related and cell wall-related processes. A subset of 75, 68 and 97 DEGs were annotated to be associated with starch, sucrose and raffinose metabolism, respectively, while 32 and 70 DEGs were associated with suberin and lignin biosynthesis, respectively. Weighted gene co-expression network analysis revealed modules and hub genes closely related to drought and rehydration. This study provides novel insights into root transcriptomic changes in response to drought dynamics in Masson pine and serves as a fundamental work for further molecular investigation on drought tolerance in conifers.

Combined analysis of transcriptome and metabolome reveals the molecular mechanism and candidate genes of Haloxylon drought tolerance

Haloxylon ammodendron and Haloxylon persicum, as typical desert plants, show strong drought tolerance and environmental adaptability. They are ideal model plants for studying the molecular mechanisms of drought tolerance. Transcriptomic and metabolomic analyses were performed to reveal the response mechanisms of H. ammodendron and H. persicum to a drought environment at the levels of transcription and physiological metabolism. The results showed that the morphological structures of H. ammodendron and H. persicum showed adaptability to drought stress. Under drought conditions, the peroxidase activity, abscisic acid content, auxin content, and gibberellin content of H. ammodendron increased, while the contents of proline and malondialdehyde decreased. The amino acid content of H. persicum was increased, while the contents of proline, malondialdehyde, auxin, and gibberellin were decreased. Under drought conditions, 12,233 and 17,953 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum , respectively, including members of multiple transcription factor families such as FAR1, AP2/ERF, C2H2, bHLH, MYB, C2C2, and WRKY that were significantly up-regulated under drought stress. In the positive ion mode, 296 and 452 differential metabolites (DEMs) were identified in H. ammodendron and H. persicum, respectively; in the negative ion mode, 252 and 354 DEMs were identified, primarily in carbohydrate and lipid metabolism. A combined transcriptome and metabolome analysis showed that drought stress promoted the glycolysis/gluconeogenesis pathways of H. ammodendron and H. persicum and increased the expression of amino acid synthesis pathways, consistent with the physiological results. In addition, transcriptome and metabolome were jointly used to analyze the expression changes of the genes/metabolites of H. ammodendron and H. persicum that were associated with drought tolerance but were regulated differently in the two plants. This study identified drought-tolerance genes and metabolites in H. ammodendron and H. persicum and has provided new ideas for studying the drought stress response of Haloxylon.