Genome-wide identification and characterization of PePYLs in poplar
We identified 12 putative PYLs from the NCBI P. euphratica genome based on 14 AtPYLs amino acid sequences and named PePYL1-PePYL12 according to sequence similarity. The characteristics of PePYLs, including gene name, gene ID, coding DNA sequence (CDS) length, theoretical isoelectric point (pI), molecular weight (MW), and protein length are presented in (Table S1).
To investigate the phylogenetic relationship between PePYLs and AtPYLs, the Neighbor-Joining method was used and the tree was constructed using MEGA 7 (Fig. 1a). The Gene Structure Display Server (GSDS v2.0) further analyzed the exon-intron structure of the PePYLs and AtPYLs gene (Fig. 1b). The conserved domain of the PePYLs was analyzed by DNAMAN software. The results showed that 12 PePYLs proteins can be divided into three subfamilies, among them, PePYL7 to PePYL11 belong to subfamily I; PePYL3 to PePYL6 belong to subfamily II; PePYL1, PePYL2, and PePYL12 belong to subfamily III. The amino acid sequences of 12 PePYLs all contained four highly conserved surface loops CL1-CL4 (Fig. 1c).
The promoter cis-acting elements analysis of PePYLs in poplar
To study the potential functions of PePYLs, the cis-acting elements on the PePYLs promoter sequence were analyzed by PlantCARE online (Table S3). The results showed that the 12 PePYLs promoter sequences mainly include five hormone response-related elements (Fig. 2): abscisic acid-responsive element, MeJA-responsive element, salicylic acid-responsive element, auxin-responsive elements and, gibberellin-responsive elements; three abiotic stress-related elements: defense and stress-responsive elements, low-temperature-responsive element and MYB binding site involved in drought-inducibility. Some PePYLs promoters contain elements related to growth and development, such as, circadian elements, which are involved in circadian control; MBSI elements, which are MYB binding sites involved in the regulation of flavonoid biosynthetic genes; HD-Zip 1 elements, which regulate the differentiation of the palisade mesophyll cells. In summary, PePYLs with different cis-acting elements may be involved in a variety of physiological processes.
Subcellular localization of PePYLs
To clarify the localization of PePYLs in cells, 35S:PePYLs-GFP (green fluorescent protein) fusion protein was transiently transfected into tobacco leaves. The results showed that: PePYL2, PePYL4, PePYL5, PePYL6, and PePYL9 are located in the cytoplasm and nucleus, the fluorescence signals of PePYL10 and PePYL12 were weak, among which PePYL10 is more obvious in the nucleus, while PePYL1, PePYL3, PePYL7, PePYL8, and PePYL11 did not detect the fluorescence signal, and its localization in cells cannot be determined (Fig. 3).
The expression patterns of PePYL6 and PePYL9
The expression of PePYL6 and PePYL9 was detected in different tissues (Fig. 4a, d). Previous studies showed that the transcription levels of PePYL6 and PePYL9 were up-regulated under ABA and mannitol treatment. To further determine the expression patterns of PePYL6 and PePYL9, we generated Arabidopsis plants with GUS driven by PePYL6 and PePYL9 promoter. The results showed that ABA and mannitol treatment enhanced β-glucuronidase (GUS) staining (Fig. 4b, e), GUS activity further confirmed GUS staining (Fig. 4c, f), PePYL6 and PePYL9 were induced by ABA and mannitol.
To further investigate the performance of the PYLs gene in response to drought stress, PePYL6 and PePYL9 were stably transformed into Arabidopsis. Col-0, pyl6 and pyl9 were used in this study. Transgenic Arabidopsis was verified by PCR and qRT-PCR, the high expression levels plants were selected for experiments (Fig. 4g-i) (Czechowski et al., 2004).
Overexpression of PePYLs increases ABA sensitivity of seed germination and root growth
As a potential ortholog of AtPYLs, the function of PePYLs in the process of seed germination and root growth regulated by ABA has been studied. The inhibition of seed germination by exogenous ABA changed with the concentration of ABA, and the higher ABA concentration, the lower seed germination rate (Fig. 5a). On the 1/2 MS medium, the Col-0, pyl6, and pyl9 plants germinated normally, while OEPePYLs plants showed delayed germination (Fig. 5b). On the medium supplemented with 0.5 µM ABA, the germination rates of Col-0, pyl6, and pyl9 plants were comparable, OEPePYLs plants were lower than those of Col-0 and mutant plants (Fig. 5c). After 24 h of growth on the medium supplemented with 1.0 µM ABA, the germination of OEPePYLs plants was about 15.35%, while the germination of Col-0 plants was 55.95%, and the difference of germination rate between Col-0 and OEPePYLs plants was significant with time (Fig. 5d). Therefore, overexpression of PePYL6 and PePYL9 increased the ABA sensitivity during seed germination, even without ABA treatment.
Root length is another phenotype for assessing ABA sensitivity under exogenous ABA treatment (Fig. 6a). The results showed the average root length of Col-0 plants decreased from 4.56 cm to 2.15 cm when grown on the 1/2 MS medium supplemented with 0, 5 and 10 µM ABA, while OEPePYLs plants root length decreased overall more than Col-0 plants, such as OEPePYL6 plants root length decreased from 4.21cm to 1.23 cm, the root length of OEPePYL9 plants decreased as much as that of OEPePYL6 plants, and when OEPePYL9 plants on the medium supplemented with 10 µM ABA, the average root length of OEPePYL9 plants was only 0.95 cm. In conclusion, overexpression of PePYL6 and PePYL9 increased the ABA sensitivity during root growth (Fig. 6b-d).
Overexpression of PePYLs improves the drought resistance of transgenic Arabidopsis
It has been reported that increasing ABA sensitivity enhances the drought resistance of plants (He et al., 2018). To investigate whether overexpression of PePYL6 and PePYL9 increases ABA sensitivity affects drought resistance of transgenic Arabidopsis. Col-0, OEPePYLs, and mutant plants were transplanted into the soil for 8 d without water. Before the drought stress, the phenotypes of Col-0, OEPePYLs, and mutant plants were not significantly different. After the drought stress, Col-0 and mutant plants showed more severe wilt, however, leaves of OEPePYLs plants did not wither, remained green (Fig. 7a). The physiological analysis showed that OEPePYLs plants had a higher net photosynthetic rate (Pn) than the Col-0 and mutant plants, mutant plants had the lowest net photosynthetic rate (Fig. 7b), the transpiration rate (Tr) of OEPePYLs transgenic plants were lower than that of Col-0 and mutant plants (Fig. 7c), which resulted in higher instantaneous water use efficiency (iWUE) (iWUE = Pn/Tr) of OEPePYLs transgenic plants under the same conditions (Fig. 7d).
The survive rate, water loss, leaf RWC, and Maximal PSII quantum yield (Fv/Fm) of all lines before and after the drought were measured. The survive rate of all lines was 95%-97% under the control condition, however, the survive rate of OEPePYLs plants was 70%-85%, higher than the Col-0 and mutant plants (Fig. 7e). The data of water loss of OEPePYLs plants showed that transgenic plants were lower than that of Col-0 and mutant plants (Fig. 7f), and the RWC and Fv/Fm values were also higher than other plants (Fig. 7g, h) under the drought stress. The RWC of Col-0 plants decreased from 82.38–53.24%, the RWC of pyl6 and pyl9 mutant plants was similar to that of Col-0 plants, while RWC of OEPePYLs plants decreased slightly. Drought-induced proline accumulation can stabilize the metabolic process of the plant under adverse conditions (Szabados and Savouré, 2010 ). After the drought stress, OEPePYLs plants had higher proline content than that of Col-0 and mutant plants (Fig. 7i). Therefore, overexpression of PePYL6 and PePYL9 enhances drought resistance of transgenic plants.
Overexpression PePYL6 and PePYL9 altered the expression of downstream genes
To study whether overexpression of PePYL6 and PePYL9 enhanced the drought tolerance of transgenic Arabidopsis by regulating the expression of downstream stress-related genes, we compared the expression of ABF2, RAB18, P5CS1, RD29A, and RD29B genes (Cutler et al., 2010) previously reported to be induced by drought and ABA in Col-0 and OEPePYLs plants under water and withholding water conditions. The data showed that the expression of all these genes in the transgenic plants was significantly higher than those of the Col-0 plants (Fig. 8), which suggested that PePYLs overexpressed altered the expression patterns of downstream stress-related genes and thus contributed to the enhancement of drought tolerance of transgenic plants.