TaPYL4, an ABA receptor gene of wheat, positively regulates plant drought adaptation through modulating the osmotic stress-associated processes

doi:10.21203/rs.3.rs-1741352/v1

Download PDF

Research Article

TaPYL4, an ABA receptor gene of wheat, positively regulates plant drought adaptation through modulating the osmotic stress-associated processes

https://doi.org/10.21203/rs.3.rs-1741352/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: Abscisic acid receptors (ABR) involve transduction of the ABA signaling in plants, impacting largely on stress-defensive physiological processes and plant osmotic stress response. In this study, we characterized TaPYL4, a gene of ABR family in T. aestivum, in mediating plant drought tolerance given scarcity of functional characterization on wheat ABR members thus far.

Results: TaPYL4 harbors nine conserved domains shared by its PYL counterparts, targeting onto plasma membrane and nucleus after endoplasmic reticulum assortment. TaPYL4 interacts with TaPP2C2 whereas the latter with TaSnRK2.1, which establish a core module of the ABA signaling pathway. TaPYL4 expression was upregulated in root and aerial tissues upon drought stress. Overexpressing TaPYL4 conferred plants improved growth traits whereas knockdown expression of target gene alleviated growth feature compared with wild type under drought treatment. The TaPYL4-enhanced drought adaptation associates gene function in positively regulating stomata movement, osmolyte biosynthesis, and root system architecture (RSA) establishment. Transgene analysis on the P5CS family genes involving proline biosynthesis indicated that TaP5CS1 exerts critical roles in promoting osmolytes accumulation in drought-challenged TaPYL4 lines. TaPIN9, a PIN-FORMED gene modulating cellular auxin translocation, was validated to function as a crucial mediator in defining RSA establishment underlying TaPYL4 regulation. Transcriptome analysis revealed that TaPYL4 controls transcription of numerous genes, which impact on physiological processes associated with ‘biological process’, ‘molecular component’, and ‘cellular process’. Moreover, the differentially expressed genes mediated by TaPYL4 were closely related to jasmonic acid (JA) signaling pathways

Conclusions: Our investigation suggested that TaPYL4 acts as a positive regulator in plant drought tolerance and a valuable target for engineering drought-tolerant cultivars in T. aestivum.

Drought stress negatively regulates the plant growth and development by inhibiting various biochemical and physiological processes [1]. On the other hand, plants have evolved a suite of effective mechanisms to acclimate them to adverse effects initiated by drought [2]. Upon perception of drought signaling, the signal transduction systems in root tissue transduce the stress cues across whole plant through distinct signaling transduction pathways [3, 4]. Thus far, a large set of signaling members functional in the plant drought signal transduction system have been documented [5–7], which contribute to the drought adaptation of plants by regulating transcription efficiency of genes in families of transcription factor (TF) [8], ion channel and transporter [9], and antioxidant enzyme [10], etc., mainly through the abscisic acid (ABA)-dependent pathways [11].

The ABA molecule modulates diverse physiological processes associated with plant growth, development and stress responses. Upon osmotic stresses, the ABA contents in root tissue are swiftly elevated due to the altered transcription of genes regulating the ABA biosynthesis [12]. The ABA molecules induced in the tissue were further transported and distributed throughout various organs mediated by corresponding cellular import transporters, finally eliciting the plant response and adaptation to distinct osmotic stressors [12].

Perception of the ABA signaling upon abiotic stresses is mediated mainly by the ABA-dependent pathway. Of which, the ABRs are functional as the critical core components. In plant species, ABA receptors are present as a family of proteins, possessing the ability to interact with their partners referred to as phosphatases in clade A PP2C. Thus far, the members of ABR family, such as pyrabactin resistance (PYR) 1, PYR1-like, and proteins in the PYL family, have been characterized in the model plant A. thaliana [13, 14, 15]. Protein structural analysis on them has defined the ternary crystal features among ABRs, the substrate ABA, and the PP2C proteins [16, 17]. These findings provide valuable information in elucidating the mechanisms as to the ABA signaling perception for osmotic stresses in plant species.

To date, the molecular processes underlying the ABA-mediated plant drought response has been explored in depth given powerful research approach and valuable mutant materials. Under osmotic stress conditions (i.e., drought), the ABRs (PYR/PYL proteins) bound to ABA molecules induced and then interact with PP2Cs, removing the inhibition of PP2C activities on the kinase proteins of SnRK2 family. The activation of the SnRK2 family members phosphorylates/ activates the transcription factors (TFs), membrane channels, and transcription of the stress-defensive genes essential in plant drought tolerance [9]. Therefore, further characterization on the module ABR/PP2C/SnRK2 of the ABA signaling pathway can deepen understanding the mechanisms of plant osmotic stress response underlying the ABA-dependent pathways.

Currently, the ABR family members have been extensively investigated in various plant species, especially in model plant A. thaliana [18]. Moreover, the components constituting the core functional module of ABA signaling pathway, such as the partners interacting with ABRs (i.e., PP2C proteins) and those with PP2C proteins (SnRK2 members) have also been characterized in Arabidopis [19]. In T. aestivum, a set of of genes encoding proteins of the PYR/PYL families have been documented [20]. However, detailed characterizations on the components of the core ABA signaling module are still to be determined. Additionally, the mechanisms underlying the ABA receptor-mediated plant drought tolerance are also largely to be elusive in wheat. In this study, we characterized TaPYL4 (TaPYL4-4A), a member of the ABR family in T. aestivum, in mediating the plant response to drought stress. Our findings suggest that TaPYL4 acts as one of critical mediators in plant drought adaptation by improving osmotic stress-associated physiological processes. This ABR gene also regulates the transcription of genes at global level functional in modulating various physiological processes mainly through the jasmonic acid signaling pathway.

Characterization of TaPYL4

TaPYL4 has a full length cDNA of 1112 bp long, encoding a 179-aa polypeptide with a molecular mass of 14.93 kDa and an isoelectric point (pI) of 5.10. The TaPYL4 protein harbors nine conserved domains same as its plant counterparts (i.e., CL1 to CL9), which are involved in the binding to ABA molecule and in the interaction with its the downstream partners, such as PP2C proteins (Fig. 1A). At nucleic acid level, TaPYL4 shares high similarities to its homologous genes distributed in various plant species (Additional file 1), suggesting its nature to be one of the PYL family members in T. aestivum. Based on an experiment to define the sub-cellular position of target-GFP fusion at cellular level, the signals derived from TaPYL4-GFP in epidermal cells of N. benthamiana were confined onto plasma membrane and nucleus (Fig. 1B, Additional file 2). These results suggested that TaPYL4 targets onto both above locations after ER assortment where exerts distinct biological roles.

The Components Constitutes Aba Signaling Module With Tapyl4

Based on yeast two-hybridization assays, the component in clade A class PP2C family interacting with TaPYL4 was identified. Results indicated that TaPYL4 specifically interacted with the PP2C member TaPP2C2 (Fig. 2A). Similar assays were performed to determine the partner to be involved in interacting with TaPP2C2. It was shown that TaSnRK2.1, a member of the SnRK2 family in T. aestivum, specifically interacted with this wheat PP2C protein mentioned above (Fig. 2B). Therefore, the results in our protein-protein interaction analysis indicated that TaPYL4 constitutes a core ABA signaling module with components of TaPP2C2 and TaSnRK2.1, namely, TaPYL4-TaPP2C2-TaSnRK2.1. This module was speculated to play crucial roles in the ABA signal pathways through transduction of signaling initiated by distinct internal or environmental cues through protein phosphorylation mechanism.

Expression Patterns of TaPYL4 under Drought Stress Conditions

The transcripts of TaPYL4 in tissues of root and aerial tissues were evaluated to characterize the expression patterns of target gene in response to drought stress. Upon drought stress treatment, the expression levels of TaPYL4 in both tissues were significantly upregulated following the intensified extent of drought stressor, showing to be gradually elevated along with the increased polyethylene glycol (PEG) in growth media (0 to 15% (w/v) PEG-6000) (Fig. 3A). Additionally, under drought stress condition (10% PEG-6000), the TaPYL4 expression in tissues displayed a temporal-dependent pattern in response to drought stressor within a 27 h regime, reaching peak level at end of the treatment. Moreover, the induced transcripts of target gene in tissues under drought stress condition were gradually recovered following a 27 h of normal recovery condition (Fig. 3B). These results together suggested that TaPYL4 sensitively responds to drought stress at transcriptional level, mediating plant drought response through its spatiotemporal mode of transcription in response to drought stress.

Growth Properties of TaPYL4 Transgenic Lines under Drought Treatment

Sen 2 and Sen 3, two transgenic lines at T3 generation with more TaPYL4 transcripts and Anti 1 and Anti 2, two lines with significant repression of target expression with one copy inserted (Additional file 3), were selected to define the gene function in mediating plant drought tolerance. Under normal growth condition, all of transgenic lines (i.e., Sen 2, Sen 3, Anti 1, and Anti 2) were comparable on phenotypes, biomass of aerial tissue and root, and root volumes with the WT plants (Figs. 4A to 4E). Under drought stress treatment, however, the transgenic lines were shown to be dramatically modified on growth traits mentioned with respect to WT. Of which, Sen 2 and Sen 3 were much more improved on phenotypes (Fig. 4A), biomass of aerial tissue (Fig. 4B), root biomass (Fig. 4C), plant biomass (Fig. 4D), and root volumes (Fig. 4E) than WT plants. In contrast, compared with WT, Anti 1 and Anti 2 drastically alleviated phenotypes (Fig. 4A), biomass of aerial tissue (Fig. 4B), root biomass (Fig. 4C), plant biomass (Fig. 4D), and root volumes (Fig. 4E) in plants treated by drought stress. The significantly modulated growth behavior was shown in the drought-challenged transgenic lines, suggesting that TaPYL4 plays an important role in mediating plant adaptation to drought stress.

Stress Response-associated Physiological Traits of TaPYL4 Transgenic Lines

A suite of stress-associated traits and photosynthetic parameters were measured in transgenic lines (Sen 2, Sen 3, Anti 1, and Anti 2) under drought stress conditions to understand the physiological processes mediated by TaPYL4. In line with the growth traits described above, all the transgenic lines were comparable on the stress-associated traits (i.e., stomata closing rate (SCR), leaf water losing rate (WSR), and contents of proline and soluble sugar) (Figs. 5A-5E), and the photosynthetic parameters (i.e., photosynthetic rate (Pn), stomatal conductance (gs), ΦPSII, and NPQ) (Figs. 5F-5I) with WT plants under normal growth condition. Under drought stress treatment, compared with WT, Sen 2 and Sen 3 displayed promoted SCR (Figs. 5A-5B), slowed WLR (Fig. 5C), increased contents of osmolytes (proline and soluble sugar) (Figs. 5D-5E), and enhanced Pn (Fig. 5F), gs (Fig. 5G), and ΦPSII (Fig. 5H), and reduced NPQ (Fig. 5I). In contrast, Anti 1 and Anti 2 drastically alleviated the stress-associated traits and the photosynthetic parameters, displaying lower SCR (Figs. 5A-5B), higher elevated leaf WLR (Fig. 5C), less contents of osmolytes (Figs. 5D-5E), and more alleviated photosynthetic function (lower Pn, gs, and ΦPSII, and higher NPQ values) than the WT plants (Figs. 5F-5I). The physiological traits associated with plant stress response were in concordance with growth features of the transgenic lines under drought treatment, suggesting that the TaPYL4-mediated drought tolerance of plants ascribes partly to the gene function in improving the stress responsive-associated physiological processes.

Expression Patterns Of P5cs And Pin-formed Family Genes

Expression patterns of the genes in delta-1-pyrroline-5-carboxylate synthetase (P5CS) family impacting proline biosynthesis and in PIN-FORMED (PIN) family involving root system architecture (RAS) establishment were investigated, with an aim to understand the molecular processes related to the TaPYL4-mediated osmolyte accumulation and RSA property under drought conditions. Among five P5CS family genes examined, TaP5CS1 was modified significantly in transcription in the drought-challenged TaPYL4 transgenic lines, displaying higher expression levels in Sen 2 and Sen 3 whereas lower ones in Anti 1 and Anti 2 compared with the WT plants (Fig. 6A). Among six genes in the PIN family, similar to TaP5CS1, TaPIN9 displayed modified transcripts in the TaPYL4 transgenic lines with respect to WT under drought treatment, with more transcripts shown in Sen 2 and Sen 3 whereas less expression levels in Anti 1 and Anti 2 than the WT plants (Fig. 6B). These genes in P5CS and PIN families (i.e., TaP5CS1 and TaPIN9) modified significantly on expression efficiency in TaPYL4 transgenic lines under drought treatment, contrasting to other genes in P5CS and PIN families examined that unaltered transcripts abundance among the transgenic and WT plants. Therefore, it is suggested that the TaPYL4-improved drought adaptation is associated with the modified transcription of distinct genes in P5CS and PIN-FORMED families that regulate osmolyte accumulation and RSA establishment.

Functions of TaP5CS1 and TaPIN9 in Mediating Plant Drought Response

TaP5CS1 and TaPIN9, two members in P5CS and PIN families that modified transcription underlying TaPYL4 regulation, were subjected to transgene analyses to address their functions in mediating plant drought stress response. Results indicated that the lines with significant knockdown expression of TaP5CS1, namely, AntiP5CS1-1 and AntiP5CS1-2 (Additional file 4), were significantly alleviated on phenotypes (Fig. 7A), proline accumulation in plants (Fig. 7B), and plant dry mass production (Fig. 7C) relative to WT under drought treatment. Likewise, the lines with drastic knockdown expression of TaPIN9 (i.e., AntiPIN9-2 and AntiPIN9-3) (Additional file 5) displayed significant modification on RSA establishment and plant drought response. Compared with WT, AntiPIN9-2 and AntiPIN9-3 drastically alleviated root growth traits and plant biomass under drought treatment, showing smaller stature of plants (Fig. 7A), less dry mass accumulation in whole plant (Fig. 7D) and in roots (Fig. 7E), and lower root volume (Fig. 7F) than the wild type. Therefore, the transgene results on TaP5CS1 and TaPIN9 validated the gene functions in promoting osmolyte accumulation and in improving RSA establishment underlying TaPYL4 regulation, respectively. These genes act as the crucial modulators in plant adaptation to drought stress via enhancement of osmolyte-regulatory capacity and improvement of RSA establishment modulated by the core ABA signaling module constituting TaPYL4 and its partners as aforementioned.

Transcriptome Profile Mediated by TaPYL4 under Drought Stress Conditions

Transcriptome profiles in the drought-challenged TaPYL4 transgenic line (Sen 2) and WT were investigated based on high throughput RNA-seq analyses to systematically understand the molecular processes underlying TaPYL4 modulation. Results indicated that in total of 3850 genes, including 2613 to be upregulated and 1237 downregulated, were shown to be differentially expressed (DE) in the transgenic line after drought stress treatment (Additional files 6–7). To be sure of reproducibility for the transcriptome results, ten DE genes including five with upregulated whereas another five with downregulated in expression were selected and subjected to qRT-PCR analysis. As expected, all of the five genes with upregulated pattern (i.e., TaWRK2, TaWRKY28, TaCML31, TaMPK18, and TaZFP1) displayed higher expression levels in Sen 2 with respect to WT, with comparable folds increased as shown in RNA-seq analysis. Likewise, the five DE genes with downregulated in expression (i.e., TaPAO, TaCA, TaCP450, TaUBI6, and TaFR1) exhibited decreased transcripts in Sen 2 relative to WT plant under drought condition (Additional file 8). Moreover, these DE genes in Anti 1 exhibited reverse expression pattern to tha shown in Sen 2 mentioned above, namely, the five upregulated DE genes were significantly decreased whereas the five downregulated DE ones were drastically elevated on expression levels compared with the WT plants (Additional file 9). These qRT-PCR results confirmed the credible for the transcriptome profile underlying modulation of TaPYL4 upon drought stress.

Based on the gene ontology (GO) analysis on the DE genes, the DE genes they were categorized into the GO terms associated with “biological process”, “molecular components”, and “cellular process”. Among these, the DE genes in the GO term “biological process” are related to processes of cation binding, metal ion binding, calcium ion binding, transcription regulator activity, oxidoreductase activity, acting as single donors with incorporation of molecular, oxygen, transcription coregulator activity, transcription corepressor activity, all-trans-beta-apo-10'-carotenal cleavage oxygenase activity, phenylalanine ammonia-lyaseactivity, ammonia-lyase activity, ferrous iron binding, and isocitrate lyase activity; the DE genes in GO term ‘molecular component’ translate proteins constituting external encapsulating structure, cell wall, and plant-type cell wall; the DE genes in GO term‘cellular process’ associate with cellular functions of aromatic compound metabolic process, organic cyclic compound biosynthetic process, aromatic compound biosynthetic process, heterocycle biosynthetic process, nucleobase-containing compound biosynthetic process, nucleicacid-templated transcription, RNA biosynthetic process, negative regulation of biological process, negative regulation of cellular process, regulation of response to stress, negative regulation of nucleobase-containing compound metabolic process, negative regulation of cellular biosynthetic process, regulation of signaling, regulation of cell communication, regulation of signal transduction, negative regulation of biosynthetic process, negative regulation of nucleic acid-templated transcription, negative regulation of RNA biosynthetic process, negative regulation of RNA metabolic process, negative regulation of macromolecule biosynthetic process, response to wounding, regulation of defense response, protein polyubiquitination, response to jasmonic acid, regulation of jasmonic acid mediated signaling pathway, cellular response to jasmonic acid stimulus, jasmonic acid mediated signaling pathway, negative regulation of mitotic sister chromatid separation, negative regulation of chromosome separation, erythrose 4-phosphate/phosphoenolpyruvate family amino acid catabolic process, negative regulation of metaphase/anaphase transition of cell cycle, negative regulation of mitotic cell cycle phase transition, negative regulation of cell cycle phase transition, hydrocarbon catabolic process, mitotic spindle checkpoint, spindle assembly checkpoint, negative regulation of chromosome segregation, glyoxylate metabolic process, terpene catabolic process, negative regulation of mitotic cell cycle, negative regulation of mitotic metaphase/anaphase transition, negative regulation of mitotic nuclear division, negative regulation of mitotic sister chromatid segregation, negative regulation of sister chromatid segregation, and spindle checkpoint (Fig. 8A). KEGG analysis revealed that the DE genes identified in the drought-challenged TaPYL4 overexpression lines suggests the a large set of the enriched biochemical pathways underlying control of TaPYL4, which were overrepresented by the biochemical pathways associated with aminoacyl-tRNA biosynthesis, glycosamino glycan degradation, biotin metabolism, lipoic acid metabolism, sphingolipid metabolism, C-type lectin receptor signaling pathway, inositol phosphate metabolism, amino sugar and nucleotide sugar metabolism, sulfur relay system, ether lipid metabolism, beta-Alanine metabolism, ubiquinone and other terpenoid-quinone biosynthesis, phosphatidylinositol signaling system, tropane, piperidine and pyridine alkaloid biosynthesis, Arachidonic acid metabolism, Phenylalanine, tyrosine and tryptophan biosynthesis, butanoate metabolism, phosphonate and phosphinate metabolism, lysine biosynthesis, tyrosine metabolism, monobactam biosynthesis, circadian rhythm-plant, fatty acid biosynthesis, ascorbate and aldarate metabolism, glutathione metabolism, phagosome, citrate cycle (TCA cycle), fructose and mannose metabolism, riboflavin metabolism, sulfur metabolism, arginine biosynthesis, thiamine metabolism, necroptosis, photosynthesis-antenna proteins, glycerophospholipid metabolism, pyruvate metabolism, folate biosynthesis, starch and sucrose metabolism, vitamin B6 metabolism, valine, leucine and isoleucine degradation, glycine, serine and threonine metabolism, carotenoid biosynthesis, pentose phosphate pathway, alanine, aspartate and glutamate metabolism, and glycolysis/gluconeogenesis, etc. (Fig. 8B). These results from the transcriptome profiles mediated by TaPYL4 suggested that this ABR gene globally modulates the transcription of quantities of genes, whose modified expression levels lead to plant drought response through regulating diverse stress responsive-associated physiological processes and biochemical pathways, such as those related to stomata movement, osmolytes biosynthesis, and RSA establishment.

As an essential phytohormone in plant species, abscisic acid (ABA) regulates diverse biological processes associated with plant growth and development as well as and responses to a large set of environmental stresses, including the osmotic stressor [12]. Thus far, genetic and biochemical analyses have identified a large set of ABR (PYR/PYL) family members across various plant species, especially the model plant Arabidopsis [21]. Moreover, a subset of the PYL proteins acting as key regulators in ABA signaling have been functionally characterized [22]. In this study, we performed molecular characterization analysis on TaPYL4 and revealed that it shares high similarities to its PYL counterparts distributed in diverse plant species. TaPYL4 also harbors nine conserved domains (CL1 to CL9) involving binding ABA and interacting with its partners (i.e., PP2Cs). Therefore, TaPYL4 acts as one of the PYL family members in T. aestivum. Our experimental analysis in determination of the subcellular localization of TaPYL4 indicated its position onto plasma membrane and nucleus after ER assortment, suggesting where it percepts the ABA signaling and exerts biological roles.

Transcripts of the genes in PYL family were detected in various plant species, such as rice [23], tomato [24, 25], cotton [26], and the rubber tree [27]. However, although the expression of PYL genes is shown in various tissues, such as seeds, tissues root and leaf, and organs flower and fruit, much more of the PYL transcripts are specified in distinct organs/tissues. For example, rice OsPYL7 and OsPYL8 are abundantly expressed in embryos, OsPYL3 and OsPYL5 primarily express in leaves, while OsPYL1 predominantly expresses in roots [28]. Moreover, upon osmotic stress, a subset of PYR/PYL/RCAR genes induces expression levels due to enhanced cellular ABA levels in plants [29, 30]. In this study, we evaluated the TaPYL4 expression patterns in root and aerial tissues upon drought stress. TaPYL4 transcripts in tissues were enhanced by drought stressor and its upregulated expression initiated by drought was restored to low levels following progression of normal recovery condition. These expression results suggested the sensitive response of TaPYL4 to external drought signaling. Thus far, characterization on cis-acting regulatory elements in the promoters of drought stress-responsive genes has indicated that CRE, which is specified by motifs CCGCGT, CACGTGT, and AAGTCAA, frequently situates in promoter region and regulates gene transcription in responding to osmotic stresses. This element leads to promoted gene transcription efficiencies by interacting with DNA binding domains harboring in proteins of the bZIP transcription factors [31]. Therefore, TaPYL4 modified expression upon drought stress associates distinct cis-acting regulatory elements in its promoter, such as CRE or other uncharacterized ones. By searching TaPYL4 promoters using online tool PLACE, we identified a set of elements involving stress responses, such as CRE, ABRE, and DRE/CRT, etc. Further characterizing these elements can provide insights into the transcription mechanisms of the PYL family members in response to drought signaling.

The PYL family members mediate extensively plant growth and osmotic stress responses. For example, transgenic lines overexpressing rice OsPYL/RCAR5 are hypersensitive to ABA during seedling stage [32]. In contrast, PYL genes OsPYL/RCAR5, OsPYL3, and OsPYL9 confer plants improved drought tolerance [33]. Elevation of AtPYL5 and AtPYL9 abundance enhances plant ABA responses and drought tolerance through interacting with a set of clade A PP2C proteins [17]. These findings together suggested that the PYL genes act as crucial modulators in plant adaptation to drought stress. In this study, we used transgene analysis to characterize the function of TaPYL4 in mediating plant drought response. The lines with TaPYL4 overexpression (i.e., Sen 2 and Sen 3) drastically improved plant growth, RSA establishment, and biomass production whereas those with target knockdown (i.e., Anti 1 and Anti 2) alleviated significantly on above growth traits with respect to wild type under drought treatment. These transgene results confirmed the essential role of TaPYL4 in regulating plant drought tolerance. Additionally, we found that TaPYL4 overexpression did not show penalties on plant dry mass production compared with wild type under normal growth conditions, suggesting its potential value in molecular breeding drought-tolerant cultivars in T. aestivum.

Stomata movement acts as one effective mechanism for plant response to osmotic stresses via an ABA-dependent pathway. Enhanced capacity of plants for water maintenance under drought condition associates promoted stomatal closure rate (SCR) initiated by induced ABA amounts, given the mediation of ion transport systems situated in plasma membrane of guard cells that modify turgor and volume of cells [34]. In this study, we investigated the stomata closing nature upon drought signaling underlying TaPYL4 regulation. Results indicated significantly modified stomata closing rate (SCR) in TaPYL4 transgenic lines compared with WT. The SCR behaviors in Sen 2 and Sen 3 were promoted whereas Anti 1 and Anti 2 slowed within a 1 h-regime of drought duration. These SCR analysis results were in concordance with the positive role of TaPYL4 in regulating water retention capacity of the drought-treated plants. Previously, anion channels in guard cells, such as those categorized into S-type and R-type, were documented to be activated upon osmotic stresses, impacting stomata movement and SCR by mediating anion efflux from guard cells, depolarizethe potassium efflux channels situated in plasma membrane, and modify guard cell osmotic potential [35]. Further characterizing the mechanisms underlying TaPTL4-mediated anion efflux in guard cells can deepen understanding the stomata movement mediated by ABRs in T. aestivum. Another issue should be mentioned although the lines overexpressing TaPYL4 (i.e., Sen 2 and Sen 3) were much faster on stomatal closing upon drought signaling than wild type in our investigation, they also possessed increased photosynthetic rate (Pn). We speculate that the improved photosynthesis mediated by TaPYL4 under drought condition is largely ascribed to its improved photosynthetic system activation (i.e, improved PSII and lowered NPQ) that compensates its adverse effects of promoting stomata movement on Pn behavior. Similar finding has been documented on PeCHYR1, a ubiquitin E3 ligase gene in P. euphratica, whose up-regulation confers plants promoted stomatal closure together with enhanced photosynthetic activity and biomass accumulation [36]. Further investigation on the relation between stomata closure and photosynthetic function mediated by TaPYL4 can provide insights into plant drought adaptation for cereal crops.

Distinct PYL members consist of core ABA signaling modules with PP2C and SnRK2 proteins to be functional in plant osmotic stress responses [37, 38]. Once interacted by ABRs upon osmotic stressor, PP2C proteins lead to released inhibition on the activity of SNF1-related protein kinase 2 (SnRK2) proteins [16], which activate the downstream regulators in ABA-dependent pathway, such as transcription factors, S-type anion channel SLAC1, and a large set of proteins functional in various groups thorough phosphorylation mechanism. These processes thereby lead to drought-acclimated stomata movement and plant stress responses [39]. In this study, we performed yeast two-hybridization assays to identify the core signaling module components aside from TaPYL4. Our results indicated that TaPP2C2, a clade A PP2C member and TaSnRK2.1, a member of the SnRK2 family, constitute a core ABA signaling module, namely, TaPYL4/TaPP2C2/TaSnRK2.1, given their specific protein-protein interactions. In this regard, it is speculated that the TaPYL4-mediates plant drought response to be accomplished via this ABA signaling module that initiate plant drought stress responses. Further functional characterization of TaPP2C2 and TaSnRK2.1 can benefit understanding of the molecular mechanisms underlying the ABA signal transduction pathways in T. aestivum.

Plant accumulation of osmolytes under osmotic stresses decreases cellular osmotic potential, enhancing cell water retention capacity and alleviating injury extent of plants initiated by stressors [40]. Proline and soluble sugar are two types of crucial osmolytes, whose induction under drought stress contribute largely to plant drought adaptation [41]. In this study, we investigated the TaPYL4-mediated modification on above two osmolytes by assessing contents of them in drought-challenged transgenic lines. As expected, we found that TaPYL4 positively regulates proline and soluble sugar accumulation; higher contents of osmolytes in Sen 2 and Sen 3 whereas lower contents of those were observed in Anti 1 and Anti 2 compared with wild type. These results suggested that TaPYL4 enhances osmolytes accumulation in plants once challenged by drought stress. Based on expression analysis, we found the transcripts of TaP5CS1, a member of Δ1-pyrroline-5-carboxylate synthetase (P5CS) family acting as rate-limiting regulator in proline biosynthesis [7, 42], significantly modified expression in transgenic lines under drought treatment. Our online search analysis revealed that the ABA responsiveness motif ABRE (ACGTG, -1274, -529, 225, 406, and 1401) and that the drought-inducible element MBS (CAACTG, -1157 and − 313) frequently situate in TaP5CS1 promoter, being consistent with its upregulated expression upon drought stress. Further transgene analysis on TaP5CS1 validated its positive role in regulating proline accumulation under drought conditions. Therefore, TaPYL4 constitutes a putative functional module with TaP5CS1, namely, TaPYL4/TaP5CS1, to contribute to plant drought adaptation through modulating osmolyte biosynthesis. To date, a subset of drought-responsive proteins including those involving osmoprotectant biosynthesis metabolisms has been document [7]. Further addressing the transcriptional mechanisms of TaP5CS1 mediated by TaPYL4 can help elucidate the osmotic-regulatory mechanism underlying PYL members in drought response in T. aestivum.

Auxin modulates widely plant physiological processes associated with morphogenesis, organogenesis, apical dominance, embryo formation, vascular differentiation, and stress responses [25]. The subcellular localization of auxin possesses polar nature in root tissue, which is mediated largely by the PIN-FORMED (PIN) proteins given their role in determining directional flow and polar transport of auxin at cellular level [43]. Distinct PIN proteins also mediate polar re-allocation of auxin from tissue shoot to root tip [44], by which to regulate lateral root differentiation [45] and root system architecture (RSA) establishment [46]. In this study, we addressed the issue whether TaPYL4-mediated RSA establishment to contribute to plant drought tolerance by investigating expression patterns of the PIN family genes in drought-challenged TaPYL4 transgenic lines. We found that TaPIN9 is modified on transcription, with more transcripts in Sen 2 and Sen 3 and less ones in Anti 1 and Anti 2 than WT plants. Online search analysis revealed that the elements associated with drought response, such as ABRE (-241, -630, -1682, and 1770) and MBS (-1472 and 942), are situated in TaPIN9 promoter region, suggesting the involvement of them in mediating gene transcription under drought stress condition. Further tangsgene analysis confirmed the positive function of TaPIN9 in regulating RSA establishment and plant drymass production under drought treatment. Therefore, our investigation suggested that TaPYL4 modulates RSA establishment by transcriptionally regulating TaPIN9 that impacts cellular auxin transportation. Previously, the clade A PP2A phosphatases were documented to mediate cellular auxin transport through a reversible phosphorylation mechanism [47]. This finding suggests the putative functional pathways constituted by PYL/PP2C/PIN to regulate auxin transport at cellular level and to impact on RSA behavior of plants challenged by drought stress.

ABA signaling initiated by osmotic stressors modifies a quantity of physiological and biochemical processes which are corporately integrated in plant drought responses [11, 48]. In this study, we performed high throughput RNA-seq analyses to identify the TaPYL4-mediated transcriptome profile under drought stress condition. We identified a large set of differentially expressed (DE) genes in Sen 2 compared with wild type after drought treatment. These DE genes are enriched into the GO terms associated with “biological process”, “molecular components”, and “cellular process” and overrepresented by KEGG pathways related to diverse biochemical pathway, suggesting the complex nature in TaPYL4-mediated plant drought response. Based on GO terms and KEGG pathways enrich by the DE genes, we provide the working model of TaPYL4 in regulating plant drought responses (Fig. 9). TaPYL4 mediates plant adaptation to drought stress through a suite of biochemical pathways, namely, ‘metabolic process (GO:0008152)’, ‘biological regulation (GO:0065007)’, ‘cellular process (GO:0009987)’, and ‘response to stimulus (GO:0050896)’. Of which, pathway ‘metabolic process’ leads to modifications of cellular biosynthesis (GO:0009058, GO:0044237, GO:0044249, and GO:0006725), which finally modifies the metabolism of aromatic compound biosynthesis; the pathway ‘biological regulation process’ impacts physiological processes associated with the regulation of biology and signaling (GO:0050789; GO:0023052; GO:0023051, and 0048583), which improves plant acclimation to stress signal (GO:0009966); the pathway ‘cellular process’ exerts roles in modulating cellular response to stimulus (GO:0051716), signal transduction (GO:0007165), which impact on plant hormone-mediated signaling (GO:0009755); and the pathway ‘response to stimulus’ is involved in modulating cellular response to stress (GO:0051716), which adjusts plant responses to a set of signals (GO:0009719, GO:0010033, GO:0070887, and GO:1901700) and enhance plant drought adaptation mainly via the jasmonic acid signaling pathway (GO:0009753, GO:0009867, and GO:2000022). These metabolism- and biosynthesis-associated biological processes modulated by the DE genes underlying TaPYL4 are comprehensively integrated in plant drought tolerance, through modulating stomata movement, osmolytes biosynthesis, RSA establishment, and photosynthesis behaviors.

Wheat ABR gene TaPYL4 shares high similarities to its homologous genes across various plant species. The expression levels of TaPYL4 in root and aerial tissues are sensitive in response to drought stress. TaPYL4 targets onto plasma membrane and nucleus where establish a putative core ABA signaling module TaPYL4/TaPP2C2/TaSnRK2.1. Overexpression of TaPYL4 confers plants significantly improved drought tolerance, due to gene function in positively regulating stomata movement, osmolytes biosynthesis, and RSA establishment. TaP5CS1, a gene in the P5CS family and TaPIN9, a member of the PIN-FORMED family, are transcriptionally regulated by TaPYL4, contribute to plant drought adaptation by promoting proline accumulation and RSA formation, respectively. TaPYL4 is also involved in the regulation of transcription efficiencies of numerous genes, which are functional in various biological processes associated with plant drought stress defensiveness mainly through the jasmonic acid signaling pathway.

Characterization Analysis on TaPYL4

TaPYL4 (GenBank accession No. XM_044507276) was significantly upregulated in expression upon drought stress in our previous RNA-seq analysis (unpublished data). This finding prompted us to characterize it in mediating the plant drought response in more detail. This ABR gene has three homeologs across chromosomes 4A (TaPYL4-4A, TraesCS4A02G114400), 4B (TaPYL4-4B, TraesCS4B02G189800), and 4D (TaPYL4-4D, TraesCS4D02G191200), which share high similarities at nucleic acid and amino acid levels each other (98.5–99.2%, Supplemental Figure S1A-S1B). Of which, TaPYL4-4A displayed strongly induced transcripts upon drought stress relative to other homeologs (Additional file 10). Therefore, we concentrated on characterization of TaPYL4-4A in this investigation (simplified by TaPYL4 hereafter). The conserved domains harbored in the TaPYL4 protein were defined base on previous literature [37]. The homologous genes of TaPYL4 in plant species were identified by sequence search analysis against the GenBank database for National Center in Biotechnology Information (NCBI). Phylogenetic relations among TaPYL4 and its plant counterparts were established using the MegAlign algorithm supplemented in the DNAStar software.

Determination Of Tapyl4 Subcellular Localization

The subcellular localization of TaPYL4 was defined based on detecting the signals of green fluorescent protein (GFP) initiated by fusion GFP-TaPYL4 at cellular level. With this purpose, the open reading frame (ORF) of TaPYL4 was amplified using RT-PCR together with gene specific primers (Additional file 11). It was then integrated into the binary vector pCAMBIA3300 in frame with the reporter gene GFP downstream of the CaMV35S promoter. The expression cassette harboring GFP-TaPYL4 was subjected to genetic transformation onto A. tumefaciens strain (EHA105) following the conventional heat-shock transformation approach, and then transformed epidermal cells of N. benthamiana through an A. tumefaciens-mediated transformation method. The GFP signals derived from the epidermal cells transformed with TaPYL4-GFP and those from control harboring empty vector (with GFP as reporter gene) were observed under fluorescent microscope after 48 h of infiltration as described previously [49]. Additionally, we used protoplast expression system of N. benthamiana to define the location of TaPYL4 at cellular level. With this aim, young leaves N. benthamiana were obtained following the procedure of previously described [50]. Transformation on isolated protoplasts using the cassettes harboring TaPYL4-GFP was conducted following a modified PEG-mediated protocol [51]. The signals in protoplasts initiated by the fusion were observed under fluorescent microscope after 24 h of infiltration with A. tumefaciens transformants.

Yeast Two-hybrid Assays

Yeast two-hybrid assays were performed to identify the partner interacting with TaPYL4 as well as that involving interaction with the latter, with an aim to identify the putative core ABA signaling module covering TaPYL4. In this context, six genes in the PP2C family and seven ones in the SnRK2 family in T. aestivum were identified from the GenBank database in NCBI. The members of PP2C family subjected to assays included TaPP2C1-1A, TaPP2C2-1A, TaPP2C3-1B, TaPP2C4-1D, TaPP2C5-1A, and TaPP2C6-3B (simplified by TaPP2C1 to TaPP2C6). The genes of SnRK2 family subjected to assays contained TaSnRK2.1-2A, TaSnRK2.2-2A, TaSnRK2.3-1A, TaSnRK2.4-3A, TaSnRK2.5-2A, TaSnRK2.6-2A, and TaSnRK2.7-1A (simplified by TaSnRK2.1-TaSnRK2.7). Information of the genes in PP2C and SnRK2 families subjected to yeast two-hybrid assays is shown in Additional file 11. During assays, TaPYL4 was firstly used as a bait whereas the proteins in PP2C family as preys separately to identify target PP2C partners interacted by the PYL member. Likewise, to determine the SnRK2 members involving the constitution of core ABA signaling module, the PP2C member identified was used as bait whereas the proteins in SnRK2 family as preys separately. Procedures for these yeast two-hybrid assays were performed following the manufacturer’s instructions [52]. Gene specific primers used for the amplification of TaPYL4 and the genes in PP2C and SnRK2 families mentioned are shown in Additional file 11.

Expression Analysis of TaPYL4 upon Drought Stress

The expression patterns of TaPYL4 upon drought stress were evaluated using Shimai 22, a cultivar to be drought-tolerant (kindly provided for research purpose by the Wheat Breeding Institute, Shijiazhuang Research Academy of Agriculture and Forestry). Briefly, the seeds were regularly germinated in a growth chamber and the young seedlings were cultured in a standard MS solution as described previously [53]. At the second-leaf stage, they were subjected to simulated drought treatment by growing in a modified MS solution containing different concentrations of polyethylene glycol 6000 (PEG-6000), namely, 0, 1, 5, 10, and 15% (w/v) which correspond to osmotic potentials of 0, -0.35, -0.81, -1.08, and − 1.27 mPa, respectively. The roots and leaves of seedlings were sampled after 27 h of the treatments. To understand the temporal effects of drought stressor on gene expression, tissues mentioned above were collected at 1, 3, 9 and 27 h after drought stress treatment (10% of PEG). Additionally, aliquots of seedlings treated by 27 h of 10% PEG were re-subjected to regular growth condition by growing in standard MS solution to define the gene response to normally recovered conditions. With this purpose, tissues mentioned were sampled at 1, 3, 9, and 27 h following the recovery treatment. Transcripts of the target gene in collected samples were evaluated based on qRT-PCR performed as described previously [49], using gene specific primers (Additional file 10). A constitutive gene Tatubulin in T. aestivum was used as an internal reference to normalize the target transcripts.

Assays for Growth Traits of TaPYL4 Transgenic Lines

A suite of transgenic lines with modified expression of TaPYL4 were generated to characterize the target gene function in mediating plant drought response. With this aim, the ORF of TaPYL4 was amplified in both sense and anti-sense orientations based on RT-PCR using gene specific primers (Additional file 11). They were then separately inserted into the NcoI/BastEII sites in binary vector pCAMBIA3301 under control of the CaMV35S promoter. The expression cassettes harboring target gene were subjected to transformation onto A. tumefaciens (strain EHA105) using the conventional heat-shock approach, which were further subjected to genetic transformation onto T. aestivum (cv. Shimai 22, a wheat cultivar kindly provided for research purpose by the Wheat Breeding Institute, Shijiazhuang Research Academy of Agriculture and Forestry) based on an A. tumefaciens–mediated transformation method as described previously [49]. The expression levels of target gene in transgenic lines, namely, six lines with target overexpression (Sen 1 to Sen 6) and five ones with knockdown expression of the target gene (Anti 1 to Anti 5), were evaluated based on qRT-PCR using gene specific primers (Additional file 11).

Sen 2 and Sen 3, two lines at T3 generation with more target transcripts and Anti 1 and Anti 2, two lines with lowered TaPYL4 expression with one copy insertion each (Additional file 3), were selected to address the gene function in mediating plant drought response. To this end, these transgenic lines together with wild type (WT) were cultured in pots filled by mixed soil (half vermiculate and another half fertile soil) under normal growth condition (supplied daily by tap water to sustain 70–75% relative soil moisture) to the second-leaf stage. Afterwards, they were then subjected to drought treatment by sustaining 55–60% of relative soil moisture, with water amount applied daily monitored by soil water potential meter (TRS-IIN, China). Three weeks after the treatments, the growth traits in transgenic lines including phenotypes, biomass in aerial tissues and roots, and root volumes were assessed. Of which, phenotypes were shown as images taken by a digital camera; biomass of plant tissues and root volumes mentioned was obtained following the conventional assay approach.

Assays of Drought Stress-associated Physiological Traits and Photosynthetic Parameters

The drought response-associated traits and photosynthetic parameters were evaluated in the transgenic lines after drought treatment to understand the TaPYL4-mediated stress responsive-associated processes. The drought response-associated traits assessed included stomata closing rate (SCR), leaf water losing rate (WLR), and osmolytes contents. Of which, SCR was determined based on recordation of stomata aperture at indicated time points (i.e., 0, 0.5, 1, and 2 h) following the drought treatment as previously [54]; WLR values were calculated according to the decreased fresh weights in detached leaves at indicated time points (0.25, 0.5, 1, 2, and 3 h) following drought stressor with respect to that at the time point 0 h (prior to treatment) as described by Imadi et al. (2015) [55]; the contents of osmolytes, including proline and soluble sugar, were assessed according to the methods reported by Liu et al. (2020) [56]. The photosynthetic parameters assayed in the drought-challenged transgenic lines included photosynthetic rate (Pn), stomatal conductance (gs), photosystem II photochemical efficiency (ΦPSII), and non-photochemical quenching coefficient (NPQ). Of which, Pn and gs were measured in the upper fully expanded leaves using a photosynthesis system (LiCOR-6200) following the manufacturers’ instructions; ΦPSII and NPQ in the representative leaf samples were evaluated as reported previously [49].

Analyses Of Expression And Function Of P5cs And Pin-formed Family Genes

To characterize the molecular processes associated with the improved osmolyte biosynthesis and the root system architecture (RSA) establishment underlying TaPYL4 modulation, five genes in delta-1-pyrroline-5-carboxylate synthetase (P5CS) family impacting on proline biosynthesis in T. aestivum, including TaP5CS1 to TaP5CS5, were identified in the NCBI GenBank database. Likewise, six genes in PIN-FORMED (PIN) family of wheat, including TaPIN1 to TaPIN3, TaPIN5, TaPIN8, and TaPIN9, were obtained. Information of the genes in P5CS and PIN families mentioned are shown in Additional file 11. The transcripts of above P5CS and PIN family genes in the transgenic lines and WT after drought treatment were evaluated based on qRT-PCR using gene primer pairs (Additional file 11).

TaP5CS1 and TaPIN9, two genes showing significantly more transcripts in Sen 2 and Sen 3 whereas less ones in Anti 1 and Anti 2 compared with wild type, were subjected to functional characterizations based on transgene analysis. With this purpose, the ORFs of TaP5CS1 and TaPIN9 were amplified in antisense orientation based on RT-PCR. They were then integrated separately into the NcoI/BstEII sites in the binary vector pCAMBIA3301 under control of the CaMV35 promoter. Genetic transformation for these genes onto A. tumefaciens strain EHA105 and further genetic transformation of them onto T. aestivum (cv. Shimai 22) were performed similarly in generation of the knockdown expression lines of TaPYL4 mentioned above. Two lines with significantly downregulated TaP5CS1 expression, namely, AntiP5CS1-1 and AntiP5CS1-2 (Additional file 11) and two ones with lowered TaPIN9 transcription, AntiPIN9-2 and AntiPIN9-3 (Additional file 11), were subjected to drought treatment as mentioned previously. To this end, the knockdown lines mentioned and the WT plants were cultured under normal growth condition to the second-leaf stage and then subjected to three weeks of drought treatment as aforementioned. After that, proline contents and biomass in lines AntiP5CS1-1 and AntiP5CS1-2 were assessed. Likewise, biomass in aerial tissues and roots and root volumes in lines AntiPIN9-2 and AntiPIN9-3 were analyzed after drought treatment. The procedures for assessing above traits were similar to those conducted in the drought-challenged TaPYL4 transgenic lines.

Transcriptome Analysis

To characterize the transcriptome profile modulated by TaPYL4 under drought condition, Sen 2 that significantly upregulated expression of TaPYL4 together with WT was subjected to high throughput RNA-seq analyses after drought treatment. Briefly, total RNA in the roots of Sen 2 and WT was separately extracted using the TRIzol reagent (Invitrogen). They were then used to construct the strand-specific RNA-seq libraries in triplicates, following the approach described previously [57]. The transcripts in the libraries were sequenced using an Illumina HiSeq 2500 sequencing platform. Raw reads generated from the libraries were further processed by removing the adaptors in reads and those to be low-quality using the tool referred to as Trimmomatic [58]. The resulting ones generated were then subjected to an alignment analysis by searching against the T. aestivum reference transcript database (Novogene Co, LTd, Beijing). Differentially expressed (DE) genes in Sen 2 after drought treatment were defined for those with modified transcripts over 2-fold variations in the transgenic line compared with WT, using the edger program that is effective in accurately calculating the raw count data in libraries [58]. During which, the raw P values were defined using a false discovery rate (FDR) less than 0.05 [39]. An online tool referred to as Plant MetGenMap (http://bioinfo.bti.cornell.edu/cgi-bin/MetGenMAP/home.cgi) was adopted to characterize the KEGG pathways of the DE genes identified, using the pearl module CPAN as described previously [60].

To validate the results derived from transcriptome analysis, we randomly selected ten DE genes, including five with upregulated expression pattern whereas another five to be downregulated in expression to be subjected to qRT-PCR analysis, using same aliquots of root tissues of lines Sen 2 and Anti 1 after drought treatment as samples. The five genes with upregulated pattern for RNA-seq analysis included TaWRK2, wall-associated receptor kinase 2 (TraesCS3B02G007300); TaWRKY28, transcription factor WRKY28 (TraesCS7B02G418400); TaCML31, calcium-binding protein CML31 (TraesCS3B02G553900); TaMPK18, mitogen-activated protein kinase kinase kinase 18 (TraesCS3B02G288100); and TaZFP1, zinc finger protein 1 (TraesCS5A02G401200). The five DE genes shown downregulated in expression for analysis included TaPAO, primary amine oxidase-like (TraesCS4A02G020900); TaCA, carbonic anhydrase (TraesCS3A02G230000); TaCP450, cytochrome P450 (TraesCS7D02G271100); TaUBI6, putative E3 ubiquitin-protein ligase SINA-like 6 (TraesCS3B02G288100); and TaFR1, fatty acyl-CoA reductase 1-like (TraesCS3B02G016500). qRT-PCR was performed similarly in detecting the transcripts of TaPYL4 as aforementioned, using gene specific primers (Additional file 11). The constitutive gene Tatubulin was used as internal standard to normalize target gene transcripts.

Statistical Analysis

Averages of gene transcripts, growth traits, plant biomass, stress-associated physiological traits, and the photosynthetic parameters were all derived from the triplicate results. Standard errors for the averages and statistical significance analyses for the traits assessed were determined based on the Statistical Analysis System software (SAS Corporation, Cory, NC, USA).

Acknowledgements

Not applicable.

Authors’ contributions

Kai Xiao conceived the research and wrote the paper. Yanyang Zhang, Yingjia Zhao, Tianjiao Li, Chengyang Ni, Le Han, and Pingping Du performed the experiments and analyzed the data. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (No. 31671686; 31872869), National Key Research and Development program, Science and Technology Planning Project of Hebei Province (No. 216Z6401G).

Availability of data and materials

The transcriptome sequencing data reported in this paper have been deposited in the National Center for Biotechnology Information (NCBI) database under project number PRJNA815379 and the SRA ID with number SRR18310810.

Ethics approval and consent to participate

The authors confirm that all the experimental methods complied with relevant institutional, national, and international guidelines and legislation. We also specify that all the permissions or licenses regarding the seed collection were obtained.

Consent for publication

Not applicable.

Competing interests

The author declare that they have no competing interests.

Fahad S, Bajwa AA, Nazir U, Anjum SA, Farooq A, Zohaib A, Sadia S, Nasim W, Adkins S, Saud S, Ihsan MZ, Alharby H, Wu C, Wang D, Huang J. Crop production under drought and heat stress: plant responses and managements. Front. Plant Sci. 2017; 8: 1147
DaCosta M, Wang Z, Huang B. Physiological adaptation of kentucky bluegrass to localized soil drying. Crop Sci. 2004; 44: 1307–1314
Zhang F, Wei Q, Shi J, Jin X, He Y, Zhang Y, Luo Q, Wang Y, Chang J, Yang G, He G. Brachypodium distachyon BdPP2CA6 interacts with BdPYLs and BdSnRK2 and positively regulates salt tolerance in transgenic Arabidopsis. Front. Plant Sci. 2017; 8: 264
Monroe JG, Powell T, Price N, Mullen JL, Howard A, Evans K, Lovell JT, McKay JK. Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function. ELife 2018; 7: e41038
Baldoni, E., Genga, A., and Cominelli, E. Plant MYB transcription factors: their role in drought response mechanisms. Int. J. Mol. Sci. 2015; 16: 15811–15851
Tan W, Zhang D, Zhou H, Zheng T, Yin Y, Lin H. Transcription factor HAT1 is a substrate of SnRK2.3 kinase and negatively regulates ABA synthesis and signaling in Arabidopsis responding to drought. PLoS Genet. 2018; 14
Zhang Z, Ali S, Zhang T, Wang W, Xie L. Identification, evolutionary and expression analysis of PYL-PP2C-SnRK2s gene families in soybean. Plants (Basel) 2020; 9: 1356
Meraj TA, Fu J, Raza MA, Zhu C, Shen Q, Xu D, Wang Q. Transcriptional factors regulate plant stress responses through mediating secondary metabolism. Genes (Basel) 2020; 11: 346
Yang M, Duan X, Wang Z, et al. Overexpression of a voltage-dependent anion-selective channel (VDAC) protein-encoding gene, MsVDAC, from Medicago sativa confers cold and drought tolerance to transgenic tobacco. Genes (Basel) 2021; 12: 1706.
Li X, Li G, Li Y, Kong X, Zhang L, Wang J, Li X, Yang Y. ABA receptor subfamily III enhances abscisic acid sensitivity and improves the drought tolerance of Arabidopsis. Int. J. Mol. Sci. 2018; 19: 1938.
Yoshida T, Mogami J, Yamaguchi-Shinozaki, K. ABA-dependent and ABA-independent signaling in response to osmotic stress in plants. Curr. Opin. Plant Biol. 2014; 21: 133–139
Gupta MK, Lenka SK, Gupta S, Rawal RK. Agonist antagonist and signaling modulators of ABA receptor for agronomic and post-harvest management. Plant Physiol. Biochem. 2020; 148: 10–25
Park SY, Fung P, Nishimura N, et al. Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science 2009; 324: 1068–1071
Zhao, Y., Zhang, Z., Gao, J., Wang, P., Hu, T., Wang, Z., Hou, Y.J., Wan, Y., Liu, W., Xie, S., Lu, T., Xue, L., Liu, Y., Macho, A.P., Tao, W.A., Bressan, R.A., and Zhu, J.K. (2018). Arabidopsis duodecuple mutant of PYL ABA receptors reveals PYL repression of ABA-independent SnRK2 activity. Cell Rep. 23, 3340–3351
Peirats-Llobet M, Rodriguez PL, Geilfus CM, et al. The role of Arabidopsis ABA receptors from the PYR/PYL/RCAR family in stomatal acclimation and closure signal integration. Nat. Plants. 2019 ; 5: 1002–1011
Schweighofer, A., Hirt, H., Meskiene, I. Plant PP2C phosphatases: emerging functions in stress signaling. Trends Plant Sci. 2004; 9: 236–243
Santiago J, Dupeux F, Betz K, Antoni R, Gonzalez-Guzman M, Rodriguez L, Márquez JA, Rodriguez PL. Structural insights into PYR/PYL/RCAR ABA receptors and PP2Cs. Plant Sci. 2012; 182: 3–11
Chen Y, Feng L, Wei N, et al. Overexpression of cotton PYL genes in Arabidopsis enhances the transgenic plant tolerance to drought stress. Plant Physiol. Biochem. 2017;115: 229–238
Xie T, Ren R, Zhang YY, Pang Y, Yan C, Gong X, He Y, Li W, Miao D, Hao Q, Deng H, Wang Z, Wu JW, Yan N. Molecular mechanism for inhibition of a critical component in the Arabidopsis thaliana abscisic acid signal transduction pathways, SnRK2.6, by protein phosphatase ABI1. J. Biol. Chem. 2012; 287: 794–802
Mega R, Abe F, Kim JS, Tsuboi Y, Tanaka K, Kobayashi H, Sakata Y, Hanada K, Tsujimoto H, Kikuchi J, Cutler SR, Okamoto M. Tuning water-use efficiency and drought tolerance in wheat using abscisic acid receptors. Nat. Plant 2019; 5: 153–159
Gonzalez-Guzman, M, Pizzio GA, Antoni R, et al. Arabidopsis PYR/PYL/RCAR receptors play a major role in quantitative regulation of stomatal aperture and transcriptional response to abscisic acid. Plant Cell, 2012; 24: 2483–2496
Zhao H, Nie K, Zhou H, Yan X, Zhan Q, Zheng Y, Song CP. ABI5 modulates seed germination via feedback regulation of the expression of the PYR/PYL/RCAR ABA receptor genes. New Phytol. 2020; 228: 596–608
Yadav SK, Santosh VV, Verma RK, Yadav P, Saroha A, Wankhede DP, Chaudhary B, Chinnusamy V. Genome-wide identification and characterization of ABA receptor PYL gene family in rice. BMC Genom. 2020; 21: 676
González-Guzmán M, Rodríguez L, Lorenzo-Orts L, Pons C, Sarrión-Perdigones A, Fernández MA, Peirats-Llobet M, Forment J, Moreno-Alvero M, Cutler SR, Albert A, Granell A, Rodríguez PL. Tomato PYR/PYL/RCAR abscisic acid receptors show high expression in root, differential sensitivity to the abscisic acid agonist quinabactin, and the capability to enhance plant drought resistance. J. Exp. Bot. 2014; 65: 4451–4464
Chen P, Sun YF, Kai WB, Liang B, Zhang YS, Zhai XW, Jiang L, Du YW, Leng P. Interactions of ABA signaling core components (SlPYLs, SlPP2Cs, and SlSnRK2s) in tomato (Solanum lycopersicon). J. Plant Physiol. 2016; 205: 67–74
Chen Y, Feng L, Wei N, Liu ZH, Hu S, Li XB. Overexpression of cotton PYL genes in Arabidopsis enhances the transgenic plant tolerance to drought stress. Plant Physiol. Biochem. 2017; 115: 229–238
Guo D, Zhou Y, Li HL, Zhu JH, Wang Y, Chen XT, Peng SQ. Identification and characterization of the abscisic acid (ABA) receptor gene family and its expression in response to hormones in the rubber tree. Sci. Rep. 2017;7: 45157
Lenka SK, Muthusamy SK, Chinnusamy V, Bansal KC. Ectopic expression of rice PYL3 enhances cold and drought tolerance in Arabidopsis thaliana. Mol. Biotechn. 2018; 60: 350–361
Antoni R, Gonzalez-Guzman M, Rodriguez L, Rodrigues A, Pizzio GA, Rodriguez PL. Selective inhibition of clade A phosphatases type 2C by PYR/PYL/RCAR abscisic acid receptors. Plant Physiol. 2012;158: 970–980
Fan W, Zhao M, Li S, Bai X, Li J, Meng H, Mu Z. Contrasting transcriptional responses of PYR1/PYL/RCAR ABA receptors to ABA or dehydration stress between maize seedling leaves and roots. BMC Plant Biol. 2016; 16: 99
Ma S, Bachan S, Porto M, Bohnert HJ, Snyder M, Dinesh-Kumar SP. Discovery of stress responsive DNA regulatory motifs in Arabidopsis. PLoS One 2012;7: e43198
Kim H, Hwang H, Hong JW, Lee YN, Ahn IP, Yoon IS, Yoo SD, Lee S, Lee SC, Kim BG. A rice orthologue of the ABA receptor, OsPYL/RCAR5, is a positive regulator of the ABA signal transduction pathway in seed germination and early seedling growth. J. Exp. Bot. 2012; 63: 1013–1024
Kim H, Lee K, Hwang H, Bhatnagar N, Kim DY, Yoon IS, Byun MO, Kim ST, Jung KH, Kim BG. Overexpression of PYL5 in rice enhances drought tolerance, inhibits growth, and modulates gene expression. J. Exp. Bot. 2014; 65: 453–464
Agurla S, Gahir S, Munemasa S, Murata Y, Raghavendra AS. Mechanism of stomatal closure in plants exposed to drought and cold stress. Adv. Exp. Med. Biol. 2018; 1081: 215–232
Azhar N, Su N, Shabala L, Shabala S. Exogenously applied 24-epibrassinolide (EBL) ameliorates detrimental effects of salinity by reducing K⁺ efflux via depolarization-activated K⁺ channels. Plant Cell Physiol. 2017; 58: 802–810
He F, Wang H, Li H, Su Y, Li S, Yang Y, Feng C, Yin W, Xia X. PeCHYR1, a ubiquitin E3 ligase from Populus euphratica, enhances drought tolerance via ABA-induced stomatal closure by ROS production in Populus. Plant Biotechnol. J. 2018; 16: 1514–1528
Bai G, Xie H, Yao H. Genome-wide identification and characterization of ABA receptor PYL/RCAR gene family reveals evolution and roles in drought stress in Nicotiana tabaccum. BAC Genom. 2019; 20: 575
Xu P, Zhang X, Su H, Liu X, Wang Y, Hong G. Genome-wide analysis of PYL-PP2C-SnRK2s family in Camellia sinensis. Bioengineered 2020;11: 103–115
Zhang J, Wang N, Miao Y, Hauser F, McCammon JA, Rappel WJ, Schroeder JI. Identification of SLAC1 anion channel residues required for CO₂/bicarbonate sensing and regulation of stomatal movements. Proc. Nat. Acad. Sci. 2018; 115: 11129–11137
Chakhchar A, Chaguer N, Ferradous A, Filali-Maltouf A, El Modafar C. Root system response in Argania spinosa plants under drought stress and recovery. Plant Sig. Behav. 2018; 13: e1489669
Fu Y, Ma H, Chen S, Gu T, Gong J. Control of proline accumulation under drought via a novel pathway comprising the histone methylase CAU1 and the transcription factor ANAC055. J. Exp. Bot. 2018; 69: 579–588
Funck D, Baumgarten L, Stift M, von Wirén N, Schönemann L. Differential contribution of P5CS isoforms to stress tolerance in Arabidopsis. Front. Plant Sci. 2020;11: 565134
Sauer M, Kleine-Vehn J. PIN-FORMED and PIN-LIKES auxin transport facilitators. Review 2019; 146: 168088
Adamowski M, Friml J. PIN-dependent auxin transport: action, regulation, and evolution. Plant Cell 2015; 27: 2032
Díaz-Sala C. A perspective on adventitious root formation in tree species. Plants (Basel) 2020; 9: 1789
Julkowska MM, Hoefsloot HC, Mol S, Feron R, de Boer GJ, Haring MA, Testerink C. Capturing Arabidopsis root architecture dynamics with ROOT-FIT reveals diversity in responses to salinity. Plant Physiol. 2014; 166: 1387–1402
Reynhout S, Janssens V. Physiologic functions of PP2A: Lessons from genetically modified mice. Biochimica ET Biophysica Acta Mol. Cell Res. 2019; 1866: 31–50
Verslues PE, Kim YS, Zhu JK, Altered ABA. Proline and hydrogen peroxide in an Arabidopsis glutamate:glyoxylate aminotransferase mutant. Plant Mol. Biol. 2007; 64: 205–217
Guo C, Zhao X, Liu X, Zhang L, Gu J, Li X, Lu W, Xiao K. Function of wheat phosphate transporter gene TaPHT2;1 in Pi translocation and plant growth regulation under replete and limited Pi supply conditions. Planta 2013; 237:1163–1178
Dovzhenko A, Bergen U, Koop HU. Thin-alginate-layer technique for protoplast culture of tobacco leaf protoplasts: shoot formation in less than two weeks. Protoplasma 1998; 204: 114–118
Yoo SD, Cho YH, Sheen J. Arabidopsis mesophyll protoplasts: a versatile cell system for transient gene expression analysis. Nat. Prot. 2007; 2: 1565–1572
Sheerin DJ. Investigation of light-regulated protein-protein interactions using yeast two-hybrid assays. Meth. Mol. Biol. 2019; 2026: 1–19
Sun Z, Ding C, Li X, Xiao K. Molecular characterization and expression analysis of TaZFP15, a C2H2-type zinc finger transcription factor gene in wheat (Triticum aestivum L.). J. Integr. Agric. 2012; 11: 31–42
Yang T, Yao S, Hao L, Zhao Y, Lu W, Xiao K. Wheat bHLH-type transcription factor gene TabHLH1 is crucial in mediating osmotic stresses tolerance through modulating largely the ABA-associated pathway. Plant Cell Rep. 2016; 35: 2309–2323
Imadi SR, Kazi AG, Ahanger MA, Gucel S, Ahmad P. Plant transcriptomics and responses to environmental stress: an overview. J. Genet. 2015; 94: 525–537
Liu L, Huang L, Lin X, et al. Hydrogen peroxide alleviates salinity-induced damage through enhancing proline accumulation in wheat seedlings. Plant Cell Rep. 2020; 39: 567–575
Zhong S, Joung JG, Zheng Y, Chen YR, Liu B, Shao Y, Xiang JZ, Fei Z, Giovannoni JJ. High-throughput Illumina strand-specific RNA sequencing library preparation. Cold Spring Harbor Protocol 2011; 8: 940–949.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinform. 2014; 30: 2114–2120
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinform. 2010; 26: 139–140
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G. GO: TermFinder-open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinform. 2004; 20: 3710–3715

No competing interests reported.

Additionalfiles15811.docx
Additional file 1: Phylogenetic relations among TaPYL4 and its homologous genes distributed in various plant speciesAdditional file 2: The signals initiated by TaPYL4-GFP in protoplast expression system of N. benthamiana observed under florescent microscopeAdditional file 3: Expression levels insertion copies of the target gene in TaPYL4 transgenic lines Additional file 4: Expression levels of the target gene detected in transgenic lines with TaP5CS1 knockdown expressionAdditional file 5: Expression levels of the target gene detected in transgenic lines with TaPIN9 knockdown expressionAdditional file 8: qRT-PCR results in roots for the differentially expressed genes with upregulated expression pattern identified based on RNA-seq analysisAdditional file 9: qRT-PCR results in roots for the differentially expressed genes with downregulated expression pattern identified based on RNA-seq analysisAdditional file 10: Alignment results among TaPYL4 and its homeologs as well as expression of TaPYL4 homeologs upon low-Pi stressAdditional file 11: PCR primers used in this study
Additionalfile6.xls
Additional file 6: Up-regulated differentially expressed (DE) genes underlying modulation of TaPYL4
Additionalfile7.xls
Additional file 7: Down-regulated differentially expressed (DE) genes underlying modulation of TaPYL4

Download PDF

Editorial decision: Major revision
11 Jul, 2022
Reviews received at journal
08 Jul, 2022
Reviews received at journal
28 Jun, 2022
Reviewers agreed at journal
23 Jun, 2022
Reviewers invited by journal
23 Jun, 2022
Editor assigned by journal
23 Jun, 2022
Editor invited by journal
23 Jun, 2022
Submission checks completed at journal
23 Jun, 2022
First submitted to journal
09 Jun, 2022

You are reading this latest preprint version

TaPYL4, an ABA receptor gene of wheat, positively regulates plant drought adaptation through modulating the osmotic stress-associated processes

Status:

Version 1

Abstract

Figures

Background

Results

The Components Constitutes Aba Signaling Module With Tapyl4

Expression Patterns Of P5cs And Pin-formed Family Genes

Discussion

Conclusions

Methods

Determination Of Tapyl4 Subcellular Localization

Yeast Two-hybrid Assays

Assays of Drought Stress-associated Physiological Traits and Photosynthetic Parameters

Analyses Of Expression And Function Of P5cs And Pin-formed Family Genes

Transcriptome Analysis

Statistical Analysis

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1