Deepika Deepika
Nikita Poddar
Shailesh Kumar
Amarjeet Singh- 1Stress Signaling Lab, National Institute of Plant Genome Research, New Delhi, India
- 2Bioinformatics Lab, National Institute of Plant Genome Research, New Delhi, India
Calcium-dependent protein kinases (CDPKs) are a major group of calcium (Ca2+) sensors in plants. CDPKs play a dual function of “Ca2+ sensor and responder.” These sensors decode the “Ca2+ signatures” generated in response to adverse growth conditions such as drought, salinity, and cold and developmental processes. However, knowledge of the CDPK family in the legume crop chickpea is missing. Here, we have identified a total of 22 CDPK genes in the chickpea genome. The phylogenetic analysis of the chickpea CDPK family with other plants revealed their evolutionary conservation. Protein homology modeling described the three-dimensional structure of chickpea CDPKs. Defined arrangements of α-helix, β-strands, and transmembrane-helix represent important structures like kinase domain, inhibitory junction domain, N and C-lobes of EF-hand motifs. Subcellular localization analysis revealed that CaCDPK proteins are localized mainly at the cytoplasm and in the nucleus. Most of the CaCDPK promoters had abiotic stress and development-related cis-regulatory elements, suggesting the functional role of CaCDPKs in abiotic stress and development-related signaling. RNA sequencing (RNA-seq) expression analysis indicated the role of the CaCDPK family in various developmental stages, including vegetative, reproductive development, senescence stages, and during seed stages of early embryogenesis, late embryogenesis, mid and late seed maturity. The real-time quantitative PCR (qRT-PCR) analysis revealed that several CaCDPK genes are specifically as well as commonly induced by drought, salt, and Abscisic acid (ABA). Overall, these findings indicate that the CDPK family is probably involved in abiotic stress responses and development in chickpeas. This study provides crucial information on the CDPK family that will be utilized in generating abiotic stress-tolerant and high-yielding chickpea varieties.
Introduction
Chickpea, the world's second most important food legume, is majorly produced by South Asia. Importantly, as the largest producer of chickpea, India contributes about 70% (5.9 million tons (mt) annually) to the world's chickpeas production (Varshney et al., 2009). Chickpea seeds are of high nutrition value as they contain 20–30% crude protein, 40% carbohydrate, 3–6% oil, and an abundance of minerals, such as calcium, magnesium, potassium, phosphorus, iron, and zinc (Gil et al., 1996; Ibrikci et al., 2003). Unfortunately, due to continuously fluctuating and harsh environmental conditions in the semi-arid tropics where chickpea is majorly cultivated, its productivity is severely affected. An estimated 40–60% of annual global chickpea yield losses are accounted for by abiotic stress factors only. While drought inflicts severe damage and accounts for almost 50% of chickpea yield loss, temperature changes and soil salinity together account for about 25% of chickpea yield loss (Varshney et al., 2014). Consequently, a wide gap is developed between chickpea demand and supply. Importantly, damage to chickpea yield dents several chickpea-producing countries economically. Drought, cold, and salinity, respectively, are reported to cost ~1.3 billion, 186, and 354 million US dollars (Ryan, 1997). Abiotic stresses negatively impact the flower set, pollen viability, pod set/abortion, and retention. As these developmental stages determine the seed number in chickpea, a negative impact on them significantly hampers chickpea yield. Thus, identification of crucial stress-related genes and their utilization in breeding programs to generate stress-tolerant chickpea varieties are urgently required.
Abiotic stresses are known to elicit a profound increase in cytosolic calcium (Ca2+) levels. The spatio-temporal accumulation of Ca2+ develops a specific “Ca2+ signature” in the form of Ca2+ spikes, waves, and oscillations. The Ca2+ signature encoded by a specific stimulus is decoded by various Ca2+ sensors, toward a specific response (Tang et al., 2020). Major Ca2+ sensors that have been identified and characterized in plants include calmodulin (CaM) and CaM-like proteins (CMLs) (Zeng et al., 2015), calcineurin B-like proteins (CBLs) (Luan, 2009), and Ca2+-dependent protein kinases (CDPKs) (Singh et al., 2017). Among these, CDPKs are one of the largest Ca2+ sensor groups in the plant kingdom and constitute a multi-gene family (Schulz et al., 2013; Xiao et al., 2017). A total of 34 CDPK have been found in the Arabidopsis thaliana genome (Cheng et al., 2002), 31 genes in the rice (Oryza sativa) genome (Ray et al., 2007), and 20 genes in the wheat (Triticum aestivum L.) genome (Li et al., 2008). CDPKs are unique enzymes with a dual function of a Ca2+ sensor and responder, attributed to their peculiar structural features. Differing from other Ca2+ sensors, CDPKs can sense, respond to, and translate Ca2+ signals into protein phosphorylation events (Poovaiah et al., 2013). Plant CDPKs are comprised of a conserved domain structure and are monomeric in nature. A variable N-terminal domain is followed by a ser/thr kinase domain and a CDPK activation domain (CAD). An auto-inhibitory region called as “inhibitory junction domain” and a calcium-binding CaM-like domain (CaM-LD) is located within CAD (Schulz et al., 2013). CaM-LD contains four elongation factors (EF) hand motifs that are responsible for Ca2+ binding. These EF hands are organized as N-terminal and C-terminal lobes (each consists of two EF hands) (Boudsocq and Sheen, 2013; Liese and Romeis, 2013). At the basal state, the C-terminal lobe shows high affinity to Ca2+ thus, it remains loaded with Ca2+ even at low Ca2+ concentration. The C-terminal lobe via interaction with auto-inhibitory junction maintains the kinase in an inactive state. The binding of Ca2+ induces a conformational change in the N-terminal lobe, which disrupts auto-inhibitory junction-kinase interaction. That removes auto-inhibition and leads to the activation of kinase (Boudsocq and Sheen, 2013). In plants, different CDPK isoforms have been found to exhibit distinct expression patterns, which possibly accounts for their functional specificity (Yang et al., 2017; Zhang et al., 2017). CDPK proteins localize in the cytosol and subcellular organelles, including the nucleus, plasma membrane, endoplasmic reticulum, tonoplast, mitochondria, and chloroplasts (Simeunovic et al., 2016). This suggests that CDPKs might target variable substrates throughout the plant cell. In plants, CDPKs have been involved in regulating important functions, including biotic and abiotic responses, hormone signaling, and development (Schulz et al., 2013). Knowledge of plants' CDPK functions has been generated majorly from the research with the model plant Arabidopsis thaliana. Information about molecular features of the CDPK family and its role in important legume crop chickpea is missing. Molecular characterization of the CDPK family will help in understanding their functions in chickpeas.
Here, we have unearthed the entire repertoire of CDPK encoding genes in the chickpea genome. Gene and domain structure analysis confirmed the authenticity and integrity of CDPKs. Phylogenetic analysis and chromosomal localization provided crucial insight into the evolution and expansion of the chickpea CDPK family. Homology modeling was used to understand the three-dimensional structure of chickpea CDPK proteins. Subcellular localization analysis showed that chickpea CDPK proteins are mainly nuclear and/or cytoplasmic. In-silico promoter analysis showed the presence of stress, hormone, and development-related cis-regulatory elements in CDPK promoters. Extensive expression analysis of the chickpea CDPK family was performed under abiotic stresses (drought, salinity, and cold) and during different developmental stages using public RNA-Seq data and qRT-PCR analysis. Expression analysis indicated the involvement of the CDPK family in abiotic stress signaling and plant development in chickpea.
Materials and Methods
Identification of CDPKs in the Chickpea Genome
The chickpea genome at National Centre for Biotechnology Information (NCBI) (Varshney et al., 2013) was explored to identify the CDPK encoding genes. CDPK protein sequences of rice and Arabidopsis thaliana were retrieved from Uniprot (Swiss-Prot), and used for BLAST homology search in the chickpea genome database. Further, the Hidden Markov Model (HMM) profile of CDPK was extracted from the Pfam (http://pfam.xfam.org/) database and was used as a query to search the chickpea database at NCBI. All obtained sequences were compiled, redundant entries were removed, and only unique entries were used in further analysis. Protein sequences of putative CDPKs were scrutinized for the presence of canonical domains using in-silico tools, such as SMART (http://smart.embl-heidelberg.de/), Interpro (https://www.ebi.ac.uk/interpro/), Prosite (https://prosite.expasy.org/) and Pfam (https://pfam.xfam.org/). Various attributes of CDPKs, such as gene ID, protein ID, CDS size, protein size, Introns, molecular weight (MW), isoelectric point (pI), and chromosomal coordinates were extracted from NCBI and ExPASy (https://web.expasy.org/compute_pi/).
Phylogenetic Analysis
Multiple Sequence Alignment (MSA) was performed using non-redundant protein sequences of CDPKs from chickpea, Arabidopsis, rice, and soybean (Glycine max), using Clustal W at default settings. A phylogenetic tree was generated in MEGA X version 10.1.8 (Pennsylvania State University, USA) by the neighbor-joining method. Bootstrap values were calculated in 1,000 replicates to determine the phylogenetic relationship among the CDPKs. The web-server iTOL (Letunic and Bork, 2021) was used to mark the different clades of CIPKs with different colors for better visualization.
Gene Structures and Domain Prediction
To investigate the gene structure of CDPKs, their CDS and the genomic sequences were extracted from NCBI. These sequences were submitted at the Gene Structure Display Server 2.0. (http://gsds.cbi.pku.edu.cn/index.php?input=ite) to generate the gene structure diagram. Identification of domains was carried out using a standalone package of InterPro Scan. Co-ordinates of essential domains and active sites were extracted and used as input in Illustrator for Biological Sequences for the visualization.
Gene Nomenclature, Chromosomal Localization, and Gene Duplication
The nomenclature of chickpea CDPK genes was done based on sequence closeness to their Arabidopsis orthologs, and phylogenetic analysis. Genes were named as CaCDPK followed by a number (1–22) corresponding to their respective Arabidopsis orthologs. The information of chromosome co-ordinates of genes was obtained from NCBI and further used to display chromosomal localization. The MCScanX software package (Wang et al., 2012; Athens, USA) was used to assess the gene duplication within the chickpea CDPK gene family. Genes located within 20 kb distance on the same chromosome were considered as tandemly duplicated genes (Feng et al., 2015).
Homology Modeling of CaCDPK Proteins
The three-dimensional (3D) structures of all CaCDPK proteins were predicted by homology modeling using the PHYRE2 web portal (http://www.sbg.bio.ic.ac.uk/phyre2). PHYRE2 uses advanced remote homology detection methods to build 3D models for protein sequences (Kelley et al., 2015). All the proteins were modeled with 100% confidence by the single highest scoring template model.
In-Silico Subcellular Localization of CaCDPK Proteins
The full-length protein sequences of CaCDPKs were used as input in the subCELlular LOcalization predictor: CELLO online tool (Yu et al., 2006) to predict their subcellular localization.
Constructs Preparation for In-Planta Subcellular Localization
The protein coding sequence (ORF excluding their stop codon) of CaCDPK5, CaCDPK16, and CaCDPK21 genes were amplified from chickpea complementary DNA (cDNA) with gene-specific primers using iProof high fidelity DNA polymerase (Bio-Rad) through PCR in a thermocycler (Applied Biosystems). The list of these primers is given in Supplementary Table 1. The ORFs after amplification cloned into gateway entry vector p-ENTR-D-TOPO (Invitrogen). Genes were subsequently mobilizing into a compatible destination vector pSITE3CA under the control of 2XCaMV35S promoter by LR recombination protocol. The authenticity of all the constructs was ensured by PCR and sequencing.
Agro-Infiltration Into Nicotiana benthamiana and Confocal Microscopy
Agrobacterium tumefaciens (GV3101::pMP90) cells were transformed with the YFP constructs of the respective genes. Transformed Agrobacterium cells were used to transfect 6-week-old N. benthamiana plant leaves. Plants were grown in a growth chamber with the following conditions: 12/12 h photoperiod, 25–26°C temperature, and 60% relative humidity for 48–72 h. Transiently transformed Nicotiana leaf discs were analyzed under Total confocal scanner (TCS) SP5 laser scanning electron microscope (Leica, Germany) to detect the florescence, according to Deepika et al. (2022).
In-Silico Promoter Analysis
The 2 kb sequence, upstream of translational start site of CDPK genes was retrieved from NCBI. This sequence was used as input in PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) for the detection of various cis-regulatory elements and motifs. Important cis-regulatory related to abiotic stresses, hormonal response, and plant development were selected.
RNA-Seq Expression Analysis in Developmental Stages and Tissues
To investigate the expression pattern of CaCDPK genes during developmental stages, RNA-Seq data was extracted from the NCBI-Sequence Read Archive (SRA) (ID:SRP121085). The raw reads obtained from SRA were processed using the FASTP tool (Shenzhen, China). An index of the reference genome was built and mapping of raw reads onto the reference genome was done using HISAT2 (Kim et al., 2019). StringTie (Pertea et al., 2015) was used to assemble the aligned sequences into potential transcripts. Transcript abundance was calculated by fragments per kilobase of transcript per million reads (FPKM) values. Expression dynamics were analyzed in 27 tissues representing different developmental stages, such as germination (radicle, plumule, and embryo), seedling (Epicotyl and primary root), vegetative (root, petiole, stem, and leaf), reproductive (Petiole, stem, nodules, root, flowers, buds, pods, immature seeds, and leaf), and senescence (immature seeds, mature seeds, seed coat, stem, petiole, root, nodules, leaf, and yellow leaf). Log2 transformed expression values were used to generate the heat-map using the MeV4 tool (Maryland, USA). RNA-Seq data for different seed developmental stages (S1–S7) in two distinct desi chickpea varieties (JGK3 and Himchana 1) was extracted from SRA number SRP072563 and SRP072564.
Plant Growth and Stress Treatment
Desi chickpea (var. ICC4958) was used for gene expression analysis. The seeds were surface sterilized and plants were grown according to Sagar et al. (2020). Ten-day-old seedlings were subjected to different stress treatments. For drought stress, water was withdrawn and seedlings were air-dried within the folds of tissue paper at 22–23°C temperature. Samples were harvested in replicates after 0 (untreated control), 1, 3, and 6 h of drought treatment. For salt stress, seedlings were kept in 150 mM sodium chloride (NaCl) solution in a beaker, and samples were collected after 0, 3, 6, and 12 h. Abscisic Acid (ABA) seedlings were kept in 100 μM (±) ABA in sterile water in a beaker under light and samples were collected after 0, 3, 6, and 12 h treatment. Same aged seedlings were kept in sterile water for control at 22–23°C.
RNA Extraction and cDNA Synthesis
A total of 100 mg tissue of the control, drought, salinity, and ABA treated root and shoot samples were used for RNA extraction using the TRIzol reagent (Ambion, Life technologies, USA) according to the manufacturer's protocol. The RNA obtained was purified to remove any genomic DNA contamination using an RNeasy Min Elute Clean-up Kit (QIAGEN, Hilden, Germany). The quantity and quality of RNA were ascertained by the ratio 1.8–2.0 for A260:A280 and 2.0–2.3 for A260:A230 using a Nano Drop ONEc (Thermo Scientific, USA) nano-spectrophotometer. Subsequently, MOPS-agarose gel electrophoresis was done to confirm the integrity of the RNA. A total of 1 μg total RNA was used to synthesize the first strand cDNA using a RevertAid first-strand cDNA synthesis kit (Thermo Scientific) according to the manufacturer's protocol.
Expression Analysis by qRT-PCR
The qRT-PCR primers for selected genes were synthesized by the PRIMER EXPRESS SOFTWARE (Applied Biosystems, USA) according to Singh and Pandey (2015). Their specificity was analyzed using the RGAP BLAST tool and melt curve analysis after a real-time PCR run. The details of all the primers are given in Supplementary Table 1. Three biological replicate samples (with three technical replicates of each biological replicate) of control and nutrient-deficient root and shoots were used to assess the expression pattern. iTaq Universal SYBR Green supermix (Bio-Rad) was used to detect the expression in Bio-rad CFX96 real-time PCR system (Bio-Rad) according to Sagar et al. (2021).
Statistical Analysis
For statistical significance, all expression and quantitative experiments were replicated three times. The data have been presented as the mean of three replicates ± SD. A two-tailed student's t-test was performed to determine the statistical significance among the replicate samples. A p < 0.05 was considered statistically significant (denoted by *), p < 0.01 (denoted by **) and p < 0.005 (denoted by ***).
Results and Discussion
Identification and Organization of CDPK Family in the Chickpea Genome
A thorough investigation of various databases resulted in the identification of 22 non-redundant CDPK encoding genes in the chickpea genome. This is consistent with previous findings, where different plant species including Arabidopsis thaliana, rice (Oryza sativa), wheat (Triticum aestivum L.), and tomato (Solanum lycopersicum) have been found to encode for about 20–30 CDPK genes (Cheng et al., 2002; Ray et al., 2007; Li et al., 2008; Hu et al., 2016; Wang et al., 2016). Gene structure analysis showed that CaCDPKs are made up of multiple exons and introns. Most CDPK genes are comprised of 6-8 introns except CaCDPK16, which is comprised of 11 introns. Domain analysis revealed that all the CaCDPKs harbored a canonical catalytic kinase domain toward the N-terminus with a typical Ser/Thr kinase active site and ATP binding site (Figure 1). In addition, four EF-hand motifs were present in all the CDPK proteins toward the C-terminus, which are crucial for Ca2+ binding. Strikingly, a crucial Ca2+ binding site was found to be missing in one of the EF hand motifs of CaCDPK6, CaCDPK7, CaCDPK8, CaCDPK13, and CaCDPK22. The absence of such an important site might impair the Ca2+ binding ability of these CDPKs. However, variations in the number of EF-hand motifs and Ca2+ binding sites have also been previously reported in other plants (Cheng et al., 2002; Asano et al., 2005; Kong et al., 2013; Zuo et al., 2013). The in-depth analysis of N-terminal sequences of CaCDPKs revealed that out of 22 proteins, 11 have a myristoylation site (Table 1). The N-terminal domain of many CDPKs in several plant species has been found to contain potential N-myristoylation and N-palmitoylation sites. For example, AtCPK2, AtCPK3, AtCPK6, AtCPK9, AtCPK13, AtCPK5, and AtCPK16 (Benetka et al., 2008; Mehlmer et al., 2010; Lu and Hrabak, 2013), N. tabacum NtCDPK2/NtCDPK3 (Witte et al., 2010) and potato (Solanum tuberosum) StCDPK4/StCDPK5 (Asai et al., 2013) have been marked with N-myristoylation. These N-terminal modifications determine the membrane targeting of plant CDPKs, as mutations in N-myristoylation or N-palmitoylation sites have been found to hamper their membrane targeting (Singh et al., 2017). CDPKs show a high degree of structural conservation in terms of gene and protein structure, as similar exon-intron and domain arrangement pattern has been observed in diverse plant species including tomato, barley (Hordeum vulgare L.), Brachypodium distachyon, and Medicago truncatula (Hu et al., 2016; Yang et al., 2017; Wen et al., 2020; Zhao et al., 2021). Chickpea CDPK proteins length was found to be in the range of 497–764 amino acids, and molecular weight varied between 56.16 and 68.74 KDa. Interestingly, CDPK protein had a highly variable isoelectric point (pI) falling in the range of 5.09 (CaCDPK22) and 8.99 (CaCDPK16) (Table 1). This indicates that the chickpea CDPK protein might function optimally in a diverse microenvironment.
Figure 1. Structural features of the chickpea calcium-dependent protein kinases (CDPK) family. (A) Exon-intron organization is shown for CaCDPK genes and gene names are mentioned at the left. The scale at the bottom represents gene length in kb. (B) Protein domain structure of the chickpea CDPK family is showing conserved protein kinase and EF-hand along with some important sites present. ATP binding site and Serine/threonine-protein kinase active site are located at the protein kinase domain.
Table 1. Summary of various features of the chickpea CDPK gene family.
Evolutionary Analysis of CDPK Family
To understand the evolution of chickpea CDPKs, phylogenetic analysis was performed with Arabidopsis, rice, soybean, and chickpea CDPK proteins. All the CDPKs from different plant species could be demarcated into four sub-clades; group I-IV(Figure 2). While, groups I and II each contained eight chickpea CDPK members, group III contained five CDPKs (CaCDPK7, 8, 10, 13, and 14), and group IV contained only a single CDPK member i.e., CaCDPK16. This phylogenetic distribution is conserved across different plant species, as group I comprises the highest number of CDPK members and group IV comprises the least members (Hu et al., 2016). This analysis suggests the evolutionary conservation of CDPKs across plant species. However, within the separate clades, chickpea CDPKs were closer to the dicot plants Arabidopsis and soybean, whereas distantly placed from monocot plant rice. Such distribution indicates the evolutionary divergence of monocot and dicot CDPKs.
Figure 2. Phylogenetic relationship among CDPKs from different plants. An un-rooted neighbor-joining phylogenetic tree was constructed from the CDPK protein sequences of Arabidopsis thaliana, rice, soybean, and chickpea. All the CDPKs could be divided into four groups; I- IV. Numbers above the nodes represent bootstrap values calculated from 1,000 replicates.
Chickpea CDPK genes were mapped on the different chromosomes using their chromosomal coordinates. Out of 22 CDPK genes, 19 could be mapped on seven out of eight chickpea chromosomes. Three genes (CaCDPK7, 13, and 17) were placed on the scaffold region and none of the genes was localized on chromosome 1. CaCDPK genes were variably distributed on chromosomes, with chromosome 4 containing the highest five genes. Four and three CDPK genes were located on chromosomes 2 and 6, respectively. Two genes each were located on chromosomes 5, 7, and 8. Whereas, only a single gene, CaCDPK8 was located on chromosome 3 (Figure 3). Gene duplication analysis revealed that three pairs of CDPK genes (CaCDPK3/6, CaCDPK9/19, and CaCDPK11/12) were segmentally duplicated. In addition, one gene pair, CDPK2/20 was found to be tandemly duplicated. Gene duplication is considered an important process for the evolution of gene families in plants (Singh et al., 2010, 2014; Sagar et al., 2021). Similar to chickpea, the CDPK family has been evolved through segmental and tandem duplication in diverse plant species, such as Arabidopsis, rice, cotton, poplar, moss (P. patens), and cabbage (Brassica rapa) (Hrabak et al., 2003; Ray et al., 2007; Zuo et al., 2013; Hamel et al., 2014; Liu et al., 2014; Wu et al., 2017). However, the extent of duplication has been found to be varied in different species. For instance, chickpea had only three segmentally duplicated pairs, whereas nine and eight gene pairs were segmentally duplicated in the rice and Arabidopsis genomes, respectively (Ray et al., 2007). This indicates that gene duplication is a conserved mechanism of the evolution and expansion of the CDPK family in plants.
Figure 3. Chromosomal locations of chickpea CDPK genes. Green color bars represent the chromosomes, the location of genes has been marked alongside. Chromosome numbers are given at the bottom of each chromosome. Except chromosome 1, the CDPK genes are distributed variably on seven chromosomes.
3D Protein Structure Analysis
To get an insight into the 3-D protein structure of chickpea CDPKs, homology-modeling was performed for all the CaCDPKs. The protein structures were obtained by comparing with related template protein pbanka-031420 (PDB id 3Q5I). All protein structures were modeled with the template protein with 100% confidence. Most proteins were modeled with a high coverage level ranging from 72 to 90%. The 3-D structures of different CaCDPKs were comprised of a variable number of α-helix, β-strands, disordered region, and transmembrane (TM) - helix (Supplementary Table 2). Among these, the most frequently occurring secondary structure was the α-helix, which warrants the stability of protein structure (Neelamathi et al., 2009). The percentage of α- helices varied from 36 to 48, whereas the β-strands contributed to 11–19% of the protein structure (Supplementary Table 2). In all the CaCDPK proteins, structural folds made up of blue and green α- helices and β-strands represent the catalytic kinase domain. Orange and red color ribbons represent the N-lobe and C-lobe, respectively with each lobe containing two EF hand motifs. Whereas, an inhibitory junction domain (JD) is represented by yellow helical ribbons (Figure 4). Thus, the presence of important characteristic domains and motifs in all the CaCDPK proteins confirms their authenticity and integrity. Similar features of the 3-D structure of CDPK proteins were obtained in different plant species including Arabidopsis, rice, maize, ginger (Zingiber officinale), and sorghum (Mittal et al., 2017; Vivek et al., 2017). Suggesting that the CDPK protein structure is highly conserved, with similar structural folds and conformational arrangements in higher plants. This also hints toward their conserved functional and structural mechanism.
Figure 4. Three-dimensional structure of chickpea CDPK proteins. A three-dimensional (3-D) structure was generated for all 22 members of the chickpea CDPK family. Each CDPK protein is made up of a variable number of α-helix, β-strands, transmembrane helix, and disordered region. Structure comprising blue and green α-helices and β-strands represents kinase domain, yellow helical ribbons represent inhibitory junction domain (JD), orange helices indicate N-lobe, and red helices indicate C-lobe with each containing two elongation factors (EF) hand motifs.
Subcellular Localization of CaCDPK Proteins
To identify the possible sites of activity of CaCDPK proteins in order to have an insight into their function, in-silico subcellular localization was performed. By considering the highest reliability score, 15 CaCDPK proteins were exclusively localized in the cytoplasm, whereas only two proteins, CaCDPK16 and 19 were exclusively localized in the nucleus (Figure 5A). Five CDPK proteins, CaCDPK2, 5, 9, 15, and 18 could be localized both in the cytoplasm and nucleus. To validate the in-silico localization pattern, fluorescence-based subcellular localization was performed in N. benthamiana for three randomly selected candidates; CaCDPK5, 16, and 21. Confocal microscopy analysis revealed that CaCDPK5 was localized in both the nucleus and cytoplasm and CaCDPK21 was localized at the cytoplasm (Figure 5B) which is consistent with in-silico prediction. Whereas, CaCDPK16, which was predicted to be nucleus localized was actually localized in the cytoplasm of Nicotiana cells. Thus, in-planta localization significantly supported and validated the in-silico localization pattern. These findings indicate that CaCDPKs might phosphorylate different substrates in the cytoplasm and in the nucleus to regulate various cellular functions in chickpea. CDPKs have been found to be localized at diverse subcellular locations in plants. Most of the CDPKs in Arabidopsis are either membrane-localized or are localized at both cytoplasm and membrane. Whereas, a few AtCDPKs are exclusively cytoplasm localized (Simeunovic et al., 2016). N-myristoylation and palmitoylation are crucial modifications of CDPK proteins that determine CDPK localization (Simeunovic et al., 2016; Zheng et al., 2019). Recently, a study on legume plant Medicago truncatula showed that the CDPKs with N-terminal acylation sites were localized at the plasma membrane whereas, those lacking the N-acylation sites were distributed in cytosol and nucleus (Zhao et al., 2021). Consistently, in our study CaCDPK2, 5, and 18 were devoid of N- myristoylation site and distributed in nucleus and cytosol. Surprisingly, several CaCDPK proteins were marked with N- myristoylation site but none of them was found to be localized at the plasma membrane. It is known that N-myristoylation is an irreversible acylation, it requires a second post-translational signal i.e., reversible palmitoylation to sustain membrane localization of a CDPK (Witte et al., 2010). Therefore, the absence of a palmitoylation site in CaCDPK proteins could be a possible reason for them not to be localized at the plasma membrane. Also, the reversibility of the second post-translational signal could lead to the ferrying of CDPK between membrane and cytosol or nucleus (Boudsocq and Sheen, 2013). Moreover, the localization of CaCDPK proteins could be coupled with a stimulus, and a particular condition may lead to a change in their subcellular location. For example, the subcellular location of Arabidopsis AtCPK10, AtCPK30, and AtCPK32 proteins was changed due to variable availability (Liu et al., 2017).
Figure 5. Subcellular localization of chickpea CDPK proteins. (A) Localization was predicted using an online tool, CELLO. The figure shows chickpea CDPK proteins localized majorly in the cytoplasm and nucleus based on the highest reliability score with black color. Proteins that were predicted to be localized both in the cytoplasm and nucleus are shown by blue color at both locations. (B) N. benthamiana cells expressing the YFP- fusion protein driven by the 2XCaMV 35S promoter. Confocal images of fluorescence expressing CaCDPK proteins are showing their variable distribution in different compartments. CaCDPK5 is localized to the cytoplasm and nucleus, whereas CaCDPK16 and CaCDPK21 are localized at the cytoplasm. Cells transformed with vector only (CaMV35S-YFP) are shown in the lowermost row. Scale bar = 20 μm.
cis-Elements in CDPK Promoters
The cis-regulatory elements in the promoter are crucial for the transcriptional regulation of a gene. Therefore, to get an insight into the transcriptional regulation of CaCDPK genes their promoters were investigated for cis-regulatory elements. A total of 16 types of cis-regulatory elements were found to be variably distributed in CaCDPK promoters (Figure 6, Supplementary Table 3). These elements majorly belong to three categories; stress-responsive (GT1-motif, LTR, WUN-motif, ARE, TC-rich repeats, Box-4, AE-box, MYB, Myb-like, MYC, chs-CMA1a, and STRE), hormone-responsive (ABRE, ERE), and plant development related (AAGAA-motif, GATA-motif) (Hughes et al., 2000; Abe et al., 2003; Yamaguchi-Shinozaki and Shinozaki, 2005; Sharma et al., 2011; Nakashima and Yamaguchi-Shinozaki, 2013; Deepika and Singh, 2021; Hou et al., 2021). Abiotic stress-related motif Box-4 (light-responsive) was present in all CaCDPK promoters, and MYB motif (drought, low temperature, salt stress-related) was found in all CaCDPK promoters except CaCDPK2. This suggests that by binding these motifs, specific TFs might regulate the expression of CaCDPKs under different abiotic stresses. Markedly, 14 out of 21 CaCDPK promoters which contain the MYB motif were also found to contain the ABRE motif (ABA-responsive element), suggesting that these genes could regulate abiotic stresses, such as drought and salinity via the ABA-dependent pathway. It has been well-understood that different developmental events and abiotic stresses are interconnected through ABA in plants (Schroeder et al., 2001; Singh et al., 2015). Especially, during the later stages of seed development, a programmed dehydration event is triggered that results in seed dormancy (Hetherington, 2001; Schroeder et al., 2001). ABA commonly regulates these developmental and abiotic stress responses (Nakashima et al., 2009; Hubbard et al., 2010), and genes, which are involved in the regulation of such responses have been found to contain ABRE elements in their promoters (Singh et al., 2012, 2013). Thus, CaCDPK promoters with elements like MYB and ABRE could be involved in regulating chickpea plant development under abiotic stress conditions. Apart from two promoters (CaCDPK4 and 16), all the CaCDPK promoters also contained one or more ethylene-responsive, ERE motifs. Ethylene is an important plant hormone that is involved in plant defense and development (Mishra et al., 2013). The presence of the ERE element in most of the CaCDPK promoters reinforces the involvement of the CDPK family in plant development, biotic and abiotic stress responses in chickpea. Consistently, CDPK family members have been found to regulate development during plant-pathogen interaction, herbivore attack, and wounding in diverse plant species (Cai et al., 2015; Zhang et al., 2015; Hettenhausen et al., 2016; Xiao et al., 2017).
Figure 6. Promoter analysis of chickpea CDPK genes. Different cis-regulatory elements in the 2 kb upstream promoter region of CaCDPK genes are illustrated by different colors in the bar chart. The X-axis represents the name of genes and Y-axis represents the number of different cis-elements in each promoter. The names of cis-regulatory elements are mentioned on the right.
Expression Pattern of CaCDPKs in Developmental Stages
To get an insight into the functional role of CDPK genes in plant development, expression analysis was undertaken during various stages of chickpea development. Expression profiles of all CaCDPK genes were generated during 27 development stages. These stages represent germination, seedling, vegetative stages, reproductive stages, and senescence. All genes except CaCDPK2 showed a differential expression pattern during multiple developmental stages (Figure 7, Supplementary Table 4) A total of 10 genes showed ubiquitous expression during most developmental stages these include, CaCDPK1, 5, 6, 8, 10, 11, 12, 13, 16, and 19. Out of these, CaCDPK8 had the strongest expression during most stages of reproductive development, and some stages of vegetative development and senescence. Notably, the expression of all these 10 CaCDPK genes was relatively lower in three stages of senescence namely, immature seeds, mature seeds, and seed coat. This suggests that these CaCDPK genes could be involved in a range of developmental processes from germination to senescence. Some of the genes, including CaCDPK7, 14, 17, and 21 showed a significant and specific expression in flowers, suggesting their role in flower development. This observation is consistent with previous findings where CDPKs were implicated in flower development. FLOWERING LOCUS T (FT) is an important component of the regulatory network of flower development and it forms a complex with its interdependent partner AtFD (also known as AtbZIP14). Arabidopsis CDPK members, AtCPK4/AtCPK6/AtCPK33 phosphorylate AtFD at T282 in-vitro in presence of Ca2+, which is crucial for the formation and function of this complex. This complex triggers the transcriptional activation of floral meristem identity genes which control floral transition (Abe et al., 2005; Kawamoto et al., 2015). Thus, genes like CaCDPK7, 14, 17, and 21 could be involved in similar regulatory networks controlling flower development in chickpea.
Figure 7. Expression profiles of CDPK genes in different developmental stages of chickpea. The heatmap represents the expression pattern of CaCDPK genes in different tissues of developmental stages, such as germination, seedling, vegetative, reproduction, and senescence. The genes are named on the left and different tissues/developmental stages are labeled at the top. Scale bar represents the normalized log_2 FPKM values.
In Arabidopsis, expression and functional analysis indicated the role of multiple CDPKs in pollen tube growth. Analysis of single and double mutants of CDPK17 and CDPK34 revealed that their redundant activity is required for pollen tube tip growth and pollen fitness (Myers et al., 2009). A shaker-type K+ inward channel mediates the K+ influx which is crucial for pollen tube growth in Arabidopsis. CPK11 and CPK24 work together to inhibit this channel, thereby suppressing pollen tube growth (Zhao et al., 2013). Similarly, CPK2 and CPK20 control pollen tube growth via regulating another ion channel SLAH3 (Gutermuth et al., 2013). Furthermore, CDPKs regulate plant development via modulating phytohormone signaling. AtCDPK28 regulates the expression of key GA biosynthesis genes (e.g., gene coding for GA3ox1) and thus, controls various facets of plant development, including petiole and stem elongation (Matschi et al., 2013). Similarly, in Nicotiana attenuata, NaCDPK4 and 5 regulating GA activity positively regulate stem elongation and plant height (Heinrich et al., 2013). In our study, CaCDPK15 was specifically and strongly down-regulated during the germination stage of plumule, vegetative stage of petiole and stem, and senescence stage of the petiole. Thus, CaCDPK15 could be a crucial regulator of the petiole, stem elongation, and plant height in chickpea.
Expression Pattern in Seed Development in Two Contrasting Chickpea Cultivars
The seed is the most important part of the chickpea plant from the yield point of view. Therefore, to understand the role of the CDPK family in seed development, expression analysis was performed in two contrasting chickpea cultivars i.e., large-seeded JGK3 and small-seeded Himchana 1. The expression pattern was analyzed during seven stages of chickpea seed development representing early-embryogenesis (S1), mid-embryogenesis (S2), late embryogenesis (S3–S4), mid-maturation (S5), and late-maturation (S6–S7) stages (Rajkumar et al., 2020). Total 17 CaCDPK genes expressed differentially in both chickpea cultivars. Out of these, nine genes (CaCDPK2, 3, 6, 8, 9, 12, 14, 16, and 20) were found to be up-regulated, whereas eight genes (CaCDPK4, 7, 11, 13, 15, 18, 19, and 22) were down-regulated in one or more seed stages in both the cultivars (Figure 8, Supplementary Table 5). Notably, most up-regulated CaCDPK genes were induced during S1–S4, while their expression was insignificant during later seed stages (S5–S7). In contrast, the most down-regulated genes had significant negative expression during later stages of seed development (S5–S7). This indicates that different sets of CaCDPK genes are involved in positively regulating early to late embryogenesis, and negatively regulating mid-late maturation phases of seed development, both in JGK3 and Himchana1 cultivars. CaCDPK3 and CaCDPK14 showed similar expression patterns (up-regulation) from the S1–S4 stages with the highest expression in the S4 stage. however, their expression level was higher in JGK3 than Himchana 1. CaCDPK2 had significant expression during S3 and S4 stages in both the cultivars. CaCDPK4 showed common as well as cultivar-specific expression as it was specifically down-regulated in the S6 stage of both JGK3 and Himchana1, and S5 and S7 stages of Himchana1 and JGK3, respectively. CaCDPK15 showed a peculiar expression pattern, it was specifically down-regulated only in S2 and S4 stages of Himchana1. CaCDPK7 could be a major negative regulator of the entire process of seed development as it was ubiquitously down-regulated during all the seed stages in both chickpea cultivars. Like chickpea, rice orthologs of CDPKs, such as OsCPK21, OsCPK23, and OsCPK31 were found to be strongly up-regulated during different seed developmental stages (Ray et al., 2007). Overexpression of seed-specific OsCPK31 results in early grain filling, and seed maturation in rice (Manimaran et al., 2015). As suggested by previous reports, CaCDPKs could phosphorylate specific targets/substrates which are important components of the seed development process in chickpea. For instance, a castor (Ricinus communis) CDPK protein, RcCDPK2 phosphorylate a sucrose synthase, RcSUS1 at Ser-11 in developing castor oil seed (COS). RcSUS1 cleaves the imported sucrose, thereby, it generated the storage end product in developing COS (Fedosejevs et al., 2016). Similarly, RcCDPK1 in vivo phosphorylates Phosphoenolpyruvate carboxylase (PEPC) interacting bacterial-type (BTPC) subunits at Ser-451, and this phosphorylation is inhibitory for BTPC activity. Thus, RcCDPK1 regulates respiratory CO2 refixation and anaplerotic photosynthate partitioning in developing COS (Ying et al., 2017).
Figure 8. Expression profiles of CDPK genes during different stages of seed development in chickpea. Expression profiles are shown for seven seed developmental stages (S1–S7) in two chickpea cultivars: large-seeded JGK3 and small-seeded Himchana1. The genes are named on the left and stages/cultivars are labeled at the top. The scale bar represents the log_2 Fold Change FPKM values.
qRT-PCR Expression Analysis Under ABA, Drought, and Salinity Stress
The onslaught of abiotic stresses, such as drought and salinity triggers the increase in cytosolic Ca2+ concentration, leading to the generation of stimulus-specific “Ca2+ signatures.” These signatures once sensed, trigger specific and overlapping downstream signaling cascades toward an adaptive response. CDPKs are a major class of Ca2+ sensors in plant cells that plays a crucial role in the plant's adaptation to stress conditions. Therefore, to understand their role in abiotic stress signaling in chickpea, qRT-PCR analysis was performed for 10 selected genes i.e., CaCDPK3, 4, 5, 7, 8, 14, 16, 18, 21, and 22 under abiotic stresses such as drought, salt stress, and phytohormone ABA. Expression profiles were generated at different time points, in root and shoot tissues separately. Based on the expression fold change value ≥1.5 w.r.t. untreated control, all 10 CaCDPK genes were found to be differentially expressed under one or more stress conditions (Figure 9, Supplementary Table 6). In roots, except CaCDPK7 and 8, eight genes express differentially under drought conditions. CaCDPK5, 14, and 21 were significantly up-regulated after 1, 3, and 6 h of drought treatment, while CaCDPK22 was up-regulated after 3h and 6h of drought treatment. CaCDPK3 and 16 were exclusively induced after 1h of drought treatment and CaCDPK4 showed significant up-regulation after 6h drought only (Figure 9). Thus, these CaCDPK genes could be categorized as early inducing, late inducing, and ubiquitously expressed. Under ABA treatment, CaCDPK4, 14, 16, 18, and 22 were found to be differentially expressed. Drought-induced CaCDPK4, 16, and 22 were also up-regulated during the mid (6 h) and late stages (12 h) of the ABA treatment. ABA is known to regulate various aspects of plant growth, including root growth and seed germination (Sagar and Singh, 2019). An overlapping expression of CaCDPK4, 16, and 22 under ABA and drought indicates that these CaCDPKs might be involved in regulating root growth via the ABA-dependent pathway, and improving water uptake under drought stress. The ABA pathway commonly regulates response to osmotic stresses, including drought and salinity stress (Singh et al., 2016). Numerous studies have revealed the co-regulation of several genes under drought and salinity stress (Ray et al., 2007; Singh et al., 2015, 2017). Surprisingly, in our study, except CaCDPK3 most drought-induced CaCDPK genes were found to be down-regulated under salt stress. Salt stress is known to be inflicted through two components; osmotic and ionic stress (Verslues et al., 2006). The osmotic component is common in both drought and salt stress, whereas the ionic component is unique to salt stress. It is possible that CaCDPKs might regulate salt stress response through its ionic component only, however, it requires detailed functional investigation. Like in roots, CaCDPK genes are expressed differentially in one or more conditions in the shoot. Seven genes (CaCDPK4, 7, 8, 14, 16, 21, and 22) showed differential expression under drought stress conditions. CaCDPK8 was down-regulated in both the stresses and under ABA treatment. Similar to the root, CaCDPK21 showed a continuous induction (1–6 h) under drought stress, and CaCDPK4 was induced after 6 h. CaCDPK22 which showed strong induction after 3 and 6 h drought in roots was found to be up-regulated after 1 h but down-regulated after 3 and 6 h drought treatment. Interestingly, CaCDPK5 which was continuously expressed in roots showed insignificant expression in shoot whereas, CaCDPK7 which had insignificant expression in roots was up-regulated after 6 h in the shoot, under drought. This suggests that some CaCDPKs preferentially express in the root or shoot tissues, and possibly regulate tissue-specific functions. Four genes, including CaCDPK7, 14, 16, and 22 were found to be up-regulated both under drought stress and ABA treatment. ABA is produced under drought conditions and it regulates the stomata closure. Drought-activated CaCDPKs might phosphorylate and regulate guard cell ion channels to regulate the stomata movement in an ABA-dependent manner, to reduce the transpiration rate and water loss under drought stress in chickpea. Previously, AtCPK3, AtCPK6, AtCPK21, and AtCPK23 have been shown to phosphorylate and regulate SLOW ANION CHANNEL (SLAC1), Ca2+ permeable channels, and SLOW ANION CHANNEL-ASSOCIATED (SLAH3) to promote ABA-dependent stomata closure (Mori et al., 2006; Geiger et al., 2010, 2011). Besides phosphorylation, CaCDPKs may also regulate the stomata movement in chickpea through interaction with different proteins, as previously shown by the interaction of AtCPK8 with CATALASE3 (CAT3) and that of OsCPK21 with 14–3–3 protein (OsGF14e) (Zou et al., 2015; Chen et al., 2017). Strikingly, CaCDPK14 which was strongly induced during most of the seed developmental stages (Figure 8), was also induced under drought and ABA treatment. Careful analysis showed that the promoter of CaCDPK14 contains an ABRE element. As discussed earlier, seed development and drought stress are interconnected, and ABA is known to simultaneously regulate these processes. Probably by regulating drought response CaCDPK14 might control seed development via the ABA-dependent pathway in chickpea. Such genes could be the potential targets for improving abiotic stress tolerance and crop yield, thus could be of great biotechnological importance.
Figure 9. Expression profiles of CDPK genes under drought, salt, and ABA treatments in chickpea. qRT-PCR analysis was performed to generate the expression profile of CaCDPK genes under drought, salt, and ABA treatment in root and shoot of desi chickpea (ICC4958). Different treatments and time points are indicated on X-axis and the relative expression level of the gene is indicated on Y-axis. Each bar represents the mean value of three replicates. Standard error among the samples is indicated by error bars. *p < 0.05, **p < 0.01 and ***p < 0.005 for treated samples w.r.t. untreated control.
Unlike roots, at least four genes, including CaCDPK3, 5, 18, and 21 were up-regulated under salt stress in the shoot. CaCDPK21 showed continuous expression after 3, 6, and 12 h of salinity treatment, CaCDPK5 was induced after 3 and 12 h whereas, CaCDPK3 and 18 expressed during later stage (12 h) of salt stress (Figure 9, Supplementary Table 6). These CaCDPKs might regulate root-shoot translocation of Na+ ion, sequestration of excess Na+ from subcellular compartments into the vacuole, and in the maintenance of Na+/K+ homeostasis in plant cells under salinity stress. CPK3 (Mehlmer et al., 2010) and CPK27 (Zhao et al., 2015) have been shown to regulate Na+ toxicity and K+ uptake, respectively, in Arabidopsis that leads to salt stress tolerance. CDPKs also control oxidative damage or ROS homeostasis to protect the plants from the deleterious effects of stresses, including drought and salinity. For instance, OsCPK12 promotes the expression of ROS scavenger gene OsAPX2/OsAPX8 and suppresses the expression of NADPH oxidase gene OsRBOHI, thereby regulating ROS homeostasis and improving tolerance to salt stress in rice (Asano et al., 2012). Overall, this analysis revealed the spatiotemporal pattern of CaCDPK genes expression in chickpea. This hints toward specific and common phosphorylation targets/substrates of CaCDPKs under abiotic stresses in the root and shoot tissues in chickpea.
Conclusions
In conclusion, a total of 22 non-redundant CDPK encoding genes are identified in the chickpea genome. Protein structure and phylogenetic analysis revealed that the CDPK family is highly conserved in terms of structure and evolution, in higher plants. Subcellular localization suggests that most CaCDPK members may target various substrates in the cytoplasm and nucleus to regulate diverse functions. The expression pattern suggests the crucial role of several CDPK members in abiotic stress tolerance and plant development in chickpea. Importantly, some candidate CaCDPKs might regulate drought stress response and seed development. Such genes could be important in determining the crop yield thus, will be of great biotechnological importance. In the future, efforts must be directed for the identification of targets/substrates of key CaCDPKs to delineate CDPK mediated signaling networks underlying abiotic stress response and plant development. Ultimately, knowledge obtained from this study will help in developing abiotic stress-tolerant and high yielding chickpea cultivars that will contribute to world food security.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
AS and SK conceptualized the study, supervised the project, and designed all the analyses and experiments. DD and NP performed the experiments. DD, NP, SK, and AS analyzed the data. AS, DD, and NP wrote the manuscript. All authors read and approved the final version of manuscript.
Funding
The authors acknowledge financial support from NIPGR core grant. SK acknowledged the BT/PR40146/BTIS/137/4/2020 project grant from the Department of Biotechnology (DBT), Government of India. AS acknowledges research grant from Science and Engineering Research Board (SERB)—Department of Science and Technology (DST), Government of India (Grant No. EEQ/2018/000106).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
DD is thankful to the council of scientific and industrial research (CSIR), India for research fellowships. The authors are thankful to Mr. Ankit, senior research fellow in AS lab for compiling the references. DBT—eLibrary Consortium (DeLCON), is acknowledged for providing e-resources.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.831265/full#supplementary-material
References
Abe, H., Urao, T., Ito, T., Seki, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2003). Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 15, 63–78. doi: 10.1105/tpc.006130
Abe, M., Kobayashi, Y., Yamamoto, S., Daimon, Y., Yamaguchi, A., Ikeda, Y., et al. (2005). FD, a bZIP protein mediating signals. Science 309, 1052–1056. doi: 10.1126/science.1115983
Asai, S., Ichikawa, T., Nomura, H., Kobayashi, M., Kamiyoshihara, Y., Mori, H., et al. (2013). The variable domain of a plant calcium-dependent protein kinase (CDPK) confers subcellular localization and substrate recognition for NADPH oxidase. J. Biol. Chem. 288, 14332–14340. doi: 10.1074/jbc.M112.448910
Asano, T., Hayashi, N., Kikuchi, S., and Ohsugi, R. (2012). CDPK-mediated abiotic stress signaling. Plant Signal. Behav. 7, 817–821. doi: 10.4161/psb.20351
Asano, T., Tanaka, N., Yang, G., Hayashi, N., and Komatsu, S. (2005). Genome-wide identification of the rice calcium-dependent protein kinase. Plant Cell Physiol. 46, 356–366.
Benetka, W., Mehlmer, N., Maurer-Stroh, S., Sammer, M., Koranda, M., Neumüller, R., et al. (2008). Experimental testing of predicted myristoylation targets involved in asymmetric cell division and calcium-dependent signalling. Cell Cycle 7, 3709–3719. doi: 10.4161/cc.7.23.7176
Boudsocq, M., and Sheen, J. (2013). CDPKs in immune and stress signaling. Trends Plant Sci. 18, 30–40. doi: 10.1016/j.tplants.2012.08.008
Cai, H., Cheng, J., Yan, Y., Xiao, Z., Li, J., Mou, S., et al. (2015). Genome-wide identification and expression analysis of calcium-dependent protein kinase and its closely related kinase genes in Capsicum annuum. Front. Plant Sci. 6, 737. doi: 10.3389/fpls.2015.00737
Chen, Y., Zhou, X., Chang, S., Chu, Z., Wang, H., Han, S., et al. (2017). Calcium-dependent protein kinase 21 phosphorylates 14-3-3 proteins in response to ABA signaling and salt stress in rice. Biochem. Biophys. Res. Commun. 493, 1450–1456. doi: 10.1016/j.bbrc.2017.09.166
Cheng, S. H., Willmann, M. R., Chen, H. C., and Sheen, J. (2002). Calcium signaling through protein kinases. The Arabidopsis calcium-dependent protein kinase gene family. Plant Physiol. 129, 469–485. doi: 10.1104/pp.005645
Deepika, D., Ankit, J. S., Mali, K. V., Sinha, A. K., and Singh, A. (2022). Molecular analysis indicates the involvement of Jasmonic acid biosynthesis pathway in low-potassium (K+) stress response and development in chickpea (Cicer arietinum). Environ. Exp. Bot. 194, 104753. doi: 10.1016/j.envexpbot.2021.104753
Deepika, D., and Singh, A. (2021). Expression dynamics indicate the role of Jasmonic acid biosynthesis pathway in regulating macronutrient (N, P and K+) deficiency tolerance in rice (Oryza sativa L.). Plant Cell Rep. 40, 1495–1512. doi: 10.1007/s00299-021-02721-5
Fedosejevs, E. T., Gerdis, S. A., Ying, S., Pyc, M., Anderson, E. M., Snedden, W. A., et al. (2016). The calcium-dependent protein kinase RcCDPK2 phosphorylates sucrose synthase at Ser11 in developing castor oil seeds. Biochem. J. 473, 3667–3682. doi: 10.1042/BCJ20160531
Feng, Z. J., He, G. H., Zheng, W. J., Lu, P. P., Chen, M., Gong, Y. M., et al. (2015). Foxtail millet NF-Y families: genome-wide survey and evolution analyses identified two functional genes important in abiotic stresses. Front. Plant Sci. 6, 1142. doi: 10.3389/fpls.2015.01142
Geiger, D., Maierhofer, T., Al-Rasheid, K. A., Scherzer, S., Mumm, P., Liese, A., et al. (2011). Stomatal closure by fast abscisic acid signaling is mediated by the guard cell anion channel SLAH3 and the receptor RCAR1. Sci. Signal. 4, ra32. doi: 10.1126/scisignal.2001346
Geiger, D., Scherzer, S., Mumm, P., Marten, I. A. P. P. L., Ache, P., Matschi, S., et al. (2010). Guard cell anion channel SLAC1 is regulated by CDPK protein kinases with distinct Ca2+ affinities. Proc. Nat. Acad. Sci U.S.A. 107, 8023–8028. doi: 10.1073/pnas.0912030107
Gil, J., Nadal, S., Luna, D., Moreno, M. T., and Haro, A. D. (1996). Variability of some physico-chemical characters in Desi and Kabuli chickpea types. J. Sci. Food Agric. 71, 179–184. doi: 10.1002/(SICI)1097-0010(199606)71:2<179::AID-JSFA566>3.0.CO;2-7
Gutermuth, T., Lassig, R., Portes, M. T., Maierhofer, T., Romeis, T., Borst, J. W., et al. (2013). Pollen tube growth regulation by free anions depends on the interaction between the anion channel SLAH3 and calcium-dependent protein kinases CPK2 and CPK20. Plant Cell 25, 4525–4543. doi: 10.1105/tpc.113.118463
Hamel, L. P., Sheen, J., and Séguin, A. (2014). Ancient signals: comparative genomics of green plant CDPKs. Trends Plant Sci. 19, 79–89. doi: 10.1016/j.tplants.2013.10.009
Heinrich, M., Hettenhausen, C., Lange, T., Wünsche, H., Fang, J., Baldwin, I. T., et al. (2013). High levels of jasmonic acid antagonize the biosynthesis of gibberellins and inhibit the growth of N icotiana attenuata stems. Plant J. 73, 591–606. doi: 10.1111/tpj.12058
Hetherington, A. M (2001). Guard cell signaling. Cell 107, 711–714. doi: 10.1016/S0092-8674(01)00606-7
Hettenhausen, C., Sun, G., He, Y., Zhuang, H., Sun, T., Qi, J., et al. (2016). Genome-wide identification of calcium-dependent protein kinases in soybean and analyses of their transcriptional responses to insect herbivory and drought stress. Sci. Rep. 6, 18973. doi: 10.1038/srep18973
Hou, F., Du, T., Qin, Z., Xu, T., Li, A., Dong, S., et al. (2021). Genome-wide in silico identification and expression analysis of beta-galactosidase family members in sweet potato [Ipomoea batatas (L.) Lam]. BMC Genomics 22, 140. doi: 10.1186/s12864-021-07436-1
Hrabak, E. M., Chan, C. W., Gribskov, M., Harper, J. F., Choi, J. H., Halford, N., et al. (2003). The Arabidopsis CDPK-SnRK superfamily of protein kinases. Plant Physiol. 132, 666–680. doi: 10.1104/pp.102.011999
Hu, Z., Lv, X., Xia, X., Zhou, J., Shi, K., Yu, J., et al. (2016). Genome-Wide identification and expression analysis of calcium-dependent protein kinase in tomato. Front. Plant Sci. 7, 469. doi: 10.3389/fpls.2016.00469
Hubbard, K. E., Nishimura, N., Hitomi, K., Getzoff, E. D., and Schroeder, J. I. (2010). Early abscisic acid signal transduction mechanisms: newly discovered components and newly emerging questions. Genes Dev. 24, 1695–1708. doi: 10.1101/gad.1953910
Hughes, J. D., Estep, P. W., Tavazoie, S., and Church, G. M. (2000). Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae. J. Mol. Biol. 296, 1205–1214. doi: 10.1006/jmbi.2000.3519
Ibrikci, H., Knewtson, S. J., and Grusak, M. A. (2003). Chickpea leaves as a vegetable green for humans: evaluation of mineral composition. J. Sci. Food Agric. 83, 945–950. doi: 10.1002/jsfa.1427
Kawamoto, N., Sasabe, M., Endo, M., Machida, Y., and Araki, T. (2015). Calcium-dependent protein kinases responsible for the phosphorylation of a bZIP transcription factor FD crucial for the florigen complex formation. Sci. Rep. 5, 8341. doi: 10.1038/srep08341
Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N., and Sternberg, M. J. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10, 845–858. doi: 10.1038/nprot.2015.053
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Kong, X., Lv, W., Jiang, S., Zhang, D., Cai, G., Pan, J., et al. (2013). Genome-wide identification and expression analysis of calcium-dependent protein kinase in maize. BMC Genomics 14, 433. doi: 10.1186/1471-2164-14-433
Letunic, I., and Bork, P. (2021). Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296. doi: 10.1093/nar/gkab301
Li, A. L., Zhu, Y. F., Tan, X. M., Wang, X., Wei, B., Guo, H. Z., et al. (2008). Evolutionary and functional study of the CDPK gene family in wheat (Triticum aestivum L.). Plant Mol. Biol. 66, 429–443. doi: 10.1007/s11103-007-9281-5
Liese, A., and Romeis, T. (2013). Biochemical regulation of in vivo function of plant calcium-dependent protein kinases (CDPK). Biochim. Biophys. Acta. 1833, 1582–1589. doi: 10.1016/j.bbamcr.2012.10.024
Liu, K. H., Niu, Y., Konishi, M., Wu, Y., Du, H., Chung, H. S., et al. (2017). Discovery of nitrate–CPK–NLP signalling in central nutrient–growth networks. Nature 545, 311–316. doi: 10.1038/nature22077
Liu, W., Li, W., He, Q., Daud, M. K., Chen, J., and Zhu, S. (2014). Genome-wide survey and expression analysis of calcium-dependent protein kinase in Gossypium raimondii. PLoS ONE 9, e98189. doi: 10.1371/journal.pone.0098189
Lu, S. X., and Hrabak, E. M. (2013). The myristoylated amino-terminus of an Arabidopsis calcium-dependent protein kinase mediates plasma membrane localization. Plant Mol. Biol. 82, 267–278. doi: 10.1007/s11103-013-0061-0
Luan, S (2009). The CBL-CIPK network in plant calcium signaling. Trends Plant Sci. 14, 37–42. doi: 10.1016/j.tplants.2008.10.005
Manimaran, P., Mangrauthia, S. K., Sundaram, R. M., and Balachandran, S. M. (2015). Constitutive expression and silencing of a novel seed specific calcium dependent protein kinase gene in rice reveals its role in grain filling. J. Plant Physiol. 174, 41–48. doi: 10.1016/j.jplph.2014.09.005
Matschi, S., Werner, S., Schulze, W. X., Legen, J., Hilger, H. H., and Romeis, T. (2013). Function of calcium-dependent protein kinase CPK 28 of Arabidopsis thaliana in plant stem elongation and vascular development. Plant J. 73, 883–896. doi: 10.1111/tpj.12090
Mehlmer, N., Wurzinger, B., Stael, S., Hofmann-Rodrigues, D., Csaszar, E., Pfister, B., et al. (2010). The Ca(2+) -dependent protein kinase CPK3 is required for MAPK-independent salt-stress acclimation in Arabidopsis. Plant J. 63, 484–498. doi: 10.1111/j.1365-313X.2010.04257.x
Mishra, M., Kanwar, P., Singh, A., Pandey, A., Kapoor, S., and Pandey, G. K. (2013). Plant omics: genome-wide analysis of ABA repressor1 (ABR1) related genes in rice during abiotic stress and development. Omics 17, 439–450. doi: 10.1089/omi.2012.0074
Mittal, S., Mallikarjuna, M. G., Rao, A. R., Jain, P. A., Dash, P. K., and Thirunavukkarasu, N. (2017). Comparative analysis of CDPK family in maize, arabidopsis, rice, and sorghum revealed potential targets for drought tolerance improvement. Front. Chem. 5, 115. doi: 10.3389/fchem.2017.00115
Mori, I. C., Murata, Y., Yang, Y., Munemasa, S., Wang, Y. F., Andreoli, S., et al. (2006). CDPKs CPK6 and CPK3 function in ABA regulation of guard cell S-type anion-and Ca2+-permeable channels and stomatal closure. PLoS Biol. 4, e327. doi: 10.1371/journal.pbio.0040327
Myers, C., Romanowsky, S. M., Barron, Y. D., Garg, S., Azuse, C. L., Curran, A., et al. (2009). Calcium-dependent protein kinases regulate polarized tip growth in pollen tubes. Plant J. 59, 528–539. doi: 10.1111/j.1365-313X.2009.03894.x
Nakashima, K., Fujita, Y., Kanamori, N., Katagiri, T., Umezawa, T., Kidokoro, S., et al. (2009). Three Arabidopsis SnRK2 protein kinases, SRK2D/SnRK2. 2, SRK2E/SnRK2. 6/OST1 and SRK2I/SnRK2. 3, involved in ABA signaling are essential for the control of seed development and dormancy. Plant Cell Physiol. 50, 1345–1363. doi: 10.1093/pcp/pcp083
Nakashima, K., and Yamaguchi-Shinozaki, K. (2013). ABA signaling in stress-response and seed development. Plant Cell Rep. 32, 959–970. doi: 10.1007/s00299-013-1418-1
Neelamathi, E., Vasumathi, E., Bagyalakshmi, S., and Kannan, R. (2009). In silico prediction of structure and functional aspects of a hypothetical protein of Neurospora crassa. J. Cell Tissue Res. 9, 1989.
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Poovaiah, B. W., Du, L., Wang, H., and Yang, T. (2013). Recent advances in calcium/calmodulin-mediated signaling with an emphasis on plant-microbe interactions. Plant Physiol. 163, 531–542. doi: 10.1104/pp.113.220780
Rajkumar, M. S., Gupta, K., Khemka, N. K., Garg, R., and Jain, M. (2020). DNA methylation reprogramming during seed development and its functional relevance in seed size/ weight determination in chickpea. Commun. Biol. 3, 340. doi: 10.1038/s42003-020-1059-1
Ray, S., Agarwal, P., Arora, R., Kapoor, S., and Tyagi, A. K. (2007). Expression analysis of calcium-dependent protein kinase gene family during reproductive development and abiotic stress conditions in rice (Oryza sativa L. ssp. indica). Mol. Genet. Genom. 278, 493–505. doi: 10.1007/s00438-007-0267-4
Ryan, J. G (1997). “A global perspective on pigeonpea and chickpea sustainable production system : Present status and future potential,” in Recent Advances in Pulses Research, eds A. N. Asthana and M. Ali (Kanpur: Indian Institute of Pulses Research), 1–32.
Sagar, S., Biswas, D. K., Chandrasekar, R., and Singh, A. (2021). Genome-wide identification, structure analysis and expression profiling of phospholipases D under hormone and abiotic stress treatment in chickpea (Cicer arietinum). Int. J. Biol. Macromol. 169, 264–273. doi: 10.1016/j.ijbiomac.2020.12.102
Sagar, S., Biswas, D. K., and Singh, A. (2020). Genomic and expression analysis indicate the involvement of phospholipase C family in abiotic stress signaling in chickpea (Cicer arietinum). Gene 753, 144797. doi: 10.1016/j.gene.2020.144797
Sagar, S., and Singh, A. (2019). “Abscisic acid, a principal regulator of plant abiotic stress responses,” in Plant Signaling Molecules, eds I. Khan, S. P. Reddy, A. Ferrante, and N. S. Khan (Woodhead Publishing), 341–373. doi: 10.1016/B978-0-12-816451-8.00021-6
Schroeder, J. I., Kwak, J. M., and Allen, G. J. (2001). Guard cell abscisic acid signalling and engineering drought hardiness in plants. Nature 410, 327–330. doi: 10.1038/35066500
Schulz, P., Herde, M., and Romeis, T. (2013). Calcium-dependent protein kinases: hubs in plant stress signaling and development. Plant Physiol. 163, 523–530. doi: 10.1104/pp.113.222539
Sharma, N., Russell, S. D., Bhalla, P. L., and Singh, M. B. (2011). Putative cis-regulatory elements in genes highly expressed in rice sperm cells. BMC Res. Notes 4, 319. doi: 10.1186/1756-0500-4-319
Simeunovic, A., Mair, A., Wurzinger, B., and Teige, M. (2016). Know where your clients are: subcellular localization and targets of calcium-dependent protein kinases. J. Exp. Bot. 67, 3855–3872. doi: 10.1093/jxb/erw157
Singh, A., Baranwal, V., Shankar, A., Kanwar, P., Ranjan, R., Yadav, S., et al. (2012). Rice phospholipase A superfamily: organization, phylogenetic and expression analysis during abiotic stresses and development. PLoS ONE 7, e30947. doi: 10.1371/journal.pone.0030947
Singh, A., Giri, J., Kapoor, S., Tyagi, A. K., and Pandey, G. K. (2010). Protein phosphatase complement in rice: genome-wide identification and transcriptional analysis under abiotic stress conditions and reproductive development. BMC Genomics 11, 435. doi: 10.1186/1471-2164-11-435
Singh, A., Jha, S. K., Bagri, J., and Pandey, G. K. (2015). ABA inducible rice protein phosphatase 2C confers ABA insensitivity and abiotic stress tolerance in Arabidopsis. PLoS ONE 10, e0125168. doi: 10.1371/journal.pone.0125168
Singh, A., Kanwar, P., Pandey, A., Tyagi, A. K., Sopory, S. K., Kapoor, S., et al. (2013). Comprehensive genomic analysis and expression profiling of phospholipase C gene family during abiotic stresses and development in rice. PLoS ONE 8, e62494. doi: 10.1371/journal.pone.0062494
Singh, A., Kanwar, P., Yadav, A. K., Mishra, M., Jha, S. K., Baranwal, V., et al. (2014). Genome-wide expressional and functional analysis of calcium transport elements during abiotic stress and development in rice. FEBS J. 281, 894–915. doi: 10.1111/febs.12656
Singh, A., Pandey, A., Srivastava, A. K., Tran, L. S. P., and Pandey, G. K. (2016). Plant protein phosphatases 2C: from genomic diversity to functional multiplicity and importance in stress management. Crit. Rev. Biotechnol. 36, 1023–1035. doi: 10.3109/07388551.2015.1083941
Singh, A., and Pandey, G. K. (2015). “Primer design using primer express® for SYBR green-based quantitative PCR,” in PCR Primer Design (New York, NY: Humana Press), 153–164. doi: 10.1007/978-1-4939-2365-6_11
Singh, A., Sagar, S., and Biswas, D. K. (2017). Calcium dependent protein kinase, a versatile player in plant stress management and development. CRC Crit. Rev. Plant Sci. 36, 336–352. doi: 10.1080/07352689.2018.1428438
Tang, R. J., Wang, C., Li, K., and Luan, S. (2020). The CBL–CIPK calcium signaling network: unified paradigm from 20 years of discoveries. Trends Plant Sci. 25, 604–617. doi: 10.1016/j.tplants.2020.01.009
Varshney, R. K., Hiremath, P. J., Lekha, P., Kashiwagi, J., Balaji, J., Deokar, A. A., et al. (2009). A comprehensive resource of drought- and salinity- responsive ESTs for gene discovery and marker development in chickpea (Cicer arietinum L.). BMC Genomics 10, 523. doi: 10.1186/1471-2164-10-523
Varshney, R. K., Song, C., Saxena, R. K., Azam, S., Yu, S., Sharpe, A. G., et al. (2013). Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat. Biotechnol. 31, 240–246. doi: 10.1038/nbt.2491
Varshney, R. K., Thudi, M., Nayak, S. N., Gaur, P. M., Kashiwagi, J., Krishnamurthy, L., et al. (2014). Genetic dissection of drought tolerance in chickpea (Cicer arietinum L.). Theor. Appl. Genet. 127, 445–462. doi: 10.1007/s00122-013-2230-6
Verslues, P. E., Agarwal, M., Katiyar-Agarwal, S., Zhu, J., and Zhu, J. K. (2006). Methods and concepts in quantifying resistance to drought, salt and freezing, abiotic stresses that affect plant water status. Plant J. 45, 523–539. doi: 10.1111/j.1365-313X.2005.02593.x
Vivek, P. J., Resmi, M. S., Sreekumar, S., Sivakumar, K. C., Tuteja, N., and Soniya, E. V. (2017). Calcium-dependent protein kinase in ginger binds with Importin-α through its junction domain for nuclear localization, and further interacts with NAC transcription factor. Front. Plant Sci. 7, 1909. doi: 10.3389/fpls.2016.01909
Wang, J. P., Xu, Y. P., Munyampundu, J. P., Liu, T. Y., and Cai, X. Z. (2016). Calcium-dependent protein kinase (CDPK) and CDPK-related kinase (CRK) gene families in tomato: genome-wide identification and functional analyses in disease resistance. Mol. Genet. Genom. 291, 661–676. doi: 10.1007/s00438-015-1137-0
Wang, Y., Tang, H., DeBarry, J. D., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49–e49. doi: 10.1093/nar/gkr1293
Wen, F., Ye, F., Xiao, Z., Liao, L., Li, T., Jia, M., et al. (2020). Genome-wide survey and expression analysis of calcium-dependent protein kinase (CDPK) in grass Brachypodium distachyon. BMC Genomics 21, 53. doi: 10.1186/s12864-020-6475-6
Witte, C. P., Keinath, N., Dubiella, U., Demoulière, R., Seal, A., and Romeis, T. (2010). Tobacco calcium-dependent protein kinases are differentially phosphorylated in vivo as part of a kinase cascade that regulates stress response. J. Biol. Chem. 285, 9740–9748. doi: 10.1074/jbc.M109.052126
Wu, P., Wang, W., Duan, W., Li, Y., and Hou, X. (2017). Comprehensive analysis of the CDPK-SnRK superfamily genes in Chinese cabbage and its evolutionary implications in plants. Front. Plant Sci. 8, 162. doi: 10.3389/fpls.2017.00162
Xiao, X. H., Yang, M., Sui, J. L., Qi, J. Y., Fang, Y. J., Hu, S. N., et al. (2017). The calcium-dependent protein kinase (CDPK) and CDPK-related kinase gene families in Hevea brasiliensis—comparison with five other plant species in structure, evolution, and expression. FEBS Open Bio 7, 4–24. doi: 10.1002/2211-5463.12163
Yamaguchi-Shinozaki, K., and Shinozaki, K. (2005). Organization of cis-acting regulatory elements in osmotic-and cold-stress-responsive promoters. Trends Plant Sci. 10, 88–94. doi: 10.1016/j.tplants.2004.12.012
Yang, Y., Wang, Q., Chen, Q., Yin, X., Qian, M., Sun, X., et al. (2017). Genome-wide survey indicates diverse physiological roles of the barley (Hordeum vulgare L.) calcium-dependent protein kinase genes. Sci. Rep. 7, 5306. doi: 10.1038/s41598-017-05646-w
Ying, S., Hill, A. T., Pyc, M., Anderson, E. M., Snedden, W. A., Mullen, R. T., et al. (2017). Regulatory phosphorylation of bacterial-type PEP carboxylase by the Ca2+-dependent protein kinase RcCDPK1 in developing castor oil seeds. Plant Physiol. 174, 1012–1027. doi: 10.1104/pp.17.00288
Yu, C. S., Chen, Y. C., Lu, C. H., and Hwang, J. K. (2006). Prediction of protein subcellular localization. Proteins 64, 643–651. doi: 10.1002/prot.21018
Zeng, H., Xu, L., Singh, A., Wang, H., Du, L., and Poovaiah, B. W. (2015). Involvement of calmodulin and calmodulin-like proteins in plant responses to abiotic stresses. Front. Plant Sci. 6, 600. doi: 10.3389/fpls.2015.00600
Zhang, H., Wei, C., Yang, X., Chen, H., Yang, Y., Mo, Y., et al. (2017). Genome-wide identification and expression analysis of calcium-dependent protein kinase and its related kinase gene families in melon (Cucumis melo L.). PLoS ONE 12, e0176352. doi: 10.1371/journal.pone.0176352
Zhang, K., Han, Y. T., Zhao, F. L., Hu, Y., Gao, Y. R., Ma, Y. F., et al. (2015). Genome-wide identification and expression analysis of the CDPK gene family in grape, Vitis spp. BMC Plant Biol. 15, 164. doi: 10.1186/s12870-015-0552-z
Zhao, L. N., Shen, L. K., Zhang, W. Z., Zhang, W., Wang, Y., and Wu, W. H. (2013). Ca2+-dependent protein kinase11 and 24 modulate the activity of the inward rectifying K+ channels in Arabidopsis pollen tubes. Plant Cell 25, 649–661. doi: 10.1105/tpc.112.103184
Zhao, P., Liu, Y., Kong, W., Ji, J., Cai, T., and Guo, Z. (2021). Genome-Wide identification and characterization of calcium-dependent protein kinase (CDPK) and CDPK-related kinase (CRK) gene families in medicago truncatula. Int. J. Mol. Sci. 22, 1044. doi: 10.3390/ijms22031044
Zhao, R., Sun, H., Zhao, N., Jing, X., Shen, X., and Chen, S. (2015). The Arabidopsis Ca2+-dependent protein kinase CPK27 is required for plant response to salt-stress. Gene 563, 203–214. doi: 10.1016/j.gene.2015.03.024
Zheng, L., Liu, P., Liu, Q., Wang, T., and Dong, J. (2019). Dynamic protein S-acylation in plants. Int. J. Mol. Sci. 20, 560. doi: 10.3390/ijms20030560
Zou, J. J., Li, X. D., Ratnasekera, D., Wang, C., Liu, W. X., Song, L. F., et al. (2015). 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 27, 1445–1460. doi: 10.1105/tpc.15.00144
Keywords: abiotic stress, CDPK, chickpea, gene expression, signaling, development
Citation: Deepika D, Poddar N, Kumar S and Singh A (2022) Molecular Characterization Reveals the Involvement of Calcium Dependent Protein Kinases in Abiotic Stress Signaling and Development in Chickpea (Cicer arietinum). Front. Plant Sci. 13:831265. doi: 10.3389/fpls.2022.831265
Received: 08 December 2021; Accepted: 11 March 2022;
Published: 12 April 2022.
Edited by:
Girdhar Kumar Pandey, University of Delhi, IndiaReviewed by:
Gaurav Zinta, Institute of Himalayan Bioresource Technology (CSIR), IndiaSung Chul Lee, Chung-Ang University, South Korea
Aryadeep Roychoudhury, St. Xavier's College (Autonomous), Kolkata, India
Copyright © 2022 Deepika, Poddar, Kumar and Singh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Amarjeet Singh, amarjeet.singh@nipgr.ac.in; orcid.org/0000-0001-9538-5406; Shailesh Kumar, shailesh@nipgr.ac.in; orcid.org/0000-0002-1872-9903
†These authors have contributed equally to this work and share first authorship