Skip to main content

Advertisement

Genome-wide identification and characterization of ABA receptor PYL/RCAR gene family reveals evolution and roles in drought stress in Nicotiana tabacum

Abstract

Background

Abscisic acid (ABA) is an important phytohormone for plant growth, development and responding to stresses such as drought, salinity, and pathogen infection. Pyrabactin Resistance 1 (PYR1)/PYR1-Like (PYL)/Regulatory Component of ABA Receptor (RCAR) (hereafter referred to as PYLs) has been identified as the ABA receptors. The PYL family members have been well studied in many plants. However, the members of PYL family have not been systematically identified at genome level in cultivated tobacco (Nicotiana tabacum) and its two ancestors. In this study, the phylogenic relationships, chromosomal distribution, gene structures, conserved motifs/regions, and expression profiles of NtPYLs were analyzed.

Results

We identified 29, 11, 16 PYLs in the genomes of allotetraploid N. tabacum, and its two diploid ancestors N. tomentosiformis and N. sylvestris, respectively. The phylogenetic analysis revealed that NtPYLs can be divided into three subfamilies, and each NtPYL has one counterpart in N. sylvestris or N. tomentosiformis. Based on microarray analysis of NtPYL transcripts, four NtPYLs (from subfamily II, III), and five NtPYLs (from subfamily I) are highlighted as potential candidates for further functional characterization in N. tabacum seed development, response to ABA, and germination, and resistance to abiotic stresses, respectively. Interestingly, the expression profiles of members in the same NtPYL subfamily showed somehow similar patterns in tissues at different developmental stages and in leaves of seedlings under drought stress, suggesting particular NtPYLs might have multiple functions in both plant development and drought stress response.

Conclusions

NtPYLs are highlighted for important functions in seed development, germination and response to ABA, and particular in drought tolerance. This work will not only shed light on the PYL family in tobacco, but also provides some valuable information for functional characterization of ABA receptors in N. tabacum.

Background

Abscisic acid (ABA) is an important phytohormone that plays crucial roles in plant growth and development, including seed dormancy and germination, leaf senescence and abscission, bud development, fruit ripening and stomatal aperture control [1,2,3,4,5,6]. ABA is also a stress indicator for plant in response to stress, such as drought, salinity, and pathogen infection [3, 7,8,9]. Moreover, many evidences show that there are crosstalks and interactions between ABA and other phytohormones on plant development and responses to environmental cues [3, 10,11,12,13].

In plants, ABA is perceived by the ABA receptors to activate ABA signaling cascade. Recently, the Pyrabactin Resistance 1 (PYR1)/PYR1-Like (PYL)/Regulatory Component of ABA Receptor (RCAR) (hereafter referred to as PYLs) were identified as ABA receptors in Arabidopsis through genetic and biochemical screens [14,15,16]. The ABA-bound PYLs interact with the clade A protein phosphatase type 2Cs (PP2Cs) to prevent PP2Cs from inhibiting of SUCROSE NONFERMENTING 1 (Snf1)-related protein kinase 2s (SnRK2s). Consequently, the liberated SnRK2s could further phosphorylate the downstream substrates to switch on ABA signaling pathway [17,18,19,20,21,22,23,24,25]. Thus, PYLs, PP2Cs, and SnRK2s consist of the core ABA signaling components [3, 26,27,28,29].

Given the important functions of ABA receptors (PYLs) as the core regulators of ABA signaling in plants, many PYLs have been identified/functional characterized in Arabidopsis [8, 14,15,16, 30,31,32,33,34,35,36], rice [37,38,39,40], tomato [41, 42], soybean [43], wheat [44], maize [45, 46], poplar [47], rubber tree [48], strawberry [49], cotton [50], and Brachypodium distachyon [51, 52]. In the Arabidopsis genome, AtPYR1 encodes a member of the cyclase subfamily of the START (star-related lipid transfer) domain superfamily, and 13 genes that are similar to AtPYR1 were named as AtPYL (PYR1-like) 1 to AtPYL13. Thus, the Arabidopsis PYL family consists of 14 members (AtPYR1 and AtPYL1–13) with a feature of START domain. The 14 AtPYLs have been classified to 3 subfamilies [14, 15, 53]. Moreover, 14 ABA receptor genes were identified from the genome of rubber tree (Hevea brasiliensis), and phylogenetic analysis demonstrated that HbPYLs can be divided into three subfamilies [48]. In tomato (Solanum lycopersicum), eight SlPYLs were isolated, and can be clustered into three subfamilies according to the phylogenetic relationship [54]. In addition, gene structure analysis of the PYLs in H. brasiliensis and Brachypodium distachyon revealed that genomic sequences of HbPYLs and BdPYLs can be divided into two clades with or without introns [48, 52].

The Static Light Scattering (SLS) and Analytical Ultracentrifugation (AUC) analysis of recombinant AtPYLs in solution revealed that AtPYR1, AtPYL1, and AtPYL2 could form homodimers, while AtPYL4, AtPYL5, AtPYL6, AtPYL8, AtPYL9 and AtPYL10 are monomers [55]. In general, these monomeric AtPYLs interact with AtPP2Cs in an ABA-independent manner, while the dimeric AtPYLs bind to AtPP2Cs in an ABA-dependent manner. The ABA-dependence of AtPYLs-AtPP2Cs interactions are determined by a conserved region which called CL2 [55]. The CL2 region forms an ABA-binding pocket with the other three highly conserved surface loops, CL1, CL3, and CL4 [56]. Therefore, the conserved CL2 region in PYLs is crucial for ABA signaling transduction in plants.

Furthermore, transcription levels of several ABA receptor genes could be regulated in different tissues and by ABA and abiotic stresses. For example, the expression profiles of Gossypium hirsutum PYLs were tissue-specific [50]. Twelve BdPYLs were identified from the monocot model plant B. distachyon genome, and expression levels of BdPYL11 are significantly down-regulated in response to ABA, NaCl, and osmotic stress [52]. Moreover, expression levels of most SlPYLs, except SlPYL1 and SlPYL8, were also down-regulated by dehydration stress in tomato leaves [54].

Cultivated tobacco (Nicotiana tabacum) is not only an important economic crop in many countries, but also widely used as a model for plant biology research. It is well known that the allotetraploid N. tabacum (2n = 48, TTSS) was originated from a hybridization event between N. tomentosiformis (2n = 24, TT) and N. sylvestris (2n = 24, SS) [57, 58]. Therefore, the gene redundancy and diversity of NtPYL family might be more complex than those of AtPYL and OsPYL family in diploid Arabidopsis and rice.

Lots of studies on the ABA receptors have been carried out in many plants, mainly in Arabidopsis [8, 14,15,16, 30, 31, 36] and rice [37,38,39,40]. However, little is known about the ABA receptor family in cultivated tobacco [59]. In this study, we identified 29 NtPYLs in the allotetraploid N. tabacum, 11 NtomPYLs and 16 NsylPYLs in the two diploid wild tobacco species N. tomentosiformis and N. sylvestris, respectively. The phylogenic relationships, chromosomal distribution, gene structures, and the conserved motifs/regions in these PYLs were analyzed. Moreover, the expression profiles of NtPYLs in tissues at different developmental stages, and in response to drought stress were further investigated. These results provide some valuable information for further functional characterization of ABA receptors in N. tabacum.

Results

Genome-wide analysis of PYLs in three Nicotiana species

To identify the PYL family members in the Nicotiana tomentosiformis (2n = 24, TT), N. sylvestris (2n = 24, SS), and N. tabacum (2n = 48, TTSS) genomes, the coding sequences and amino acid sequences of 14 Arabidopsis AtPYLs were applied as queries to search against the NCBI database and China Tobacco Genome Database (V2.0) in China Tobacco Gene Research Centre at Zhengzhou Tobacco Research Institute. Based on amino acid sequence similarity to AtPYLs, a total of 11, 16, and 29 PYL genes were retrieved from the genomes of two progenitor diploid species N. tomentosiformis and N. sylvestris, and the descendant tetraploid specie N. tabacum, respectively (Additional files 1, 2, 3: Tables S1, S2, S3). As showed in Table 1, the gene and open reading frame (ORF) lengths of NtPYLs range from 465 to 3879 bp and 465 to 660 bp, respectively. For the deduced NtPYLs, the protein lengths vary from 154 to 219 amino acids, the molecular weight (MW) from 16.89 to 24.23 kDa, and the Isoelectric Point (pI) from 4.50 to 8.41. Notably, the NtPYL13 and NtPYL22 are the shortest proteins (154 aa) among the NtPYLs.

Table 1 Basic information of ABA receptor gene family in Nicotiana tabacum

The basic information of the NtomPYLs and NsylPYLs, including the gene ID, lengths of gene and ORF, MW, and pI are listed in Table 2 and Table 3, respectively. For the NtomPYLs, the gene lengths vary from 561 to 3275 bp, and ORF lengths range from 534 to 660 bp. For the deduced NtomPYLs, the protein lengths range from 177 to 219 amino acids, MW from 20.23 to 24.92 kDa, and pI from 4.73 to 7.97 (Table 2). For the NsylPYLs, the gene lengths vary from 630 to 3888 bp, and ORF lengths range from 528 to 657 bp. For the deduced NsylPYLs, the protein lengths range from 175 to 218 amino acids, MW from 19.59 to 24.11 kDa, and pI from 4.97 to 8.66 (Table 3). Analysis of the physical properties in the deduced Nicotiana PYL members showed that these PYLs are highly conserved. Most of the deduced Nicotiana PYLs have the similar amino acid lengths, MW, and theoretical pI. Majority of the putative Nicotiana PYLs possess 175–219 amino acids. MW and pI of the deduced Nicotiana PYLs range from 16.89 to 25.01 kDa, and 4.50 to 8.66, respectively.

Table 2 Basic information of ABA receptor gene family in Nicotiana tomentosiformis
Table 3 Basic information of ABA receptor gene family in Nicotiana sylvestris

Phylogenetic analysis of PYL family in Arabidopsis, soybean and N. tabacum

To characterize the phylogenetic relationship among PYLs from Arabidopsis, soybean, and cultivated tobacco, an unrooted neighbor-joining (NJ) tree was constructed using the MEGA software from the alignment of 29 NtPYLs in N. tabacum, 14 AtPYLs in Arabidopsis thaliana, and 23 GmPYLs in Glycine max (Fig. 1). The phylogenetic analysis of amino acid sequences from the deduced NtPYLs, PYLs in Arabidopsis and soybean revealed that the 29 NtPYLs could be grouped with their orthologous PYLs from Arabidopsis and soybean (Fig. 1). Thus, the NtPYLs were renamed as NtPYL1 to NtPYL29 according to their sequence similarities to AtPYLs.

Fig. 1
figure1

Phylogenetic analysis of PYL proteins from Arabidopsis, soybean and cultivated tobacco. A total of 14 AtPYLs from Arabidopsis (Arabidopsis thaliana), 23 GmPYLs from soybean (Glycine max), and 29 NtPYLs from cultivated tobacco (Nicotiana tabacum) were used to generate the unrooted neighbor-joining (NJ) tree with 1000 bootstrap replicates. The PYL proteins are classified into 3 subfamilies (marked as I, II, III), and distinguished by different colors: AtPYLs labeled in yellow, GmPYLs labeled in blue, and NtPYLs labeled in red

According to the created phylogenetic tree, these PYLs could be classified to 3 subfamilies. Subfamily I include the NtPYL1–18, AtPYR1, AtPYL1 to AtPYL6, and GmPYL1 to GmPYL14. Subfamily II contains NtPYL19 to NtPYL27, AtPYL7 to AtPYL10, and GmPYL15 to GmPYL20. Subfamily III consists of NtPYL28 and NtPYL29, AtPYL11 to AtPYL13, and GmPYL21 to GmPYL23 (Fig. 1). There are 18 NtPYLs, 7 AtPYLs, and 14 GmPYLs in PYL subfamily I; 9 NtPYLs, 4 AtPYLs, and 6 GmPYLs in PYL subfamily II; 2 NtPYLs, 3 AtPYLs, and 3 GmPYLs in PYL subfamily III. Therefore, PYL subfamily I and III have the largest and minimum number of members among the Arabidopsis, soybean, and cultivated tobacco, respectively.

Chromosomal distributions of NtPYLs

The localizations of the PYLs in the chromosomes of N. tabacum were further determined. The information of a physical maps (including the length and number of each chromosome, and gene locus) in each N. tabacum chromosome were obtained from the China Tobacco Genome Database (V2.0). A simplified physical map which shows the location of NtPYLs in the N. tabacum chromosomes was drawn by the software MapGene2Chromosome (Fig. 2). In general, NtPYLs were unevenly distributed in the N. tabacum chromosomes. For example, 24 NtPYLs were distributed in the 12 chromosomes, and four NtPYLs (NtPYL2, NtPYL23, NtPYL26, and NtPYL28) were located on the chromosome 3. Each of chromosome 16 and 18 contains 3 NtPYLs (NtPYL8, NtPYL11, NtPYL12 and NtPYL4, NtPYL5, NtPYL25, respectively), while four chromosomes (chromosome 4, 6, 23, and 24) harbor 2 NtPYLs, and only 1 NtPYL were located on the each of 4 chromosomes (chromosome 7, 9, 14, and 15). However, there are 5 NtPYLs that were located on the scaffolds could not be distributed on the N. tabacum chromosomes. In addition, NtomPYLs and NsylPYLs could only be located on the scaffolds but not on the chromosomes of N. tomentosiformis and N. sylvestris, due to the lacking information of physical maps for the N. tomentosiformis and N. sylvestris.

Fig. 2
figure2

The chromosomal location of the ABA receptor genes (NtPYLs) on the Nicotiana tabacum chromosomes. Chromosome size is indicated by its relative length. The scale on the left is in megabases (Mb). The bars on the chromosomes indicate the positions of the ABA receptor genes. The figure was generated and modified using the MapGene2Chrom program

Interestingly, according to the analysis of amino acid sequences similarity for PYLs among the N. tabacum, N. tomentosiformis, and N. sylvestris, there should be a tandem gene duplication (NtPYL11 and NtPYL12, both derive from Ntom0128520) and gene duplication (NtPYL4 and NtPYL5, both derive from Nsyl0370740) event on chromosome 16 and 18, respectively (Fig. 2, Table 4).

Table 4 ABA receptor gene family in Nicotiana tabacum and its putative ancestors N. sylvestris and N. tomentosiformis

Phylogenetic analysis of PYLs in three Nicotiana species

Cultivated tobacco (N. tabacum) is a tetraploid crop, and the genome of N. tabacum (TTSS) is likely derived from two different genomes of two wild diploid tobacco (T genome from N. tomentosiformis and S genome from N. sylvestris, respectively) [57, 58]. To analyze the phylogenetic relationship among PYLs from N. tomentosiformis, N. sylvestris, and N.tabacum, an unrooted neighbor-joining tree was constructed via the MEGA software by comparing 11 NtomPYLs, 16 NsylPYLs, and 29 NtPYLs (Fig. 3, left panel). The PYLs in these tobacco species could be divided into 3 subfamilies and further classified to 13 groups, and each NtPYL has a putative orthologous gene in either N. sylvestris or N. tomentosiformis (Fig. 3, Table 4).

Fig. 3
figure3

Phylogenetic relationship and gene structure of PYL family in three Nicotiana species. The unrooted neighbor-joining (NJ) evolutionary tree was constructed using MEGA 7.0 with 1000 bootstrap replicates based on the deduced full-length amino acid sequences of the PYLs in Nicotiana tabacum, N. tomentosiformis, and N. sylvestris (left panel). Exon-intron analyses of the identical Nicotiana PYL genes were performed with GSDS 2.0 (right panel). Exons and introns are represented by black rectangle and black lines, respectively. The lengths of the exons and introns for each Nicotiana PYL gene are shown proportionally

Notably, NtPYL4, NtPYL5 in group 2, and NtPYL22, NtPYL23 in group 10 were all derived from N. sylvestris, while NtPYL11, NtPYL12 in group 5 are originated from N. tomentosiformis. Namely, NtPYL4 and NtPYL5 in group 2 are derived from Nsyl0370740, and NtPYL11 and NtPYL12 in group 5 are derived from Ntom0128520, indicating that NtPYL4 and NtPYL5 might have the same origin from N. sylvestris, and NtPYL11 and NtPYL12 might share the same origin from N. tomentosiformis (Fig. 3, Table 4).

Gene structure of PYLs in three Nicotiana species

To better understand the Nicotiana PYLs’ structural features, organization of intron/exon was detected via alignment of genomic DNA and ORF sequences of Nicotiana PYL family according to the phylogenetic relationship of Nicotiana PYLs (Fig. 3, right panel). Fifteen of the 29 NtPYLs (NtPYL6–18, NtPYL28, NtPYL29) are intronless, six NtPYLs (NtPYL1–5, NtPYL22) have one intron, and 8 NtPYLs (NtPYL19–21, NtPYL23–27) have two introns. For the NtPYL subfamily I (NtPYL1 to NtPYL18), five NtPYLs (NtPYL1 to NtPYL5) have one intron, and other 13 members (NtPYL6 to NtPYL18) does not have intron. Most members in the NtPYL subfamily II (NtPYL19–27) have two introns, except NtPYL22 has only one intron. All members in the NtPYL subfamily III (NtPYL27 and NtPYL28) have no intron (Fig. 3, right panel).

Notably, only NtPYL22 in the NtPYL subfamily II has one intron, and amino acid length of NtPYL22 (154 aa) is the shortest among the ABA receptors of N. tabacum (Table 1). Through comparing the length of exons and introns of members (NtPYL19 to NtPYL27) in the NtPYL subfamily II, the first exon and intron might had lost in NtPYL22 gene.

Conserved motifs and CL2 regions/loops of NtPYLs

Amino acid alignment analysis revealed that all the identified NtPYLs share a highly similar helix-grip organization with three α-helices separated by 7 β-sheets, and several conserved CL regions/loops (Fig. 4a), which have been well characterized in the PYR/PYL/RCAR ABA receptor gene family [14, 15].

Fig. 4
figure4

Alignment and conserved motifs of NtPYLs. Amino acid sequence alignment of the 29 NtPYLs and AtPYL2 was performance by ClustalW. a Secondary structural elements are indicated above the primary sequence. Helices and sheets/strands are shown as black helices and arrows, respectively. The four conserved ABA receptor region CL1-CL4 are indicated with red lines. The conserved motifs analysis of the NtPYLs based on their phylogenetic relationship were identified using MEME software. b In left panel, the members of each subfamily are indicated with the same color and different NtPYL subfamilies are represented by the Roman numeral I-III in the phylogenetic tree. In right panel, grey lines represent non-conserved sequences, and colored boxes numbered at the bottom indicate different motifs. The length of motifs in each NtPYL protein is shown proportionally

The conserved motifs in 29 NtPYLs were predicted via MEME and Pfam software. In total, four motifs were identified in NtPYLs. One motif (motif 4) was identified only in 17 NtPYLs from subfamily I (NtPYL1 to NtPYL18) (Fig. 4b). Most of the 29 NtPYLs share the same motifs (motif 1, 2, and 3) suggesting their conserved functions in the NtPYLs. Among the 29 NtPYLs, only 2 NtPYLs (NtPYL13 from subfamily I and NtPYL22 from subfamily II) have two motifs (motif 1 and 3). Notably, NtPYL13 and NtPYL22 are the shortest receptors (154 aa) in the NtPYL family (Table 1). Moreover, most of the members in subfamily I, except NtPYL13, have the motif 4 at the N-terminal, indicating that most members in NtPYL subfamily I might carry out special biological functions.

Previous study showed that the conserved CL2 region/loop in PYLs is important for the interaction between PYLs and PP2Cs in an ABA dependent or independent manner in Arabidopsis and soybean [15, 43]. To identify CL2 region in the NtPYLs, alignment of AtPYLs, GmPYLs and NtPYLs was investigated. The conserved CL2 regions including 10 amino acid residues in the NtPYLs show a certain extent of similarity and polymorphism to those in AtPYLs, and GmPYLs (Fig. 5). The No. 3 and 4 residues in CL2 region are the amino acids that are critical for the monomeric or dimeric status, ABA dependence of PYL-PP2C interactions, and activities of AtPYLs [55]. In Arabidopsis, the combination of two key amino acid residues are VI and VV, VK and LK, VV and LV in AtPYL subfamily I, II and III, respectively. In Glycine max, the combination of two key amino acid residues are VI and VV, VK, IT and VT in GmPYL subfamily I, II and III, respectively. In Nicotiana tabacum, the combination of two key amino acid residues are VI and VV, VK and VR, LV in NtPYL subfamily I, II and III, respectively (Fig. 5).

Fig. 5
figure5

Alignment of the conserved CL2 regions in ABA receptors (PYLs) in Arabidopsis, soybean, and cultivated tobacco. Multiple amino acid sequence alignment of the full-length PYLs from Arabidopsis, soybean, and cultivated tobacco was carried out by ClustalW algorithm. Black line indicates the conserved CL2 loop/region in ABA receptor PYLs; asterisk indicates the position of the two key amino acids

Common patterns of two key amino acid residues in the CL2 region of AtPYLs, GmPYLs and NtPYLs could be illustrated as followed: VI and VV for subfamily I, VK for subfamily II, and LV is the common pattern in subfamily III of AtPYLs and NtPYLs. The patterns of two key amino acid residues in the CL2 region in AtPYLs, GmPYLs, and NtPYLs are conserved in subfamily I and II among these three plant species, but distinguished in subfamilies III, indicating that the members in different PYL subfamilies might have similar functions in different plants.

Expression profiles of NtPYLs in tissues at different developmental stages

To investigate the expression patterns of NtPYLs in different tissues and developmental stages, the relative expression levels of NtPYLs was analyzed in N. tabacum dry seeds, germination seeds, cotyledons, leaves and roots from two, four, six, and ten true leaves stages, and flowers at squaring stage by Microarray analysis (Fig. 6, Additional file 5: Table S5).

Fig. 6
figure6

The expression profile of 29 NtPYLs in tissues at different developmental stages. The relative transcript abundances of 29 NtPYLs were examined via microarray and visualized as a heatmap. The gene expression profiles of NtPYLs in 12 different samples, including dry seeds, germination seeds, cotyledons, leaves from two-true leaf stage (labeled as two true leaf_leaf), roots from two-true leaf stage (two true leaf_root), leaves from four-true leaf stage (four true leaf_leaf), roots from four-true leaf stage (four true leaf_root), leaves from six-true leaf stage (six true leaf_leaf), roots from six-true leaf stage (six true leaf_root), leaves from 10-true leaf stage (ten ture leaf_leaf), roots from ten-true leaf stage ten ture leaf_root), and flowers at squaring stage (squaring stage_flower). The X axis is the samples in tissues at different developmental stages. The color scale represents Log2 expression values

The gene expression levels in NtPYL subfamily I have a broader range from 2.83 (NtPYL3) to 8.33 (NtPYL18) in dry seeds. When compared with the expression levels in the dry seeds, expression levels of NtPYL3, NtPYL4, and NtPYL5 were constantly higher in germination seeds, cotyledons, leaves and roots from two, four, and six true leaves stages, and leaves in the ten true leaves stage; expression levels of NtPYL10, NtPYL11, and NtPYL12 were constantly higher in cotyledons, leaves, and roots in two true leaves stage; expression levels of NtPYL15 and NtPYL16 were higher from germination seeds to roots in the two true leaves stage; expression levels of NtPYL6 and NtPYL18 were constantly lower from four true leaves stage to flowering stage. Notably, the NtPYL7 and NtPYL17 had lower expression levels in leaves comparison with roots in the four, six, and ten true leaves stage.

In NtPYL subfamily II, the expression levels of most genes (8 of 9 genes) ranged from 5 to 8, except the lower expression level of NtPYL23 (3.3) in dry seeds. Compared with the expression levels in the dry seeds, expression level of NtPYL23 showed constantly higher in germination seeds till flowers, and expression level of NtPYL26 and NtPYL27 have constantly higher in cotyledons till leaves in ten true leaves stage. However, NtPYL21 and NtPYL22 showed lower expression levels in leaves of four, six, and ten true leaves stage and flowers than those in dry seeds. Notably, expression levels of NtPYL21 and NtPYL22 are lower in leaves compared with roots at four, six, and ten true leaf stages.

For the two members of NtPYL family III, NtPYL29 (5.13) showed higher expression level than that of NtPYL28 (3.84) in dry seeds. Interestingly, expression levels of both NtPYL28 and NtPYL29 are constantly lower in all the developmental stages after dry seeds stage compared with those in the dry seeds stage (Fig. 6, Additional file 5: Table S5).

Expression profiles of NtPYLs in response to drought stress

To understand the possible function of NtPYLs in plant response to drought stress, we analyzed the expression profiles of NtPYLs in the tobacco seedlings after drought treatment for indicated time. The expression levels and patterns of all NtPYLs were detected by the Microarray. Notably, three pairs of gene, NtPYL4 and NtPYL5, NtPYL11 and NtPYL12, and NtPYL26 and NtPYL27, have identical expression levels and patterns (Fig. 7, Additional file 6: Table S6). The NtPYL4 and 5 are derived from Nsy10370740, while NtPYL11 and 12 are derived from Ntom0128520 (Fig. 3, Table 4). The amino acid sequence similarity between NtPYL26 and 27 is 98.87% (Additional file 4: Table S4). Therefore, the microarray experiment could not distinguish the expression profiles for each of three gene pairs.

Fig. 7
figure7

The expression profile of 29 NtPYLs in response to drought treatment. The relative transcript abundances of 29 NtPYLs were examined via microarray and visualized as a heatmap. The 6–7 weeks old Nicotiana tabacum seedings grown in the soil were treated by air dehydration. The X axis is represented for the indicated time after drought treatment. The color scale represents Log2 expression values

In addition, the expression value of each NtPYL in CK was assigned as 1, and the expression ratio for the treatment at indicated time/CK was calculated for the expression pattern of each NtPYLs (Additional file 6: Table S6). The ratio between treatment and CK that is more than 1.2 or less than 0.8 was recognized as up- or down-regulation, respectively. For the NtPYLs in subfamily I, expression levels of 4 NtPYLs (NtPYL6, NtPYL15, NtPYL16 and NtPYL18) are down-regulated after dehydration treatment for 2 h, expression levels of 7 NtPYLs (NtPYL7, NtPYL8, NtPYL9, NtPYL10, NtPYL11, NtPYL12 and NtPYL17) are up-regulated after dehydration treatment for 0.5 h and down-regulated thereafter, and expression levels of 7 NtPYLs (NtPYL1 to NtPYL5, NtPYL13 and NtPYL14) were changed slightly. For the NtPYLs in subfamily II, expression levels of 2 NtPYLs (NtPYL24, 25) are up-regulated after dehydration treatment for 2 h. For the subfamily III, the expression level of NtPYL29 is slightly up-regulated after dehydration treatment for 8 h (Fig. 7, Additional file 6: Table S6). Notably, the expression levels of most NtPYLs (11 NtPYLs in 18 NtPYLs) in subfamily I are decreased after dehydration treatment for 2 h, which is consistent with the previous study in tomato [54].

To further confirm the microarray data, gene expression of NtPYLs were performed by quantitative real-time PCR. The specific primers of NtPYLs were used for qRT-PCR validation (Additional file 7: Table S7). Consisted with microarray data, gene expressions of NtPYL6, NtPYL16 and NtPYL18 were down-regulated in response to dehydration treatment, and expression levels of NtPYL9, NtPYL11 and NtPYL17 were up-regulated for 0.5 h after dehydration treatment. While gene expressions of NtPYL1, NtPYL3, NtPYL4 and NtPYL14 were not significantly changed under dehydration treatment (Fig. 8). In addition, transcript levels of NtPYL24 and NtPYL25 were up-regulated after dehydration treatment, which is consisted with microarray data (Fig. 8).

Fig. 8
figure8

Gene expressions of 20 NtPYLs in response to dehydration treatment. The relative transcript levels of 20 NtPYLs were examined via quantitative real-time PCR. The 6–7 weeks old Nicotiana tabacum seedings grown in the soil were treated by air dehydration

Discussion

Abscisic acid (ABA) is an important phytohormone for plant growth, development, and response to many environmental stresses, particularly in drought stress [1,2,3, 35]. Recently, genetic and biochemical studies revealed that AtPYLs are ABA receptors in Arabidopsis [14,15,16], and the ABA receptor PYLs are the key regulators to perceive ABA for initiating ABA signaling [1, 3, 4, 8, 9, 21, 22, 26,27,28, 60, 61]. PYLs have been characterized in many plants, including Arabidopsis [8, 14,15,16, 30,31,32,33,34,35,36], soybean [43], tomato [41, 42], wheat [44], maize [45, 46], poplar [47], rubber tree [48], cotton [50], and rice [37,38,39,40]. However, knowledge about systematical identification and characterization of Nicotiana PYL family is very limited at the genome level [59] . In this study, we identified 29, 11, and 16 PYLs based on the similarity to AtPYLs in the genomes of N. tabacum, N. tomentosiformis and N. sylvestris, respectively. The phylogenic relationships, gene structures, chromosomal location, conserved motifs/regions, and expression profiles of NtPYL family were analyzed.

According to the coding and amino acid sequences of 14 AtPYLs, 11 NtomPYLs, 16 NsylPYLs, and 29 NtPYLs were retrieved from the genomes of the tetraploid specie N. tabacum and its two progenitor diploid species N. tomentosiformis, and N. sylvestris, respectively. NtPYLs, NtomPYLs, and NsylPYLs and their putative proteins exhibited similar physical properties. In addition, amino acid sequence alignment analysis showed that all the 29 NtPYLs share highly similar structure characterized by several α-helices, β-sheets and conserved CL loops, similar with the AtPYL2, a functional member of the ABA receptor family in Arabidopsis [62]. Therefore, these NtomPYLs, NsylPYLs, and NtPYLs are the ABA receptors in N. tomentosiformis, N. sylvestris, and N. tabacum, respectively. In cotton, PYL families had been identified recently at genome level in two ancestral diploid species Gossypium raimodii and G. arboretum, and two tetraploid species G. barbadense and G. hirsutum derived from G. raimodii and G. arboretum, respectively. Analysis of the physical properties of the 20 GrPYLs, 21 GaPYLs, 39 GrPYLs, and 40 GaPYLs revealed that the Gossypium PYL families share similar ORF and amino acid lengths, MW and pI [50]. Compared with the published plant PYL families [14, 15, 38, 41, 48, 50], we found that the Nicotiana PYLs had similar physical properties regarding to the amino acid lengths, ORF and MW. Notably, the number of PYL members (29) in N. tabacum is larger than those in the tetraploid soybean (23) [43] and other diploid plant species, such as Arabidopsis (14) [14,15,16], and tomato (15) [41], but lesser than that in the tetraploid cotton (39 or 40) [50], indicating that tetraploid plant species that harbor large PYL families might have more complex ABA responses in plant development and respond to environmental stress.

In Arabidopsis and soybean, PYLs could be classified into 3 subfamilies, and each subfamily might have diverse function [14, 15, 43]. The phylogenetic analysis of amino acid sequences of NtPYLs, AtPYLs and GmPYLs showed that NtPYLs could also be grouped into 3 subfamilies. Further analysis also identified 3 subfamilies among three Nicotiana species (N. tomentosiformis, N. sylvestris, and N. tabacum). Amino acid sequence alignment of Nicotiana PYLs suggested that these PYLs could be further divided into 13 groups (Table 4). The phylogenetic relationships of Nicotiana PYLs revealed that each of the 29 NtPYLs had an orthologous ancestor in either N. tomentosiformis or N. sylvestris. Notably, 16 NsylPYLs and 11 NtomPYLs have been identified in N. sylvestris and N. tomentosiformis, respectively. These results suggested that both S-genome from N. sylvestris (SS) and T-genome from N. tomentosiformis (TT) had contributed the PYLs for the genome of N. tabacum (TTSS), and S-genome might contribute more PYLs than that of T-genome during evolution [57, 58].

It was showed that NtPYL4 and NtPYL5 are located in the different positions of chromosome 18 in the N. tabacum genome, while NtPYL11 and NtPYL12 are almost located in the identical position of N. tabacum chromosome 16. These observations suggest that there might be an event of gene duplication for Nsyl0370740, the putative ancestor of NtPYL4 or NtPYL5, on N. tabacum chromosome 18. Given the physical locations of NtPYL11 and NtPYL12 were almost identical on N. tabacum chromosome 16, a tandem gene duplication event of Ntom0128520 might happen during the evolution. However, errors during sequencing, assembling, and/or annotation might result in different annotations for the same gene.

Gene structures of Nicotiana PYLs could also provide useful information to understand the evolution relationship of PYL family from N. sylvestris and N. tomentosiformis to N. tabacum. For example, since NtPYL13, NtPYL22 is the shortest Nicotiana PYLs (154 aa) identified in this study (Fig. 3, Fig. 4 a), alignment analysis of the gene structures and conserved motifs of NtPYL13, NtPYL22, and their putative ancestors revealed that both NtPYL13 and NtPYL22 are lacking the N-terminal motifs when compared with their putative ancestor Ntom0177100 and Nsyl0129950, respectively (Fig. 4 b). These results indicated that NtPYL13 had lost the 5′-region of Ntom0177100, and NtPYL22 had lost the first exon and intron of Nsyl0129950 during the evolution.

The CL2 loops/regions, particularly the combination of two key amino acid residues in the CL2 region, in AtPYLs determine the ABA-dependence of PYL-PP2C interactions in Arabidopsis [55, 56]. Recent studies on the mutation and degradation of PYLs for PYL-PP2C interactions, together with the downstream substrates of SnRKs, suggest the potential to engineer the core PYLs-PP2Cs-SnRKs components of ABA signaling [30, 32, 63,64,65]. For the combination of two key amino acid residues in the CL2 region, NtPYLs share the same pattern with those in the AtPYLs, GmPYLs: VI and VV for PYL subfamily I, and VK for PYL subfamily II. AtPYL13, NtPYL28, and NtPYL29 of the PYL subfamily III share the pattern of LV. The conserved combination of two key amino acid residues in the CL2 region in N. tabacum, Arabidopsis and soybean, suggest that the CL2 loop/region play vital roles in PYL-PP2C interactions in plants. Whether N. tabacum have the similar ABA dependence of PYL-PP2C interactions as those in Arabidopsis deserves further study using yeast two-hybrid and bimolecular fluorescence complementation (BiFC) assays.

Expression profiles analysis in many plants, such as Arabidopsis [5, 35, 36], rice [37, 38, 40], tomato [41, 54], maize [46], cotton [50], and rubber tree [48], revealed that the majority of PYLs are expressed in various tissues including seeds, root, leaves, flowers and fruits. For example, most PYLs were expressed in all tissues of rice, and OsPYL7/8 were highly expressed in embryos, OsPYL3 and OsPYL5 were primarily expressed in leaves, while OsPYL1 was predominantly expressed in roots [38]. In solanaceous tomato, most of the PYL genes (including two members of subfamily I, five members of subfamily II, and two members of subfamily III) were highly expressed in roots, while one member (Sl1g095700) from PYL subfamily III showed the highest expression levels than those of other PYLs in leaves. Notably, one gene (Sl6g061180) from PYL subfamily I and most of the members in PYL subfamily II and III exhibited high expression levels during tomato fruit development. Interestingly, transcripts of many fruit-expressed SlPYLs gradually increased to the highest level and declined afterwards. Sun et al. also reported that transcripts of SlPYL1, SlPYL2 (two tomato orthologue of AtPYL7 and AtPYL9), SlPYL3 (a tomato orthologue of AtPYL8 and AtPYL10), and SlPYL6 (a tomato orthologue of AtPYL4) within the SlPYL family fluctuated during fruit development and ripening in tomato [54]. Additionally, expression profile analysis of PYLs in a tetraploid cotton (G. hirsutum) demonstrated that 22 GhPYLs were preferentially expressed in flowers, 10 GhPYLs were dominantly expressed in roots, and 3 GhPYLs were highly expressed in the fiber [50]. As a tetraploid solanaceous specie, N. tabacum might have similar tissue- and developmental stage-specific expression of PYLs as those reported in tomato and cotton.

In this study, expression profiles of NtPYLs in different tissues and developmental stages were analyzed by chip-Microarray. Transcripts of 29 NtPYLs could be detected in all the tissue/developmental stage checked (Fig. 6, Additional file 5: Table S5). Notably, NtPYL28 and NtPYL29 from subfamily III showed a distinct expression pattern during the development: transcripts of NtPYL28 and NtPYL29 reached the highest level in dry seeds, reduced dramatically in germination seeds, and remain almost constant low level afterwards. The highly accumulation of NtPYL28 and NtPYL29 transcripts particularly in dry seeds, suggest that NtPYL28 and NtPYL29 might play important roles in N. tabacum seed dormancy. Kim et al. reported that overexpressing OsPYL/RCAR5 (a rice orthologue of AtPYL8) resulted in hypersensitive to ABA during seed germination and early seedling growth [37]. Overexpressing OsPYL3 (a rice orthologue of AtPYL8) and OsPYL9 (a rice orthologue of AtPYL3) confer ABA hypersensitivity during seed germination in rice [38]. These studies in rice indicated some OsPYLs (OsPYL3, 5, and 9) are involved in seed germination and early seedling growth in response to ABA treatment. Recently, AtPYL8 and AtPYL9 have been shown to play important roles in regulating lateral root growth in the presence of ABA [5, 36]. Given NtPYL19 and NtPYL20 are cluster closely with AtPYL7, AtPYL9, AtPYL8 and AtPYL10 in phylogenetic analysis (Fig. 1), together with the previous studies on functions of the orthologue of AtPYL8 and AtPYL10 in tomato fruit development, and unique expression patterns of NtPYL28 and NtPYL29 in dry seeds during developmental process, the functions of NtPYL19, 20, 28, 29 in seed development, germination, response to ABA deserve further study.

Since the phytohormne ABA play crucial roles in plant responses to drought, we analyzed the expression levels of NtPYLs in the whole seedling after dehydration treatment using chip-Microarray (Fig. 7, Additional file 6: Table S6). In general, the transcripts of members in NtPYL subfamily II and III remain constant in the control seedlings (CK, without drought stress) and during the drought treatment, indicating that NtPYL19–27, NtPYL28, and NtPYL29 might not involve in the N. tabacum seedling responses to drought stress. Interestingly, for members in NtPYL subfamily I, the transcriptional responses to drought treatment could be divided into three types: (1) NtPYL1–5, 8, 9, 13, and 14 showed constant expression levels; (2) NtPYL7, 10–12, and 17 exhibit quick induction (0.5 and 1 h after drought treatment) with a gradually decrease in the later time points (2, 4, and 8 h after drought treatment); (3) NtPYL6, 15, 16, and 18 have a constant transcript accumulation in the early time point (0.5 and 1 h after drought treatment) with a clear decline in the later time points (2, 4, and 8 h after drought treatment). The distinct expression patterns of members in NtPYL subfamily I suggest that some NtPYLs (such as NtPYL7, 10–12, and 17) mainly function in the early response to drought stress, while other NtPYLs (such as NtPYL6, 15, 16, and 18) might play roles in the later response to dehydration.

Overexpression OsPYL/RCAR5 (a rice orthologue of AtPYL8), OsPYL3 (a rice orthologue of AtPYL8), and OsPYL9 (a rice orthologue of AtPYL3) conferred improved drought stress tolerance in rice [37, 38, 40]. Moreover, in leaves of tomato seedlings subjected to dehydration, transcript accumulation of 6 SlPYLs (SlPYL2–7) were deduced, and expression levels of 3 SlPYLs (SlPYL4, 6, and 7) were recovered after re-watering [54].

The unexpected results in Arabidopsis showed that expression levels of AtPYLs were down-regulated by stress [12]. Overexpression of AtPYL5 and AtPYL9 resulted in enhanced ABA responses and drought resistance in Arabidopsis through PYL-mediated inhibition of clade A PP2Cs [16, 34, 35]. Importantly, Gonzalez-Guzman et al. revealed that overexpression of tomato monomeric-type, but not dimeric ABA receptors in Arabidopsis confers enhanced resistance to drought stress [41]. In this study, NtPYL6–12 in NtPYL subfamily I were clustered closely with AtPYL4 and AtPYL5 (Fig. 1), which are monomeric-type ABA receptors [55]. Together with the NtPYLs that showed down-regulated expression pattern in NtPYL subfamily I, NtPYL6, 7, 10, 11, 12 will be the key candidate ABA receptors for functional identification in N. tabacum response to drought, salt, and osmotic stresses.

Notably, a comprehensive comparison comparing transcriptional profiles of the core ABA signaling components under osmotic/dehydration stress or ABA treatment between roots and leaves of maize (Zea mays) seedlings grew in hydroponic culture revealed that, after treating roots with ABA, the expression of ZmPYLs homologous to monomeric-type AtPYLs were reduced, whereas those of ZmPYLs homologous to dimeric-type AtPYLs were increased in maize primary roots. Surprisingly, the opposite pattern was observed in the leaves of the same experiments [46]. This interesting organ-specific expression patterns for ABA receptor genes between roots and leaves of maize in response to ABA suggest that there might be a contrast transcriptional patterns for monomeric- and dimeric-type NtPYLs in roots and leaves under abiotic stresses and ABA treatment, which deserves further study in the future study.

Interestingly, the expression patterns of many members in each NtPYL subfamily are similar in different tissues/developmental stages and in seedlings under drought stress, suggesting these NtPYLs could be regulated by the developmental signals and dehydration treatment as well. Characterization of these NtPYLs is required for understanding their functions during developmental process and in response to stresses.

In the current study, our study identified 29 putative NtPYLs in the allotetraploid cultivated tobacco (N. tabacum), 11 NtomPYLs and 16 NsylPYLs in the two diploid ancestral tobacco species, N. tomentosiformis and N. sylvestris, respectively. We further investigated the physical properties, phylogenetic relationship, protein motifs, and gene structures of Nicotiana PYLs. The Nicotiana PYLs could also be divided to 3 subfamilies, in consistent with the results from studies in other plants. Phylogenetic, gene structure, and protein motif analysis revealed NtPYLs were originated from NtomPYLs and NsylPYLs, NtPYL22 had lost the first exon and intron compared with its origin in N. sylvestris0129950, and NtPYL13 was derived from N-terminal truncation of N. tomentosiformis0177100 during evolution. Moreover, analysis of NtPYLs expression profiles via chip-Microarray assay indicated that each NtPYL subfamily might have diverse functions in different tissues/developmental stages and in response to abiotic stresses. Furthermore, four NtPYLs (NtPYL 19, 20, 28, and 29) in NtPYL subfamily II and III were suggested to play important roles in seed development, germination, and response to ABA. Finally, five NtPYLs (NtPYL6, 7, 10, 11, 12) in NtPYL subfamily I were highlighted as potential candidates for further characterizing their functions in N. tabacum resistance to abiotic stresses. Taken together, the results from genome-wide identification of Niacotiana PYLs will provide some insights on understanding the roles of PYLs in ABA signaling and in response to abiotic stresses in tetraploid plants, which will facilitate the improvement of crop resistance to drought stress.

Conclusions

ABA receptors (PYLs) play central roles in ABA signaling and plant response to many environmental stresses. Although lots of PYL genes and family have been identified in many plant species, the information of Nicotiana PYL family is still missing. Here we conducted a genome-wide identification and expression analysis of the PYL family in N. tabacum. A total of 29, 11, 16 Nicotiana PYL genes were identified in the genome of N. tabacum and its two ancestors N. tomentosiformis and N. sylvestris, respectively. Phylogenetic and gene structure analysis revealed that Nicotiana PYL could be divided into 3 subfamilies and 13 groups. Furthermore, each NtPYL might have a putative orthologous gene in either N. sylvestris or N. tomentosiformis, in consistent with the evolutionary origin of N. tabacum. The microarray-based analysis of NtPYLs expression profiles in tissues at different developmental stages, and in response to drought stress revealed that members in different NtPYL subfamily might paly specific roles in N. tabacum growth, developmental, and drought stress responses. Interestingly, the expression profiles of members in the same NtPYL subfamily showed somehow similar patterns in tissues at different developmental stages and in leaves of seedlings under drought stress, suggesting particular NtPYLs might have multiple functions in both plant development and drought stress response. Importantly, four NtPYLs (NtPYL 19, NtPYL20, NtPYL28, and NtPYL29) are highlighted for the potential functions in seed development, germination and response to ABA. Moreover, five NtPYLs (NtPYL6, NtPYL7, NtPYL10, NtPYL11, NtPYL12) might play important roles in response to abiotic stresses, particularly in drought. Taken together, these results will facilitate further functional characterization of NtPYLs in plant development and in response to abiotic and biotic stresses in tetraploid plants.

Methods

Plant materials and growth conditions

The Nicotiana tabacum L. cv. Honghuadajingyuan seeds were obtained from Yunnan Academy of Tobacco Agricultural Sciences (Yunnan, China) [66]. Surface-sterilized seeds were directly sowed into the soil in pots. The tobacco young seedlings were grown in the plant growth chamber with a 16-h-light/8-h-dark photoperiod under continuous white light (75 mol m− 2 s− 1) at 28 °C-day/ 23 °C-night. All plants were kept well-watered after sowing.

For expression profiling of ABA receptor genes in response to drought stress, the plants were grown for 7–8 weeks with 6–7 leaves. The plants were moved out from the pots carefully without disturbing the root, and the surface soil was washed out softly. Then the plants were put on the bench for air drying which termed as drought stress treatment. The whole seedlings were collected at the indicated time after treatment, and immediately frozen in liquid nitrogen for RNA extraction for microarray assay. Five biological replicates were used for sample harvesting at each indicated time of the treatment.

For expression profiling of ABA receptor genes in different developmental stages and tissues in tobacco, tobacco plants were kept in the growing condition mentioned above. Samples were harvested from 12 tissues at different developmental stages, including dry seed, germination seeds, cotyledon, leaves and roots from two, four, six, and ten true-leaf stage, respectively, and flowers at squaring stage. Samples were immediately frozen in liquid nitrogen for RNA extraction for microarray assay. Three to five biological replicates were used for sample harvesting.

Analyses of phylogenetic, conserved motif, isoelectric point prediction, gene structure, and chromosome localization

The Nicotiana PYL gene sequences were retrieved from NCBI (https://www.ncbi.nlm.nih.gov/) and China tobacco genome database V2.0 (data not published). The genomic DNA, open reading frame, and deduced protein sequences of Nicotiana PYL family are provided in the (Additional files 1, 2, 3: Tables S1, S2, S3), and had been submitted to NCBI database (currently waiting for assigning the accession numbers).

The sequences of GmPYLs and AtPYLs were retrieved from the NCBI GeneBank database. The sequences of NtPYLs, AtPYLs, and GmPYLs were aligned using ClustalW [67], and an unrooted phylogenetic tree was generated using MEGA 7.0 software (http://www.megasoftware.net) by the neighbor-joining method with 1000 replicates of bootstrap analysis. For analyzing secondary structure of NtPYLs, the above alignment results were further treated by ESPript 3.0 software with default parameter settings (http://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi).

Protein motifs were predicted by motif elicitation program MEME (http://meme-suite.org/tools/meme). The conserved motifs were further queried in the InterPro database (http://www.ebi.ac.uk/interpro). The isoelectric point and molecular weight of deduced NtPYL proteins were predicted by ProtParam tool (http://web.expasy.org/protparam/).

For gene structure analysis, the open reading frame sequence of each NtPYL was aligned with its corresponding genomic DNA sequence to identify the gene structure by GSDS 2.0 (http://gsds.cbi.pku.edu.cn/). The chromosome localization of NtPYLs were mapped with localization MapGene2Chromosome V2 (http://mg2c.iask.in/mg2c_v2.0/).

Transcriptomic microarray analysis

Total RNA was extracted with SuperPure Plantpoly RNA Kit (GeneAnswer, China). All RNA samples were treated with RNase-free DNase I (GeneAnswer) and analyzed for integrity on a Bioanalyzer 2100 (Agilent technologies, USA). About 33.3 ng total RNA were used for amplification with WT Amplification Kit (Affymetrix, Thermo Fisher Scientific, USA). 5.5 μg of the amplified product were fragmented by uracil-DNA glycosylase and apurinic/apyrimidinic endonuclease 1 (Affymetrix, Thermo Fisher Scientific, USA).

The fragmented cDNA was labeled by terminal deoxynucleotidyltransferase using the DNA Labeling Reagent (Affymetrix, Thermo Fisher Scientific, USA) that was covalently linked to biotin. The resulting labeled cDNAs (5.2 μg) were dissolved in 160 μl of hybridization mix solution, then denatured at 99 °C for 5 min. The mixed hybridization buffer was loaded into a microarray, and then the both septa were covered by round labels to prevent leaks and evaporation.

An Affymetrix custom Tobacco Genome Array with feature Size 5 μm was used. Eighty thousand six hundred fifty-two tobacco genes were covered within this array. Tobacco L25, EF1-alpha, Ntubc2, PP2A genes were used as housekeeping genes. RMA method provided by the R package, affy package, was used to conduct background correction, normalization, probe-specific background correction, probe summarization and to convert probe level data to expression values.

The hybridizations were performed in a hybridization oven (Affymetrix, Thermo Fisher Scientific, USA) at 45 °C for 16 h. After hybridization, microarrays were washed by Fluidics Station 450 with wash buffer A and B (Affymetrix, Thermo Fisher Scientific, USA). Three biological replicates were used in the Microarrays assay. The expression levels of members in NtPYL family in several tissues at different developmental stages and in response to drought stress were documented in Additional file 5: Table S5 and Additional file 6: Table S6, respectively.

qRT-PCR validation of chip-microarray

For RT-qPCR validation of expression pattern of ABA receptor genes in response to drought stress, the plants were grown for 7–8 weeks with 6–7 leaves. The plants were moved out from the pots carefully without disturbing the root, and the surface soil was washed out softly. Then the plants were put on the bench for air drying which termed as drought stress treatment. The whole shoots were collected at the indicated time after treatment, and immediately frozen in liquid nitrogen for RNA extraction. Three to five biological replicates were used for sample harvesting at each indicated time of the treatment.

RNA was extracted from three to five biological replicates of the whole shoots collected at the indicated time after drought treatment using the Qiagen RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) following the manufacturer instructions.

2 μg of total RNA in a 20 μl reaction was converted to cDNA with a SuperScript III Reverse Transcriptase (Invitrogen, USA) by manufacturer instructions on a Eppendorf Mastercycler thermocycler (Eppendorf AG, Germany) with the following conditions: 25 °C for 5 min, 50 °C for 60 min, 70 °C for 15 min, followed by a hold at 4 °C until use in RT-qPCR reaction. 60 μl deionized water was added into 20 μl cDNA, and 1 μl of diluted cDNA mixture was used as the input for qPCR reaction. qPCR reactions were made with a SuperReal PreMix Plus SYBR Green Kit (TIANGEN Biotech, China) following manufacturer instructions in a 20 μl volume. The specific primers of NtPYLs were used for qRT-PCR validation (Additional file 7: Table S7).

qPCR was done on an Applied Biosystems™ QuantStudio™ 6 Flex Real-Time PCR System (ThemoFisher Scientific, USA) with the following cycling conditions: 95 °C for 15 min, followed by 40 cycles of 95 °C for 10 s, 60 °C for 20 s, and 72 °C for 32 s. Melt curve conditions were 95 °C for 15 s, 60 °C for 1 min, 95 °C for 15 s. All samples had only one melt temperature peak. Fold change between experimental samples (with drought treatment) and control samples (without drought treatment) was calculated by the 2-ΔΔCT method using 26S as a reference gene. CT values represent the average of three technical replicates.

Statistical analysis

The presented values are the means ± SE of three individual experiments with three replicated measurements. An analysis of variance (ANOVA) was used to compare significant differences based on significance levels of P < 0.05 and P < 0.01.

Availability of data and materials

The original data that support the findings of this study are available from National Tobacco Gene Research Centre at Zhengzhou Tobacco Research Institute but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of National Tobacco Gene Research Centre at Zhengzhou Tobacco Research Institute.

Abbreviations

ABA:

Abscisic acid

MW:

Molecular weight

NJ:

Neighbor-joining

ORF:

Open reading frame

pI:

isoelectric point

PYL:

PYR1-Like

PYR:

Pyrabactin Resistance

RCAR:

Regulatory Component of ABA Receptor

References

  1. 1.

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

  2. 2.

    Nambara E, Marion-Poll A. Abscisic acid biosynthesis and catabolism. Annu Rev Plant Biol. 2005;56:165–85.

  3. 3.

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

  4. 4.

    Lumba S, Cutler S, McCourt P. Plant nuclear hormone receptors: a role for small molecules in protein-protein interactions. Annu Rev Cell Dev Biol. 2010;26:445–69.

  5. 5.

    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. Sci Rep. 2016;6:27177.

  6. 6.

    Weng JK, Ye M, Li B, Noel JP. Co-evolution of hormone metabolism and signaling networks expands plant adaptive plasticity. Cell. 2016;166(4):881–93.

  7. 7.

    Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53:247–73.

  8. 8.

    Yang Z, Liu J, Tischer SV, Christmann A, Windisch W, Schnyder H, et al. Leveraging abscisic acid receptors for efficient water use in Arabidopsis. Pro Natl Acad Sci. 2016;113(24):201601954.

  9. 9.

    Zhao Y, Gao J, Im KJ, Chen K, Bressan RA, Zhu JK. Control of plant water use by ABA induction of senescence and dormancy: an overlooked lesson from evolution. Plant Cell Physiol. 2017;58(8):1319–27.

  10. 10.

    Chinnusamy V, Schumaker K, Zhu JK. Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. J Exp Bot. 2004;55(395):225–36.

  11. 11.

    Aleman F, Yazaki J, Lee M, Takahashi Y, Kim AY, Li Z, et al. 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:28941.

  12. 12.

    Chan Z. Expression profiling of ABA pathway transcripts indicates crosstalk between abiotic and biotic stress responses in Arabidopsis. Genomics. 2012;100(2):110–5.

  13. 13.

    Wang P, Zhao Y, Li Z, Hsu CC, Liu X, Fu L, et al. Reciprocal regulation of the TOR kinase and ABA receptor balances plant growth and stress response. Mol Cell. 2018;69(1):100–12 e6.

  14. 14.

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

  15. 15.

    Park SY, Fung P, Nishimura N, Jensen DR, Fujii H, Zhao Y, et al. Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science. 2009;324(5930):1068–71.

  16. 16.

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

  17. 17.

    Leung J, Bouvier-Durand M, Morris PC, Guerrier D, Chefdor F, Giraudat J. Arabidopsis ABA response gene ABI1: features of a calcium-modulated protein phosphatase. Science. 1994;264(5164):1448–52.

  18. 18.

    Santiago J, Dupeux F, Round A, Antoni R, Park SY, Jamin M, et al. The abscisic acid receptor PYR1 in complex with abscisic acid. Nature. 2009;462(7273):665–8.

  19. 19.

    Melcher K, Ng LM, Zhou XE, Soon FF, Xu Y, Suino-Powell KM, et al. A gate-latch-lock mechanism for hormone signalling by abscisic acid receptors. Nature. 2009;462(7273):602–8.

  20. 20.

    Melcher K, Xu Y, Ng LM, Zhou XE, Soon FF, Chinnusamy V, et al. Identification and mechanism of ABA receptor antagonism. Nat Struct Mol Biol. 2010;17(9):1102–8.

  21. 21.

    Melcher K, Zhou XE, Xu HE. Thirsty plants and beyond: structural mechanisms of abscisic acid perception and signaling. Curr Opin Struct Biol. 2010;20(6):722–9.

  22. 22.

    Miyakawa T, Fujita Y, Yamaguchi-Shinozaki K, Tanokura M. Structure and function of abscisic acid receptors. Trends Plant Sci. 2013;18(5):259–66.

  23. 23.

    Miyazono K, Miyakawa T, Sawano Y, Kubota K, Kang HJ, Asano A, et al. Structural basis of abscisic acid signalling. Nature. 2009;462(7273):609–14.

  24. 24.

    Nishimura N, Hitomi K, Arvai AS, Rambo RP, Hitomi C, Cutler SR, et al. Structural mechanism of abscisic acid binding and signaling by dimeric PYR1. Science. 2009;326(5958):1373–9.

  25. 25.

    Peterson FC, Burgie ES, Park SY, Jensen DR, Weiner JJ, Bingman CA, et al. Structural basis for selective activation of ABA receptors. Nat Struct Mol Biol. 2010;17(9):1109–13.

  26. 26.

    Hauser F, Waadt R, Schroeder JI. Evolution of abscisic acid synthesis and signaling mechanisms. Curr Biol. 2011;21(9):R346–55.

  27. 27.

    Hubbard KE, Nishimura N, Hitomi K, Getzoff ED, Schroeder JI. Early abscisic acid signal transduction mechanisms: newly discovered components and newly emerging questions. Genes Dev. 2010;24(16):1695–708.

  28. 28.

    Joshi-Saha A, Valon C, Leung J. Abscisic acid signal off the STARting block. Mol Plant. 2011;4(4):562–80.

  29. 29.

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

  30. 30.

    Sang-Youl P, Peterson FC, Assaf M, Jin Y, Volkman BF, Cutler SR. Agrochemical control of plant water use using engineered abscisic acid receptors. Nature. 2015;520(7548):545–8.

  31. 31.

    Izabela S, Klaus R, Michal K, Simone D, Yue M, Arthur K, et al. Closely related receptor complexes differ in their ABA selectivity and sensitivity. Plant J. 2010;61(1):25–35.

  32. 32.

    Tischer SV, Wunschel C, Papacek M, Kleigrewe K, Hofmann T, Christmann A, et al. Combinatorial interaction network of abscisic acid receptors and coreceptors from Arabidopsis thaliana. Pro Natl Acad Sci. 2017;114(38):10280–5.

  33. 33.

    Yizhou W, Zhong-Hua C, Ben Z, Adrian H, Blatt MR. PYR/PYL/RCAR abscisic acid receptors regulate K+ and cl- channels through reactive oxygen species-mediated activation of Ca2+ channels at the plasma membrane of intact Arabidopsis guard cells. Plant Physiol. 2013;163(2):566–77.

  34. 34.

    Zhao J, Zhao L, Zhang M, Zafar SA, Fang J, Li M, et al. Arabidopsis E3 ubiquitin ligases PUB22 and PUB23 negatively regulate drought tolerance by targeting ABA receptor PYL9 for degradation. Int J Mol Sci. 2017;18(9):E1841.

  35. 35.

    Zhao Y, Chan Z, Gao J, Xing L, Cao M, Yu C, et al. ABA receptor PYL9 promotes drought resistance and leaf senescence. Pro Natl Acad Sci. 2016;113(7):1949–54.

  36. 36.

    Zhao Y, Xing L, Wang X, Hou YJ, Gao J, Wang P, et al. The ABA receptor PYL8 promotes lateral root growth by enhancing MYB77-dependent transcription of auxin-responsive genes. Sci Signal. 2014;7(328):ra53.

  37. 37.

    Kim H, Hwang H, Hong J-W, Lee Y-N, Ahn IP, Yoon IS, et al. 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(2):1013–24.

  38. 38.

    Tian X, Wang Z, Li X, Lv T, Liu H, Wang L, et al. Characterization and functional analysis of Pyrabactin resistance-like abscisic acid receptor family in Rice. Rice (N Y). 2015;8(1):28.

  39. 39.

    He Y, Hao Q, Li W, Yan C, Yan N, Yin P. Identification and characterization of ABA receptors in Oryza sativa. PLoS One. 2014;9(4):e95246.

  40. 40.

    Kim H, Lee K, Hwang H, Bhatnagar N, Kim DY, Yoon IS, et al. Overexpression of PYL5 in rice enhances drought tolerance, inhibits growth, and modulates gene expression. J Exp Bot. 2014;65(2):453–64.

  41. 41.

    Mou W, Li D, Luo Z, Mao L, Ying T. Transcriptomic analysis reveals possible influences of ABA on secondary metabolism of pigments, flavonoids and antioxidants in tomato fruit during ripening. PLoS One. 2015;10(6):e0129598.

  42. 42.

    Gonzãlez-Guzmãn M, Rodrìguez L, Lorenzo-Orts L, Pons C, Sarrión Perdigones A, Fernãndez MA, et al. 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(15):4451–64.

  43. 43.

    Bai G, Yang D-H, Zhao Y, Ha S, Yang F, Ma J, et al. Interactions between soybean ABA receptors and type 2C protein phosphatases. Plant Mol Biol. 2013;83(6):651–64.

  44. 44.

    Gordon CS, Rajagopalan N, Risseeuw EP, Surpin M, Ball FJ, Barber CJ, et al. Characterization of Triticum aestivum abscisic acid receptors and a possible role for these in mediating Fusairum head blight susceptibility in wheat. PLoS One. 2016;11(10):e0164996.

  45. 45.

    Wang YG, Fu FL, Yu HQ, Hu T, Zhang YY, Tao Y, et al. Interaction network of core ABA signaling components in maize. Plant Mol Biol. 2018;96(3):245–63.

  46. 46.

    Fan W, Zhao M, Li S, Xue B, Jia L, Meng H, et al. 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(1):99.

  47. 47.

    Yu J, Yang L, Liu X, Tang R, Wang Y, Ge H, et al. Overexpression of poplar Pyrabactin resistance-like abscisic acid receptors promotes abscisic acid sensitivity and drought resistance in transgenic Arabidopsis. PLoS One. 2016;11(12):e0168040.

  48. 48.

    Guo D, Zhou Y, Li HL, Zhu JH, Wang Y, Chen XT, et al. 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.

  49. 49.

    Chai YM, Jia HF, Li CL, Dong QH, Shen YY. FaPYR1 is involved in strawberry fruit ripening. J Exp Bot. 2011;62(14):5079–89.

  50. 50.

    Zhang G, Lu T, Miao W, Sun L, Tian M, Wang J, et al. Genome-wide identification of ABA receptor PYL family and expression analysis of PYLs in response to ABA and osmotic stress in Gossypium. PeerJ. 2017;5:e4126.

  51. 51.

    Pri-Tal O, Shaar-Moshe L, Wiseglass G, Peleg Z, Mosquna A. Non-redundant functions of the dimeric ABA receptor BdPYL1 in the grass Brachypodium. Plant J. 2017;92(5):774–86.

  52. 52.

    Zhang F, Wei Q, Shi J, Jin X, He Y, Zhang Y, et al. Brachypodium distachyon BdPP2CA6 interacts with BdPYLs and BdSnRK2 and positively regulates salt tolerance in transgenic Arabidopsis. Front Plant Sci. 2017;8:264.

  53. 53.

    Iyer LM, Koonin EV, Aravind L. Adaptations of the helix­grip fold for ligand binding and catalysis in the START domain superfamily. Proteins. 2001;43(2):134–44.

  54. 54.

    Sun L, Wang YP, Chen P, Ren J, Ji K, Li Q, et al. Transcriptional regulation of SlPYL, SlPP2C, and SlSnRK2 gene families encoding ABA signal core components during tomato fruit development and drought stress. J Exp Bot. 2011;62(15):5659–69.

  55. 55.

    Hao Q, Yin P, Li W, Wang L, Yan C, Lin Z, et al. The molecular basis of ABA-independent inhibition of PP2Cs by a subclass of PYL proteins. Mol Cell. 2011;42(5):662–72.

  56. 56.

    Yin P, Fan H, Hao Q, Yuan X, Wu D, Pang Y, et al. Structural insights into the mechanism of abscisic acid signaling by PYL proteins. Nat Struct Mol Biol. 2009;16(12):1230–7.

  57. 57.

    Lim KY, Matyasek R, Kovarik A, Leitch AR. Genome evolution in allotetraploid Nicotiana. Biol J Linn Soc. 2004;82(4):599–606.

  58. 58.

    Renny-Byfield S, Chester M, Kovarik A, Le Comber SC, Grandbastien MA, Deloger M, et al. Next generation sequencing reveals genome downsizing in allotetraploid Nicotiana tabacum, predominantly through the elimination of paternally derived repetitive DNAs. Mol Biol Evol. 2011;28(10):2843–54.

  59. 59.

    Lackman P, Gonzalez-Guzman M, Tilleman S, Carqueijeiro I, Perez AC, Moses T, et al. Jasmonate signaling involves the abscisic acid receptor PYL4 to regulate metabolic reprogramming in Arabidopsis and tobacco. Pro Natl Acad Sci. 2011;108(14):5891–6.

  60. 60.

    Li W, Wang L, Sheng X, Yan C, Zhou R, Hang J, et al. Molecular basis for the selective and ABA-independent inhibition of PP2CA by PYL13. Cell Res. 2013;23(12):1369.

  61. 61.

    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.

  62. 62.

    Xiaoqiu Y, Ping Y, Qi H, Chuangye Y, Jiawei W, Nieng Y. Single amino acid alteration between valine and isoleucine determines the distinct pyrabactin selectivity by PYL1 and PYL2. J Biol Chem. 2010;285(37):28953.

  63. 63.

    Assaf M, Peterson FC, Sang-Youl P, Jorge LJ, Volkman BF, Cutler SR. Potent and selective activation of abscisic acid receptors in vivo by mutational stabilization of their agonist-bound conformation. Pro Natl Acad Sci. 2011;108(51):20838–43.

  64. 64.

    Umezawa T, Sugiyama N, Takahashi F, Anderson JC, Ishihama Y, Peck SC, Shinozaki K, et al. Genetics and phosphoproteomics reveal a protein phosphorylation network in the abscisic acid signaling pathway in Arabidopsis thaliana. Sci Signal. 2013;6(270):rs8.

  65. 65.

    Chen HH, Qu L, Xu ZH, Zhu JK, Xue HW. EL1-like casein kinases suppress ABA signaling and responses by phosphorylating and destabilizing the ABA receptors PYR/PYLs in Arabidopsis. Mol Plant. 2018;11(5):706–19.

  66. 66.

    Yang JK, Tong ZJ, Fang DH, Chen XJ, Zhang KQ, Xiao BG. Transcriptomic profile of tobacco in response to Phytophthora nicotianae infection. Sci Rep. 2017;7(1):401.

  67. 67.

    Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.

Download references

Acknowledgments

We thank Dr. Peijian Cao (National Tobacco Gene Research Centre at Zhengzhou Tobacco Research Institute) for retrieving Nicotiana PYL sequences from the China Tobacco Genome Database (CTGD) and Chip-Microarray experiments.

Funding

This work was funded by the National Natural Science Foundation of China (grant number 31760072 to G. Bai, and grant number 31860413 to H. Xie) and the Yunnan Academy of Tobacco Agricultural Sciences (grant numbers YNTC-2016YN22 to H. Xie, YNTC-2016YN24, YNTC-2015YN02, YNTC-2018530000241002, and YNTC-2019530000241003 to D.-H. Yang).

Author information

Conceptualization, GB, HX, YL and D-HY; Methodology, GB, HX, FL; Software, GB, HX, FL; Validation, XC, YZ and BX; Formal Analysis, GB, HX, HY, and D-HY; Writing-Original Draft Preparation, GB, HX, HY, and D-HY; Writing-Review & Editing, GB, HX, YL and D-HY; Supervision, JY, YL and D-HY. All authors read and approved the final manuscript.

Correspondence to Yongping Li or Da-Hai Yang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the result.

Additional information

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Genomic DNA sequences of NtPYLs, NtomPYLs, NsylPYLs. (XLSX 34 kb)

Additional file 2:

Table S2. Coding sequences of NtPYLs, NtomPYLs, NsylPYLs. (XLSX 19 kb)

Additional file 3:

Table S3. Deduced amino acid sequences of NtPYLs, NtomPYLs, NsylPYLs. (XLSX 14 kb)

Additional file 4:

Table S4. Matrix of amino acid sequence similarity of NtPYLs. (XLSX 16 kb)

Additional file 5:

Table S5. Expression levels of NtPYLs in tissues at different developmental stages. (XLSX 37 kb)

Additional file 6:

Table S6. Expression levels of NtPYLs in response to drought treatment. (XLSX 22 kb)

Additional file 7:

Table S7. List of primers used for qRT-PCR validation. (XLSX 10 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • ABA receptor
  • Drought stress
  • Gene family
  • Gene expression
  • Tetraploid
  • Nicotiana tabacum
  • Nicotiana tomentosiformis
  • Nicotiana sylvestris