Skip to main content

Overexpression VaPYL9 improves cold tolerance in tomato by regulating key genes in hormone signaling and antioxidant enzyme



Abscisic acid (ABA) has been reported in controlling plant growth and development, and particularly dominates a role in resistance to abiotic stress. The Pyrabactin Resistance1/PYR1-Like /Regulatory Components of ABA receptors (PYR1/PYL/RCAR) gene family, of which the PYL9 is a positive regulator related to stress response in ABA signaling transduction. Although the family has been identified in grape, detailed VaPYL9 function in cold stress remains unknown.


In order to explore the cold tolerance mechanism in grape, VaPYL9 was cloned from Vitis amurensis. The subcellular localization showed that VaPYL9 was mainly expressed in the plasma membrane. Yeast two-hybrid (Y2H) showed VaPCMT might be a potential interaction protein of VaPYL9. Through the overexpression of VaPYL9 in tomatoes, results indicated transgenic plants had higher antioxidant enzyme activities and proline content, lower malondialdehyde (MDA) and H2O2 content, and improving the ability to scavenge reactive oxygen species than wild-type (WT). Additionally, ABA content and the ratio of ABA/IAA kept a higher level than WT. Quantitative real-time PCR (qRT-PCR) showed that VaPYL9, SlNCED3, SlABI5, and antioxidant enzyme genes (POD, SOD, CAT) were up-regulated in transgenic tomatoes. Transcriptome sequencing (RNA-seq) found that VaPYL9 overexpression caused the upregulation of key genes PYR/PYLPYL4MAPK17/18, and WRKY in transgenic tomatoes under cold stress.


Overexpression VaPYL9 enhances cold resistance of transgenic tomatoes mediated by improving antioxidant enzymes activity, reducing membrane damages, and regulating key genes in plant hormones signaling and antioxidant enzymes.

Peer Review reports


As sessile organism, plant can’t move randomly when they are subjected to various adverse environments such as soil salinity, drought, and extreme temperatures. Therefore, they must cope with these abiotic stress factors to adapt to changing environments by various mechanisms [1]. Adverse environmental factors affect plant growth that ultimately causes losses in crop quality and yield all over the world [2, 3]. Although plant has formed a complex response regulation networking to combat stress challenges, it has certain limitations. Therefore, transgenic engineering and transcriptome sequencing are employed to improve adversity stress and enhance plant tolerance [2, 4, 5]. Cold stress is one of the primary environmental factors for plant growth, and cold stress responses and the generation of cold tolerance plants have received much attention in recent years [6]. A full understanding of the molecular mechanisms of cold tolerance will be conducive to plant breeding of cold tolerance species and provide cold tolerance candidate genes for the following researches.

ABA is an important phytohormone that regulates plant growth and development including cell elongation, seed dormancy, leaf senescence, fruit abscission, and stomatal closure [7, 8]. In addition, ABA is also an important stress response phytohormone, and it’s called as a ‘stress hormone’. Many researches have verified that ABA functions as an endogenous messenger in abiotic stress responses, especially in drought, salinity and cold stresses, and plants keep a high ABA accumulation [9,10,11,12,13]. The PYR/PYL/RCAR proteins are currently considered as the typical ABA receptors according to structural, biochemical, and genetic evidences implemented by two independent research groups [10, 11, 14, 15]. The PYR/PYL/RCAR proteins belong to Bet v 1-like superfamily members and include special STAR-RELATED LIPID-TRANSFER (START) related to metabolism, and hydrophobic compounds such as lipids, hormones, and antibiotics [16, 17]. A number of studies found that PYR/PYL/RCAR which was bound to ABA negatively regulated type 2C protein phosphates (PP2C) in ABA signaling transduction, in turn which negatively regulated SNF1-related protein kinase 2 (SnRK2). These three core components formed a double negative regulatory system in ABA signal transduction. Under normal growth conditions, PP2C inactivated SnRK2 through dephosphorylation, and ABA signal was non-activation. When plants were subjected to adverse environmental conditions, ABA production maintained a high level, ABA bound to PYR/PYL/RCAR proteins and promoted interaction function with PP2C to inhibit PP2C phosphatase activity, relieving the inhibition of PP2C on the protein kinase SnRK2, and SnRK2 was activated, and then phosphorylated downstream transcription factors or membrane proteins to arouse ABA signaling reaction [18, 19].

There existed 14 PYL gene members in Arabidopsis thaliana, except PYL13, other gene sequences and structures were highly conserved, and could bind to ABA inhibiting PP2C enzyme activity [11, 20]. Nine of the 14 PYR/PYL/RCAR family members were the most robustly co-purified proteins with ABI1. The data supported that ABI1 interacted with multiple members of the PYR/PYL/RCAR family in Arabidopsis [21]. The knockout of AtPYL8 reduced the sensitivity of ABA on primary root growth and lateral root formation and the overexpression AtPYL9 promoted the lateral root elongation in the presence of ABA. AtPYL9 was found to interact with MYB77 and MYB44 in vivo to regulate auxin-responsive genes affecting roots growth [22, 23]. Additionally, overexpression SlPYL9 accelerated tomato fruits ripening [24]. PYR/PYL receptors regulated stomatal aperture highlights the relevance of the PYR/PYL pathway to cope with drought stress [14]. PYR/PYL/RCAR ABA receptors were close to abiotic stresses, RCAR12 and RCAR13 were confirmed to play positive roles in regulating extreme temperature, including cold and high temperature in Arabidopsis [25]. ABA receptors, RCAR11, RCAR12, RCAR13 and RCAR14 showed a hypersensitivity for ABA, increased genes expression related stress responses, as well as enhanced drought tolerance [26]. Under drought and low temperature treatments, the ectopic expression of OsPYL3 could enhance the cold and drought resistance of Arabidopsis thaliana, and OsPYL10 overexpression had the potential ability to improve both drought and cold stress tolerance of rice [27, 28]. It has been reported that RCAR5/PYL11 was involved in response to cold stress by delaying seed germination and controlling stomatal closure via ABA-dependent and ABA-independent pathways, respectively [29]. Constitutive expression of PtPYRL1 and PtPYRL5 significantly enhanced the resistance to drought, osmotic, and cold stresses by positively regulating ABA signaling in Populus L. [30]. Overexpression AtPYL5 promoted drought tolerance [31, 32]. Exogenous hormones induced AtPYL8 expression to improve drought stress [33]. MdPYL9 improved tolerance to drought stress in transgenic apple [34]. The expression level of GhPYL9-11A was higher in drought-tolerant cotton cultivars than in drought-sensitive cottons under drought treatment [35]. These functional achievements of PYR/PYL/RCAR showed that researches on stresses tolerance of PYR/PYL/RCAR mainly focused on drought. Although cold stress of PYL genes has been studied in several species, its knowledge was less than drought.

To date, a total of eight VvPYL genes were identified from the grape genome database, which could respond to cold, salt, and PEG stresses. The researcher successfully cloned six PYL genes from the wild Chinese V. yeshanensis (‘Yanshan-1’), and overexpression VyPYL9 in Arabidopsis enhanced ABA sensitivity, and improved drought tolerance under ABA and drought treatments [36, 37]. VvPYL1 modulated ABA signaling pathway by inhibiting PP2C expression [38]. VlPYL1 from ‘Kyoho’ promoted anthocyanin accumulation in grape berry skin [39]. There were most widely studies on CBF and COR genes related to cold stress in grape [40, 41]. However, knowledge of PYLs related to cold stress is relatively narrow in grape. Therefore, in this study, VaPYL9 was cloned from Vitis Amurensis with strong cold tolerance, and it was ectopically expressed in tomatoes. Results showed that the VaPYL9 overexpression was induced by cold stress. The overexpression of VaPYL9 in tomatoes contributed to an increasing of the ability to scavenge reactive oxygen species and reduce membrane damage, and keep a high ABA level in transgenic tomatoes. Through RNA-seq, the result showed VaPYL9 bound to ABA leading to a series of reactions in metabolism and biosynthesis, and mainly focused on MAPK and plant hormones signaling transduction paths under cold stress. In all, these findings indicate that VaPYL9 dominates a positive role in plant response to cold stress and may be an important candidate gene for molecular breeding for the improvement of stress tolerance in grape.


Detailed information and evolutionary analysis of VaPYL9

In order to have a basic understanding of VaPYL9 (GSVIVG01027078001), the physicochemical property was analyzed in supplementary Table S1. The online prediction found that isoelectric point (PI) was 6.38, the molecular weight was 20.05 KD, grand average of hydropathicity (GRAVY) was -0.193, instability index(II) was 37.65, and the aliphatic index(AI) was 95.67. The subcellular localization prediction showed that VaPYL9 was expressed mainly in cytoplasm. VaPYL9 protein was stable because instability index was less than 40. Protein was hydrophilic amino acid when GRAVY was negative value. Therefore, VaPYL9 protein had strong hydrophilicity. Phylogenetic tree analysis showed that PYL genes from grape, tomato, and Arabidopsis were divided into three subgroups, and VaPYL9 was clustered in III subgroup (Fig. 1). VaPYL9 and most tomato genes were clustered in a close evolutionary branch. Multiple sequences alignment showed that PYL genes had a polyketide cyclase/dehydrase and lipid transport domain (polyketide-cyc2, Pfam accession: PF10604) and highly conservative amino acid sequences (Supplementary Fig. S1). The finding indicated VaPYL9 might have a similar function to PYL genes of tomatoes.

Fig. 1
figure 1

Evolution analysis of PYL genes from grape, Arabidopsis thaliana, and tomato. The phylogenetic tree was constructed using PYL sequences of grape, tomato, and A. thaliana in MEGA X, which was clustered into three subgroups. Different background colors represent different subgroups. VaPYL9 marked with red ellipse was in the III class. The green circle, blue circle, and pink circle represent grapes, A. thaliana, and tomato PYL genes, respectively 

Transient expression of VaPYL9 in tobacco

Based on the above online prediction, VaPYL9 was located mainly in the cytoplasm. To further verify the prediction result, the pART-CAM-EGFP-VaPYL9 vector was constructed. The pART-CAM-EGFP-VaPYL9 and pART-CAM-EGFP empty construct plasmids were transformed into Agrobacterium tumefaciens strain GV3101 (Supplementary Fig. S2), and then were injected into tobacco leaves. The physical investigation found that VaPYL9 was not only located in plasma membrane and nucleus but also in cytoplasm (Fig. 2a).

Fig. 2
figure 2

(a). Subcellular localization of VaPYL9. Four-weekend tobaccos were infected with bacterial liquid. After 24 h, infected leaves were observed and photoed in cell membrane and cytoplasm. (b). Yeast two -hybrid (Y2H) assay. Combinations of pGADT7-T with pGBKT7-53 and pGBKT7-Lam were served as positive and negative controls, respectively. VaPYL9 and VaPCMT were commonly transformed into Y2H gold. The result indicated that there was an interaction relationship between VaPYL9 and VaPCMT Proteins 

The cloning of VaPYL9 and obtaining of transgenic tomatoes

To specifically investigate cold tolerance of the VaPYL9, the Open Reading Framework (ORF) with 591 amino acids was inserted into overexpression vector pCAMBIA1301. The presence of VaPYL9 cDNA fragment was verified by PCR and sequenced in Escherichia coalition (DH5α) and Agrobacterium tumefaciens strain GV3101 (Supplementary Fig. S3). Through a Agrobacterium tumefaciens medicated method, VaPYL9 fragment was successfully introduced into tomatoes and identified by PCR (Supplementary Fig S4). Three stable and homozygous T3 lines (#2, #6, #9) were screened from transgenic tomato lines and named OE1, OE2, and OE3, respectively.

The interaction proteins screening of VaPYL9

To preliminary investigate interaction proteins of VaPYL9 from AD library-related cold stress, pGBKT7-VaPYL9 was constructed and introduced into Y2H gold (Supplementary Fig. S5a). Self-activation and toxicity detection were detected in DDQ and DDQ (x-α-gal) plates. Bacterial plaque didn’t become blue in DDQ and DDQ (x-α-gal) plates. This indicated that VaPYL9 had no self-activation and toxicity (Supplementary Fig. S6a). Through the mating hybridization method, interaction proteins of VaPYL9 were preyed from AD library related-cold stress. The presence of potential interaction proteins was verified by PCR using AD vector primers. Finally, VaPCMT (protein-L-isoaspartate O-methyltransferase 1) was screened by blast in NCBI. pGADT7-VaPCMT was commonly constructed and introduced into Y2H gold (Supplementary Fig. S5b). Using the same procedure, self-activation and toxicity were detected, indicating that VaPCMT didn’t produce self-activation and toxicity (Supplementary Fig. S6b). pGBKT7-VaPYL9 and pGADT7-VaPCMT proteins were co-transformation into Y2H gold in QDO and QDO(x-α-gal) plates (Supplementary Fig. S7). Compared to negative and positive controls, bacterial plaque became gradually blue in DDQ and DDQ (x-α-gal) plates after 2 days. VaPCMT was a potential interaction protein of VaPYL9. Comprehensive above results, the findings showed that VaPYL9 and VaPCMT might be a potential interaction function in vivo (Fig. 2b).

Phenotype observation and degree of membrane damage in tomatoes under cold stress

To evaluate the role of VaPYL9 in response to cold stress, WT and OE lines were exposed to a long-term cold environment. qRT-PCR showed that relative expression of VaPYL9 was upregulated when grapes were induced by cold stress, and up to a maximum in 48 h (Fig. 3a). Compared with WT tomatoes, transgenic tomatoes were highly expressed in cold stress (Fig. 3b). Through phenotype observation in WT and OE lines under 4 ℃ cold stress for 0 h, 48 h, results indicated that WT and OE lines showed wilting difference but no significant in appearance, and WT showed cold damage and leaves coiled and wilted more heavily than OE lines (Fig. 3c). Some indexes involved in membrane damages were determined. In general, REL is frequently associated with abiotic stress, the membrane permeability tends to the increasing of REL when plants are exposed to abiotic stress. With the prolongation of a cold duration time, REL showed an increasing trend, indicating the degree of membrane damage was increasing in WT and OE lines. From 12 h, relative electrolyte leakage (REL) was quickly increasing. REL of OE plants was lower than WT under cold stress (Fig. 3d). The result showed VaPYL9 overexpression enhanced membrane protection feature.

Fig. 3
figure 3

Phenotype investigation and degree of membrane damage between WT and OE. Six-week-old tomatoes were used to detect physiological indexes related to cold stress. (a). Expression level of VaPYL9 in grape under 4 ℃ stress for 0 h, 6 h, 12 h, 24 h, 48 h. (b). Expression level of VaPYL9 in WT and transgenic tomatoes under cold stress and nonstress. (c). Tomatoes phenotype was observed in WT and OE under 4 ℃ for 0 h, 48 h. (d). The determination of relative electrolyte leakage in WT and OE plants under 4 ℃ for 0 h, 6 h, 12 h, 24 h, 48 h. (e). The DAB and NBT staining of WT and OE plants under 4 ℃ for 0 h, 12 h, 48 h. (f). The content of H2O2 was determined in WT and OE plants under 4 ℃ for 0 h, 6 h, 12 h, 24 h, 48 h. (g). The content of MDA was determined in WT and OE plants under 4 ℃ for 0 h, 6 h, 12 h, 24 h, 48 h. (h). The content of proline was determined in WT and OE plants under 4 ℃ for 0 h, 6 h, 12 h, 24 h, 48 h. The bar represents the value of standard error (p < 0.01). Lowercase letters represent the significance level. Means and SE values were calculated from at least three independent experiments

DAB and NBT staining were used to evaluate H2O2 and O−2 content at 0 h, 12 h, 48 h. Leaves staining showed there was no significant colour difference in WT and OE lines under non-stress. H2O2 and O−2 didn’t accumulate and kept a homeostasis in leaves. Under cold stress, the degree and area of leaves staining in WT were deeper and expanded than OE lines, respectively (Fig. 3e). WT plants hold a high accumulation of H2O2 and O−2 under cold stress for 48 h. As was shown in Fig. 3f, through the determination of H2O2, result found H2O2 content of OE was lower than WT at all time points, and OE was lowest at 48 h. The data supported that OE plants had prominent scavenging effects against superoxide anion radicals. VaPYL9 expression was related to oxidative stress responses. Malondialdehyde (MDA) content showed an increasing trend in WT and OE with the time of cold stress increasing, but in all time points MDA of OE was lower than WT lines. Low-level MDA indicated that lipid peroxidation of membrane was relatively low (Fig. 3g). There were almost no differences in proline content between OE and WT under non-stress. Proline content of OE was higher than WT and reached a maximum at 48 h. High-level proline protected plant cells from cold stress (Fig. 3h). Taken together, the overexpression VaPYL9 improved the ability to scavenge reactive oxygen species (ROS), and reduce the degree of membrane lipid peroxidation and cell dehydration property to protect tomatoes from cold stress.

The changing of antioxidant enzymes activity and endogenous hormones level under cold stress

Catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) play an important role in maintaining the homeostasis of ROS metabolism and antioxidation defense system. The results showed that there was no significant difference in enzyme activities of OE and WT at 0 h, while CAT and POD activities of WT and OE lines showed an increasing trend in all cold time points. Interestingly, VaPYL9 expression tomatoes exhibited higher activities of SOD, CAT, and POD enzymes compared to WT under cold stress. These three enzyme activities reached a maximum under cold stress for 48 h (Fig. 4a). In all, the results indicated that VaPYL9 expression caused an increase of antioxidant enzymes activity, reducing ROS damage for tomatoes.

figure 4

The changing of antioxidant enzyme activities and endogenous hormones in WT and OE plants. Six-weekend tomatoes were used as materials to detect enzyme activities and endogenous hormones. (a). POD, SOD, and CAT enzymes activities were determined in WT and OE plants under 4 ℃ for 0 h, 6 h, 12 h, 24 h, 48 h. b. (b). The content of IAA and ABA was determined, and ABA/IAA was calculated in WT and OE plants under 4 ℃ for 0 h, 6 h, 12 h, 24 h, 48 h. The bar represents the value of standard error (p < 0.01). Lowercase letters represent the significance level. Means and SE values were calculated from at least three independent experiments

Endogenous hormones were determined by High Performance Liquid Chromatography (HPLC). The content of ABA was increasing while IAA was decreasing with the prolongation of cold duration time under cold stress. VaPYL9 expression tomatoes exhibited higher content of ABA compared to WT under cold stress, and ABA content reached a maximum under cold stress for 48 h. On the contrary, VaPYL9 overexpression tomatoes exhibited lower content of IAA compared to WT under cold stress. ABA/IAA of OE was gradually increasing compared to WT with cold duration time increasing (Fig. 4b). The findings supported that VaPYL9 overexpression resulted in higher ABA content and ABA/IAA, promoting cold tolerance for plants.

qRT-PCR analysis

qRT-PCR was employed to perform expression level measurement under cold stress for 0 h, 48 h. The expressions of SOD (NM_001311084.1), CAT (NM_001247257), and POD (NM_001302921.2) key genes involved in antioxidant enzymes were significantly upregulated under cold stress. High expresssion levels of antioxidant enzymes led to the increase of SOD, CAT, and POD activities (Fig. 5a). Compared to WT, SlNCED3 (Solyc07g056570.1), a key rate-limiting enzyme gene in the process of ABA synthesizes, and SlABI5 (XM_010327616.3), a key gene in ABA signaling transduction was up-regulated in OE lines than WT under cold stress (Fig. 5b). Altogether, these findings indicated that VaPYL9 gene was highly expressed in transgenic plants when plants were subjected to cold stimulation, and VaPYL9 gene could positively respond to cold stress.

Fig. 5
figure 5

qRT-PCR analysis. (a). Expression level detection of key enzyme genes related to POD, SOD, and CAT genes under 4 ℃ for 0 h, 48 h. (b). Expression level detection of NCED3 and ABI5 genes related to ABA path under 4 ℃ for 0 h, 48 h. (c). Expression level detection of MAPK17/18, GsSRK, MYC2, PR5, PYL4, PYR/PYL, WRKY key genes in RNA-seq. The bar represents the value of standard error (p < 0.01). Lowercase letters represent significance level. Means and SE values were calculated from at least three independent experiments

Transcriptome sequencing analysis

To investigate the role of VaPYL9 gene in regulating cold tolerance, the RNA-Seq analysis was performed using the six-week-old transgenic and WT tomato seedlings incubated at 4 ℃ for 48 h and identified differentially expressed genes (DEGs). Basic information of sequencing data was listed in Table 1, a total of 44 Gb clean data was gained, GC content of every sample was up to over 40%, and base percentage Q30 of every sample was over 94.43%. According to the statistics of the comparison results, the comparison efficiency between reads of each sample and reference genome ranged from 95.00% to 95.97%. A total of 1270 DEGs were found in WT vs OE lines subset based on the criteria fold change > 2 and P < 0.05, of which 740 DEGs (VaPYL9-activated genes) were up-regulated, 530 DEGs were down-regulated (VaPYL9-repressed genes) (Fig. 6a, b). GO (Gene Ontology) showed that key functional annotations were associated with transcription factor activity, oxidoreductase activity, signaling receptor activity, abscisic acid binding, peroxidase activity, and protein serine/threonine kinase activity (Fig. 6c). KEGG analysis found that VaPYL9 overexpression mediated MAPK and plant signaling transduction paths to regulate cold tolerance. There were 34 DEGs in MAPK path, of which 25 DEGs were up-regulated in OE compared to WT (Fig. 6d). MAPK activated the upstream ABA signaling to function in cold defense, MAPK17/18 was highly expressed in adapting to stress (Fig. 6d). There were 43 DEGs in plant signaling transduction path, of which 29 DEGs were up-regulated in OE compared to WT (Fig. 6e). Key PYR/PYL and MYC2 genes were up-regulated in OE compared to WT. Transcription factors mainly were involved in WRKY, and were up-regulated in OE lines than WT (Fig. 6f). Signaling receptor activity expression was higher than WT, mainly involved in PYL4, PR5 and GsSRK (Fig. 6g). Through qRT-PCR verification, MAPK17/18, PYR/PYL, MYC2, WRKY, PYL4, PR5 and GsSRK, some key genes involved in MAPK, plant signaling transduction, signaling receptor activity, and transcription factors activity were up-regulated in transgenic under cold stress (Fig. 5c). In summary, these data demonstrated that VaPYL9 played a positive role in cold stress tolerance in tomatoes mediated by MAPK and plant signaling transduction paths, improving transcription factors and signaling receptor activity (Fig. 6h).

Table 1 Detailed information of RNA-seq data
Fig. 6
figure 6

The analysis of transcriptome sequencing in WT and OE plants under 4 ℃ for 48 h. (a) The heatmap of total differential expression genes. Rectangles in red and blue represent up and down-regulated expression. (b) The number of DEGs was shown by the volcano plot. Dots in red, green and red represent up, down and normal regulated genes, respectively. The number of the dots indicates the number of DEG. (c) The number of DEGs in main GO enrichments. (d) The expression level of key differential expression genes in MAPK path. Rectangles in red and blue represent up and down-regulated expression. (e) Expression level of differential expression genes in plant hormone signaling path. Rectangles in red and blue represent up and down-regulated expression. (f) Key differential expression genes expression in transcription factors activity. Rectangles in red and blue represent up and down-regulated expression. (g) Key differential expression genes expression in signaling receptor activity. Rectangles in red and blue represent up and down-regulated expression. (h) VaPYL9 overexpression participated in MAPK and plant hormone signaling KEGG path, and improved the activity of transcription factors and signaling receptors


ABA regulates cold tolerance in plant

Low temperature is a serious threat to crop yield stability, plant growth, quality improvement, and plant geography distribution [42]. To adapt to cold stress, a variety of biochemical and physiological processes have formed in plant. Plants modulate multiple morphological and physiological processes to maintain cellular and metabolic homeostasis, including changes in phytohormone homeostasis [43]. ABA is an important stress hormone that adjusts plant stress such as extreme temperature, salt, drought [1]. ABA has been demonstrated to be involved in the cold stress response through the regulation of a series of specific stress-responsive genes, and ABA can coordinate several phytohormones that integrated with cold signaling to regulate cold stress [13]. In recent years, an increasing number of studies have attached the importance of studying the mechanisms of ABA in cold tolerance. Exogenous ABA can improve the content of proline and soluble sugar, enhance water retention, reduce membrane lipid peroxidation, protect membrane integrity, and reduce photosynthetic characteristics [44, 45]. In addition, exogenous ABA application altered the expression of or related-cold and ABA genes, improving resistance to cold in Bermudagrass, maize, and cucumber [46,47,48]. Here, through previous RNA-seq, ABA receptor VaPYL9 was selected from DEGs related to cold stress to further verify the function of cold tolerance.

The ABA receptors PYLs were involved in cold stress

The function of PYL proteins on ABA-mediated stress responses has been well studied [18]. ABA receptors OsPYL3, OsPYL10, and AtPYL11 enhanced the cold tolerance of plants [27,28,29]. Vitis Amurensis Rupr, which has strong cold hardiness can withstand -35 ℃ [49]. In this study, VaPYL9 gene from Vitis Amurensis Rupr was cloned, evolutionary analysis showed that PYLs could be classified into three subfamilies (Fig. 1). The result was consistent with clusters in Brassica napusin and rice [50, 51]. The alignment of amino acid sequences revealed that these PYLs were highly conserved (Fig. 1b), especially amino acid residues that involved in ABA binding and PP2C interaction [52]. VaPYL9 had a polyketide-cyc2 domain related to polyketide cyclase/dehydrase and lipid transport as the same result in the previous study [53]. Most AtPYLs were mainly localized to the cytoplasm and the nucleus, and transduced ABA signals to activate a series of responses to various environmental factors [10, 11]. The same location result was found in this study, through predication online and physical transient transformation in tobacco, which showed that the VaPYL9 was also localized in the cytoplasm and nucleus (Fig. 2a). When plants were subjected to abiotic stress, ABA was synthesized and bound to PYL proteins, and the ABA-PYL complex interacts with PP2Cs to initiate ABA signal transduction [31, 54]. It has been reported that PYLs interacted with specific PP2C members in an ABA-dependent or ABA-independent manner [20]. In the present study, through screening of the Y2H method, VaPCMT was identified to be a potential protein with VaPYL9 commonly to mediate cold tolerance (Fig. 2b), and PCMT was related to oxidative stress in the previous study [55].

Overexpression of VaPYL9 improved cold tolerance mediated by scavenging reactive oxygen species, reducing membrane damage in tomato

Many previous evidences have shown that PYLs play crucial roles in ABA-mediated responses to cold stress [27,28,29]. Previous research reported that AtPYL9, the homolog sequence of VaPYL9 in Arabidopsis, enhanced ABA sensitivity and drought stress when expressed to high levels in transgenic Arabidopsis [56]. But PYL9 function was seldom involved in cold stress compared to drought. To preliminarily explore the function of VaPYL9 in response to cold in grape, overexpressing-VaPYL9 tomatoes lines were generated to evaluate the response of these transgenic plants under cold stress.

To have a good understanding of mechanisms associated with VaPYL9-induced cold tolerance, changes in several physiological and biochemical characteristics which are known to be associated with response to cold stress were determined. When plants are exposed to abiotic stress leading to an increasing of ROS production, which can cause oxidative damage to cellular components, and the intracellular ROS level is considered as an important indicator of cell death in response to stress conditions [57]. Here, results showed that REL, H2O2 and MDA accumulation in WT plants was higher and proline was lower than in the transgenic lines after cold stress (Fig. 3d, e, f, g, h), suggesting that VaPYL9 expression is related to oxidative stress responses. To protect cells and membrane integrity from effects of excessive ROS, plants have evolved a series of complex antioxidant defense mechanisms to keep the homeostasis of ROS metabolism, such as SOD, CAT, and POD [58,59,60]. The results of the present study showed that VaPYL9 overexpression, transgenic tomatoes exhibited higher activities of SOD, CAT, and POD enzymes compared to WT under cold stress (Fig. 4a), suggesting that the reducing ROS in transgenic tomatoes may result from increasing ROS metabolism by these enzymes. qRT-PCR showed SOD, CAT and POD gene of OE lines was up-regulated under cold stress compared to WT (Fig. 5a). Coincidentally with plant stress, overexpression of GhPYL9-11A, PtPYL1, and PtPYL5 genes also caused an increasing in activities of antioxidant enzymes and the ability to scavenge ROS [30, 35, 61]. SlNCED3, a key rate-limiting enzyme gene in the process of ABA synthesizes, and SlABI5, a key gene in ABA signaling transduction was up-regulated in OE lines under cold stress compared to WT (Fig. 5b). Consistence with the result that NCED3 promoted ABA biosynthesis under stress and ABI5 was related to multiply hormones adapting to various stresses [62, 63]. With the prolongation of a cold duration time, the ratio of ABA/ IAA was increased in OE lines than in WT (Fig. 4b). The result was similar to the previous study, the ratio of ABA/ IAA was increased in sugarcane leaves under cold stress [64]. The result showed that the increasing activities of antioxidant enzymes, reducing membrane damage, improving the ability to scavenge ROS, and high ABA content constitute an important physiological mechanism in plants driven by VaPYL9 overexpression under cold stress.

VaPYL9 improved cold tolerance mediated by MAPK pathway and multiply hormones coordination

Additionally, RNA-seq found a total of 1270 DEGs were found in WT vs OE lines subset, including 740 DEGs (VaPYL9-activated genes) that were up-regulated, and 530 DEGs were down-regulated (VaPYL9-repressed genes) (Fig. 6a, b). Through the KEGG map, the result showed overexpression VaPYL9 activated MAPK pathway and ethylene, and ABA commonly promoting adaption to cold stress in plant (Fig. 6d, e). MAPK was induced by ABA during abiotic stress. A combination of MAP3Ks, MAP3K17/18, MAP2K MKK3 and MAPKs MPK1/2/7/14 MAPK cascade related to PYR/PYL/RCAR mediated ABA signaling pathway, MPK1 and MPK2 may also phosphorize many ABA effector proteins [65, 66]. Similarly, PYR/PYL and MAP3K17/18 were up-regulated in transgenic plants under cold stress (Fig. 6d). PYL6 responded to JA by interacting with MYC2. In the presence of ABA, the interaction of the two become a strength between ABA and JA signaling pathways [67]. The ERF1 binding to response element to regulate abiotic stress and mediated ROS signaling to response salt stress [68, 69]. Coincidentally, in this study MYC2 indirectly regulated ERF1 to connect JA and ethylene enhancing defending stress, the two were up-regulated in overexpression VaPYL9 induced by cold stress (Fig. 6d, e). In all, VaPYL9 response to cold stress probably thereby MAPK and hormone signaling paths.


In summary, a working model of the role of VaPYL9 in the response of tomatoes to cold stress was proposed (Fig. 7). VaPYL9 overexpression induces the high expression of SlNCED3 and SlABI5 leading to ABA accumulation. VaPYL9 might be active in MAPK path of the downstream of ABA signaling response and plant hormones signaling by binding to ABA. The overexpression of VaPYL9 enhanced cold sensitivity in tomatoes by reducing membrane damage and improving activities of enzymatic antioxidant defense systems. RNA-seq indicated that MAPK and other hormones also might be involved in improving cold response. Our works suggest that the VaPYL9 may be a good candidate gene for breeding cold-tolerant crops.

Fig. 7
figure 7

A systematic model of the role of VaPYL9 in the response of tomatoes under 25 ℃, 4 ℃. Green and red arrows represent down and up-regulated expression levels in WT vs OE, respectively. The size of the red areas represents the degree of leaves damage

Materials and methods

Plant materials and growth conditions

One-year-old grapevines (Vitis amurensis Rupr. var. ‘Zuoshan 1’) derived from the vineyard located in Gansu Agricultural University (Lanzhou, Gansu, China), were used to extract RNA. Tobacco (Nicotiana benthamiana), tomato (Solanum lycopersicum cv. ‘Micro-Tom’) seeds were purchased from Nanjing Fengshuo Horticulture Ltd., Company (Nanjing, Jiangsu, China). Tomato seeds were soaked using sterilized distilled water for 3 h, surface sterilized in 75% (v/v) ethanol for 30 s, and then treated with 2% (v/v) NaClO for 10 min, followed by four washes with sterilized distilled water. Seeds were planted on 1/2 Murashige and Skoog (MS) medium supplemented with 30 g/L sucrose and 6.5 g/L agar, pH was 5.8–6.0. Tomatoes with two cotyledons expansion about six days were used for genetic transformation. Transgenic tomato seeds (Transgenic generation 3, T3) and WT seeds were soaked in a conical flask with warm water for 30 min, and then placed on a shaking table until most seeds showed white sprouts. Seeds with white sprouts and tobacco seeds were planted in plastic plots with seedling substrate containing a mixture of soil and vermiculite (3:1), and cultured in growth chambers at 28 °C, 8 h, the light at 135 μmol·m−2·s−1/20 °C, 16 h, dark growth condition. The six-week-old T3 transgenic and WT tomato seedlings were used for the following cold assay. Four-weeks old tobacco seedlings were used for transient transformation.

sDetailed information and evolutionary analysis of VaPYL9

The physicochemical property was analyzed by ProtParam in EXPASY ( including theoretical isoelectric point (pI), molecular weight, hydrophilicity, instability index, and grand average of hydropathicity (GRAVY). The protein is stable when the instability index is less than 40. Protein belongs to a neutral amino acid when GRAVY is in 0.5 ~ -0.5 (GRAVY < 0, strong hydrophilicity; GRAVY > 0, hydrophobicity). WoLF PSORT ( was employed to predict the subcellular localization of the VaPYL9. The protein sequences of PYL genes were downloaded from Phytozomev12.1: Home (, including 6 grape genes (Vitis vinifera L.), 16 Arabidopsis genes (Arabidopsis thaliana), and 16 tomato genes (Solanum Lycopersicum). Using these protein sequences, phylogenetic tree was constructed in MEGAX (Pennsylvania State University, State College, PA, USA) with the Neighbor-joining (NJ) algorithm setting 1000 bootstrap replicates [70]. Jalview was employed to perform multiple sequence alignment, and the muscle method with default parameter was adopted [71].

Transient expression in tobacco

Four-weeks-old tobacco leaves were used to perform transient transformation for observing VaPYL9 expression location in plant cell. The first strand cDNA of grape phloem was employed as a template to clone VaPYL9 gene. The purified cDNA fragment was combined with pART-CAM-EGFP vector to form a fusion protein vector. The fusion plasmid with pART-CAM-VaPYL9-EGFP was introduced into Agrobacterium tumefaciens strain GV3101. pART-CAM-VaPYL9-EGFP Agrobacterium tumefaciens liquid was injected into tobacco leaves back, and then was incubated in the dark at 26 °C for 12 h. After dark incubation, injected tobaccos were placed in the light at 26 °C for 24 h. At last, green fluorescence was motivated in tobaccos and tobaccos leaves were snapped by confocal microscopy using laser scanning confocal microscope (Olympus FV1000 viewer) [72].

Y2H protein binding assay

The CDS sequence of VaPYL9 was cloned into pGBKT7 and fused to the DNA-binding domain to generate a bait construct. Primer sequences used in these assays are listed in Supplementary Table S1. The recombinant plasmid was transformed into Y2H Gold yeast strain, which was cultured on double dropout (DDO) medium (SD/ − Leu/ − Trp) for 3 days at 30 ℃ to detect self-activation and toxicity. PCR was used to verify the presence of VaPYL9 in the Y2H Gold strain. The Mating hybridization method was adopted to screen potential interaction proteins from the AD library related to cold treatment which was constructed and preserved in Physiology and Biotechnology Laboratory (Gansu Agricultural University, China). Screened potential protein sequences were cloned into pGADT7 and fused to activation domain sequences to generate prey construct. Self-activation and toxicity detection were detected using the same method. The two constructs recombinant plasmids were co-transformed into Y2H Gold strain yeast cells and then cultured on a double dropout (DDO) medium (SD/ − Leu/ − Trp) at 30 ℃ for 3 days. After Co-transformation, the colonies were transferred to a quadruple dropout (QDO) medium (SD/ − Leu/ − Trp/ − Ade/ − His) with 10 mM 5-Bromo-4-chloro-3-indolyl-α-D-pyrangalactoside (X-α-Gal) to examine potential physical binding and interaction. Combinations of pGADT7-T with pGBKT7-53 and pGBKT7-Lam were served as positive and negative controls, respectively. All plates were incubated at 30 °C for 4 days and then photographed.

Cloning of VaPYL9 and obtaining of transgenic tomatoes

Based on our project's previous transcriptome sequencing on Amur grape phloem under cold stress, VaPYL9 was screened and identified as a key gene from grape in cold stress. To further verify VaPYL9 gene function, the full-length VaPYL9 cDNA was cloned using grape phloem cDNA as a template into the overexpression vector pCAMBIA1301 digested by NcoI and BstEII enzymes, and Escherichia coalition (DH5α) carried recombinant plasmid was sequenced by Sangon Biotech (Shanghai, China) Co. Ltd. After successfully sequenced, the recombinant plasmid was transformed into Agrobacterium tumefaciens strain GV3101. Agrobacterium tumefaciens strain GV3101 with overexpression pCAMBIA1301-VaPYL9 were transformed into WT tomato plants by the Agrobacterium-mediated leaf disk transformation method as described previously [73]. Gained sterile tomato cotyledons were cut as infection materials. Cut cotyledons were planted in MS (30 g/L sucrose + 6.5 g/L + 1 mg/L KT) to preincubate in dark condition for 2 days. Prepared infection bacteria solution using 1/2 MS suspension was used to immerse preincubated cotyledons for 5 min. Then infected cotyledons were coincubated for 3 days in dark condition. After cocultivation, they were transferred to generate callus and adventitious buds in MS (30 g/L sucrose + 6.5 g/L agar + 0.2 mg/L IAA + 2 mg/L ZT + 350 mg/L Timentin + 30 g/L hygromycin (Hyg). When the length of adventitious buds was about to 3 cm, cutting adventitious buds were transferred into MS (30 g/L sucrose + 6.5 g/L agar + 0.05 mg/L IAA + 350 mg/L Timentin + 30 g/L Hyg) to form roots. Transgenic plant lines were screened for 30 mg/L Hyg resistance and identified by the special primer according to instruction of FTL MIX (MF848) kit (Mei5 Biotechnology, Beijing, China Co., Ltd). Homozygous transgenic plants were selected in T1 generation derived from T0, and were stably maintained up to T3 generations. Ultimately, homozygous and stable T3 lines were gained.

Cold assay

Three stable T3 transgenic lines were selected from previously gained tomatoes. By cooling plant culture chamber up to 4 °C step by step, some growth well and consistent six-week-old tomato plants including transgenic plants and WT plants were selected and placed in the culture chamber keeping 4 °C for 48 h. Tomato and grape leaves were collected in 0 h, 6 h, 12 h, 24 h, 48 h, quickly frozen in liquid nitrogen and stored at − 80 °C for qRT-PCR and determination of physiological indexes related to cold tolerance. Three biological replicates were applied in this experiment.

Determination of H2O2, MDA, proline contents and antioxidant enzyme activities

H2O2 and malondialdehyde (MDA) content and catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activities were performed as described in the instruction of the commercial ELISA kit purchased from Suzhou Comin Biotechnology (Suzhou, Jiangsu, China) and reaction mixtures were measured in Multimode Reader (Spark® Multifunctional microplate detector, TECAN, Pudong, Shanghai). The proline contents were determined as previously described [74].

Determination of relative electrolyte leakage and content of IAA and ABA

The relative electrolyte leakage (REL) was measured using a model DDS307A device (Shanghai Leica Instrument) and the ratio of C1 (value of Leaf before boiling) to C2 (value of Leaf after boiling) was calculated as previously described [75]. 20 round leaves were gained using punches, and placed in glass test tubes with 15 mL ultrapure water. After 12 h, C1 was detected by a conductivity instrument. Then boiling for 15 min, C1 was detected. C1/C2% was viewed as relative electrical conductivity value. The content of IAA and ABA was determined by High Performance Liquid Chromatography (HPLC) mentioned in previous experiments [76]. Transgenic and WT Tomatoes leaves were ground and extracted using 5 mL 80% (V/V) chromatographic methanol at 4 ℃ for 12 h. During extraction, samples were constantly shaken to make hormone sufficiently dissolve. After centrifuged at 8000 rpm, 4 ℃,10 min, supernatants were transferred in new tubes. Repeat steps above once, supernatants were transferred in the same tube and constant volume was up to 10 mL. 2 mL supernatants were evaporated to dryness using a rotary evaporator. After evaporated to dryness, residue samples were thawed using 2 mL 50% (V/V) chromatographic methanol. HPLC condition was listed: Symmetry C18 chromatographic column (4.6 mm* 250 mm、5 μm), flowing phase (10% chromatographic methanol and 90%, 0.1% phosphoric acid), flow velocity (1 mL/min), sample volume (10 μL), determination waves length (254 nm), the temperature column (30 ℃). Hormones were detected in American Waters Acquity Arc high performance liquid chromatography.

The DAB and NBT staining

3, 3, 9-Diaminobenzidine (DAB) and Nitrotetrazolium Blue Chloride (NBT) staining methods were performed as previously described [76, 77]. 10 mM 2-(N-Morpholino) ethanesulfonic acid (MES) was dissolved in distilled water and PH was 5.5. 1 mg/L DAB (Sigma-Aldrich) was dissolved in MES Buffer Solution and PH was 5.5. 0.5 mg/L NBT (Sigma-Aldrich) was dissolved in distilled water and PH was 5.8. Prepared DAB and NBT stain solution were stored at 4 ℃ and kept in a dark place. The destaining solution was ethanol: lactic acid: glycerin (3:1:1). Detached Leaves of tomatoes were immersed in culture dishes containing DAB and NBT staining solution for 12 h. The stained leaves were boiled with 95% ethanol for 10 min and leaves were decolorized with a destaining solution for 2 h. Decolorized leaves were placed in culture dishes with distilled water and photoed in EPSON Scan of root scanner (Regent company, Canada).

RNA extraction and qRT-PCR

Total RNA was extracted from the leaves of grape and transgenic and WT tomato plants through RNAprep Pure Plant Kit (TIANGEN, Beijing, China). RNA was reversed as cDNA using the Prime Reverse Transcriptase Kit (Takara Bio, Shiga, Japan). Reversed cDNA was stored at − 20 °C and used as templates to amplify genes and perform qRT-PCR. All primers were designed on an online website and synthesized by Sangon Biotech (Shanghai) Co. Ltd (Supplementary Table S2). qRT-PCR was performed through the Light Cycler® 96 real-time PCR system (Roche, Switzerland). The VvNAPDH (GenBank accession no.CB973647) and SlActin (GenBank accession U60480) gene were used as internal controls. The 20 µL reaction solution contained 2 µL cDNA, 2 µL primers, 10 µL SYBR Green Master Mix Reagent (Takara Bio, Shiga, Japan), and 6 µL ddH2O. PCR amplification conditions were set as 50 cycles of 95 °C for 30 s, 95 °C for 10 s, 58 °C for 30 s, 72 °C for 20 s. The relative expressions of the genes were obtained from three biological replicates and calculated using the 2Ct method [78].

RNA-seq library preparation and sequencing in tomatoes under cold stress

Based on comprehensive data analysis in tomatoes after cold stress, transgenic and WT tomato leaves at 4 ℃ for 48 h were screened for RNA-seq. RNA extraction was conducted as described above. RNA-seq was commissioned by Biomarker Technologies Co., Ltd. (Beijing, China). Three biological replicates of each sample were used for RNA-seq analysis. RNA quality and concentration were determined by agarose gel electrophoresis using NanoPhotometer spectrophotometer (Implant, Germany) and Agilent 2100 BioAnalyzer (Agilent Technologies, USA). RNA samples with inspection qualified were used to construct the library. Eukaryotic mRNA was enriched by magnetic beads with Oligo (dT). The mRNA was randomly interrupted by Fragmentation Buffer. Using mRNA as a template, the first cDNA strand was synthesized using random hexamers, and then the second cDNA strand was synthesized by adding buffer, dNTPs, RNase H and DNA polymerase I. The cDNA strand was purified using AMPure XP Beads. The purified double-stranded cDNA was then terminally-repaired, added Poly-A and linked to a sequenced adaptor. AMPure XP Beads were then made segment size selection. Finally, the cDNA library was obtained by PCR enrichment. The libraries were sequenced on an Illumina the Novaseq6000 platform. FPKM value of each gene was calculated to estimate gene expression levels [79]. Gene expression difference was estimated using the DESeq2 [80]. The differentially expressed genes (DEGs) were screened by setting false discovery rate (FDR) < 0.05 and |Log2 FC (fold change)|≥ 2 as thresholds screening measurement. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis of DEGs was implemented by the cluster Profiler [81, 82].

Statistical analysis

All data were determined in three independent biological replicates for each experiment. Data analyses were done by using one-way analysis of variance (ANOVA) in SPSS 22.0 (SPSS, Inc., Chicago, IL, USA). Duncan’s multiple range tests were employed to test significant differences (P < 0.01). Origin 9.0 was used to draw diagrams.

Availability of Data and materials

All generated data in this study are included in supplementary information files. RNA-seq raw data have been deposited in NCBI ( under BioProject PRJNA830676.



Abscisic acid


Pyrabactin Resistance1/PYR1-Like /Regulatory Components of ABA receptors


Yeast two-hybrid




Quantitative real-time PCR


Transcriptome sequencing




Type 2C protein phosphates


SNF1-related protein kinase 2


Murashige and Skoog


Transgenic generation 3


Isoelectric point


Grand average of hydropathicity




Double dropout


Quadruple dropout










Superoxide dismutase




Relative electrolyte leakage


High Performance Liquid Chromatography


3, 3, 9-Diaminobenzidine


Nitrotetrazolium Blue Chloride


2-(N-Morpholino) ethanesulfonic acid


Differentially expressed genes


False discovery rate




Gene Ontology


Kyoto encyclopedia of genes and genomes


Protein-L-isoaspartate O-methyltransferase 1


  1. Zhu J. Abiotic stress signaling and responses in plants. Cell. 2016;167(2):313–24.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. Qin F, Shinozaki K, Yamaguchi-Shinozaki K. Achievements and challenges in understanding plant abiotic stress responses and tolerance. Plant Cell Physiol. 2011;52(9):1569–82.

    CAS  Article  PubMed  Google Scholar 

  3. Bailey-Serres J, Lee SC, Brinton E. Waterproofing crops: effective flooding survival strategies. Plant Physiol. 2012;160(4):1698–709.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. Gokul A, Niekerk LA, Carelse MF, Keyster M. Transgenic technology for efficient abiotic stress tolerance in plants (Chapter 5). In: Kiran U, Abdin M, Kamaluddin, editors. Transgenic Technology Based Value Addition in Plant Biotechnology. Amsterdam: Elsevier; 2020. P. 95-122.

  5. Lata C, Shivhare R. Engineering cereal crops for enhanced abiotic stress tolerance. P Indian Natl Sci Ac. 2021;87(1):6–83.

    Article  Google Scholar 

  6. Guan Y, Liu S, Wu W, Hong K, Shi J. Genome-wide identification and cold stress-induced expression analysis of the CBF gene family in liriodendron chinense. J Forestry Res. 2021;32:2531–43.

    CAS  Article  Google Scholar 

  7. Koornneef M, Bentsink L, Hilhorst H. Seed dormancy and germination. Curr Opin Plant Biol. 2002;5(1):33–6.

    CAS  Article  PubMed  Google Scholar 

  8. Parent B, Hachez C, Redondo E, Simonneau T, Chaumont F, Tardieu F. Drought and abscisic acid effects on aquaporin content translate into changes in hydraulic conductivity and leaf growth rate: a trans-scale approach. Plant Physiol. 2009;149(4):2000–12.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Fujii H, Chinnusamy V, Rodrigues A, Rubio S, Antoni R, Park SY, Zhu JK. In vitro reconstitution of an abscisic acid signalling pathway. Nature. 2009;462(7273):660–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. Ma Y, Szostkiewicz I, Korte A, Moes D, Yang Y, Christmann A. Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science. 2009;324(5930):1064–8.

    CAS  Article  PubMed  Google Scholar 

  11. Park SY, Fung P, Nishimura N, Jensen DR, Fujii H, Zhao Y, Lumba S, Santiago J, Rodrigues A, Chow TF, Alfred SE, Bonetta D, Finkelstein R, Provart NJ, Desveaux D, Rodriguez PL, McCourt P, Zhu JK, Schroeder JI, Volkman BF, Cutler SR. Abscisic acid inhibits type 2c protein phosphatases via the PYR/PYL family of START proteins. Science. 2009;324(5930):1068–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. Umezawa T, Nakashima K, Miyakawa T, Kuromori T, Tanokura M, Shinozaki K, Yamaguchi-Shinozaki K. Molecular basis of the core regulatory network in ABA responses: sensing signaling and transport. Plant Cell Physiol. 2010;51(11):1821–39.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. Shi Y, Yang S. ABA regulation of the cold stress response in plants (Chapter 17). In: Zhang DP, editor. Abscisic acid: metabolism transport and signaling. New York: Springer; 2014. p. 337–363.

  14. Gonzalez-Guzman M, Pizzio GA, Antoni R, Vera-Sirera F, Merilo E, Bassel GW, Fernandez MA, Holdsworth MJ, Perez-Amador MA, Kollist H, Rodriguez PL. 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(6):2483–96.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. Zhao Y, Zhang Z, Gao J, Wang P, Hu T, Wang Z, Liu W, Xie X, Lu T, Xue L, Liu Y, Macho AP, Tao WA, Bressan RA, Zhu J. Arabidopsis duodecuple mutant of PYL ABA receptors reveals PYL repression of ABA-independen SnRK2 activity. Cell Rep. 2018;23(11):3340–51.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. Radauer C, Lackner P, Breiteneder H. The Bet v 1 fold: An ancient versatile scaffold for binding of large hydrophobic ligands. BMC Evo Biol. 2008;8(1):286.

    CAS  Article  Google Scholar 

  17. Klingler JP, Batelli G, Zhu JK. ABA receptors: the START of a new paradigm in phytohormone signalling. J Exp Bot. 2010;61(12):3199–210.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. Raghavendra AS, Gonugunta VK, Christmann A, Grill E. ABA perception and signalling. Trends Plant Sci. 2010;15(7):395–401.

    CAS  Article  PubMed  Google Scholar 

  19. Cutler SR, Rodriguez PL, Finkelstein RR, Abrams SR. Abscisic acid: emergence of a core signaling network. Annu Rev Plant Biol. 2010;61:651–79.

    CAS  Article  PubMed  Google Scholar 

  20. Yuan X, Yin P, Hao Q, Yan C, Wang J, Yan N. Single amino acid alteration between valine and isoleucine determines the distinct pyrabactin selectivity by PYL1 and PYL2. J Biol Chem. 2010;285(37):28953–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. Nishimura N, Sarkeshik A, Nito K, Park S-Y, Wang A, Carvalho PC, Lee S, Caddell DF, Cutler SR, Chory J, Yates JR. SchroederJI, PYR/PYL/RCAR family members are major in vivo ABI1 protein phosphatase 2C-interacting proteins in Arabidopsis. Plant J. 2009;61(2):290–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Antoni R, Gonzalez-Guzman M, Rodriguez L, Peirats-Llobet M, Rodriguez PL. PYRABACTIN RESISTANCE1-LIKE8 plays an important role in the regulation of abscisic acid signaling in root. Plant Physiol. 2013;161:931–41.

    CAS  Article  PubMed  Google Scholar 

  23. Xing L, Zhao Y, Gao J, Xiang C, Zhu JK. The ABA receptor PYL9 together with PYL8 plays an important role in regulating lateral root growth. Scientific Rep. 2016;6(1):27177.

    CAS  Article  Google Scholar 

  24. Kai W, Wang J, Liang B, Fu Y, Zheng Y, Zhang W, Leng P. Role of an ABA receptor SlPYL9 in tomato fruit ripening. J Ex Bot. 2019;70(21):6305–19.

    CAS  Article  Google Scholar 

  25. Zhang Q, Kong X, Yu Q, Ding Y, Li X, Yang Y. Responses of PYR/PYL/RCAR ABA receptors to vontrasting stresses, heat and cold in Arabidopsis. Plant Signaling&Beha. 2019;4(12):1670596.

    CAS  Article  Google Scholar 

  26. 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(7):1938.

    CAS  Article  PubMed Central  Google Scholar 

  27. Lenka SK, Muthusamy SK, Chinnusamy V, Bansal KC. Ectopic expression of rice PYL3 enhances cold and drought tolerance in Arabidopsis thaliana. Molecular Biotech. 2018;60(5):350–61.

    CAS  Article  Google Scholar 

  28. Verma RK, Santosh Kumar VV, Yadav SK, Pushkar S, Rao MV, Chinnusamy V. Overexpression of ABA Receptor PYL10 gene confers drought and cold tolerance to Indica rice. Front Plant Sci. 2019;10:1488–1488.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Lim CW, Lee SC. ABA-Dependent and ABA-Independent functions of RCAR5/PYL11 in response to cold stress. Front Plant Sci. 2020;11:620.

    Article  Google Scholar 

  30. Yu J, Ge H, Wang X, Tang R, Wang Y, Zhao F, Lan W, Luan S, Yang L. Overexpression of Pyrabactin Resistance-Like Abscisic Acid receptors enhances drought, osmotic, and cold tolerance in transgenic poplars. Front Plant Sci. 2017;8:1752.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Santiago J, Rodrigues A, Saez A, Rubio S, Antoni R, Dupeux F, Rodriguez PL. Modulation of drought resistance by the abscisic acid receptor PYL5 through inhibition of clade A PP2Cs. Plant J. 2009;60(4):575–88.

    CAS  Article  PubMed  Google Scholar 

  32. Quan W, Hu Y, Mu Z, Shi H, Chan Z. Overexpression of AtPYL5 under the control of guard cell specific promoter improves drought stress tolerance in Arabidopsis. Plant Physiol Bioch. 2018;129:150–7.

    CAS  Article  Google Scholar 

  33. Lim CW, Baek W, Han SW, Lee SC. Arabidopsis PYL8 plays an important role for ABA signaling and drought stress responses. Plant Pathol J. 2013;29(4):471–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Yang J, Wang M, Zhou S, Xu B, Chen P, Ma F, Mao K. The ABA receptor gene MdPYL9 confers tolerance to drought stress in transgenic apple (Malus domestica). Environ Exp Bot. 2021;194:104695.

  35. Liang C, Liu Y, Li Y, Meng Z, Yan R. Zhu Tao, Wang Y, Kang S, Ali Abid M, Malik W, Sun G, Guo S, Zhang Rui, Activation of ABA Receptors gene GhPYL9-11A is positively correlated with cotton drought tolerance in transgenic Arabidopsis. Front Plant Sci. 2017;8:1453.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Liu J, Zhao FL, Guo Y, Fan X, Wang Y, Wen YQ. The ABA receptor-like gene VyPYL9 from drought-resistance wild grapevine confers drought tolerance and ABA hypersensitivity in Arabidopsis. Plant Cell Tiss Org. 2019;138:543–58.

    CAS  Article  Google Scholar 

  37. Zhao Y, Qi G, Ren F, Wang Y, Wang P, Wu X. Analysis of PYL genes and their potential relevance to stress tolerance and berry ripening in grape. J Amer Soc Hort Sci. 2020;145(5):308–17.

    CAS  Article  Google Scholar 

  38. Li G, Xin H, Zheng XF, Li S, Hu Z. Identification of the abscisic acid receptor VvPYL1 in Vitis vinifera. Plant Biol. 2011;14(1):244–8.

    CAS  Article  PubMed  Google Scholar 

  39. Gao Z, Li Q, Li J, Chen Y, Luo M, Li H, Ma C. Characterization of the ABA receptor VlPYL1 that regulates anthocyanin accumulation in grape berry skin. Front Plant Sci. 2018;9:592.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Siddiqua MK. Functional analysis of grape CBF genes. [Ph.D. Thesis.], Oran: University of Guelph. 2007.

  41. Chang D, Yang Q, Wang B, Lu H, Feng Y. COR-like gene is involved in induced-expression response to multiple abiotic stresses in grape vine (vitis amurensis) tissues. Agr Sci China. 2014;5(7):604–10.

    CAS  Article  Google Scholar 

  42. Sudesh KY. Cold stress tolerance mechanisms in plants. A review Agron. 2010;30(3):515–27.

    Article  Google Scholar 

  43. Theocharis A, Clement C, Barka EA. Physiological and molecular changes in plants grown at low temperatures. Planta. 2012;235(6):1091–105.

    CAS  Article  PubMed  Google Scholar 

  44. Zhou BY, Guo ZF. Effect of ABA and its biosynthesis inhibitor on chilling resistance and anti-oxidant enzymes activity. Acta Prata Sci. 2005;14(6):94–6.

    Google Scholar 

  45. He HY, Xue L, Tian LP, Chen YL. Effect of low temperature stress on the chlorophyll contents and chlorophyll fluorescence parameters in muskmelon seedling leaves. Northern Hort. 2008;4:121–7.

    Google Scholar 

  46. Huang X, Shi H, Hu Z, Liu A, Amombo E, Chen L, Fu J. ABA is involved in regulation of cold stress response in bermudagrass. Front Plant Sci. 2017;8:1613.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Battal P, Erez ME, Turker M, Berber I. “Molecular and Physiological Changes in Maize (Zea mays) Induced by Exogenous NAA, ABA and MeJa during Cold Stress. Ann Bot Fenn. 2008;45(3):173–85.

    Article  Google Scholar 

  48. Zhang Y, Yu H, Yang X, Li Q, Ling J, Wang H, Jiang W. CsWRKY46, a WRKY transcription factor from cucumber, confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA-dependent manner. Plant Physiol Bioch. 2016;108:478–87.

    CAS  Article  Google Scholar 

  49. Fennell A. Freezing tolerance and injury in grapevines. J Crop Improv. 2004;10(1–2):201–35.

    Article  Google Scholar 

  50. Di F, Jian H, Wang T, Chen X, Ding Y, Du H, Lu K, Li J, Liu L. Genome-Wide analysis of the pyl gene family and identification of pyl genes that respond to abiotic stress in Brassica napus. Genes. 2018;9(3):156.

    CAS  Article  PubMed Central  Google Scholar 

  51. Yadav SK, Santosh Kumar VV, Verma RK, Yadav P, Saroha A, Wankhede DP, Chinnusamy V. Genome-wide identification and characterization of ABA receptor PYL gene family in rice. BMC Genomics. 2020;21(1):676.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. Yin P, Fan H, Hao Q, Yuan X, Wu D, Pang Y, Yan C, Li W, Wang J, Yan N. Structural insights into the mechanism of abscisic acid signaling by PYL proteins. Nat Struct Mol Biol. 2009;16(12):1230–6.

    CAS  Article  Google Scholar 

  53. Tsujishita Y, Hurley J. Structure and lipid transport mechanism of a StAR-related domain. Nat Struct Mol Biol. 2000;7:408–14.

    CAS  Article  Google Scholar 

  54. Saavedra X, Modrego A, Rodríguez D, Gonzálezgarcía MP, Sanz L, Nicolás G, Lorenzo O. The nuclear interactor PYL8/RCAR3 of Fagus sylvatica FsPP2C1 is a positive regulator of abscisic acid signaling in seeds and stress. Plant Physiol. 2010;152(1):133–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. Cimmino A, Capasso R, Muller F, Sambri I, Masella L, Raimo M, Ingrosso D. Protein isoaspartate methyltransferase prevents apoptosis induced by oxidative stress in endothelial cells: role of BCL-XL deamidation and methylation. PLoS ONE. 2008;3(9):e3258.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. Zhao Y, Chan Z, Gao J, Xing L, Cao M, Yu C, Shi H, Zhu Y, Gong Y, Mu Z, Wang H, Deng X, Bressan RA, Zhu JK. ABA receptor PYL9 promotes drought resistance and leaf senescence. P Natl A Sci India B. 2016;113(7):1949–54.

    CAS  Article  Google Scholar 

  57. Yildirim K, Kaya Z. Gene regulation network behind drought escape, avoidance and tolerance strategies in black poplar (Populus nigra L). Plant Physiol Biochem. 2017;115:183–99.

    CAS  Article  PubMed  Google Scholar 

  58. Apel K, Hirt H. Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004;55:373–99.

    CAS  Article  PubMed  Google Scholar 

  59. Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.

    CAS  Article  PubMed  Google Scholar 

  60. Sucu S, Yagci A, Yildirim K. Changes in morphological, physiological traits and enzyme activity of grafted and ungrafted grapevine rootstocks under drought stress. Erwerbs-obstbau. 2018;60(2):127–36.

    Article  Google Scholar 

  61. Yu JL, Yang L, Liu XB, Tang RJ, Wang Y, Ge HM, Wu MT, Zhang J, Zhao FG, Luan S, Lan WZ. Overexpression of poplar pyrabactin resistance-like abscisic acid receptors promotes abscisic acid sensitivity and drought resistance in transgenic Arabidopsis. PLoS ONE. 2016;11(12):17.

    CAS  Article  Google Scholar 

  62. Kalladan R, Lasky JR, Sharma S, Kumar MN, Juenger TE, Des Marais DL, Verslues PE. Natural variation in 9-Cis-Epoxycartenoid dioxygenase 3 and ABA accumulation. Plant Physiol. 2019;179(4):1620–31.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  63. Skubacz A, Daszkowska-Golec A, Szarejko I. The role and regulation of ABI5 (ABA-Insensitive 5) in plant development, abiotic stress responses and phytohormone crosstalk. Front Plant Sci. 2016;7:1884.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Huang X, Chen MH, Yang LT, Li YR, Wu JM. Effects of exogenous abscisic acid on cell membrane and endogenous hormone contents in leaves of sugarcane seedlings under cold stress. Sugar Tech. 2014;17(1):59–64.

    CAS  Article  Google Scholar 

  65. De Zelicourt A, Colcombet J, Hirt H. The Role of MAPK modules and ABA during abiotic stress signaling. Trends Plant Sci. 2016;21(8):677–85.

    CAS  Article  PubMed  Google Scholar 

  66. Liu Y. Roles of mitogen-activated protein kinase cascades in ABA signaling. Plant Cell Rep. 2011;31(1):1–12.

    CAS  Article  PubMed  Google Scholar 

  67. Aleman F, Yazaki J, Lee M, Takahashi Y, Kim AY, Li Z, Kinoshita KE. JR Schroeder JI, An ABA-increased interaction of the PYL6 ABA receptor with MYC2 Transcription Factor: A putative link of ABA and JA signaling. Sci Rep. 2016;6(1):28941.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. Cheng MC, Liao PM, Kuo WW, Lin TP. The Arabidopsis ETHYLENE RESPONSE FACTOR1 regulates abiotic stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals. Plant Physiol. 2013;162(3):1566–82.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. Schmidt R, Mieulet D, Hubberten HM, Obata T, Hoefgen R, Fernie AR, Fisahn J, San Segundo B, Guiderdoni E, Schippers JHM, Mueller-Roeber B. SALT-RESPONSIVE ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell. 2013;25(6):2115–31.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  70. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGAX: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. Waterhouse AM, Procter JB, Martin D, Clamp M, Barton GJ. Jalview version 2: a multiple sequence alignment and analysis workbench. Bioinformatics. 2009;25(9):1189–91.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  72. Cheng W, Zhang H, Zhou X, Liu H, Liu Y, Li J, Han S, Wang Y. Subcellular localization of rice hexokinase (oshxk) family members in the mesophyll protoplasts of tobacco. Biol Plantarum. 2011;55(1):173–7.

    CAS  Article  Google Scholar 

  73. Hou YL, Meng K, Han Y, Ban QY, Wang B, Suo JT. The persimmon 9-lipoxygenase gene DkLOX3 plays positive roles in both promoting senescence and enhancing tolerance to abiotic stress. Front Plant Sci. 2015;6:1073.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Doneva D, Pál M, Brankova L, Szalai G, Tajti J, Khalil R. Ivanovska, Beti; Velikova V, Misheva S, Janda T, Peeva V, The effects of putrescine pre-treatment on osmotic stress responses in drought-tolerant and drought-sensitive wheat seedlings. Physiol Plantarum. 2021;171:200–16.

    CAS  Article  Google Scholar 

  75. Ma X, Chen C, Yang M, Dong X, Lv W, Meng Q. Cold-regulated protein (SlCOR413IM1) confers chilling stress tolerance in tomato plants. Plant Physiol Biochem. 2018;124:29–39.

    CAS  Article  PubMed  Google Scholar 

  76. Watanabe N, Lam E. BAX inhibitor-1 modulates endoplasmic reticulum stress-mediated programmed cell death in Arabidopsis. J Biol Chem. 2008;283(6):3200–10.

    CAS  Article  PubMed  Google Scholar 

  77. Bournonville CFG, Diaz-Ricci JC. Quantitative determination of superoxide in plant leaves using a modified NBT Staining. Method. 2011;22(3):268–71.

    CAS  Article  Google Scholar 

  78. Livak K, Schmittgen T. Analysis of relative gene expression data using realtime quantitative and the 2Ct method. Methods. 2001;25(4):402–8.

    CAS  Article  PubMed  Google Scholar 

  79. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  80. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  81. Yu G, Wang LG, Han Y, He QY. Cluster Profiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  82. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(1):545–51.

    CAS  Article  Google Scholar 

Download references


We thank Prof. Jianming Xie (Gansu Agricultural University, China) for providing the root scanner to photo.


This work was supported by the FuXi Foundation of Gansu Agricultural University (No. Ganfx-03J02) and Youth Innovation and Entrepreneurship Talent Project of Longyuan (2018LYQN01).

Author information

Authors and Affiliations



The authors JM, BHC, and GJN contributed to the study conception and design. Material preparation, data collection and analysis were performed by GJN, GPL, YML, WFM, SXL, HMG, LLG. The first draft of the manuscript was written by GJN. All the authors read and approved the manuscript.

Corresponding author

Correspondence to Juan Mao.

Ethics declarations

Ethics approval and consent to participants

Plant materials (Vitis amurensis Rupr. var. ‘Zuoshan 1’) used in this study were from vineyard of Gansu Agricultural University. Micro-Tom tomato seeds purchased from Nanjing Fengshuo Horiculture Ltd., Company (Nanjing, Jiangsu, China). The laboratory experiments were conducted under logical legislation and permissions. All methods were carried out in accordance with relevant guidelines and regulations in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supplementary Table S1. Detailed bioinformatic information of VaPYL9 gene in grape.

Addirional file 2:

Supplementary Table S2. All quantitative real-time primers for the expression level.

Additional file 3:

Supplementary Fig S1. Multiply sequence alignment of PYL gene members among grape, tomato, Arabidopsis thaliana. Underline part was conserved domain (polyketide-cyc2) feature in PYL sequences. The length of VaPYL9 domain sequence was 1-178.

Additional file 4:

Supplementary Fig S2. The construction of VaPYL9 transient expression vector. The presence of VaPYL9 were verified including amplified cDNA (a), introduced into Escherichia coalition (b) and Agrobacterium tumefaciens strain GV3101 (c) by 1 % agarose gel eletrophoresis, respectively.

Additional file 5:

Supplementary Fig S3. The construct of VaPYL9 overexpression vector. The presence of VaPYL9 were verified including amplified cDNA (a), introduced into Escherichia coalition (b) and Agrobacterium tumefaciens strain GV3101 (c) by 1 % agarose gel eletrophoresis, respectively.

Additional file 6:

Supplementary Fig S4. The identification of transgenic tomatoes. Through Agrobacterium-mediated leaf disk transformation method, transgenic tomatoes were gained. Through screening, stable and homozygous lines (#2, #6, #9) were identified and named OE1, OE2, OE3 by 1 % agarose gel electrophoresis. These stable lines were used for the following experiment.

Additional file 7:

Supplementary Fig S5. The construction pGBKT7-VaPYL9 and pGADT7-VaPCMT. The presence of VaPYL9 (a) and VaPCMT (b) were verified including amplified cDNA, introduced into Y2H gold by 1 % agarose gel electrophoresis, respectively.

Additional file 8:

Supplementary Fig S6. Self-activation toxicity test of bait protein VaPYL9 and prey protein VaPCMT. Dental plaque didn’t show blue. They were verified that there were no self-activation toxicity.

Additional file 9:

Supplementary Fig S7. Co-transformation of VaPYL9 and VaPCMT. VaPYL9 and VaPCMT were together introduced into Y2H gold and its presence was identified in Y2H gold by 1 % agarose gel electrophoresis.

Additional file 10:

Supplementary Excel file. Annotation information of DEGs in heatmap.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nai, G., Liang, G., Ma, W. et al. Overexpression VaPYL9 improves cold tolerance in tomato by regulating key genes in hormone signaling and antioxidant enzyme. BMC Plant Biol 22, 344 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Abscisic acid
  • Overexpression VaPYL9
  • Cold stress
  • RNA-seq
  • Expression analysis