Genome-wide identification and expression analysis of the cucumber PP2C gene family
BMC Genomics volume 23, Article number: 563 (2022)
Type 2C protein phosphatase (PP2C) is a negative regulator of ABA signaling pathway, which plays important roles in stress signal transduction in plants. However, little research on the PP2C genes family of cucumber (Cucumis sativus L.), as an important economic vegetable, has been conducted.
This study conducted a genome-wide investigation of the CsPP2C gene family. Through bioinformatics analysis, 56 CsPP2C genes were identified in cucumber. Based on phylogenetic analysis, the PP2C genes of cucumber and Arabidopsis were divided into 13 groups. Gene structure and conserved motif analysis showed that CsPP2C genes in the same group had similar gene structure and conserved domains. Collinearity analysis showed that segmental duplication events played a key role in the expansion of the cucumber PP2C genes family. In addition, the expression of CsPP2Cs under different abiotic treatments was analyzed by qRT-PCR. The results reveal that CsPP2C family genes showed different expression patterns under ABA, drought, salt, and cold treatment, and that CsPP2C3, 11–17, 23, 45, 54 and 55 responded significantly to the four stresses. By predicting the cis-elements in the promoter, we found that all CsPP2C members contained ABA response elements and drought response elements. Additionally, the expression patterns of CsPP2C genes were specific in different tissues.
The results of this study provide a reference for the genome-wide identification of the PP2C gene family in other species and provide a basis for future studies on the function of PP2C genes in cucumber.
Plants have developed multiple response mechanisms in response to salt, drought and cold stress, among which the regulation of related gene expression and protein modification are two important pathways. Protein reversible phosphorylation, a protein modification process that regulates multiple physiological responses in plants, is catalyzed by protein kinases (PKs) and phosphatases (PPs) [1, 2]. PKs primarily phosphorylate serine (Ser), threonine (Thr) and tyrosine (Tyr), while PPs can reverse this function by eliminating phosphate groups . Based on substrate specificity, protein phosphatases (PPs) can be divided into three groups: Ser/Thr phosphatases (STPs), protein Tyr phosphatases (PTPs), and dual-specificity phosphatases (DSPTPs). Moreover, PTPs are further classified into phosphoprotein metal phosphatases (PPM) and phosphoprotein phosphatases (PPP) based on crystal structure, amino acid sequence and specific response to inhibitors . The PPP family contains a variety of protein phosphatases, including PP1, PP4, PP5, PP6, PP7, PP2A and PP2B [5, 6]. At present, very little is known about the PPP family. However, we know that PP2A changes cell morphology in the PPP family during elongation, affecting root growth and root hair growth . PP1 is a highly conserved and ubiquitous phosphatase in eukaryotes . Nine members of the PP1 gene family have been identified in Arabidopsis [9, 10]. In addition, a PP1 gene up-regulated by biological stress was found in Phaseolus vulgaris .
Type 2C protein phosphatases (PP2C) are Mg2+/Mn2+ dependent serine/threonine phosphatases that are closely related to the PPP family, but have no sequence homology, thus forming a unique group belonging to the PPM family. The PP2C gene family is evolutionarily conserved, significantly regulates stress signaling pathways, and is primarily present in bacteria, fungi, plants and animals . Plant PP2C family proteins have unique structures that contain a conserved catalytic region at the N- or C-terminus and an unconserved extension region of varying length at the other end . The diversity of their structures indicates that PP2C family proteins have different functions in the signal transduction mechanism [14, 15]. A total of 80 PP2C genes were identified in Arabidopsis, which were divided into 13 subsets (A-J) according to their evolutionary relationships. Six of the nine members of subgroup A (ABI1, ABI2, AHG1, AHG3/ ATPP2CA, HAB1 and HAB2) negatively regulate ABA signaling, and the remaining three members (HAI1, HAI2/AIP1 and HAI3) respond differently to stress than the other six members [16,17,18]. PP2C of subgroup A inactivated PP2C through dephosphorylation, inhibited ABA receptor (PYR/PYL/RCRA) activity, and negatively regulated ABA signaling . The identification of two PP2C mutants, abi1 and abi2, provides favorable support for their negative regulatory function of ABA signaling . Group B participated in the mitogen-activated protein kinase (MAPK) signaling pathway and catalyzed the dephosphorylation of MAPKs, leading to their inactivation [21, 22]. Group C is involved in the regulation of flower development . Group D members can respond to salt and alkali stress . Group E PP2C is involved in the regulation of stomatal signaling, for example, as a response to dehydration, PP2C interacts with MPK3 and MPK6 to induce stomatal closure and inhibit water loss due to transpiration [25, 26]. Group F PP2C is involved in the induction of bacterial stress response . Currently, the PP2C gene family has been identified and studied in a variety of plants other than Arabidopsis, such as rice , maize , and Brachypodium distachyon . Rice OsPP108 belonging to subgroup A could be highly expressed under ABA, drought, and salt stress; when overexpressed in Arabidopsis, transgenic Arabidopsis showed stronger stress resistance . Overexpression of the maize ZmPP2C gene in Arabidopsis reduced the sensitivity to ABA and the tolerance to drought and salt treatment . FaABI1 in strawberry had negative regulatory effects on fruit ripening . The expression level of the PP2C gene regulated by endogenous ABA in Arabidopsis was significantly up-regulated under exogenous ABA and stress . In conclusion, the PP2C gene family plays a critical role in plant development and environmental stress.
Although the PP2C gene family has been extensively studied in other species, there are few reports on the PP2C gene family in cucumber. Therefore, identifying the PP2C gene family in cucumber and analyze its expression under stress is of great significance. In this study, 56 PP2C genes were identified at the genome-wide level in cucumber by bioinformatics methods. The physicochemical properties, protein secondary structure, chromosome location, phylogenetic analysis, gene structure, conserved motifs and cis-acting elements, and relative expression of 56 PP2C genes were systematically analyzed, providing reference for the function of abundant PP2C gene families. This study will lay the foundation for understanding the potential molecular mechanism of PP2C in stress signal transduction.
Identification and basic information of the CsPP2C gene family
In this study, we determined that 56 putative CsPP2C were present in the cucumber genome through BLASTp using 80 AtPP2C protein sequences as references. From the analysis of their physical and chemical properties (Table 1), the 56 CsPP2C genes identified encode proteins varying from 233 to 813 amino acids in length, isoelectric point (PI) value ranging from 4.5 to 9.61, and the molecular weight ranging from 30 kDa to 90 kDa. The total average hydrophobic index of the 56 CsPP2C gene family members were all less than zero and therefore encoded hydrophilic proteins. The subcellular localization prediction indicated that most of the CsPP2C proteins might be located in the nucleus, chloroplast, or cytoplasm, while only CsPP2C48 might be located in the endoplasmic reticulum and CsPP2C50 in the cytoskeleton.
Chromosome distribution and collinearity analysis of the PP2C gene family in cucumber
To obtain information on the position of CsPP2C genes on the chromosome, we used TBtools to map the chromosomal location (Fig. 1). A total of 56 PP2C genes were anchored to corresponding chromosomes and designated as CsPP2C1–CsPP2C56 according to their order on the chromosomes, among which chromosomes 3 and 6 were more distributed, and chromosome 5 was the least distributed, with only 3 PP2C genes. Closely related genes located within a distance of less than 200 kb on the same chromosome were defined as tandem duplications, otherwise they were defined as segmental duplications . To further understand the expansion mechanism of the CsPP2Cs, we examined segmental and tandem duplications within the cucumber genome. Our results showed that the PP2C gene family had no tandem duplication gene pairs, but that there were seven fragment repeat gene pairs (Fig. 2 a). In the seven pairs of collinear relationships, CsPP2C49 was paired with CsPP2C15 and CsPP2C12, respectively, while the others were one-to-one paired.
In addition, we also detected homologous PP2C gene pairs between cucumber and Arabidopsis. There were 59 collinear gene pairs between 41 CsPP2Cs and 48 AtPP2Cs (Fig. 2 b). The maximum number of homologous genes in cucumber was 11 pairs on chromosome 3, while the minimum number was three pairs on chromosome 5. According to this result, we speculated that cucumber and Arabidopsis may have high homology.
Analysis of dN/ds values between collinear gene pairs
To further investigate the divergence and selection in duplication of PP2C genes, the non-synonymous substitution rate (dN), synonymous substitution rate (ds), and dN/ds values were evaluated for the homologous gene pairs among cucumber and Arabidopsis (Table S1). Where dN/dS > 1 is the positive selection, dN/dS = 1 is the neutral selection, and 0 < dN/dS < 1 is the purifying selection . The dN/ds value of all cucumber gene pairs was less than 1. Similarly, the dN/ds value of all collinear gene pairs in cucumber and Arabidopsis was less than 1. These data suggest that these genes were primarily under purifying selection during evolution and could help to maintain the basic function of this gene.
Phylogenetic analysis of CsPP2C genes
To investigate the phylogenetic relationships of PP2C genes between cucumber and Arabidopsis, we used the maximum likelihood method to construct a phylogenetic tree based on 80 PP2C genes in Arabidopsis and 56 in cucumber (Fig. 3). The phylogenetic analyses indicated that each subfamily includes PP2C proteins from cucumber and Arabidopsis, and that the genes of cucumber and Arabidopsis tended to form independent branches in each subgroup, that is, cucumber genes clustered together, and Arabidopsis genes clustered together. The 56 CsPP2C proteins were divided into 13 subgroups (A-L), while CsPP2C1, 21, and 11 were not clustered with any other group. This was similar to the grouping of PP2C in Arabidopsis and rice. Each group included seven, four, four, nine, nine, five, three, four, three, three, zero, zero, and two CsPP2C genes. Aside from subgroup J and K (only the AtPP2C gene), the distribution of PP2Cs in cucumber and Arabidopsis subgroups was similar. This suggests that the PP2C gene family may have evolved from a common ancestor.
Gene structural and conserved domain analyses of CsPP2Cs
Since the pattern diversity of exon/intron structure and protein domain plays an important role in the evolution of gene families, we studied the exon/intron structure patterns of CsPP2C genes and conserved domain based on their phylogenetic relationships (Fig. 4 a). Studies on exon/intron structure showed that most members of the same subfamily have similar exon/intron numbers but differ in length (Fig. 4 b). The structure of the CsPP2C gene in each group was mostly similar, but there were differences in the exon/intron arrangement of some genes. For example, in group F, the 3′-end of CsPP2C50 has the longest non-coding region and the number of exons in group F2 (8 exons) was almost twice that in group F1 (4–5 exons). In addition, CsPP2C35 in group I has 10 exons, and its longest gene fragment is more than 12 kb long. CsPP2C40 in group G had no non-coding regions, and only two exons, the lowest number of exons in all PP2C genes. CsPP2C14, 40, 3, 24, 37, 15, and 51 have no upstream and downstream gene sequences, and most were located in group D. This indicated that the CsPP2C gene was relatively conservative in its evolution, ensuring the integrity of its gene structure so that there is little change in its function.
To identify common motifs among different groups of CsPP2C proteins, we used the MEME motif search tool to identify 10 conserved motifs (Table 2). As shown in Fig. 4. c, proteins in the same group exhibited similar motif distribution patterns. Motifs 1, 2 (except CsPP2C3, 23), 3 (except CsPP2C51), 4, 6, and 7 were found in all CsPP2C genes. In addition to the common motif, there were specific motifs in each group. For example, motif 8 was not present in group C, but was present in all other groups. Motif 5 was present in group C and D, but not in the other groups. Motif 9 was not found in group C, D, and H, while it was found in all other groups. According to these results, the CsPP2C genes in the same subgroup had similar conserved motif composition and distribution, suggesting that the CsPP2C members in the same cluster likely share similar functions.
Cis-element analysis of the CsPP2Cs promoter in cucumber
Abundant responsive regulatory elements were found in the promoter regions of CsPP2Cs through PlantCARE analysis (Fig. 5). The cis-elements screened were divided into two categories. The first type of element was hormone response, such as TCA-element (salicylic acid response element), the ABRE (ABA response element), the TGA-element (auxin response element), the CGTCA-motif and TGACG-motif (MeJA-responsiveness response element), and the P-box (gibberellin response element), among others. The second type of element was stress response, such as LTR, MYC, TC-rich repeats and MBS. The abscisic acid ABA-responsive (ABRE) elements were identified abundantly in the promoter regions of CsPP2Cs, among which CsPP2C2 and CsPP2C3 contained 12 ABREs, which was the largest number, followed by CsPP2C7, CsPP2C20, CsPP2C26, CsPP2C40, and CsPP2C56; it was the most abundant element in the promoter region. The second was the MYC element, which was not present in only CsPP2C12–13, CsPP2C32–34, CsPP2C42–43, CsPP2C45, 46, and 48, and was most abundant in CsPP2C1. This suggests that most CsPP2C genes may respond to various abiotic stresses.
Tissue-specific expression profiles of CsPP2C genes
To better understand the role of CsPP2C genes in cucumber growth and development, the temporal and spatial expression patterns of CsPP2C genes was analyzed using RNA-Seq data of different tissues of the cucumber cultivar Chinese long ‘9930’ (Fig. 6). Only CsPP2C11, 41, 5, 33, and 50 had low expression in all 10 tissues. On the contrary, the expression levels of other CsPP2C genes were high in the fertilized ovaries, male, female, and leaf but low in other organs, such as CsPP2C12, 51, 15, 46, 37, 31, 22, and 47. Moreover, the expression level of CsPP2C49 was medium in males but low in other tissues. Similarly, the expression level of CsPP2C53 was medium in females and males but low in other tissues. In general, most cucumber CsPP2C genes showed similar expression patterns in different tissues.
Changes in relative expression of the CsPP2C gene under osmotic stress and ABA treatment
Several members of group A PP2Cs have been shown to function as negative regulators of the ABA signaling pathway in Arabidopsis. The expression of seven PP2C genes in Arabidopsis was suppressed by ABA treatment, two of which are members of subfamily D. The PP2C genes in different subfamilies might play different functional roles in distinct signaling pathways. Therefore, the expression of 56 CsPP2C genes under ABA, salt, drought, and cold treatment was analyzed by qRT-PCR. Under ABA treatment (Fig. 7 a), we observed that the expression of CsPP2C1–3, 13, 22, 29, 39, 45, 54, and 55 continually increased with the extension of ABA treatment time. Among them, the expression level of CsPP2C3 was the highest at 24 h, 20 times that at 0 h. On the contrary, the expression levels of CsPP2C4, 5, 9, 33–38, 41–43, and 46–53 had decreasing trends. The expression levels of CsPP2C6, 8, 11, 14–17, 19, 20, 23, 27–29, and 40 were increased at 6 and 12 h but decreased at 24 h; among them, CsPP2C14–17 was up-regulated most significantly, and although their expression decreased at 24 h, was still higher than that at 0 h. Compared with 0 h, the expression of CsPP2C19, 24 was up-regulated at 6 h, the expression of CsPP2C10, 12, 24, 30, and 32 was up-regulated at 12 h, while the expression of CsPP2C25 and 31 was almost unchanged. As shown in Fig. 7 b, under 10% PEG treatment, CsPP2C1, 3, 11–17, 23, 45, and 55 were significantly up-regulated, among which CsPP2C3, CsPP2C11–17 and CsPP2C45 were up-regulated most noticeably and were more than eight times higher than that at 0 h after 6 h treatment, especially CsPP2C45, which was more than 16 times higher than that at 0 h after 24 h treatment. The relative expression levels of other genes such as CsPP2C9, 18, 21, 34–36, and 46–53 were down-regulated. It is worth noting that after 200 mmol/L NaCl treatment (Fig. 7 c), the expression of many genes increased significantly after 12 and 24 h treatment and were higher than other treatments, especially CsPP2C1, CsPP2C6–8, CsPP2C11, CsPP2C12, CsPP2C17, and CsPP2C20–32. The highest was more than 15 times higher than that at 0 h after 24 h treatment (CsPP2C23). The relative expression of other genes was similar to that of drought treatment. Interestingly, the expression of almost all CsPP2C genes under cold stress (Fig. 7 d) was lower than that of other treatments, except for the genes whose expression was down-regulated, such as CsPP2C1, 3, 11–17, 23, 45, and 54–55. Compared with that at 0 h, the up-regulation of these genes was not high, and CsPP2C14 had the highest up-regulation range at 12 h, more than eight times that of 0 h. Under ABA, drought, and salt treatments, the expression level of CsPP2C3, 11–17, 23, 45, 54, and 55 was significantly up-regulated. This may play an important role in abiotic stress. These results provide a basis for the functional study of the PP2C gene in the future.
Plant PP2Cs play a role in plant development as well as in drought, salt, alkali, and fungal pathogen stress resistance [24, 35] and are important PKs in the ABA signaling pathway . In recent years, they have been identified in several plants, such as Arabidopsis , rice , maize , and wheat . In this study, we comprehensively analyzed CsPP2C genes in cucumber, including genome-wide identification, chromosome location, collinear relationship, gene structure, conserved motifs, and expression patterns. Fifty-six cucumber PP2C genes were identified by homology comparison. Compared with Arabidopsis (80) , wheat (95) , maize (97) [16, 38], rice (78)  and B. distachyon (86) , the amount of PP2C in cucumber was much less. This indicated that the increase and expansion of PP2C genes are diverse in different species and may also be related to their adaptation to complex environments or to the few chromosomes and small genome of cucumber (2n = 2x = 14) .
Chromosome mapping (Fig. 1) found that the PP2C genes located on Chr 3, 4, and 6 were the most distributed, with the number of 13, 9 and 11, respectively. The CsPP2C gene does not cluster on the chromosome, and we did not find tandem duplicate gene pairs. Collinearity results (Fig. 2 a) showed that seven homologous gene pairs existed among CsPP2Cs, and there were 59 pairs of collinear genes in Arabidopsis and cucumber (Fig. 2 b), indicating that CsPP2C genes expand primarily through segmental duplication of the chromosome. Compared with tandem repeats, segmental duplication is more conducive to the maintenance of gene function in the process of gene replication . In addition, the amplification of AtPP2C genes in Arabidopsis is also through fragment replication , which was consistent with our results. Furthermore, the dN/dS values of these gene replication pairs were all less than 1, indicating that they underwent purification selection, and their high conservatism in the evolution process was predicted (Table S1).
According to the phylogenetic tree, the PP2C genes of cucumber and Arabidopsis were divided into 13 groups (Fig. 3). This grouping is the same as that in Arabidopsis, rice, and wheat [16, 37]. Phylogenetic tree analysis showed that most subfamilies contained both cucumber and Arabidopsis genes. Cucumber and Arabidopsis had a similar number of members in the same subfamily, and the members of the two species tended to cluster separately. In addition, we analyzed the gene structure and conserved motifs of CsPP2C genes. The exon/intron structure of genes is an important marker of the evolutionary relationship among gene family members . Our results indicated that cucumber PP2C genes in the same group had similar exon/intron structures (Fig. 4 b). While the distribution pattern of exon/intron structures was similar to most genes in the same group, there were some exceptions, which may be due to a variety of reasons . Unlike other genes, CsPP2C14, 40, 3, 24, 37, 15 and 51 have no non-coding region. In addition, all CsPP2C genes contain different numbers of exons and introns, of which CsPP2C35 has the longest coding sequence and CsPP2C21 has the largest number of exons. A total of 10 conserved motifs were identified in the amino acid sequence of CsPP2C genes. The results (Fig. 4 c) showed that CsPP2C genes in the same group showed similar motif distribution, and motifs 1, 2 (except CsPP2C3, 23), 3 (except CsPP2C51), 4, 6 and 7 existed in all cucumber PP2C genes. This motif pattern was closely related to the catalytic core domain of the PP2C protein .
Understanding the subcellular location of proteins (Table 1) may provide us with necessary information to infer the biological function of proteins. CsPP2Cs were primarily located in the nucleus, cytoplasm and chloroplast, thus it is speculated that they are related to photosynthesis, respiration, and cell growth and development. At the same time, the CsPP2C gene also showed specific expression in different tissues (Fig. 6). The expression of all CsPP2C genes in female, leaf, male, and fertilized ovaries was higher (except CsPP2C5, 33, and 41), but lower in other tissues. In cotton, the majority of the GhPP2CAs was predominantly expressed in flowers . In Arabidopsis and rice, most AtPP2C genes (84%) and OsPP2C genes (72%) are expressed in more than two tissues . In maize, most ZmPP2Cs had a very broad expression spectrum; only ZmPP2C51, ZmPP2C87, ZmPP2C86, ZmPP2C80, ZmPP2C58 and ZmPP2C41 were low or not expressed in different tissues . This shows that most PP2C gene family members in different plants play a role in multiple processes of plant growth and development, and that only a small number of members act in tissue-specific biological processes.
Cis-acting elements are important regulators of hormone response in plant development and resistance to various stresses . In this study, ABRE elements were most abundant in all CsPP2C gene promoter regions, followed by MYC elements (Fig. 5). qRT-PCR results (Fig. 7) showed that CsPP2C1, 3, 6, 8, 11–17, 22, 23, 45, 54, and 55 responded significantly to ABA and drought treatment; CsPP2C1, 3, 6, 8, 11, 14–17, 54 and 55 contain many MYC elements; and CsPP2C3, 13, 15–17, 22, and 23 contain many ABA response elements. These results indicated that the number of homeopathic elements in the promoter region of CsPP2C genes may be related to its response to stress.
Studies on Arabidopsis and rice suggest that family A PP2C plays an important role in plant responses to abiotic stress, especially in ABA signaling . Overexpression of group E PP2C gene, AtPP2CF1, increased plant biomass in Arabidopsis . AtPP2CG1 in group G is an ABA-dependent positive regulator of salt tolerance . In addition, the expression patterns of PP2C genes in maize , rice  and B. distachyon  under various stresses have been examined. In B.distachyon, BdPP2C70 from subgroup D, BdPP2C13 from subgroup F and BdPP2C32 from subgroup G exhibited strongly increased expression levels in response to ABA and abiotic treatments . In rice, the expression levels of all OsPP2C genes in subgroup A increased after ABA and salt treatment . Our study showed that (Fig. 7) the relative expressions of CsPP2C3 from subgroup A, CsPP2C14, CsPP2C13, from subgroup B, CsPP2C12, CsPP2C15, CsPP2C54 from subgroup D, CsPP2C17, CsPP2C45 from subgroup E, CsPP2C55 from subgroup G, and CsPP2C23 from subgroup I were up-regulated under the four treatments. In addition, the expression of many CsPP2C genes showed a similar trend under four different treatments, which was consistent with research results for Medicago truncatula . In cucumber, subgroup A CsPP2Cs includes seven members (CsPP2C2, CsPP2C3, CsPP2C7, CsPP2C20, CsPP2C26, CsPP2C41, and CsPP2C48). The qRT-PCR results suggested that only CsPP2C3 was highly induced by exogenous ABA treatment, and that CsPP2C41 and CsPP2C48 showed a downward trend, while other genes were up-regulated but not noticeable relative to 0 h. In eight members of B.distachyon subgroup A, BdPP2C27 and BdPP2C34 were also insensitive to ABA treatment. These results showed that not all CsPP2C genes homologous to group A members of Arabidopsis responded to ABA treatment as well as CsPP2C genes in other groups, and that there may be different expression patterns of PP2C genes in different species. Under cold stress, the relative expression level of most CsPP2C genes was lower (Fig. 7d), which was not significant, compared with that of other treatments. Many genes are silenced or reduced in expression. This indicated that the PP2C gene of cucumber had a weak response to cold stress. In conclusion, CsPP2C genes in different groups have different responses to different stresses, and there were CsPP2CAs genes that did not respond to ABA treatment. The expression of CsPP2C3, 11–17, 23, 45, 54, and 55 was significantly up-regulated under the four treatments, suggesting that these genes from different subgroups may play an important role in cucumber resistance to abiotic stress. These results provide a reference for the study of PP2C under different stresses; their functions need to be further explored.
In this study, the whole genome of the cucumber CsPP2C gene family was identified, and its expression level was analyzed. Fifty-six CsPP2C genes were highly similar in gene structure and had conserved motifs. Collinearity and selection pressure indicated that CsPP2C genes were amplified by fragment replication and underwent purifying selection during evolution, ensuring the stability of their functions. In addition, qRT-PCR results showed that in subgroup A, only CsPP2C3 responded significantly to ABA and other treatments, other genes were also up-regulated but not significant, while the expressions of CsPP2C41 and CsPP2C48 were down-regulated. In addition, the members of other subgroups also have genes that respond significantly to abiotic stress, such as CsPP2C11–17, CsPP2C23, 45, 54, and 55. These genes may play an important role in cucumber growth and development. This study provided relevant information for follow-up study of PP2C gene function in cucumber.
Plant materials and growth conditions
The germinated seeds (Cucumis sativus L, “L306” cultivar) were treated by 5% NaClO for 10 min and washed with deionized water 3–5 times. Afterwards, the cucumber seeds were put into an artificial climate incubator to promote germination. The cultivation conditions were as follows: relative humidity 80%, temperature 25 °C/18 °C (day/night), with a light intensity during the day of 250 μmol m− 2 s− 1. After 4 days of culture, cucumber seedlings were soaked in Yamazaki nutrient solution for hydroponics. When cucumber seedlings grew to three true leaves, four treatments were set, T1: 100 μmol/L ABA; T2: 10% PEG; T3: 200 mmol/L NaCl; T4: cold stress. Cucumber leaves after 0, 6, 12 and 24 h of different treatments were collected, and the samples were put into liquid nitrogen and immediately stored in a refrigerator at − 80 °C.
Identification of PP2C genes in cucumber
The protein sequences of 80 PP2C genes in Arabidopsis were downloaded from the Arabidopsis information resource website (https://www.arabidopsis.org/). Then, cucumber genome data (http://cucurbitgenomics.org/) was used for BLASTp. Next, screening of candidate genes was performed, and the retrieval threshold was set as E-value < E− 10. Redundant results were removed manually. At the same time, the protein sequences of the Arabidopsis PP2C gene were compared in TBtools to further screen candidate genes. Then, the Pfam (http://pfam.xfam.org/search#tabview=tab1) and NCBI-CDD databases (https://www.ncbi.nlm.nih.gov/cdd/) were used for domain identification of candidate gene sequence structures. Candidate genes that did not contain the specific domain of the PP2C gene (registration number: PF00481) were manually removed .
Sequence analysis and basic information of the cucumber PP2C gene family
The subcellular localization of the PP2C protein in cucumber was predicted by Wolf PSORT (https://wolfpsort.hgc.jp/) . The number of amino acids and isoelectric genes were identified using ExPASy (websitehttps://web.expasy.org/protparam/) . Additionally, the secondary structure of the cucumber PP2C protein were predicted through the online web site (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl? page = npsa_sopma. html).
Analysis of chromosome location and collinearity analysis
The Chinese long ‘9930’ v3 gff3 file was downloaded from the cucumber genome database, and the chromosome location and length information of 56 PP2C genes were screened for and classified as CsPP2Cs according to their distribution on chromosomes. MG2C (http://mg2c.iask.in/mg2c_v2.0/) was used to map the chromosomal location. Using MCScanX search for homology, the protein-coding genes from the cucumber genome were compared against themselves and those from Arabidopsis genomes using BLASTp, and the retrieval threshold was set as E-value < E− 5. Others were not modified by default parameters. Whole-genome BLASTp results were used to compute collinear blocks for all possible pairs of chromosomes and scaffolds . Subsequently, TBtools was used to highlight the identified PYL collinear pairs and their collinear pairs with Arabidopsis .
Selective pressure analysis of the cucumber PP2C genes family
Construction of phylogenetic tree
The phylogenetic tree of the PP2C gene family of cucumber and Arabidopsis thaliana was constructed using MEGA 7 with muscle-sequence alignment methods selected and the Bootstrap value set as 1000 . Finally, a stable minimum neighbor tree was selected to represent their evolutionary relationship and the constructed tree was beautified using iTOL (https://itol.embl.de/).
Analysis of gene exon-intron structures and protein conserved motifs
The structure of the PP2C gene in cucumber was analyzed by GSDS (http://gsds.gao-lab.org/) . The conserved motif of the cucumber PP2C gene was analyzed by MEME (http://meme-suite.org/tools/meme) . The maximum motif number was set as 10, and the remaining parameters were set as default values.
Analysis of cis-acting elements in CsPP2C gene promoters
The first 2000 bp upstream sequence of 56 identified PP2C gene initiation codons (ATG) was extracted using TBtools. The 2000 bp sequence was then submitted to PlantCARE (http://bioinformatics.psb.ugent.be/ webtools/plantcare/html/) for cis role element forecast analysis .
Tissue expression analysis of PP2C genes
To study the gene-specific expression of CsPP2C genes in different tissues of cucumber, the accession number PRJNA80169 was used from cucumber genome data (http://cucurbitgenomics.org/) to obtain cucumber RNA samples from different tissue and organ (tendril-base, tendril, root, leaf, stems, ovary-unfertilized, ovary-fertilized, ovary) RNA-Seq data . Finally, the log2 method was used for data conversion and TBtools was used to draw the expression heat map of the CsPP2C genes.
RNA extraction and real-time PCR (qRT-PCR)
Total RNA was isolated from collected samples using a Plant RNA Extraction kit (Tiangen, China). The complementary DNA was synthesized using the fastking cDNA dispersing RT supermaxs kit (Tiangen, China) with 2 μL RNA as the template. The coding sequences of CsPP2C genes were input into the homepage of Shanghai biology company (Shanghai, China) for online primer design (Table S2). The SYBR Green kit (Tiangen, China) was used for fluorescence quantitative analysis. The volume of the reaction system was 20 μL, containing 2 μL cDNA solution, 10 μL 2*SuperReal PreMix Plus, 0.6 μL of 10 μM forward and reverse primers, 0.4 μL 50*ROX Reference Dye, and 6.4 μL of distilled deionized water. Next, qRT-PCR was performed using the LightCycler® 480 II real-time fluorescence quantitative PCR instrument. The amplification program conditions were as follows: 95 °C for 15 min, 40 cycles of 95 °C for 10 s and 60 °C for 30 s. Each sample was replicated three times. The relative expressions levels of the PP2C genes were calculated using the 2-∆∆CT method . SPSS 20.0 was used to analyze the relative expressions and Origin 9.0 was used to complete the histogram of relative expressions.
Availability of data and materials
Wang L, Liang S, Lu Y-T. Characterization, physical location and expression of the genes encoding calcium/calmodulin-dependent protein kinases in maize (Zea mays L.). Planta. 2001;213(4):556–64.
Kong X, Lv W, Jiang S, Zhang D, Cai G, Pan J, et al. Genome-wide identification and expression analysis of calcium-dependent protein kinase in maize. BMC Genomics. 2013;14(1):1–15.
Li F, Fan G, Wang K, Sun F, Yuan Y, Song G, et al. Genome sequence of the cultivated cotton Gossypium arboreum. Nat Genet. 2014;46(6):567–72.
Cohen P. The structure and regulation of protein phosphatases. Annu Rev Biochem. 1989;58(1):453–508.
Luan S. Protein phosphatases in plants. Annu Rev Plant Biol. 2003;54(1):63–92.
MacKintosh C, Coggins J, Cohen P. Plant protein phosphatases. Subcellular distribution, detection of protein phosphatase 2C and identification of protein phosphatase 2A as the major quinate dehydrogenase phosphatase. Biochem J. 1991;273(3):733–8.
Zhou H-W, Nussbaumer C, Chao Y, DeLong A. Disparate roles for the regulatory a subunit isoforms in Arabidopsis protein phosphatase 2A. Plant Cell. 2004;16(3):709–22.
Smith RD, Walker JC. Expression of multiple type 1 phosphoprotein phosphatases in Arabidopsis thaliana. Plant Mol Biol. 1993;21(2):307–16.
Takemiya A, Kinoshita T, Asanuma M, Shimazaki K-I. Protein phosphatase 1 positively regulates stomatal opening in response to blue light in Vicia faba. Proc Natl Acad Sci. 2006;103(36):13549–54.
Lin Q, Li J, Smith RD, Walker JC. Molecular cloning and chromosomal mapping of type one serine/threonine protein phosphatases in Arabidopsis thaliana. Plant Mol Biol. 1998;37(3):471–81.
Zimmerlin A, Jupe SC, Bolwell GP. Molecular cloning of the cDNA encoding a stress-inducible protein phosphatase 1 (PP1) catalytic subunit from French bean (Phaseolus vulgaris L.). Plant Mol Biol. 1995;28(3):363–8.
Cao J, Jiang M, Li P, Chu Z. Genome-wide identification and evolutionary analyses of the PP2C gene family with their expression profiling in response to multiple stresses in Brachypodium distachyon. BMC Genomics. 2016;17(1):1–17.
Song S-K, Hofhuis H, Lee MM, Clark SE. Key divisions in the early Arabidopsis embryo require POL and PLL1 phosphatases to establish the root stem cell organizer and vascular axis. Dev Cell. 2008;15(1):98–109.
Schweighofer A, Hirt H, Meskiene I. Plant PP2C phosphatases: emerging functions in stress signaling. Trends Plant Sci. 2004;9(5):236–43.
Meskiene I, Baudouin E, Schweighofer A, Liwosz A, Jonak C, Rodriguez PL, et al. Stress-induced protein phosphatase 2C is a negative regulator of a mitogen-activated protein kinase. J Biol Chem. 2003;278(21):18945–52.
Xue T, Wang D, Zhang S, Ehlting J, Ni F, Jakab S, et al. Genome-wide and expression analysis of protein phosphatase 2C in rice and Arabidopsis. BMC Genomics. 2008;9(1):1–21.
Schweighofer A, Kazanaviciute V, Scheikl E, Teige M, Doczi R, Hirt H, et al. The PP2C-type phosphatase AP2C1, which negatively regulates MPK4 and MPK6, modulates innate immunity, jasmonic acid, and ethylene levels in Arabidopsis. Plant Cell. 2007;19(7):2213–24.
Gagne JM, Clark SE. The Arabidopsis stem cell factor POLTERGEIST is membrane localized and phospholipid stimulated. Plant Cell. 2010;22(3):729–43.
Chan Z. Expression profiling of ABA pathway transcripts indicates crosstalk between abiotic and biotic stress responses in Arabidopsis. Genomics. 2012;100(2):110–5.
Merlot S, Gosti F, Guerrier D, Vavasseur A, Giraudat J. The ABI1 and ABI2 protein phosphatases 2C act in a negative feedback regulatory loop of the abscisic acid signalling pathway. Plant J. 2001;25(3):295–303.
Boudsocq M, Barbier-Brygoo H, Lauriere C. Identification of nine sucrose nonfermenting 1-related protein kinases 2 activated by hyperosmotic and saline stresses in Arabidopsis thaliana. J Biol Chem. 2004;279(40):41758–66.
Smékalová V, Doskočilová A, Komis G, Šamaj J. Crosstalk between secondary messengers, hormones and MAPK modules during abiotic stress signalling in plants. Biotechnol Adv. 2014;32(1):2–11.
Wu P, Wang W, Li Y, Hou X. Divergent evolutionary patterns of the MAPK cascade genes in Brassica rapa and plant phylogenetics. Hortic Res. 2017;4:17079.
Chen C, Yu Y, Ding X, Liu B, Duanmu H, Zhu D, et al. Genome-wide analysis and expression profiling of PP2C clade D under saline and alkali stresses in wild soybean and Arabidopsis. Protoplasma. 2018;255(2):643–54.
Haider MS, Kurjogi MM, Khalil-Ur-Rehman M, Fiaz M, Pervaiz T, Jiu S, et al. Grapevine immune signaling network in response to drought stress as revealed by transcriptomic analysis. Plant Physiol Biochem. 2017;121:187–95.
Galbiati M, Simoni L, Pavesi G, Cominelli E, Francia P, Vavasseur A, et al. Gene trap lines identify Arabidopsis genes expressed in stomatal guard cells. Plant J. 2008;53(5):750–62.
Lee MW, Jelenska J, Greenberg JT. Arabidopsis proteins important for modulating defense responses to pseudomonas syringae that secrete HopW1-1. Plant J. 2008;54(3):452–65.
Singh A, Giri J, Kapoor S, Tyagi AK, Pandey GK. Protein phosphatase complement in rice: genome-wide identification and transcriptional analysis under abiotic stress conditions and reproductive development. BMC Genomics. 2010;11(1):1–18.
Fan K, Yuan S, Chen J, Chen Y, Li Z, Lin W, et al. Molecular evolution and lineage-specific expansion of the PP2C family in Zea mays. Planta. 2019;250(5):1521–38.
Singh A, Jha SK, Bagri J, Pandey GK. ABA inducible rice protein phosphatase 2C confers ABA insensitivity and abiotic stress tolerance in Arabidopsis. PLoS One. 2015;10(4):e0125168.
Liu L, Hu X, Song J, Zong X, Li D, Li D. Over-expression of a Zea mays L. protein phosphatase 2C gene (ZmPP2C) in Arabidopsis thaliana decreases tolerance to salt and drought. J Plant Physiol. 2009;166(5):531–42.
Jia H-f, Lu D, Sun J-H, Li C-I, Xing Y, Qin L, et al. Type 2C protein phosphatase ABI1 is a negative regulator of strawberry fruit ripening. J Exp Bot. 2013;64(6):1677–87.
Cheung J, Estivill X, Khaja R, MacDonald JR, Lau K, Tsui L-C, et al. Genome-wide detection of segmental duplications and potential assembly errors in the human genome sequence. Genome Biol. 2003;4(4):1–10.
Yadav CB, Bonthala VS, Muthamilarasan M, Pandey G, Khan Y, Prasad M. Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database. DNA Res. 2015;22(1):79–90.
Akimoto-Tomiyama C, Tanabe S, Kajiwara H, Minami E, Ochiai H. Loss of chloroplast-localized protein phosphatase 2Cs in Arabidopsis thaliana leads to enhancement of plant immunity and resistance to Xanthomonas campestris pv. Campestris infection. Mol Plant Pathol. 2018;19(5):1184–95.
Park S-Y, 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.
Yu X, Han J, Wang E, Xiao J, Hu R, Yang G, et al. Genome-wide identification and homoeologous expression analysis of PP2C genes in wheat (Triticum aestivum L.). Front Genet. 2019;10:561.
Wei K, Pan S. Maize protein phosphatase gene family: identification and molecular characterization. BMC Genomics. 2014;15(1):1–20.
Yang L, Koo DH, Li Y, Zhang X, Luan F, Havey MJ, et al. Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly. Plant J. 2012;71(6):895–906.
Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;290(5494):1151–5.
Cannon SB, Mitra A, Baumgarten A, Young ND, May G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004;4(1):1–21.
Long M, Betrán E, Thornton K, Wang W. The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003;4(11):865–75.
Rogozin IB, Sverdlov AV, Babenko VN, Koonin EV. Analysis of evolution of exon-intron structure of eukaryotic genes. Brief Bioinform. 2005;6(2):118–34.
Shi Y. Serine/threonine phosphatases: mechanism through structure. Cell. 2009;139(3):468–84.
Lu T, Zhang G, Wang Y, He S, Sun L, Hao F. Genome-wide characterization and expression analysis of PP2CA family members in response to ABA and osmotic stress in Gossypium. PeerJ. 2019;7:e7105.
Nishimura N, Sarkeshik A, Nito K, Park SY, Wang A, Carvalho PC, et al. PYR/PYL/RCAR family members are major in-vivo ABI1 protein phosphatase 2C-interacting proteins in Arabidopsis. Plant J. 2010;61(2):290–9.
Bhaskara GB, Nguyen TT, Verslues PE. Unique drought resistance functions of the highly ABA-induced clade a protein phosphatase 2Cs. Plant Physiol. 2012;160(1):379–95.
Sugimoto H, Kondo S, Tanaka T, Imamura C, Muramoto N, Hattori E, et al. Overexpression of a novel Arabidopsis PP2C isoform, AtPP2CF1, enhances plant biomass production by increasing inflorescence stem growth. J Exp Bot. 2014;65(18):5385–400.
Liu X, Zhu Y, Zhai H, Cai H, Ji W, Luo X, et al. AtPP2CG1, a protein phosphatase 2C, positively regulates salt tolerance of Arabidopsis in abscisic acid-dependent manner. Biochem Biophys Res Commun. 2012;422(4):710–5.
Yang Q, Liu K, Niu X, Wang Q, Wan Y, Yang F, et al. Genome-wide identification of PP2C genes and their expression profiling in response to drought and cold stresses in Medicago truncatula. Sci Rep. 2018;8(1):1–14.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39(suppl_2):W29–37.
Xiong E, Zheng C, Wu X, Wang W. Protein subcellular location: the gap between prediction and experimentation. Plant Mol Biol Report. 2016;34(1):52–61.
Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, De Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40(W1):W597–603.
Wang Y, Tang H, DeBarry JD, Tan X, Li J, Wang X, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49–9.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.
Goldman N, Yang Z. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994;11(5):725–36.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Hu B, Jin J, Guo A-Y, Zhang H, Luo J, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015;31(8):1296–7.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(suppl_2):W202–8.
Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.
Li Z, Zhang Z, Yan P, Huang S, Fei Z, Lin K. RNA-Seq improves annotation of protein-coding genes in the cucumber genome. BMC Genomics. 2011;12(1):1–11.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.
This experiment thanks the support and help of the Key Laboratory of Arid Habitat Crop Science, China.
This research was supported by Director Foundation of Key Laboratory of Arid Habitat Crop Science, China (GHSJ-2020-Z3).
All methods were carried out in accordance with relevant guidelines and regulations.
Ethics approval and consent to participate
The experimental research and field studies on plants, including the collection of plant material, were carried out in accordance with relevant institutional, national, and international guidelines and legislation.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, G., Zhang, Z., Luo, S. et al. Genome-wide identification and expression analysis of the cucumber PP2C gene family. BMC Genomics 23, 563 (2022). https://doi.org/10.1186/s12864-022-08734-y
- Protein phosphatase 2C
- Stress response