The Transcriptomic Responses of Pinus massoniana to Drought Stress

Effects of PmaIAA27 and PmaARF15 genes on drought stress tolerance in pinus massoniana

BACKGROUND

Auxin plays an important role in plant resistance to abiotic stress. The modulation of gene expression by Auxin response factors (ARFs) and the inhibition of auxin/indole-3-acetic acid (Aux/IAA) proteins play crucial regulatory roles in plant auxin signal transduction. However, whether the stress resistance of Masson pine (Pinus massoniana), as a representative pioneer species, is related to Aux/IAA and ARF genes has not been thoroughly studied and explored.

RESULTS

The present study provides preliminary evidence for the regulatory role of the PmaIAA27 gene in abiotic stress response in Masson pine. We investigated the effects of drought and hormone treatments on Masson pine by examining the expression patterns of PmaIAA27 and PmaARF15 genes. Subsequently, we conducted gene cloning, functional testing using transgenic tobacco, and explored gene interactions. Exogenous auxin irrigation significantly downregulated the expression of PmaIAA27 while upregulating PmaARF15 in Masson pine seedlings. Moreover, transgenic tobacco with the PmaIAA27 gene exhibited a significant decrease in auxin content compared to control plants, accompanied by an increase in proline content - a known indicator of plant drought resistance. These findings suggest that overexpression of the PmaIAA27 gene may enhance drought resistance in Masson pine. To further investigate the interaction between PmaIAA27 and PmaARF15 genes, we performed bioinformatics analysis and yeast two-hybrid experiments which revealed interactions between PB1 structural region of PmaARF15 and PmaIAA27.

CONCLUSION

The present study provides new insights into the regulatory functions of Aux/IAA and ARF genes in Masson pine. Overexpression of PmaIAA gene may have negative effects on the growth of Masson pine, but may improve the drought resistance. Therefore, this study has great application prospects.

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.

Warming appears as the main risk of non-adaptedness for western Mediterranean relict fir forests under expected climate change scenarios

Circum-Mediterranean firs are considered among the most drought-sensitive species to climate change. Understanding the genetic basis of trees' adaptive capacity and intra-specific variability to drought avoidance is mandatory to define conservation measures, thus potentially preventing their extinction. We focus here on Abies pinsapo and Abies marocana, both relict tree species, endemic from south Spain and north Morocco, respectively. A total of 607 samples were collected from eight nuclei: six from Spanish fir and two from Moroccan fir. A genotyping by sequencing technique called double digestion restriction site-associated DNA sequencing (ddRAD-seq) was performed to obtain a genetic matrix based on single-nucleotide polymorphisms (SNPs). This matrix was utilized to study the genetic structure of A. pinsapo populations and to carry out selection signature studies. In order to understand how Spanish fir and Moroccan fir cope with climate change, genotype-environment associations (GEAs) were identified. Further, the vulnerability of these species to climate variations was estimated by the risk of non-adaptedness (RONA). The filtering of the de novo assembly of A. pinsapo provided 3,982 SNPs from 504 out of 509 trees sequenced. Principal component analysis (PCA) genetically separated Grazalema from the rest of the Spanish populations. However, FST values showed significant differences among the sampling points. We found 51 loci potentially under selection. Homolog sequences were found for some proteins related to abiotic stress response, such as dehydration-responsive element binding transcription factor, regulation of abscisic acid signaling, and methylation pathway. A total of 15 associations with 11 different loci were observed in the GEA studies, with the maximum temperature of the warmest month being the variable with the highest number of associated loci. This temperature sensitivity was also supported by the risk of non-adaptedness, which yielded a higher risk for both A. pinsapo and A. marocana under the high emission scenario (Representative Concentration Pathway (RCP) 8.5). This study sheds light on the response to climate change of these two endemic species.

Trichoderma longibrachiatum Inoculation Improves Drought Resistance and Growth of Pinus massoniana Seedlings through Regulating Physiological Responses and Soil Microbial Community

Drought stress poses a serious threat to Pinus massoniana seedling growth in southern China. Trichoderma species, as beneficial microorganisms, have been widely used in agriculture to enhance plant growth and drought tolerance, but the interaction mechanisms remain unclear. To investigate the effect of drought-resistant Trichoderma longibrachiatum inoculation on P. massoniana growth under drought stress, the plant physiological indicators and rhizosphere microbiome diversity were measured to identify Trichoderma-activated mechanisms. Trichoderma longibrachiatum inoculation significantly promoted P. massoniana growth under drought treatment, and enhanced nitrogen, phosphorus, and potassium absorption compared with those of non-inoculated seedlings. Trichoderma longibrachiatum treatment alleviated the damage to cell membranes and needle tissue structure, and significantly increased antioxidant enzyme activities, osmotic substance contents, and photosynthesis in P. massoniana in response to drought stress. Soil nutrient contents, activities of sucrase, phosphatase, and urease as well as the relative abundances of the dominant genera Burkholderia, Rhodanobacter, and Trichoderma were elevated in the rhizosphere soil of P. massoniana inoculated with T. longibrachiatum under drought stress. A network analysis showed that certain crucial dominant taxa driven by T. longibrachiatum inoculation, including Penicillium, Trichoderma, Simplicillium, Saitozyma, Burkholderia, Bradyrhizobium, Sinomonas, and Mycobacterium, had more correlations with other microorganisms in the soil. Trichoderma longibrachiatum enhanced P. massoniana seedling growth under drought stress by regulating physiological responses and soil microbial community.

Exogenous calcium regulates the growth and development of Pinus massoniana detecting by physiological, proteomic, and calcium-related genes expression analysis

Pinus massoniana is an important industrial crop tree species commonly used for timber and wood pulp for papermaking, rosin, and turpentine. This study investigated the effects of exogenous calcium (Ca) on P. massoniana seedling growth, development, and various biological processes and revealed the underlying molecular mechanisms. The results showed that Ca deficiency led to severe inhibition of seedling growth and development, whereas adequate exogenous Ca markedly improved growth and development. Many physiological processes were regulated by exogenous Ca. The underlying mechanisms involved diverse Ca-influenced biological processes and metabolic pathways. Calcium deficiency inhibited or impaired these pathways and processes, whereas sufficient exogenous Ca improved and benefited these cellular events by regulating several related enzymes and proteins. High levels of exogenous Ca facilitated photosynthesis and material metabolism. Adequate exogenous Ca supply relieved oxidative stress that occurred at low Ca levels. Enhanced cell wall formation, consolidation, and cell division also played a role in exogenous Ca-improved P. massoniana seedling growth and development. Calcium ion homeostasis and Ca signal transduction-related gene expression were also activated at high exogenous Ca levels. Our study facilitates the elucidation of the potential regulatory role of Ca in P. massoniana physiology and biology and is of guiding significance in Pinaceae plant forestry.

Differences in heat tolerance, water use efficiency and growth among Douglas-fir families and varieties evidenced by GWAS and common garden studies

Severe and frequent heat and drought events challenge the survival and development of long-generation trees. In this study, we investigated the genomic basis of heat tolerance, water use efficiency and growth by performing genome-wide association studies in coastal Douglas-fir (Pseudotsuga menziesii) and intervarietal (menziesii × glauca) hybrid seedlings. GWAS results identified 32 candidate genes involved in primary and secondary metabolism, abiotic stress and signaling, among other functions. Water use efficiency (inferred from carbon isotope discrimination), photosynthetic capacity (inferred from %N), height and heat tolerance (inferred from electrolyte leakage in a heat stress experiment) were significantly different among Douglas-fir families and varieties. High-elevation seed sources had increased water use efficiency, which could be a result of higher photosynthetic capacity. Similarly, families with greater heat tolerance also had higher water use efficiency and slower growth, suggesting a conservative growth strategy. Intervarietal hybrids showed increased heat tolerance (lower electrolyte leakage at 50 and 55 °C) and higher water use efficiency compared with coastal families, suggesting that hybridization might be a source of pre-adapted alleles to warming climates and should be considered for large-scale reforestation projects under increasingly arid conditions.

Comparative Transcriptome Analysis of Tolerant and Sensitive Genotypes of Common Bean (Phaseolus vulgaris L.) in Response to Terminal Drought Stress

We conducted a genome-wide transcriptomic analysis of three drought tolerant and sensitive genotypes of common bean to examine their transcriptional responses to terminal drought stress. We then conducted pairwise comparisons between the root and leaf transcriptomes from the resulting tissue based on combined transcriptomic data from the tolerant and sensitive genotypes. Our transcriptomic data revealed that 491 (6.4%) DEGs (differentially expressed genes) were upregulated in tolerant genotypes, whereas they were downregulated in sensitive genotypes; likewise, 396 (5.1%) DEGs upregulated in sensitive genotypes were downregulated in tolerant genotypes. Several transcription factors, heat shock proteins, and chaperones were identified in the study. Several DEGs in drought DB (data Base) overlapped between genotypes. The GO (gene ontology) terms for biological processes showed upregulation of DEGs in tolerant genotypes for sulfate and drug transmembrane transport when compared to sensitive genotypes. A GO term for cellular components enriched with upregulated DEGs for the apoplast in tolerant genotypes. These results substantiated the temporal pattern of root growth (elongation and initiation of root growth), and ABA-mediated drought response in tolerant genotypes. KEGG (kyoto encyclopedia of genes and genomes) analysis revealed an upregulation of MAPK (mitogen activated protein kinase) signaling pathways and plant hormone signaling pathways in tolerant genotypes. As a result of this study, it will be possible to uncover the molecular mechanisms of drought tolerance in response to terminal drought stress in the field. Further, genome-wide transcriptomic analysis of both tolerant and sensitive genotypes will assist us in identifying potential genes that may contribute to improving drought tolerance in the common bean.

Leaf transcriptome analysis of Medicago ruthenica revealed its response and adaptive strategy to drought and drought recovery

BACKGROUND

Drought is one of the main causes of losses in forage crop yield and animal production. Medicago ruthenica (L.) cv. Zhilixing is a high-yielding alfalfa cultivar also known for its high tolerance to drought. We analyzed the transcriptome profile of this cultivar throughout drought stress and recovery and we were able to describe its phased response through the expression profiles of overlapping gene networks and drought-specific genes.

RESULTS

The ABA and auxin signal transduction pathways are overlapping pathways in response to drought and drought recovery in forage crops. Medicago ruthenica (L.) cv. Zhilixing adopts different strategies at different degrees of drought stress. On the 9th day of drought, transcriptional regulations related to osmoregulation are enhanced mainly through increased activities of carbohydrate and amino acid metabolism, while photosynthetic activities were reduced to slow down growth. With drought prolonging, on the 12th day of drought, the synthesis of proline and other stored organic substances was suppressed in general. After recovery, Medicago ruthenica synthesizes flavonoids through the flavonoid biosynthesis pathway to remove accumulated ROS and repair the oxidative damage from water stress. In addition, the regulation of circadian rhythm seems to accelerate the drought recovery process.

CONCLUSIONS

Medicago ruthenica adapts to drought by regulating the osmoregulatory system and photosynthesis, which appears to involve the ABA and auxin signaling pathways as key regulators. Furthermore, the synthesis of flavonoids and the regulation of the circadian rhythm can accelerate the recovery process. These results enriched our knowledge of molecular responses to drought and drought recovery in Medicago ruthenica and provide useful information for the development of new legume forage grass varieties with improved adaptability to drought stress.

Identification of WRKY transcription factor family genes in Pinus massoniana Lamb. and their expression patterns and functions in response to drought stress

BACKGROUND

Pinus massoniana Lamb. is the timber species with the widest distribution and the largest afforestation area in China, providing a large amount of timber, turpentine and ecological products. Seasonal drought caused by climate warming severely constrains the quality and growth of P. massoniana forests. WRKY transcription factors play an important role in plant responses to abiotic stress. In this study, the molecular mechanisms by which P. massoniana responds to drought stress were analysed based on the P. massoniana WRKY (PmWRKY) family of genes.

RESULTS

Forty-three PmWRKYs are divided into three major families, 7 sub-families, and the conserved motifs are essentially the same. Among these 43 PmWRKYs express under drought stress but with different expression patterns in response to stress. PmWRKYs respond to drought stress induced by exogenous hormones of SA, ABA, and MeJA. The expression of PmWRKY6, PmWRKY10, and PmWRKY30 up-regulate in different families and tissues under drought stress, while PmWRKY22 down-regulate. Transgenetic tobaccos of PmWRKY31 are with lower malondialdehyde (MDA) content and higher proline (Pro) content than wild type (WT) tobaccos. In transgenic tobaccos of PmWRKY31, expression levels of related genes significantly improve, and drought tolerance enhance.

CONCLUSIONS

This study analysed the molecular biological characteristics of PmWRKYs and investigated the expression patterns and functions of PmWRKYs in response to drought stress in P. massoniana. The results of this study provide a basis for in-depth research of the molecular functions of PmWRKYs in response to drought stress.

Acclimation Strategy of Masson Pine (Pinus massoniana) by Limiting Flavonoid and Terpenoid Production under Low Light and Drought

Low light and drought often limit the growth and performance of Masson pines (Pinus massoniana) in the subtropical forest ecosystem of China. We speculated that stress-induced defensive secondary metabolites, such as flavonoids and terpenoids, might influence the growth of Masson pines, considering the existence of tradeoffs between growth and defense. However, the mechanisms of Masson pines responsive to low light and drought at the levels of these two metabolites remain unclear. In the present work, the compositions of flavonoids and terpenoids, as well as their biosynthetic pathways, were revealed through metabolome and transcriptome analyses, respectively, coupled with a study on carbon allocation using a ¹³CO₂-pulse-labeling experiment in two-year-old seedlings under low light (LL), drought (DR), and their combined stress (DL) compared to a control (CK). A total of 35 flavonoids and derivatives (LL vs. CK: 18; DR vs. CK: 20; and DL vs. CK: 18), as well as 29 terpenoids and derivatives (LL vs. CK: 23; DR vs. CK: 13; and DL vs. CK: 7), were differentially identified in the leaves. Surprisingly, most of them were decreased under all three stress regimes. At the transcriptomic level, most or all of the detected DEGs (differentially expressed genes) involved in the biosynthetic pathways of flavonoids and terpenoids were downregulated in phloem and xylem under stress treatments. This indicated that stress treatments limited the production of flavonoids and terpenoids. The reduction in the ¹³C allocation to stems might suggest that it is necessary for maintaining the growth of Masson pine seedlings at the whole-plant level by attenuating energetic resources to the biosynthetic pathways of flavonoids and terpenoids when facing the occurrence of adverse environments. Our results provide new insight into understanding the acclimation strategy of Masson pines or other conifers in adverse environments.

Transcriptome Approach Reveals the Response Mechanism of Heimia myrtifolia (Lythraceae, Myrtales) to Drought Stress

Drought is a major environmental condition that inhibits the development and cultivation of Heimia myrtifolia. The molecular processes of drought resistance in H. myrtifolia remain unknown, which has limited its application. In our study, transcriptome analyzes were compared across three treatment groups (CK, T1, and T2), to investigate the molecular mechanism of drought resistance. Plant leaves wilted and drooped as the duration of drought stress increased. The relative water content of the leaves declined dramatically, and relative electrolyte leakage rose progressively. Using an RNA-Seq approach, a total of 62,015 unigenes with an average length of 1730 bp were found, with 86.61% of them annotated to seven databases, and 14,272 differentially expressed genes (DEGs) were identified in drought stress. GO and KEGG enrichment analyzes of the DEGs revealed significantly enriched KEGG pathways, including photosynthesis, photosynthetic antenna proteins, plant hormone signal transduction, glutathione metabolism, and ascorbate and aldarate metabolism. Abscisic acid signal transduction was the most prevalent in the plant hormone signal transduction pathway, and other plant hormone signal transductions were also involved in the drought stress response. The transcription factors (including MYB, NAC, WRKY, and bHLH) and related differential genes on significantly enriched pathways all played important roles in the drought process, such as photosynthesis-related genes and antioxidant enzyme genes. In conclusion, this study will provide several genetic resources for further investigation of the molecular processes that will be beneficial to H. myrtifolia cultivation and breeding.

Regulatory Role of Circadian Clocks on ABA Production and Signaling, Stomatal Responses, and Water-Use Efficiency under Water-Deficit Conditions

Plants deploy molecular, physiological, and anatomical adaptations to cope with long-term water-deficit exposure, and some of these processes are controlled by circadian clocks. Circadian clocks are endogenous timekeepers that autonomously modulate biological systems over the course of the day–night cycle. Plants’ responses to water deficiency vary with the time of the day. Opening and closing of stomata, which control water loss from plants, have diurnal responses based on the humidity level in the rhizosphere and the air surrounding the leaves. Abscisic acid (ABA), the main phytohormone modulating the stomatal response to water availability, is regulated by circadian clocks. The molecular mechanism of the plant’s circadian clock for regulating stress responses is composed not only of transcriptional but also posttranscriptional regulatory networks. Despite the importance of regulatory impact of circadian clock systems on ABA production and signaling, which is reflected in stomatal responses and as a consequence influences the drought tolerance response of the plants, the interrelationship between circadian clock, ABA homeostasis, and signaling and water-deficit responses has to date not been clearly described. In this review, we hypothesized that the circadian clock through ABA directs plants to modulate their responses and feedback mechanisms to ensure survival and to enhance their fitness under drought conditions. Different regulatory pathways and challenges in circadian-based rhythms and the possible adaptive advantage through them are also discussed.

Analysis of the transcriptome of the needles and bark of Pinus radiata induced by bark stripping and methyl jasmonate

BACKGROUND

Plants are attacked by diverse insect and mammalian herbivores and respond with different physical and chemical defences. Transcriptional changes underlie these phenotypic changes. Simulated herbivory has been used to study the transcriptional and other early regulation events of these plant responses. In this study, constitutive and induced transcriptional responses to artificial bark stripping are compared in the needles and the bark of Pinus radiata to the responses from application of the plant stressor, methyl jasmonate. The time progression of the responses was assessed over a 4-week period.

RESULTS

Of the 6312 unique transcripts studied, 86.6% were differentially expressed between the needles and the bark prior to treatment. The most abundant constitutive transcripts were related to defence and photosynthesis and their expression did not differ between the needles and the bark. While no differential expression of transcripts were detected in the needles following bark stripping, in the bark this treatment caused an up-regulation and down-regulation of genes associated with primary and secondary metabolism. Methyl jasmonate treatment caused differential expression of transcripts in both the bark and the needles, with individual genes related to primary metabolism more responsive than those associated with secondary metabolism. The up-regulation of genes related to sugar break-down and the repression of genes related with photosynthesis, following both treatments was consistent with the strong down-regulation of sugars that has been observed in the same population. Relative to the control, the treatments caused a differential expression of genes involved in signalling, photosynthesis, carbohydrate and lipid metabolism as well as defence and water stress. However, non-overlapping transcripts were detected between the needles and the bark, between treatments and at different times of assessment. Methyl jasmonate induced more transcriptional responses in the bark than bark stripping, although the peak of expression following both treatments was detected 7 days post treatment application. The effects of bark stripping were localised, and no systemic changes were detected in the needles.

CONCLUSION

There are constitutive and induced differences in the needle and bark transcriptome of Pinus radiata. Some expression responses to bark stripping may differ from other biotic and abiotic stresses, which contributes to the understanding of plant molecular responses to diverse stresses. Whether the gene expression changes are heritable and how they differ between resistant and susceptible families identified in earlier studies needs further investigation.

Physiological, proteomic, and metabolomic analysis provide insights into Bacillus sp.-mediated salt tolerance in wheat

Herein, the inoculation with strain wp-6 promoted the growth of wheat seedlings by improving the energy production and conversion of wheat seedlings and alleviating salt stress. Soil salinization decreases crop productivity due to high toxicity of sodium ions to plants. Plant growth-promoting rhizobacteria (PGPR) have been demonstrated to alleviate salinity stress. However, the mechanism of PGPR in improving plant salt tolerance remains unclear. In this study, physiological analysis, proteomics, and metabolomics were applied to investigate the changes in wheat seedlings under salt stress (150 mM NaCl), both with and without plant root inoculation with wp-6 (Bacillus sp.). Under salt stress, root inoculation with strain wp-6 increased plant biomass (57%) and root length (25%). The Na⁺ content was reduced, while the K⁺ content and K⁺/Na⁺ ratio were increased. The content of malondialdehyde was decreased by 31.94% after inoculation of wp-6 under salt stress, while the content of proline, soluble sugar, and soluble protein were increased by 7.48%, 12.34%, and 4.12%, respectively. The peroxidase, catalase, and superoxide dismutase activities were increased after inoculation of wp-6 under salt stress. Galactose metabolism, phenylalanine metabolism, caffeine metabolism, ubiquinone and other terpenoid-quinone biosynthesis, and glutathione metabolism might play an important role in promoting the growth of salt-stressed wheat seedlings after the inoculation with wp-6. Interaction analysis of differentially expressed proteins and metabolites found that energy production and transformation-related proteins and six metabolites (D-arginine, palmitoleic acid, chlorophyllide b, rutin, pheophorbide a, and vanillylamine) were mainly involved in the growth of wheat seedlings after the inoculation with wp-6 under salt stress. Furthermore, correlation analysis found that inoculation with wp-6 promotes the growth of salt-stressed wheat seedlings mainly through regulating amino acid metabolism and porphyrin and chlorophyll metabolism. This study provides an eco-friendly method to increase agricultural productivity and paves a way to sustainable agriculture.

Candidate regulators and target genes of drought stress in needles and roots of Norway spruce

Drought stress impacts seedling establishment, survival and whole-plant productivity. Molecular responses to drought stress have been most extensively studied in herbaceous species, mostly considering only aboveground tissues. Coniferous tree species dominate boreal forests, which are predicted to be exposed to more frequent and acute drought as a result of ongoing climate change. The associated impact at all stages of the forest tree life cycle is expected to have large-scale ecological and economic impacts. However, the molecular response to drought has not been comprehensively profiled for coniferous species. We assayed the physiological and transcriptional response of Picea abies (L.) H. Karst seedling needles and roots after exposure to mild and severe drought. Shoots and needles showed an extensive reversible plasticity for physiological measures indicative of drought-response mechanisms, including changes in stomatal conductance (gs), shoot water potential and abscisic acid (ABA). In both tissues, the most commonly observed expression profiles in response to drought were highly correlated with the ABA levels. Still, root and needle transcriptional responses contrasted, with extensive root-specific down-regulation of growth. Comparison between previously characterized Arabidopsis thaliana L. drought-response genes and P. abies revealed both conservation and divergence of transcriptional response to drought. In P. abies, transcription factors belonging to the ABA responsive element(ABRE) binding/ABRE binding factors ABA-dependent pathway had a more limited role. These results highlight the importance of profiling both above- and belowground tissues, and provide a comprehensive framework to advance the understanding of the drought response of P. abies. The results demonstrate that a short-term, severe drought induces severe physiological responses coupled to extensive transcriptome modulation and highlight the susceptibility of Norway spruce seedlings to such drought events.

Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae)

BACKGROUND

Elucidating the effects of geography and selection on genetic variation is critical for understanding the relative importance of adaptation in driving differentiation and identifying the environmental factors underlying its occurrence. Adaptive genetic variation is common in tree species, especially widely distributed long-lived species. Pseudotaxus chienii can occupy diverse habitats with environmental heterogeneity and thus provides an ideal material for investigating the process of population adaptive evolution. Here, we characterize genetic and expression variation patterns and investigate adaptive genetic variation in P. chienii populations.

RESULTS

We generated population transcriptome data and identified 13,545 single nucleotide polymorphisms (SNPs) in 5037 unigenes across 108 individuals from 10 populations. We observed lower nucleotide diversity (π = 0.000701) among the 10 populations than observed in other gymnosperms. Significant negative correlations between expression diversity and nucleotide diversity in eight populations suggest that when the species adapts to the surrounding environment, gene expression and nucleotide diversity have a reciprocal relationship. Genetic structure analyses indicated that each distribution region contains a distinct genetic group, with high genetic differentiation among them due to geographical isolation and local adaptation. We used FST outlier, redundancy analysis, and latent factor mixed model methods to detect molecular signatures of local adaptation. We identified 244 associations between 164 outlier SNPs and 17 environmental variables. The mean temperature of the coldest quarter, soil Fe and Cu contents, precipitation of the driest month, and altitude were identified as the most important determinants of adaptive genetic variation. Most candidate unigenes with outlier signatures were related to abiotic and biotic stress responses, and the monoterpenoid biosynthesis and ubiquitin-mediated proteolysis KEGG pathways were significantly enriched in certain populations and deserve further attention in other long-lived trees.

CONCLUSIONS

Despite the strong population structure in P. chienii, genomic data revealed signatures of divergent selection associated with environmental variables. Our research provides SNPs, candidate unigenes, and biological pathways related to environmental variables to facilitate elucidation of the genetic variation in P. chienii in relation to environmental adaptation. Our study provides a promising tool for population genomic analyses and insights into the molecular basis of local adaptation.

Transcriptome Analysis of Needle and Root of Pinus Massoniana in Response to Continuous Drought Stress

Pinus massoniana Lamb. is an important coniferous tree species in ecological environment construction and sustainable forestry development. The function of gene gradual change and coexpression modules of needle and root parts of P. massoniana under continuous drought stress is unclear. The physiological and transcriptional expression profiles of P. massoniana seedlings from 1a half-sibling progeny during drought stress were measured and analyzed. As a result, under continuous drought conditions, needle peroxidase (POD) activity and proline content continued to increase. The malondialdehyde (MDA) content in roots continuously increased, and the root activity continuously decreased. The needles of P. massoniana seedlings may respond to drought mainly through regulating abscisic acid (ABA) and jasmonic acid (JA) hormone-related pathways. Roots may provide plant growth through fatty acid β-oxidative decomposition, and peroxisomes may contribute to the production of ROS, resulting in the upregulation of the antioxidant defense system. P. massoniana roots and needles may implement the same antioxidant mechanism through the glutathione metabolic pathway. This study provides basic data for identifying the drought response mechanisms of the needles and roots of P. massoniana.

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.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

Role of Biofilm Formation by Bacillus pumilus HR10 in Biocontrol against Pine Seedling Damping-Off Disease Caused by Rhizoctonia solani

The biocontrol process mediated by plant growth-promoting rhizobacteria (PGPR) relies on multiple mechanisms. Biofilm formation plays an important role in the ability of PGPR to control plant diseases. Bacillus pumilus HR10, one such PGPR, promotes the growth of Pinus thunbergii. This study showed that the wild-type strain B. pumilus HR10 produces a stable and mature biofilm in vitro. Biofilm-deficient mutants of B. pumilus HR10 with different phenotypes were screened by mutagenesis. The contents of extracellular polysaccharides (EPS) and proteins produced by the mutant strains were significantly reduced, and the biofilms of the mutants were weakened to varying degrees. The swarming abilities of the wild-type and mutant strains were positively correlated with biofilm formation. A colonization assay demonstrated that B. pumilus HR10 could colonize the roots of Pinus massoniana seedlings in a large population and persist, while biofilm-deficient mutants showed weak colonization ability. Furthermore, a biocontrol assay showed that biocontrol efficacy of the mutants was reduced to a certain degree. We determined the inhibitory activity of B. pumilus HR10 and its ability to induce systemic resistance against Rhizoctonia solani of plants. The synthesis of lipopeptide antibiotics is probably involved in biofilm formation by B. pumilus HR10. These observations not only provide a reference for further research about the coordinated action between biofilm formation and the multiple biocontrol mechanisms of B. pumilus HR10 but also improve the understanding of the regulatory pathway of biofilm formation by B. pumilus HR10.

Evaluation of Different Algorithms for Estimating the Growing Stock Volume of Pinus massoniana Plantations Using Spectral and Spatial Information from a SPOT6 Image

Precise growing stock volume (GSV) estimation is essential for monitoring forest carbon dynamics, determining forest productivity, assessing ecosystem forest services, and evaluating forest quality. We evaluated four machine learning methods: classification and regression trees (CART), support vector machines (SVM), artificial neural networks (ANN), and random forests (RF), for their reliability in the estimation of the GSV of Pinus massoniana plantations in China’s northern subtropical regions, using remote sensing data. For all four methods, models were generated using data derived from a SPOT6 image, namely the spectral vegetation indices (SVIs), texture parameters, or both. In addition, the effects of varying the size of the moving window on estimation precision were investigated. RF almost always yielded the greatest precision independently of the choice of input. ANN had the best performance when SVIs were used alone to estimate GSV. When using texture indices alone with window sizes of 3 × 5 × 5 or 9 × 9, RF achieved the best results. For CART, SVM, and RF, R2 decreased as the moving window size increased: the highest R2 values were achieved with 3 × 3 or 5 × 5 windows. When using textural parameters together with SVIs as the model input, RF achieved the highest precision, followed by SVM and CART. Models using both SVI and textural parameters as inputs had better estimating precision than those using spectral data alone but did not appreciably outperform those using textural parameters alone.

Physiological traits determining high adaptation potential of sweet briar (Rosa rubiginosa L.) at early stage of growth to dry lands

Little is known about mechanisms of sweet briar adaptation to dry habitats. The species is highly invasive and displaces native plants from dry lands of the southern hemisphere. This study evaluates physiological basis of Rosa rubiginosa L. adaptation to soil drought. We performed a pot soil drought experiment and assessed water relations, water use efficiency, gas exchange and photosynthetic apparatus activity. The study also measured the content of chlorophyll, soluble carbohydrates and proline and analyzed plant biomass growth. We hypothesized that the drought stress induced an effective mechanism enabling adaptation of young sweet briar roses to soil water deficit. The study identified several adaptation mechanisms of R. rubiginosa allowing the plant to survive soil drought. These included limiting transpiration and stomatal conductance, increasing the level of soluble sugars, reducing chlorophyll content, accumulating CO₂ in intercellular spaces, and increasing the quantum yield of electron transport from Q_A− to the PSI end electron acceptors. As a result, young sweet briar roses limited water loss and photoinhibition damage to the photosynthetic apparatus, which translated into consumption of soluble sugars for growth purposes. This study showed that photosynthesis optimization and increased activity of the photosynthetic apparatus made it possible to avoid photoinhibition and to effectively use water and sugars to maintain growth during water stress. This mechanism is probably responsible for the invasive nature of R. rubiginosa and its huge potential to displace native plant species from dry habitats of the southern hemisphere.

Transcriptional Analysis of Masson Pine (Pinus massoniana) under High CO2 Stress

To explore the molecular mechanism of the response of Masson pine (Pinus massoniana), the main coniferous tree in southern China, to high CO₂ stress, transcriptome sequencing was carried out to analyze the genome-wide responses of annual seedlings under different durations (0 h, 6 h, 12 h and 24 h) of high CO₂ stress. The results showed that a total of 3080/1908, 3110/2115 and 2684/1483 genes were up-/down-regulated after 6 h, 12 h and 24 h of treatment, respectively, compared with control check group (CK, 0 h). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that most of these differentially expressed genes (DEGs) were enriched in energy metabolism, carbohydrate synthesis, cell wall precursor synthesis and hormone regulation pathways. For energy metabolism, the expression of most genes involved in photosynthesis (including the light reaction and Calvin cycle) was generally inhibited, while the expression of genes related glycolysis, the tricarboxylic acid (TCA) cycle and PPP pathway was up-regulated. In addition, the increase in the CO₂ concentration induced the up-regulation of gene expression in the sucrose synthesis pathway. Among all starch synthesis genes, GBSS (granule-bound starch synthase) had the highest expression level. On the other hand, during the synthesis of hemicellulose and pectin (cell wall precursor substances), the expression levels of GMD (GDP-mannose 4,6-dehydratase), MGP (Mannose-1-phosphate guanylyl transferase) and RHM (Rhamnose biosynthetic enzyme) were the highest, suggesting that the synthesis of the raw materials hemicellulose and pectin in Masson pine under stress were mainly supplied by GDP-Man, GDP-Fuc and UDP-Rha. Finally, stress inhibited gene expression in the ABA (Abscisic Acid) synthesis pathway and induced gene expression in the GA (Gibberellin), SA (Salicylic acid), BR(Brassinolide) and MeJA (Methyl Jasmonate) pathways. Stomatal switches were regulated by hormonal interactions. This experiment elaborated on the response and molecular mechanism of Masson pine to CO₂ stress and aided in screening carbon sequestration genes for the corresponding molecular research of Masson pine in the future.

Transcriptome Analysis of Elm (Ulmus pumila) Fruit to Identify Phytonutrients Associated Genes and Pathways

Plant fruit is an important source of natural active phytonutrients that are profitable for human health. Elm (Ulmus pumila) fruit is considered as natural plant food in China that is rich in nutrients. In the present study, high-throughput RNA sequencing was performed in U. pumila edible fruits and leaves and 11,386 unigenes were filtered as dysregulated genes in fruit samples, including 5231 up- and 6155 downregulated genes. Hundreds of pathways were predicted to participate in seed development and phytonutrient biosynthesis in U. pumila by GO, MapMan, and KEGG enrichment analysis, including “seed maturation”, “glycine, serine, and threonine metabolism” and “phenylpropanoid biosynthesis”. ABA-mediated glucose response-related ethylene-activated signaling pathway (e.g., ABI4) were supposed to associate with elm fruit development; unsaturated fatty acids pathway (e.g., ACX2 and SAD) were predicted to participate in determination of fatty acid composition in elm fruit; flavonoid and coumarins biosynthesis (e.g., CYP98A3 and CCoAOMT1) were demonstrated to correlate with the bioactivity of elm fruits in human cancer and inflammation resistance. To provide more information about fruit developmental status, the qRT-PCR analysis for key genes of “phenylpropanoid biosynthesis” and “alpha-Linolenic acid metabolism” were conducted in samples of young fruits, ripe fruit, old fruit, and leaves. Two biosynthetic pathways for unsaturated fatty acid and Jasmonic acid (JA) were deduced to be involved in fruit development in U. pumila and the phenylpropanoid glycoside, syringin, was speculated to accumulate in the early development stages of elm fruit. Our transcriptome data supports molecular clues for seed development and biologically active substances in elm fruits.

Comparative analysis of abscisic acid levels and expression of abscisic acid-related genes in Scots pine and Norway spruce seedlings under water deficit

Abscisic acid (ABA) is one of the main participants in the regulation of plant responses to water deficiency. Knowledge of the ABA signal transduction pathways in gymnosperms is rather limited, especially in comparison with those in angiosperms. Seedlings of Scots pine and Norway spruce are known for their contrasting behaviour strategies under water deficit. To characterize the possible role of ABA in these differences, ABA dynamics were investigated under conditions of water deficit in seedlings of these two species. The content of ABA and its catabolites was followed in the roots and needles of seedlings of Pinus sylvestris and Picea abies under conditions of polyethylene glycol (PEG)-induced water deficiency (-0.15 and -0.5 MPa) for 10 days. The expression of the main genes for ABA-biosynthetic enzymes was also analysed. ABA showed more pronounced stress-dependent dynamics in pine roots than in spruce roots, whereas in needles, the response was greater for spruce than pine. The ABA increase during drought was mainly due to de novo synthesis and the shift in the balance between ABA synthesis and catabolism towards synthesis. The ABA-glucosyl ester did not serve as a reserve for the release of free ABA under water deficiency. The expression levels of the main ABA biosynthetic genes showed a weak or no correlation with changes in ABA content under water stress, i.e., the ABA content in the seedlings of both species was not directly linked to the transcript levels of the main ABA biosynthetic genes. Less-pronounced stress-induced changes in ABA in pine needles than in spruce needles may be related to pine seedlings having a less conservative strategy of growth and maintenance of water balance under water deficit.

Physiological and Transcriptome Analysis of a Yellow-Green Leaf Mutant in Birch (Betula platyphylla × B. Pendula)

Chlorophyll (Chl)-deficient mutants are ideal materials for the study of Chl biosynthesis, chloroplast development, and photosynthesis. Although the genes encoding key enzymes related to Chl biosynthesis have been well-characterized in herbaceous plants, rice (Oryza sativa L.), Arabidopsis (Arabidopsis thaliana), and maize (Zea mays L.), yellow-green leaf mutants have not yet been fully studied in tree species. In this work, we explored the molecular mechanism of the leaf color formation in a yellow-green leaf mutant (yl). We investigated the differentially expressed genes (DEGs) between yl and control plants (wild type birch (WT) and BpCCR1 overexpression line 11, (C11)) by transcriptome sequencing. Approximately 1163 genes (874 down-regulated and 289 up-regulated) and 930 genes (755 down-regulated and 175 up-regulated) were found to be differentially expressed in yl compared with WT and C11, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for DEGs revealed that photosynthesis antenna proteins represent the most significant enriched pathway. The expressions of photosynthesis antenna proteins are crucial to the leaf color formation in yl. We also found that Chl accumulate, leaf anatomical structure, photosynthesis, and growth were affected in yl. Taken together, our results not only provide the difference of phenomenal, physiological, and gene expression characteristics in leaves between yl mutant and control plants, but also provide a new insight into the mutation underlying the chlorotic leaf phenotype in birch.