- Research article
- Open Access
Functional characterization and regulatory mechanism of wheat CPK34 kinase in response to drought stress
BMC Genomics volume 21, Article number: 577 (2020)
Drought is one of the most adverse environmental factors limiting crop productions and it is important to identify key genetic determinants for food safety. Calcium-dependent protein kinases (CPKs) are known to be involved in plant growth, development, and environmental stresses. However, biological functions and regulatory mechanisms of many plant CPKs have not been explored. In our previous study, abundance of the wheat CPK34 (TaCPK34) protein was remarkably upregulated in wheat plants suffering from drought stress, inferring that it could be involved in this stress. Therefore, here we further detected its function and mechanism in response to drought stress.
Transcripts of the TaCPK34 gene were significantly induced after PEG-stimulated water deficiency (20% PEG6000) or 100 μM abscisic acid (ABA) treatments. The TaCPK34 gene was transiently silenced in wheat genome by using barley stripe mosaic virus-induced silencing (BSMV-VIGS) method. After 14 days of drought stress, the transiently TaCPK34-silenced wheat seedlings showed more sensitivity compared with control, and the plant biomasses and relative water contents significantly decreased, whereas soluble sugar and MDA contents increased. The iTRAQ-based quantitative proteomics was employed to measure the protein expression profiles in leaves of the transiently TaCPK34-silenced wheat plants after drought stress. There were 6103 proteins identified, of these, 51 proteins exhibited significantly altered abundance, they were involved in diverse function. And sequence analysis on the promoters of genes, which encoded the above identified proteins, indicated that some promoters harbored some ABA-responsive elements. We determined the interactions between TaCPK34 and three identified proteins by using bimolecular fluorescent complementation (BiFC) method and our data indicated that TaCPK34directly interacted with the glutathione S-transferase 1 and prx113, respectively.
Our study suggested that the TaCPK34 gene played positive roles in wheat response to drought stress through directly or indirectly regulating the expression of ABA-dependent manner genes, which were encoding identified proteins from iTRAQ-based quantitative proteomics. And it could be used as one potential gene to develop crop cultivars with improved drought tolerance.
Because of the ever-increasing global population and the climate change, the food safety becomes a major issue because it needs an annual increase (approximately 2.4%) in food yields, whereas yields of the major crops wheat, rice, and maize are presently increased merely only at 0.9, 1 and 1.6%, respectively . Abiotic stresses including drought, cold, heat, and salt are major adverse environmental stresses that greatly limit crop yields and quality. Therefore, it is pivotal to understand the molecular mechanisms underlying stress responses,identify important genetic genes and develop crop cultivars with enhanced stress tolerance.
Plants respond to environmental stresses by regulating the expression of numerous stress-responsive genes through complicated signaling pathways. Ca2+ has been recognized as a major second messenger in these pathways, and can perceive intracellular changes in Ca2+ concentrations and translate them into specific phosphorylation events to initiate downstream signaling processes . In plants, several classes of Ca2+ sensors: calcium-dependent protein kinases (CPKs or CDPKs), calmodulins (CaMs) and CaM-like proteins and calcineurin B-like proteins (CBLs) have been reported to involve in the perception and transmission of intracellular Ca2+ concentrations . CaMs and CBLs do not have enzymatic activities and therefore do not directly transmit the Ca2+ signals, whereas CPKs can sense Ca2+ signals and directly mediate a variety of cellular responses .
CPKs, a large multigene family of Ca2+-regulated Ser/Thr protein kinases, contain four-domain structure consisting of a variable N-terminal variable domain, a kinase or catalytic domain, a CPK activation domain, and a CaM-like domain that is composed of four EF hand motifs . CPKs are involved in biotic and abiotic signal transduction pathways through recognizing the Ca2+ signatures [2, 6]. The numbers within this kinase family vary among different plant species . For instance, 34, 31, 30, and 20 CPKs have been identified members in Arabidopsis thaliana, rice, Brachypodium distachyon, and bread wheat, respectively [8,9,10]. Different CPK members could have differential roles in response to environmental stresses and even within the same stress. For example, AtCPK6 and AtCPK10 in Arabidopsis and OsCPK13 in rice have been found to positively regulate the drought tolerance, respectively [11,12,13], whereas AtCPK23 negatively regulates drought and salt tolerance .
Growing evidence has demonstrated the crucial functions of CPKs in plant growth, development, and responses to environmental stresses, whereas the biological functions and regulatory mechanisms of many plant CPKs remain unclear [9, 15]. Bread wheat (Triticum aestivum L., hexaploid species, AABBDD) is one of the most important global cereal crops and offers a rich source of carbohydrates, proteins, minerals, and other nutrients essential for the human diet . Our previous study identified genome-scale protein expression profiles of wheat grains after anthesis. Among the identified proteins, a wheat CPK protein (CPK34) sharply increased as grain development proceeded, inferring that it could be involved in grain filling and environmental stress tolerance . To our knowledge, however, the function and mechanism of TaCPK34 in response to environmental stresses has not been examined [9, 15]. The objectives of this study were (1) to determine the role of TaCPK34 in wheat response to drought stress by using the BSMV-VIGS method and (2) to investigate the potential molecular mechanisms by which TaCPK34 influences drought tolerance by using the isobaric tagging for relative and absolute quantification (iTRAQ) proteomic analysis.
Identification of TaCPK34
During the wheat grain filling period (from 22 April to 31 May) in 2015, our lab performed a quantitative proteomics to characterize the genome-scale grain protein expression profiles of the bread wheat cultivar Zhoumai 18 at different grain filling stages (5, 10, 15, 20, 25, 30, and 35 days after anthesis) under field conditions . During the late stages of wheat grain filling period (from 6 to 31 May, 2015, corresponding to 15 ~ 35 days after anthesis), there was only 25.5 mm rainfall, which was enough for wheat growth and caused drought stress in this growth season (Table S1). From these proteomic data, one protein that matched to TuCPK34 (GenBank accession no. EMS58754) were identified from Triticum urartu (diploid, AA, the progenitor of the A subgenome of hexaploid bread wheat), and its sequence of amino acid alignment showed 87% similarity to Chinese Spring (CS) (Figure S1), suggesting that TaCPK34 was the homologue of TuCPK34 in Triticum aestivum (Figure S1). Compared to the protein abundance of TaCPK34 at 5 days after anthesis (control), its abundance in wheat grains at the subsequent sampling timepoints (15, 20, 30, and 35 days after anthesis) was more than 16-fold higher (Figure S2).
Characterization and subcellular localization of TaCPK34
The cDNA sequence of TaCPK34 was isolated and sequenced (Fig. 1a), its predicated protein had 518 amino acids (Fig. 1a), and contained the three typical domains of plant CPKs: a variable N-terminal domain, a serine/threonine kinase domain, and a CPK activation domain that comprised an inhibitory junction domain and a C-terminal Ca2+-binding domain (Fig. 1b). Phylogenetic tree of TaCPK34 protein was constructed and grouped into four branches (groups I-IV) (Fig. 2), which followed CPKs’ phylogenetic tree in rice and Arabidopsis thaliana [9, 15], using CPKs of the identified wheat (3 CPKs), rice (11 CPKs) and Arabidopsis thaliana (17 CPKs) (Table S2). The TaCPK34 protein shared high similarities (≥95%) with AtCPK32, AtCPK1, and AtCPK10 of group III (Fig. 2).
Expression patterns of TaCPK34 gene in responses to water deficiency and ABA
PEG-stimulated water deficiency (20% PEG6000) and 100 μM abscisic acid (ABA) were applied to wheat seedlings to monitor the transcriptional responses of the TaCPK34 gene. Under water deficiency condition, transcripts of the TaCPK34 gene increased quickly, reached the highest level at 3 h, and then decreased quickly to low levels at 6 h (Fig. 3a). After a 100 μM ABA treatment, the transcripts of the TaCPK34 gene quickly increased, remarkably higher than that in control seedlings (0 h) at all sampling timepoints (1, 3, and 6 h) (Fig. 3b).
Role of TaCPK34 in response to water deficiency
Barley stripe mosaic virus-induced silencing (BSMV-VIGS) method was used to verify the function of TaCPK34 under drought stress. A 219 bp-length conserved cDNA fragment (+ 505 to + 723, the translation start site “ATG” is + 1) of TaCPK34 gene, which shared high identity (100%) with the other two copies (Figure S3), was used to construct its silencing vector (BSMV-VIGS-TaCPK34), and BSMV-VIGS-GFP vector was used as the control. The silencing vectors of BSMV-VIGS-TaCPK34 and BSMV-VIGS-GFP were separately inoculated into the leaves of wheat seedlings. At 6 days after inoculation, BSMV-VIGS-induced chlorosis emerged on the leaves of both BSMV-VIGS-TaCPK34- or BSMV-VIGS-GFP-inoculated wheat seedlings (Figure S4A). Transcripts of the TaCPK34 gene were remarkably inhibited in leaves of BSMV-VIGS-TaCPK34-inoculated wheat seedlings with chlorotic symptoms (Figure S4B).
At 6 days after inoculation, inoculated wheat seedlings with chlorotic symptoms were well watered, and then suffered from water loss stress (drought) for subsequent 14 days. At this timepoint, the BSMV-VIGS-TaCPK34-inoculated wheat plants exhibited visual drought-induced damage symptoms such as leaf rolling and wilting, whereas the BSMV-VIGS-GFP-inoculated wheat seedlings remained vigorous (Fig. 4a). Coincident with these qualitative phenotypic changes, plant biomasses and relative water contents in BSMV-VIGS-TaCPK34-inoculated wheat plants were significantly lower than those in BSMV-VIGS-GFP-inoculated wheat seedlings (Fig. 4b and c). However, the contents of soluble sugar and MDA were remarkably higher in the former than those in the latter (Fig. 4d and e).
Molecular mechanism of TaCPK34-mediated drought tolerance
The iTRAQ-based quantitative proteomics was employed to identify differences in the leaf protein profiles between the BSMV-VIGS-TaCPK34- and the BSMV-VIGS-GFP-inoculated wheat seedlings at 14 days after drought stress. A total of 25,017 peptides were identified from 295,102 mass spectra (Table S3). And these corresponded to 6, 103 expressed proteins, which were detected in all four biological replicates in the above treatments (Table S4). Our mass spectrometry proteomics data have been deposited (ProteomeXchange identifier PXD008567). Ratios were used to assess the fold changes in the abundance of proteins between the BSMV-VIGS-TaCPK34- and BSMV-VIGS-GFP-inoculated wheat plants. A volcano plot illustrated the asymmetry between up-regulated (red) and down-regulated (blue) differentially expressed proteins, and the correlation of labeled samples showed good reproducibility between the four biological replicates (Fig. 5a). Fold changes of ≥1.50- or ≤ 0.67-fold were considered to be significant in the abundance of identified proteins, and 51 proteins with significantly altered abundance (P < 0.05 with 95% confidence, FDR < 0.05) were identified between BSMV-VIGS-TaCPK34- and BSMV-VIGS-GFP-inoculated wheat seedlings suffering from the above drought stress (Table 1).
On the basis of their biological and molecular functions, the differentially expressed proteins were classified into two groups of known (32, 62.7%) and unknown (19, 37.3%) functions. The function-known proteins were further divided into six categories: stress and defense (12), signal transduction (2), carbohydrate metabolism (6), photosynthesis (2), protein metabolism (4), and transportation (6) (Fig. 5b).
The promoters of the genes encoding the above function-known proteins were retrieved from the International Wheat Genome Sequencing Consortium (IWGSC) database (RefSeq 1.0 version) from Ensemblplants (http://plants.ensembl.org/index.html). Their sequences were analyzed by using PlantCARE software, and our analysis indicated that the promoter sequences of these genes encoding 19 identified proteins harbored ABA-responsive elements. These genes encoded the dehydrin DHN2, LEA3 protein, low-temperature-induced 65 kDa protein-like isoform, barwin-like, 11 kDa late embryogenesis abundant protein-like, LEA 3, germin-like protein 4–1, Prx113, glutathione S-transferase 1, peroxidase, probable xyloglucan endotransglucosylase, serine/threonine-protein kinase CTR1, putative acyl-activating enzyme 19 isoform X1, aspartyl protease AED3-like, non-specific lipid-transfer protein 4.1-like, non-specific lipid-transfer protein 3-like isoform X1, type 1 non specific lipid transfer protein precursor, short-chain dehydrogenase TIC 32, and type 1 non specific lipid transfer protein precursor (Figure S5).
Three stress and defense-related enzymes/proteins (glutathione S-transferase 1, prx113, and peroxidase), which have been reported to directly function in biotic and abiotic stress tolerance [18, 19], were selected for BiFC assay to determine whether they interacted directly with TaCPK34. In this study, TaCPK34 was fused to the N-terminal fragment of YFP, while glutathione S-transferase 1 (GST1), prx113, or peroxidase were fused to the C-terminal fragment of YFP, respectively. If the two tested proteins associated with one another, a complete fluorescent YFP should be detected. Our results showed that YFP fluorescence was observed at the plasma membrane or in nucleus when TaCPK34 was injected together with YC-GST1 or YC-prx113 (Fig. 6), whereas not with YC-peroxidase (data not indicated). These results indicated that TaCPK34 could directly interact with GST1 in the plasma membrane, and prx113 at the plasma membrane and in the nucleus.
TaCPK34 could be a novel and drought-tolerant plant CPK protein
Evolutionary analyses have revealed that gene duplication (e.g., tandem duplication, segmental duplication) has played an important role in the expansion and evolution of plant CPKs [20,21,22]. Thus, plant CPKs constitute a multigene family, and there are differences in the numbers and functions of CPK family members in different plant species [9, 15]. Only 20 CPK encoding genes (TaCPK1 ~ 20) have been identified in the very large and complex genome of bread wheat [9, 23]. In the post-genomic era, recently available genomic sequencing data for bread wheat and its relatives [24,25,26,27,28], have helped to explore more members of crucial CPKs. In our previous proteomic data, one peptide was matched to TuCPK34 , and its orthologous TaCPK34 gene was isolated in this study (Fig. 1a). TaCPK34 protein shared low identities (≤ 50%) to other wheat CPKs and belonged to different groups (Fig. 2), suggesting that TaCPK34 might be a novel member of wheat CPKs.
To cope with environmental stress, plants have evolved sophisticated regulatory mechanisms and undergone a series of morphological, physiological, developmental, and cellular changes. At the genetic level, these changes are well controlled byenvironment-tolerant genes, and therefore, the identification of environment-tolerant genes is the first step toward developing the environment-tolerant cultivars . Of the reported 20 wheat CPKs, only 4 (CPK1, 6, 9 and 18) have been found to be drought-responsive, but their drought tolerance has not been tested . In our previous study, the abundance of TaCPK34 protein was upregulated during the late stages of wheat grain filling period (from 6 to 31 May, 2015, corresponding to 15 ~ 35 days after anthesis) under natural field condition (There was only 25.5 mm rainfall, from 6 to 31 May, 2015, thus drought stress occurred), in which wheat plants suffered from various stresses including drought . Transcripts of the TaCPK34 gene were significantly induced by the PEG-stimulated drought stress (Fig. 3a). These changes in the expression patterns at the translational and transcriptional levels implied the requirement of TaCPK34 in wheat plant for drought adaptation. BSMV-VIGS-TaCPK34- and BSMV-VIGS-GFP-inoculated wheat plants displayed a drought sensitive phenotype with inhibited growth and physiological parameters (Fig. 4). These results suggested that TaCPK34 is a wheat drought-tolerant gene, and could be used for one gene resource for breeding wheat cultivars with drought tolerance.
TaCPK34 could regulate expression of the target genes in response to drought stress
Stress stimuli, as well as growth and development, induce a series of cellular events, firstly changes in Ca2+ concentration to generate specific Ca2+ signatures and subsequently activate Ca2+ sensors, then activated Ca2+ sensors, which decoded the Ca2+ signatures and induced downstream reactions, including changes in gene expression [6, 15]. However, there are functional specificities among multiple Ca2+ sensors due to their differential structures, which allow them to sense specific Ca2+ signatures or targeting specific and unique genes/proteins [9, 30]. CPKs are major Ca2+ sensors in plant cells, and play important roles in plants adaptation to stress stimuli by detecting the stress signals and inducing the expression of downstream genes. From the view of genetic engineering, therefore, they could be crucial for conferring stress tolerance in crop plants. To our knowledge, however, only few plant CPKs have functionally been confirmed to mediate drought tolerance in plants and in two model plants, these include Arabidopsis thaliana CPK4, 6, 8, 10, 11, and 21 [4, 11, 13, 31, 32], and Oryza sativa CPK1, 4, 9, and 10 [1, 33,34,35]. Of these CPKs, TaCPK34 had high similarities to AtCPK10, and they were grouped in the same group (III) (Fig. 2), suggesting that they could have similar mechanisms. AtCPK10 has been found to interact with HSP1 and participate the drought tolerance through the modulation of ABA- and Ca2+-mediated stomatal movements . OsCPK4 enhances the drought tolerance through reducing membrane lipid peroxidation in rice, possibly by upregulating a number of genes associated with the lipid metabolism . However, molecular mechanisms underlying their function in drought tolerance are largely unknown. Here, potential molecular mechanisms of TaCPK34-mediated drought tolerance in wheat were explored by using the iTRAQ-based proteomic approach.
Our proteomic experiments identified 51 proteins with significantly altered abundance (Fold changes of ≥1.50- or ≤ 0.67, P < 0.05) in BSMV-VIGS-TaCPK34-inoculated wheat plants (Table 1). These identified proteins were grouped into two types of known and unknown functions. A large proportion (19/51, 37.3%) of the unknown function proteins were not annotated (Fig. 5b, Table 1), possibly because the genome of bread wheat is huge and complex, and many of its genes/proteins have not been characterized. This has slowed research progress in wheat compared with its homologous species such as rice and maize, whose genomes are relatively small and simple . The identified proteins with known function were involved with stress defense, signal transduction, protein metabolism, carbohydrate metabolism, etc. (Table 1). It has been reported that several plant CPKs function in drought stress responses through direct interaction with some proteins. In Arabidopsis, CPK8 has been found to serve as a positive regulator in plant responses to drought stress through the direct interaction with CAT3 protein to regulate its activity . More evidences have indicated that many CPKs indirectly control the expression of a series of stress-responsive genes through various regulators [9, 14]. OsCPK14 has been identified to interact with OsDi19–4, one member of the drought-induced 19 (Di19) transcription factor family, whose overexpressing Arabidopsis plants display enhanced drought tolerance. Di19 bound to the promoters of OsASPG1 and OsNAC18 genes, and thus, OsCPK14 indirectly regulated the expression of these two genes through Di19 . Among our identified proteins of known function, many proteins or their isoforms have been reported to function in drought stress and they included dehydrin , LEAs , peroxidases , GST1 , germin-like proteins , low-temperature-induced 65 kDa protein . These implied that TaCPK34 could directly or indirectly regulate the expression of the genes encoding these proteins. Our BiFC assay showed that TaCPK34 directly interacted with GST1 and prx113 (Fig. 6), while not with peroxidase, suggesting that it could directly regulate the expression of the formers, while indirectly control the latter. The direct or indirect relations between TaCPK34 and other identified proteins need to be determined.
ABA is the most critical phytohormone in plant drought stress responses, and several plant CPKs are involved in the mechanisms of drought stress response through an ABA-dependent manner . The Arabidopsis plants overexpressing VaCPK20 gene showed increased drought tolerance, and enhanced the transcripts of COR47, NHX1, KIN1, and ABF3 genes, which were induced by ABA and drought stresses . In this study, there were ABA-responsive elements (ACGT) of the genes encoding 19 identified proteins (Figure S5). Some of these genes such as dehydrins , LEAs , prx113, and peroxidase , serine/threonine-protein kinase , and non specific lipid transfer protein , are known to be ABA-responsive, inferring that they could be ABA-dependent manner. Based on these results, we propose that TaCPK34 functions in a signal transduction pathway that is involved in wheat response to drought stress (Fig. 7).
Our data showed that transcripts of the TaCPK34 gene were significantly induced by PEG-stimulated drought stress or ABA application, and TaCPK34-silenced wheat plants exhibited more sensitive to drought stress, suggesting that TaCPK34 has a positive regulatory role in response to drought stress. Further proteomics analysis and BiFC assay suggested that TaCPK34 protein mediated the wheat drought tolerance by directly regulate the expression of the target genes.
Plant materis and treatment designing
A bread wheat cultivar Zhoumai 18 was grown at the Experimental Farm of Agricultural Faculty at Henan Agricultural University (34°N, 113°E, and 52 m elevation). Some meteorological parameters during different sampling timepoints are listed in Table S1. The grains of Zhoumai 18 wheat cultivar were sterilized with 0.01% HgCl2 and fully washed with distilled water. The sterilized grains were grown hydroponically in full-strength Hoagland solution at glass dishes (15-cm diameter) with 50 seedlings each . Germinated grains were transferred into an FPG-300C-30D illumination incubator (Ningbo Laifu Technology Co., Ltd., China) with a 16-h light photoperiod (06:00 ~ 20:00 and 20:00 ~ 06:00, day/night), 25 °C/15 °C day/night, light intensity of 250 μmol m2s− 1, and relative humidity of 60/75% (day/night). Our previous study found that, two-week-old wheat seedlings with two fully expanded leaves were in the autotrophic stage and were sensitive to abiotic stresses . Two-week-old wheat seedlings were subjected to 20% polyethylene glycol (PEG6000)-stimulated water deficiency, or 100 μM abscisic acid (ABA). The uppermost fully expanded leaves were sampled at 0, 1, 3, and 6 h after treatments.
Phylogenetic tree construction
ScanProsite (https://prosite.expasy.org/prosite.html) was employed to predicate the conserved domains and motifs of TaCPK34 protein. MEGA 7.0 software was used to construct the phylogenetic tree of wheat, Arabidopsis and rice CPKs proteins.
RNA isolation and PCR or qPCR analysis
RNA extraction and the first strand cDNA synthesis were described in our previous study . Based on the TuCPK34 sequence, a pair of primers was designed to amplify the coding sequence (CDS) of the TaCPK34 gene using the cDNA from Zhoumai 18. The primer sequences and amplified product are provided in Table S5. Transcripts of the TaCPK34 gene were measured by using quantitative real-time PCR (qPCR) method following reaction conditions: 95 °C for 5 min (1 cycle), 95 °C for 10 s, 58 ~ 60 °C for 15 s, and finally 72 °C for 20 s (40 cycles). Relative transcripts abundance of the TaCPK34 gene was calculated using the 2−ΔΔCt method with the wheat β-actin gene as the internal control. qPCR primers are listed in Table S5. qPCR was performed using three independent biological replicates with at least three plants in each.
Bread wheat is an allohexaploid species, and therefore, the majority of wheat genes have three copies and these three copies share over 95% sequence identity across coding regions . Because of functional complementation among these copies, it is difficult to explore gene function through single copy-silencing in this species . In the present study, a 219 bp-length conserved cDNA fragment (+ 505 to + 723, the translation start site “ATG” is + 1) from the identified TaCPK34 copy was used to construct its silencing vector (BSMV-VIGS-TaCPK34). Consequently, in this study, all three copies of the TaCPK34 gene could be silenced simultaneously with the BSMV-VIGS-TaCPK34 vector, and this helped alleviate functional complementation among the three TaCPK34 copies. BSMV-derived GFP-silenced wheat plants were used as the control . Grains of bread wheat cultivar Zhoumai 18 were sown in soil (nutritive soil: vermiculite with 3:2 ratio), and wheat seedlings with two fully expanded leaves and almost identical heights (12 ± 0.2 cm) were inoculated with BSMV-VIGS-TaCPK34 or BSMV-VIGS-GFP virus on the second fully expanded leaves with 10 μL transcript mixtures each leaf, respectively. The BSMV-VIGS method (e.g., viral vector construction, viral RNA transcription, viral inoculation, and TaCPK34-silenced wheat plant measurement) were performed as described in our previous studies . For drought stress, these inoculated seedlings were initially well watered, and then were cultivated in the above illumination incubator without watering for 14 days. Plant fresh and dry weights and relative water contents after water stress were calculated as described previously . And the contents of soluble sugar or malondialdehyde (MDA) in the sampled leaves, which were composed of the inoculated leaves (the second leaves of BSMV-inoculated wheat plants) or the younger fully expanded leaves (the third or fourth leaves of BSMV-inoculated wheat plants) with silencing phenotype after BSMV-inoculated, were determined spectrophotometrically with a UNICO spectrophotometer [UV-2600, UNICO Instruments Co., Ltd., Shanghai, China] as previously described [56, 57]. Experiments were independently repeated three times, and each replicate contained at least five wheat plants.
Total proteins from leaves of BSMV-VIGS-TaCPK34-, or BSMV-VIGS-GFP-inoculated wheat seedlings were extracted at 14 days after water stress using the trichloroacetic acid/acetone method. Four biological replicates were performed and each replicate with at least six plants. iTRAQ proteomics assay wasperformedas described in our previous study . Samples were labeled with 113–116 for BSMV-VIGS-TaCPK34-inoculated wheat seedlings, and 117–119 and 121 for BSMV-VIGS-GFP-inoculated wheat seedlings. The MS/MS spectra were searched against the NCBInr database (https://www.ncbi.nlm.nih.gov/taxonomy/?term=Triticeae, version 10/06/2018, 192,314 entries) to identify the proteins, which were isolated by using MASCOT software (Matrix Science, London, UK; version 2.2). The quantified peptides were calculated by using the ProGroup algorithm (AB SCIEX). Then the identified proteins were screen using several criteria (i) the cRAP database (ftp://ftp.thegpm.org/fasta/cRAP) was used for filtering the identified protein species; (ii) the confidence and the false discovery rate (FDR) of proteins were ≥ 95% and < 0.05, respectively; (iii) the proteins were commonly identified in four biological replicates; and (iv) multiple testing corrections (t < 0.05) were used to adjust P values for controlling the FDR, and ratios were further analyzed by using Benjamini and Hochberg (BH) test. Proteins whose abundances differed significantly by fold-changes of ≥1.5 or ≤ 0.67 between BSMV-VIGS-TaCPK34- and BSMV-VIGS-GFP-inoculated wheat plants were selected and considered to be used for further study.
The interactions between TaCPK34 and the identified proteins were tested by using the BiFC method as described previously . The CDS regions of the genes encoding wheat GST1, Prx113, and peroxidase identified in our proteomic data, were amplified and ligated into the BiFC vector, which contains a yellow fluorescent protein in carboxy-terminal part, to construct a YFP-C-genes (YC-gene). Likewise, the TaCPK34 gene was also inserted into the amino-terminal part of the yellow fluorescent protein (YFPN) to form YN-TaCPK34. The proteins of HAL3 and Di19, which were confirmed to interact with plant CPKs in the cytoplasm and nucleus [38, 58], respectively, were used as positive controls in this study. Then, the transient assays of interactions between TaCPK34 and other proteins were performed by using the epidermal leaf cells of tobacco, and the fluorescence of YFP was observed with confocal laser scanning microscopy (LSM710, Karl Zeiss, Jena, Germany).
All experiments were independently repeated at least three times and each sampling was analyzed separately. SPSS 20.0 software was used for statistical analyses and statistically significant differences were measured by using Duncan’s tests at P < 0.05 level.
Availability of data and materials
The mass spectrometry proteomics data used in this study have been deposited in ProteomeXchange Consortium with the dataset identifier PXD008567 (http://www.ebi.ac.uk/pride/archive/projects/PXD008567) and FTP Download: ftp://ftp.pride.ebi.ac.uk/pride/data/archive/2020/08/PXD008567. Others data used and analyzed during the current study are available in the manuscript and its additional files.
Barley stripe mosaic virus-induced silencing
CaM-like proteins and calcineurin B-like proteins
Calcium-dependent protein kinases
Isobaric tags for relative and absolute quantitation
International Wheat Genome Sequencing Consortium
Quantitative real-time PCR
Ray DK, Mueller ND, West PC, Foley A. Yield trends are insufficient to double global crop production by 2050. PLoS One. 2013;8:e66428.
Schulz P, Herde M, Romeis T. Calcium-dependent protein kinases: hubs in plant stress signaling and development. Plant Physiol. 2013;163:523–30.
Batistič O, Kudla J. Analysis of calcium signaling pathways in plants. BBA-Gen. 1820;2012:1283–93.
Zou JJ, Li XD, Ratnasekera D, Wang C, Liu WX, Song LF, Zhang WZ, Wu WH. Arabidopsis CALCIUM-DEPENDENT PROTEIN KINASE8 and CATALASE3 function in abscisic acid-mediated signaling and H2O2 homeostasis in stomatal guard cells under drought stress. Plant Cell. 2015;27:1445–60.
Liese A, Romeis T. Biochemical regulation of in vivo function of plant calcium-dependent protein kinases (CDPK). BBA-Mol Cell Res. 1833;2013:1582–9.
Romeis T, Herde M. From local to global: CDPKs in systemic defense signaling upon microbial and herbivore attack. Curr Opin Plant Biol. 2014;20:1–10.
Valmonte GR, Arthur K, Higgins CM, MacDiarmid RM. Calcium-dependent protein kinases in plants: evolution, expression and function. Plant Cell Physiol. 2014;55:551–69.
Simeunovic A, Mair A, Wurzinger B, Teige M. Know where your clients are: subcellular localization and targets of calcium-dependent protein kinases. J Exp Bot. 2016;67:3855–72.
Atif RM, Shahid L, Waqas M, Ali B, Rashid MAR, Azeem F, Nawaz MA, Wani SH, Chung G. Insights on calcium-dependent protein kinases (CPKs) signaling for abiotic stress tolerance in plants. Int J Mol Sci. 2019;20:5298.
Wen F, Ye F, Xiao Z, Liao L, Li T, Jia M, Liu X, Wu X. Genome-wide survey and expression analysis of calcium-dependent protein kinase (CDPK) in grass Brachypodium distachyon. BMC Genomics. 2020;21:53.
Xu J, Tian YS, Peng RH, Xiong AS, Zhu B, Jin XF, Gao F, Fu XY, Hou XL, Yao QH. AtCPK6, a functionally redundant and positive regulator involved in salt/drought stress tolerance in Arabidopsis. Planta. 2010;231:1251–60.
Saijo Y, Hata S, Kyozuka J, Shimamoto K, Izui K. Over-expression of a single Ca2+-dependent protein kinase confers both cold and salt/drought tolerance on rice plants. Plant J. 2000;23:319–27.
Zou JJ, Wei FJ, Wang C, Wu JJ, Ratnasekera D, Liu WX, Wu WH. Arabidopsis calcium-dependent protein kinase CPK10 functions in abscisic acid- and Ca2+-mediated stomatal regulation in response to drought stress. Plant Physiol. 2010;154:1232–43.
Ma SY, Wu WH. AtCPK23 functions in Arabidopsis responses to drought and salt stresses. Plant Mol Biol. 2007;65:511–8.
Singh A, Sagar S, Biswas DK. Calcium dependent protein kinase, a versatile player in plant stress management and development. Crit Rev Plant Sci. 2017;36:336–52.
Shewry PR, Hey SJ. The contribution of wheat to human diet and health. Food Energy Secur. 2015;4:178–202.
Xu MJ, Gao T, Wang PF, Li GZ, Kang GZ. Function of calcium dependent protein kinase 34 in grain starch synthesis in common wheat (Triticum aestivum L.). J Agr Sci Technol Chin. 2019;21:26–33 (in Chinese with English abstract).
Choi HW, Hwang BK. The pepper extracellular peroxidase CaPO2 is required for salt, drought and oxidative stress tolerance as well as resistance to fungal pathogens. Planta. 2012;235:1369–82.
Xu J, Xing XJ, Tian YS, Peng RH, Xue Y, Zhao W, Yao QH. Transgenic Arabidopsis plants expressing tomato glutathione S-transferase showed enhanced resistance to salt and drought stress. PLoS One. 2015;10:e0136960.
Chen F, Fasoli M, Tornielli GB, Dal Santo S, Pezzotti M, Zhang L, Cai B, Cheng ZM. The evolutionary history and diverse physiological roles of the grapevine calcium-dependent protein kinase gene family. PLoS One. 2013;8:e80818.
Hamel LP, Sheen J, Séguin A. Ancient signals: comparative genomics of green plant CDPKs. Trends Plant Sci. 2014;19:79–89.
Liu W, Li W, He Q, Daud MK, Chen J, Zhu S. Genome-wide survey and expression analysis of calcium-dependent protein kinase in Gossypium raimondii. PLoS One. 2014;9:e98189.
Feldman M, Levy AA, Fahima T, Korol A. Genomic asymmetry in allopolyploid plants: wheat as a model. J Exp Bot. 2012;63:5045–59.
Li AL, Zhu YF, Tan XM, Wang X, Wei B, Guo HZ, Zhang ZL, Chen XB, Zhao GY, Kong XY. Evolutionary and functional study of the CDPK gene family in wheat (Triticum aestivum L.). Plant Mol Biol. 2008;66:429–43.
Jia J, Zhao S, Kong X, Li Y, Zhao G, He W, Appels R, Pfeifer M, Tao Y, Zhang X, et al. Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation. Nature. 2013;496:91–5.
Ling HQ, Zhao S, Liu D, Wang J, Sun H, Zhang C, Fan H, Li D, Dong L, Tao Y, et al. Draft genome of the wheat A-genome progenitor Triticum urartu. Nature. 2013;496:87–90.
International Wheat Genome Sequencing Consortium (IWGSC). A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science. 2014;345:1251788.
International Wheat Genome Sequencing Consortium (IWGSC); IWGSC RefSeq principal investigators, Appels R, Eversole K, Feuillet C, Keller B, Rogers J, Stein N, Pozniak CJ, et al. IWGSC whole-genome assembly principal investigators. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018;361:661.
Heuer S, Gaxiola R, Schilling R, Herrera-Estrella L, López-Arredondo D, Wissuwa M, Delhaize E, Rouached H. Improving phosphorus use efficiency: a complex trait with emerging opportunities. Plant J. 2017;90:868–85.
Li A, Wang X, Leseberg CH, Jia J, Mao L. Biotic and abiotic stress responses through calcium-dependent protein kinase (CDPK) signaling in wheat (Triticum aestivum L.). Plant Signal Behav. 2008;3:654–6.
Boudsocq M, Droillard MJ, Regad L, Laurière C. Characterization of Arabidopsis calcium-dependent protein kinases: activated or not by calcium? Biochem J. 2012;447:291–9.
Franz S, Ehlert B, Liese A, Kurth J, Cazalé AC, Romeis T. Calcium-dependent protein kinase CPK21 functions in abiotic stress response in Arabidopsis thaliana. Mol Plant. 2011;4:83–96.
Ho SL, Huang LF, Lu CA, He SL, Wang CC, Yu SP, Chen J, Yu SM. Sugar starvation-and GA-inducible calcium-dependent protein kinase 1 feedback regulates GA biosynthesis and activates a 14-3-3 protein to confer drought tolerance in rice seedlings. Plant Mol Biol. 2013;81:347–61.
Wei S, Hu W, Deng X, Zhang Y, Liu X, Zhao X, Luo Q, Jin Z, Li Y, Zhou S. A rice calcium-dependent protein kinase OsCPK9 positively regulates drought stress tolerance and spikelet fertility. BMC Plant Biol. 2014;14:133.
Bundó M, Coca M. Calcium-dependent protein kinase OsCPK10 mediates both drought tolerance and blast disease resistance in rice plants. J Exp Bot. 2017;68:erx145.
Campo S, Baldrich P, Messeguer J, Lalanne E, Coca M, Segundo BS. Overexpression of a calcium-dependent protein kinase confers salt and drought tolerance in rice by preventing membrane lipid peroxidation. Plant Physiol. 2014;165:688–704.
Brenchley R, Spannagl M, Pfeifer M, Barker GL, D'Amore R, Allen AM, McKenzie N, Kramer M, Kerhornou A, Bolser D, Kay S, et al. Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature. 2012;491:705–10.
Wang LL, Yu CC, Chen C, He CL, Zhu YG, Huang WC. Identification of rice Di19 family reveals OsDi19-4 involved in drought resistance. Plant Cell Rep. 2014;33:2047–62.
Kumar M, Lee SC, Kim JY, Kim SJ, Aye SS, Kim SR. Over-expression of dehydrin gene, OsDhn1, improves drought and salt stress tolerance through scavenging of reactive oxygen species in rice (Oryza sativa L.). J. Plant Biol. 2014;57:383–93.
Olvera-Carrillo Y, Campos F, Reyes JL, Garciarrubio A, Covarrubias AA. Functional analysis of the group 4 late embryogenesis abundant proteins reveals their relevance in the adaptive response during water deficit in Arabidopsis. Plant Physiol. 2010;154:373–90.
Li YG, Zhang DY, Li WN, Mallano AI, Zhang YH, Wang T, Lu M, Qin ZW, Li WB. Expression study of soybean germin-like gene family reveals a role of GLP7 gene in various abiotic stress tolerances. Can J Plant Sci. 2016;96:296–304.
Shi H, Chen Y, Qian Y, Chan Z. Low temperature-induced 30 (LTI30) positively regulates drought stress resistance in Arabidopsis: effect on abscisic acid sensitivity and hydrogen peroxide accumulation. Front Plant Sci. 2015;6:893.
Cieśsla A, Mituła F, Misztal L, Fedorowicz-Strońska O, Janicka S, Tajdel-Zielińska M, Marczak M, Janicki M, Ludwików A, Sadowski J. A role for barley calcium-dependent protein kinase CPK2a in the response to drought. Front Plant Sci. 2016;7:1550.
Vivek PJ, Tuteja N, Soniya EV. CDPK1 from ginger promotes salinity and drought stress tolerance without yield penalty by improving growth and photosynthesis in Nicotiana tabacum. PLoS One. 2013;8:e76392.
Allagulova CR, Gimalov FR, Shakirova FM, Vakhitov VA. The plant dehydrins: structure and putative functions. Biochemistry-Moscow. 2003;68:945–51.
Bies-Ethève N, Gaubier-Comella P, Debures A, Lasserre E, Jobet E, Raynal M, Cooke R, Delseny M. Inventory, evolution and expression profiling diversity of the LEA (late embryogenesis abundant) protein gene family in Arabidopsis thaliana. Plant Mol Biol. 2008;67:107–24.
Washio K, Ishikawa K. Organ-specific and hormone-dependent expression of genes for serine carboxypeptidases during development and following germination of rice grains. Plant Physiol. 1994;105:1275–80.
Rapala-Kozik M, Wolak N, Kujda M, Banas AK. The upregulation of thiamine (vitamin B1) biosynthesis in Arabidopsis thaliana seedlings under salt and osmotic stress conditions is mediated by abscisic acid at the early stages of this stress response. BMC Plant Biol. 2012;12:2.
Yubero-Serrano EM, Moyano E, Medina-Escobar N, Muñoz-Blanco J, Caballero JL. Identification of a strawberry gene encoding a non-specific lipid transfer protein that responds to ABA, wounding and cold stress. J Exp Bot. 2003;54:1865–77.
Hoagland DR, Arnon DI. The water-culture method for growing plants without soil. Calif Agric Exp Sta Circ. 1950;347:357–9.
Wei L, Wang L, Yang Y, Liu G, Wu Y, Guo T, Kang G. Abscisic acid increases leaf starch content of polyethylene glycol-treated wheat seedlings by temporally increasing transcripts of genes encoding starch synthesis enzymes. Acta Physiol Plant. 2015;37:206.
Li G, Wu Y, Liu G, Xiao X, Wang P, Gao T, Xu M, Han Q, Wang Y, Guo T, Kang G. Large-scale proteomics combined with transgenic experiments demonstrates an important role of jasmonic acid in potassium deficiency response in wheat and rice. Mol Cell Proteomics. 2017;16:1889–905.
Tufan HA, Stefanato FL, McGrann GRD, MacCormack R, Boyd LA. The barley stripe mosaic virus system used for virus-induced gene silencing in cereals differentially affects susceptibility to fungal pathogens in wheat. J Plant Physiol. 2011;168:990–4.
Ma HZ, Liu GQ, Li CW, Kang GZ, Guo TC. Identification of the TaBTF3 gene in wheat (Triticum aestivum L.) and the effect of its silencing on wheat chloroplast, mitochondria and mesophyll cell development. Biochem Biophys Res Commun. 2012;426:608–14.
Kang G, Li G, Xu W, Peng X, Han Q, Zhu Y, Guo T. Proteomics reveals the effects of salicylic acid on growth and tolerance to subsequent drought stress in wheat. J Proteome Res. 2012;11:6066–79.
Güitman MR, Arnozis PA, Barneix AJ. Effect of source-sink relations and nitrogen nutrition on senescence and N remobilization in the flag leaf of wheat. Physiol Plant. 1991;82:278–84.
Demiral T, Türkan İ. Comparative lipid peroxidation antioxidant defense systems and proline content in roots of two rice cultivars differing in salt tolerance. Environ Exp Bot. 2005;53:247–57.
Su L, Shan JX, Gao JP, Lin HX. OsHAL3, a blue light-responsive protein, interacts with the floral regulator Hd1 to activate flowering in rice. Mol Plant. 2016;9:233–44.
We thank Dr. Changqing Du for helping to determine the characteristics of TaCPK34. We also thank Shanghai Hoogen Biotechnology Co., Ltd. for assisting in sequencing and bioinformatics analysis.
This work was financially supported by the National Natural Science Foundation of China (U1704110 and 31701345). And the funding bodies were not involved in the design of the study, in the collection, analysis and interpretation of data or in writing the manuscript.
Ethics approval and consent to participate
No ethics approval was required for the collection or use of animals in this study.
Consent for publication
The authors declare they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Amino acid sequence alignment between TuCPK34 (EMS58754) in Triticum urartu and TaCPK34 in Chinese Spring (CS-TaCPK34). Supplemental Figure S2. Abundance of TaCPK34 protein in wheat grains after anthesis. Supplemental Figure S3. The location of silenced gene-specific fragment in full-length of the TaCPK34 gene. Supplemental Figure S4. BSMV-VIGS-mediated silencing of the TaCPK34 gene. Supplemental Figure S5. The promoter sequences of the genes encoding the identified proteins with known functions.
Meteorological parameters (rainfall, temperature, etc.) after anthesis in the experimental field.
The selected CPKs from Arabidopsis thaliana, rice and bread wheat for the phylogenetic tree analysis.
All peptides of the identified proteins in leaves of BSMV-VIGS-TaCPK34- or BSMV-VIGS-GFP-inoculated wheat seedlings suffering from drought for 14 days.
All identified proteins in leaves of BSMV-VIGS-TaCPK34- or BSMV-VIGS-GFP-inoculated wheat seedlings suffering from drought for 14 days in all four biological replicates.
Primer sequences used in this study.
About this article
Cite this article
Li, GZ., Li, HX., Xu, MJ. et al. Functional characterization and regulatory mechanism of wheat CPK34 kinase in response to drought stress. BMC Genomics 21, 577 (2020). https://doi.org/10.1186/s12864-020-06985-1
- Calcium-dependent protein kinase
- Triticum aestivum L.