The microtubule cytoskeleton acts as a sensor for stress response signaling in plants

Dehydration-responsive cytoskeleton proteome of rice reveals reprograming of key molecular pathways to mediate metabolic adaptation and cell survival

The plant cytoskeletal proteins play a key role that control cytoskeleton dynamics, contributing to crucial biological processes such as cell wall morphogenesis, stomatal conductance and abscisic acid accumulation in repercussion to water-deficit stress or dehydration. Yet, it is still completely unknown which specific biochemical processes and regulatory mechanisms the cytoskeleton uses to drive dehydration tolerance. To better understand the role of cytoskeleton, we developed the dehydration-responsive cytoskeletal proteome map of a resilient rice cultivar. Initially, four-week-old rice plants were exposed to progressive dehydration, and the magnitude of dehydration-induced compensatory physiological responses was monitored in terms of physicochemical indices. The organelle fractionation in conjunction with label-free quantitative proteome analysis led to the identification of 955 dehydration-responsive cytoskeletal proteins (DRCPs). To our knowledge, this is the first report of a stress-responsive plant cytoskeletal proteome, representing the largest inventory of cytoskeleton and cytoskeleton-associated proteins. The DRCPs were apparently involved in a wide array of intra-cellular molecules transportation, organelles positioning, cytoskeleton organization followed by different metabolic processes including amino acid metabolism. These findings presented open a unique view on global regulation of plant cytoskeletal proteome is intimately linked to cellular metabolic rewiring of adaptive responses, and potentially confer dehydration tolerance, especially in rice, and other crop species, in general.

Characterization of Soybean Events with Enhanced Expression of the Microtubule-Associated Protein 65-1 (MAP65-1)

Microtubule-associated protein 65-1 (MAP65-1) protein plays an essential role in plant cellular dynamics through impacting stabilization of the cytoskeleton by serving as a crosslinker of microtubules. The role of MAP65-1 in plants has been associated with phenotypic outcomes in response to various environmental stresses. The Arabidopsis MAP65-1 (AtMAP65-1) is a known virulence target of plant bacterial pathogens and is thus a component of plant immunity. Soybean events were generated that carry transgenic alleles for both AtMAP65-1 and GmMAP65-1, the soybean AtMAP65-1 homolog, under control of cauliflower mosaic virus 35S promoter. Both AtMAP65-1 and GmMAP65-1 transgenic soybeans are more resistant to challenges by the soybean bacterial pathogen Pseudomonas syringae pv. glycinea and the oomycete pathogen Phytophthora sojae, but not the soybean cyst nematode, Heterodera glycines. Soybean plants expressing AtMAP65-1 and GmMAP65-1 also display a tolerance to the herbicide oryzalin, which has a mode of action to destabilize microtubules. In addition, GmMAP65-1-expressing soybean plants show reduced cytosol ion leakage under freezing conditions, hinting that ectopic expression of GmMAP65-1 may enhance cold tolerance in soybean. Taken together, overexpression of AtMAP65-1 and GmMAP65-1 confers tolerance of soybean plants to various biotic and abiotic stresses. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A multi-level approach reveals key physiological and molecular traits in the response of two rice genotypes subjected to water deficit at the reproductive stage

Rice is more vulnerable to drought than maize, wheat, and sorghum because its water requirements remain high throughout the rice life cycle. The effects of drought vary depending on the timing, intensity, and duration of the events, as well as on the rice genotype and developmental stage. It can affect all levels of organization, from genes to the cells, tissues, and/or organs. In this study, a moderate water deficit was applied to two contrasting rice genotypes, IAC 25 and CIRAD 409, during their reproductive stage. Multi-level transcriptomic, metabolomic, physiological, and morphological analyses were performed to investigate the complex traits involved in their response to drought. Weighted gene network correlation analysis was used to identify the specific molecular mechanisms regulated by each genotype, and the correlations between gene networks and phenotypic traits. A holistic analysis of all the data provided a deeper understanding of the specific mechanisms regulated by each genotype, and enabled the identification of gene markers. Under non-limiting water conditions, CIRAD 409 had a denser shoot, but shoot growth was slower despite better photosynthetic performance. Under water deficit, CIRAD 409 was weakly affected regardless of the plant level analyzed. In contrast, IAC 25 had reduced growth and reproductive development. It regulated transcriptomic and metabolic activities at a high level, and activated a complex gene regulatory network involved in growth-limiting processes. By comparing two contrasting genotypes, the present study identified the regulation of some fundamental processes and gene markers, that drive rice development, and influence its response to water deficit, in particular, the importance of the biosynthetic and regulatory pathways for cell wall metabolism. These key processes determine the biological and mechanical properties of the cell wall and thus influence plant development, organ expansion, and turgor maintenance under water deficit. Our results also question the genericity of the antagonism between morphogenesis and organogenesis observed in the two genotypes.

The role and risks of selective adaptation in extreme coral habitats

The alarming rate of climate change demands new management strategies to protect coral reefs. Environments such as mangrove lagoons, characterized by extreme variations in multiple abiotic factors, are viewed as potential sources of stress-tolerant corals for strategies such as assisted evolution and coral propagation. However, biological trade-offs for adaptation to such extremes are poorly known. Here, we investigate the reef-building coral Porites lutea thriving in both mangrove and reef sites and show that stress-tolerance comes with compromises in genetic and energetic mechanisms and skeletal characteristics. We observe reduced genetic diversity and gene expression variability in mangrove corals, a disadvantage under future harsher selective pressure. We find reduced density, thickness and higher porosity in coral skeletons from mangroves, symptoms of metabolic energy redirection to stress response functions. These findings demonstrate the need for caution when utilizing stress-tolerant corals in human interventions, as current survival in extremes may compromise future competitive fitness.

CRK41 Modulates Microtubule Depolymerization in Response to Salt Stress in Arabidopsis

The pivotal role of cysteine-rich receptor-like kinases (CRKs) in modulating growth, development, and responses to stress has been widely acknowledged in Arabidopsis. However, the function and regulation of CRK41 has remained unclear. In this study, we demonstrate that CRK41 is critical for modulating microtubule depolymerization in response to salt stress. The crk41 mutant exhibited increased tolerance, while overexpression of CRK41 led to hypersensitivity to salt. Further analysis revealed that CRK41 interacts directly with the MAP kinase3 (MPK3), but not with MPK6. Inactivation of either MPK3 or MPK6 could abrogate the salt tolerance of the crk41 mutant. Upon NaCl treatment, microtubule depolymerization was heightened in the crk41 mutant, yet alleviated in the crk41mpk3 and crk41mpk6 double mutants, indicating that CRK41 suppresses MAPK-mediated microtubule depolymerizations. Collectively, these results reveal that CRK41 plays a crucial role in regulating microtubule depolymerization triggered by salt stress through coordination with MPK3/MPK6 signalling pathways, which are key factors in maintaining microtubule stability and conferring salt stress resistance in plants.

Evolutionary Relationships and Divergence of Filamin Gene Family Involved in Development and Stress in Cotton (Gossypium hirsutum L.)

Filamin protein is characterized by an N-terminal actin-binding domain that is followed by 24 Ig (immunoglobulin)-like repeats, which act as hubs for interactions with a variety of proteins. In humans, this family has been found to be involved in cancer cell invasion and metastasis and can be involved in a variety of growth signal transduction processes, but it is less studied in plants. Therefore, in this study, 54 Filamin gene family members from 23 plant species were investigated and divided into two subfamilies: FLMN and GEX2. Subcellular localization showed that most of the Filamin gene family members were located in the cell membrane. A total of 47 Filamin gene pairs were identified, most of which were whole-genome copies. Through the analyses of cis-acting elements, expression patterns and quantitative fluorescence, it was found that GH_ A02G0519 and GH_ D02G0539 are mainly expressed in the reproductive organs of upland cotton, and their interacting proteins are also related to the fertilization process, whereas GH_A02G0216 and GH_D02G0235 were related to stress. Thus, it is speculated that two genes of the GEX2 subfamily (GH_A02G0519 and GH_D02G0539) may be involved in the reproductive development of cotton and may affect the fertilization process of cotton. This study provides a theoretical basis for the further study of the cotton Filamin gene family.

Abscisic acid induces the expression of AsKIN during the recovery period of garlic cryopreservation

Abscisic acid induced the expression of AsKIN during the recovery period of garlic cryopreservation. AsKIN was identified as a gene involved in cold and osmotic stress resistance. Cryopreservation has been proven to be effective in removing viruses from garlic. However, oxidative damage in cryopreservation has a significant impact on the survival after preservation. Abscisic acid (ABA) has been shown to reduce oxidative stress and promote the survival after cryopreservation. However, it is not clear which genes play important roles in this process. In this study, we added ABA to the dehydration step and analyzed the transcriptomic divergences between the ABA-treated group and the control group in three cryogenic steps (dehydration, unloading and recovery). By short time-series expression miner (STEM) analysis and weighted gene co-expression network analysis (WGCNA), the recovery step was identified as the period of significant changes in gene expression levels in cryopreservation. The addition of ABA promoted the upregulated expression of microtubule-related genes in the recovery step. We further identified AsKIN as a hub gene in the recovery step and verified its function. The results showed that overexpression of AsKIN enhanced the tolerance of Arabidopsis to oxidative stress in cryopreservation, influenced the expression of genes in response to cold and osmotic stress and promoted plant growth after stress. The AsKIN gene is likely to be involved in the plant response to cold stress and osmotic stress. These results reveal the molecular mechanisms of ABA in cryopreservation and elucidate the potential biological functions of the kinesin-14 subfamily.

Integrative proteomic and physiological analyses of the molecular response to dessication-stress in Auricularia fibrillifera

Drought stress is one of the main factors influencing the growth and development of an organism. Auricularia fibrillifera has strong dessication resistance. In A. fibrillifera under dessication-stress, the melanin content of fruiting bodies elevated significantly by >10-fold compared with the control. Folate content also increased sharply but decreased significantly after rehydration, and amino acid and biotin levels increased by 40.11 and 22.14%, respectively. In proteomic analysis, 1,572 and 21 differentially abundant proteins (DAPs) were identified under dessication-stress and rehydration, respectively. A large number of DAPs were annotated in "amino acid metabolism," "carbohydrate metabolism," and "translation" pathways, and the DAPs related to osmotic regulation and antioxidant enzymes were significantly increased in abundance. Transcriptome-proteome association analysis showed that most DAPs (30) were annotated in the "biosynthesis of antibiotics" pathway. DAPs and corresponding differentially expressed genes were all up-regulated in the "biotin biosynthesis" pathway and associated with "folate biosynthesis" and "phenylalanine, tyrosine, and tryptophan biosynthesis." In the analysis of protein-protein interactions, the DAPs annotated in the "phenylalanine, tyrosine, and tryptophan biosynthesis" pathway had the strongest interactions with other DAPs. These enriched pathways could enhance amino acid, folate, biotin, and melanin levels during desiccation stress, which is consistent with the physiological data (amino acid, folate, biotin, and melanin contents). In addition, many DAPs related to the cytoskeleton were significantly increased in abundance under dessication-stress. Physiological and transcriptome data were in agreement with proteomic results. This work provides valuable insight into the dessication-tolerant mechanisms of A. fibrillifera.

Chromosome-level and graphic genomes provide insights into metabolism of bioactive metabolites and cold-adaption of Pueraria lobata var. montana

Pueraria lobata var. montana (P. montana) belongs to the genus Pueraria and originated in Asia. Compared with its sister P. thomsonii, P. montana has stronger growth vigour and cold-adaption but contains less bioactive metabolites such as puerarin. To promote the investigation of metabolic regulation and genetic improvement of Pueraria, the present study reports a chromosome-level genome of P. montana with length of 978.59 Mb and scaffold N50 of 80.18 Mb. Comparative genomics analysis showed that P. montana possesses smaller genome size than that of P. thomsonii owing to less repeat sequences and duplicated genes. A total of 6,548 and 4,675 variety-specific gene families were identified in P. montana and P. thomsonii, respectively. The identified variety-specific and expanded/contracted gene families related to biosynthesis of bioactive metabolites and microtubules are likely the causes for the different characteristics of metabolism and cold-adaption of P. montana and P. thomsonii. Moreover, a graphic genome was constructed based on 11 P. montana accessions. Total 92 structural variants were identified and most of which are related to stimulus-response. In conclusion, the chromosome-level and graphic genomes of P. montana will not only facilitate the studies of evolution and metabolic regulation, but also promote the breeding of Pueraria.

Spatial variation in high temperature-regulated gene expression predicts evolution of plasticity with climate change in the scarlet monkeyflower

A major way that organisms can adapt to changing environmental conditions is by evolving increased or decreased phenotypic plasticity. In the face of current global warming, more attention is being paid to the role of plasticity in maintaining fitness as abiotic conditions change over time. However, given that temporal data can be challenging to acquire, a major question is whether evolution in plasticity across space can predict adaptive plasticity across time. In growth chambers simulating two thermal regimes, we generated transcriptome data for western North American scarlet monkeyflowers (Mimulus cardinalis) collected from different latitudes and years (2010 and 2017) to test hypotheses about how plasticity in gene expression is responding to increases in temperature, and if this pattern is consistent across time and space. Supporting the genetic compensation hypothesis, individuals whose progenitors were collected from the warmer-origin northern 2017 descendant cohort showed lower thermal plasticity in gene expression than their cooler-origin northern 2010 ancestors. This was largely due to a change in response at the warmer (40°C) rather than cooler (20°C) treatment. A similar pattern of reduced plasticity, largely due to a change in response at 40°C, was also found for the cooler-origin northern versus the warmer-origin southern population from 2017. Our results demonstrate that reduced phenotypic plasticity can evolve with warming and that spatial and temporal changes in plasticity predict one another.

Silver Nanoparticles Alter Microtubule Arrangement, Dynamics and Stress Phytohormone Levels

The superior properties of silver nanoparticles (AgNPs) has resulted in their broad utilization worldwide, but also the risk of irreversible environment infestation. The plant cuticle and cell wall can trap a large part of the nanoparticles and thus protect the internal cell structures, where the cytoskeleton, for example, reacts very quickly to the threat, and defense signaling is subsequently triggered. We therefore used not only wild-type Arabidopsis seedlings, but also the glabra 1 mutant, which has a different composition of the cuticle. Both lines had GFP-labeled microtubules (MTs), allowing us to observe their arrangement. To quantify MT dynamics, we developed a new microscopic method based on the FRAP technique. The number and growth rate of MTs decreased significantly after AgNPs, similarly in both lines. However, the layer above the plasma membrane thickened significantly in wild-type plants. The levels of three major stress phytohormone derivatives—jasmonic, abscisic, and salicylic acids—after AgNP (with concomitant Ag⁺) treatment increased significantly (particularly in mutant plants) and to some extent resembled the plant response after mechanical stress. The profile of phytohormones helped us to estimate the mechanism of response to AgNPs and also to understand the broader physiological context of the observed changes in MT structure and dynamics.

Microtubule organization changes severely after mannitol and n-butanol treatments inducing microspore embryogenesis in bread wheat

BACKGROUND

A mannitol stress treatment and a subsequent application of n-butanol, known as a microtubule-disrupting agent, enhance microspore embryogenesis (ME) induction and plant regeneration in bread wheat. To characterize changes in cortical (CMT) and endoplasmic (EMT) microtubules organization and dynamics, associated with ME induction treatments, immunocytochemistry studies complemented by confocal laser scanning microscopy (CLSM) were accomplished. This technique has allowed us to perform advanced 3- and 4D studies of MT architecture. The degree of MT fragmentation was examined by the relative fluorescence intensity quantification.

RESULTS

In uni-nucleated mannitol-treated microspores, severe CMT and EMT fragmentation occurs, although a complex network of short EMT bundles protected the nucleus. Additional treatment with n-butanol resulted in further depolymerization of both CMT and EMT, simultaneously with the formation of MT aggregates in the perinuclear region. Some aggregates resembled a preprophase band. In addition, a portion of the microspores progressed to the first mitotic division during the treatments. Bi-nucleate pollen-like structures showed a high MT depolymerization after mannitol treatment and numerous EMT bundles around the vegetative and generative nuclei after n-butanol. Interestingly, bi-nucleate symmetric structures showed prominent stabilization of EMT.

CONCLUSIONS

Fragmentation and stabilization of microtubules induced by mannitol- and n-butanol lead to new configurations essential for the induction of microspore embryogenesis in bread wheat. These results provide robust insight into MT dynamics during EM induction and open avenues to address newly targeted treatments to induce ME in recalcitrant species.

QTL dissection and mining of candidate genes for Ascochyta fabae and Orobanche crenata resistance in faba bean (Vicia faba L.)

BACKGROUND

Ascochyta blight caused by Ascochyta fabae Speg. and broomrape (Orobanche crenata) are among the economically most significant pathogens of faba bean. Several QTLs conferring resistance against the two pathogens have been identified and validated in different genetic backgrounds. The aim of this study was to saturate the most stable QTLs for ascochyta and broomrape resistance in two Recombinant Inbred Line (RIL) populations, 29H x Vf136 and Vf6 x Vf136, to identify candidate genes conferring resistance against these two pathogens.

RESULTS

We exploited the synteny between faba bean and the model species Medicago truncatula by selecting a set of 219 genes encoding putative WRKY transcription factors and defense related proteins falling within the target QTL intervals, for genotyping and marker saturation in the two RIL populations. Seventy and 50 of the candidate genes could be mapped in 29H x Vf136 and Vf6 x Vf136, respectively. Besides the strong reduction of the QTL intervals, the mapping process allowed replacing previous dominant and pedigree-specific RAPD flanking markers with robust and transferrable SNP markers, revealing promising candidates for resistance against the two pathogens.

CONCLUSIONS

Although further efforts in association mapping and expression studies will be required to corroborate the candidate genes for resistance, the fine-mapping approach proposed here increases the genetic resolution of relevant QTL regions and paves the way for an efficient deployment of useful alleles for faba bean ascochyta and broomrape resistance through marker-assisted breeding.

Gene Expression Profiles Suggest a Better Cold Acclimation of Polyploids in the Alpine Species Ranunculus kuepferi (Ranunculaceae)

Alpine habitats are shaped by harsh abiotic conditions and cold climates. Temperature stress can affect phenotypic plasticity, reproduction, and epigenetic profiles, which may affect acclimation and adaptation. Distribution patterns suggest that polyploidy seems to be advantageous under cold conditions. Nevertheless, whether temperature stress can induce gene expression changes in different cytotypes, and how the response is initialized through gene set pathways and epigenetic control remain vague for non-model plants. The perennial alpine plant Ranunculus kuepferi was used to investigate the effect of cold stress on gene expression profiles. Diploid and autotetraploid individuals were exposed to cold and warm conditions in climate growth chambers and analyzed via transcriptome sequencing and qRT-PCR. Overall, cold stress changed gene expression profiles of both cytotypes and induced cold acclimation. Diploids changed more gene set pathways than tetraploids, and suppressed pathways involved in ion/cation homeostasis. Tetraploids mostly activated gene set pathways related to cell wall and plasma membrane. An epigenetic background for gene regulation in response to temperature conditions is indicated. Results suggest that perennial alpine plants can respond to temperature extremes via altered gene expression. Tetraploids are better acclimated to cold conditions, enabling them to colonize colder climatic areas in the Alps.

Microtubule Dynamics Plays a Vital Role in Plant Adaptation and Tolerance to Salt Stress

Although recent studies suggest that the plant cytoskeleton is associated with plant stress responses, such as salt, cold, and drought, the molecular mechanism underlying microtubule function in plant salt stress response remains unclear. We performed a comparative proteomic analysis between control suspension-cultured cells (A0) and salt-adapted cells (A120) established from Arabidopsis root callus to investigate plant adaptation mechanisms to long-term salt stress. We identified 50 differentially expressed proteins (45 up- and 5 down-regulated proteins) in A120 cells compared with A0 cells. Gene ontology enrichment and protein network analyses indicated that differentially expressed proteins in A120 cells were strongly associated with cell structure-associated clusters, including cytoskeleton and cell wall biogenesis. Gene expression analysis revealed that expressions of cytoskeleton-related genes, such as FBA8, TUB3, TUB4, TUB7, TUB9, and ACT7, and a cell wall biogenesis-related gene, CCoAOMT1, were induced in salt-adapted A120 cells. Moreover, the loss-of-function mutant of Arabidopsis TUB9 gene, tub9, showed a hypersensitive phenotype to salt stress. Consistent overexpression of Arabidopsis TUB9 gene in rice transgenic plants enhanced tolerance to salt stress. Our results suggest that microtubules play crucial roles in plant adaptation and tolerance to salt stress. The modulation of microtubule-related gene expression can be an effective strategy for developing salt-tolerant crops.

A Transcriptomic Approach to Understanding the Combined Impacts of Supra-Optimal Temperatures and CO2 Revealed Different Responses in the Polyploid Coffea arabica and Its Diploid Progenitor C. canephora

Understanding the effect of extreme temperatures and elevated air (CO₂) is crucial for mitigating the impacts of the coffee industry. In this work, leaf transcriptomic changes were evaluated in the diploid C. canephora and its polyploid C. arabica, grown at 25 °C and at two supra-optimal temperatures (37 °C, 42 °C), under ambient (aCO₂) or elevated air CO₂ (eCO₂). Both species expressed fewer genes as temperature rose, although a high number of differentially expressed genes (DEGs) were observed, especially at 42 °C. An enrichment analysis revealed that the two species reacted differently to the high temperatures but with an overall up-regulation of the photosynthetic machinery until 37 °C. Although eCO₂ helped to release stress, 42 °C had a severe impact on both species. A total of 667 photosynthetic and biochemical related-DEGs were altered with high temperatures and eCO₂, which may be used as key probe genes in future studies. This was mostly felt in C. arabica, where genes related to ribulose-bisphosphate carboxylase (RuBisCO) activity, chlorophyll a-b binding, and the reaction centres of photosystems I and II were down-regulated, especially under 42°C, regardless of CO₂. Transcriptomic changes showed that both species were strongly affected by the highest temperature, although they can endure higher temperatures (37 °C) than previously assumed.

Whole-Transcriptome RNA Sequencing Reveals the Global Molecular Responses and CeRNA Regulatory Network of mRNAs, lncRNAs, miRNAs and circRNAs in Response to Salt Stress in Sugar Beet (Beta vulgaris)

Sugar beet is an important sugar-yielding crop with some tolerance to salt, but the mechanistic basis of this tolerance is not known. In the present study, we have used whole-transcriptome RNA-seq and degradome sequencing in response to salt stress to uncover differentially expressed (DE) mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in both leaves and roots. A competitive endogenous RNA (ceRNA) network was constructed with the predicted DE pairs, which revealed regulatory roles under salt stress. A functional analysis suggests that ceRNAs are implicated in copper redistribution, plasma membrane permeability, glycometabolism and energy metabolism, NAC transcription factor and the phosphoinositol signaling system. Overall, we conducted for the first time a full transcriptomic analysis of sugar beet under salt stress that involves a potential ceRNA network, thus providing a basis to study the potential functions of lncRNAs/circRNAs.

Salt Stress-Induced Structural Changes Are Mitigated in Transgenic Tomato Plants Over-Expressing Superoxide Dismutase

Various abiotic stresses cause the appearance of reactive oxygen species (ROS) in plant cells, which seriously damage the cellular structures. The engineering of transgenic plants with higher production of ROS-scavenging enzyme in plant cells could protect the integrity of such a fine intracellular structure as the cytoskeleton and each cellular compartment. We analyzed the morphological changes in root tip cells caused by the application of iso-osmotic NaCl and Na2SO4 solutions to tomato plants harboring an introduced superoxide dismutase gene. To study the roots of tomato plants cultivar Belyi Naliv (WT) and FeSOD-transgenic line, we examined the distribution of ROS and enzyme-linked immunosorbent detection of α-tubulin. In addition, longitudinal sections of the root apexes were compared. Transmission electronic microscopy of atypical cytoskeleton structures was also performed. The differences in the microtubules cortical network between WT and transgenic plants without salt stress were detected. The differences were found in the cortical network of microtubules between WT and transgenic plants in the absence of salt stress. While an ordered microtubule network was revealed in the root cells of WT tomato, no such degree of ordering was detected in transgenic line cells. The signs of microtubule disorganization in root cells of WT plants were manifested under the NaCl treatment. On the contrary, the cytoskeleton structural organization in the transgenic line cells was more ordered. Similar changes, including the cortical microtubules disorganization, possibly associated with the formation of atypical tubulin polymers as a response to salt stress caused by Na2SO4 treatment, were also observed. Changes in cell size, due to both vacuolization and impaired cell expansion in columella zone and cap initials, were responsible for the root tip tissue modification.

Transcriptome and Co-Expression Network Analyses Identify the Molecular Signatures Underlying Drought Resistance in Yellowhorn

Drought is a key factor that limits plant growth and yield. Yellowhorn is an important and vigorously promoted oil tree in China. It can survive under certain drought conditions, but a lack of water severely restricts its growth and results in yield losses in arid and semi-arid areas. Therefore, it is important to identify the key pathways and genes to understand the mechanisms of its drought resistance. In this study, we evaluated drought resistance in four types of yellowhorn, and obtained 2669 and 2451 differentially expressed genes (DEGs) via the transcriptome analysis of the comparison of water-saving/water-consuming and fast-growing/slow-growing yellowhorn, respectively, under long-term drought conditions. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of DEGs showed the key biological processes and metabolic pathways involved in drought resistance, which demonstrated that there are both the same and different biological processes involved in regulating water use efficiency (WUE) and growth in response to drought stress. Furthermore, the network analysis indicated hub genes (especially seven protein kinases) and potential co-expressed gene clusters in a greenyellow module associated with WUE and a blue module associated with growth. These identified hub genes and key biological processes can significantly enhance our knowledge about the molecular mechanisms of drought resistance in yellowhorn.