- Research article
- Open Access
Genome-wide identification of Calcineurin B-Like (CBL) gene family of plants reveals novel conserved motifs and evolutionary aspects in calcium signaling events
BMC Plant Biologyvolume 15, Article number: 189 (2015)
Calcium ions, the most versatile secondary messenger found in plants, are involved in the regulation of diverse arrays of plant growth and development, as well as biotic and abiotic stress responses. The calcineurin B-like proteins are one of the most important genes that act as calcium sensors.
In this study, we identified calcineurin B-like gene family members from 38 different plant species and assigned a unique nomenclature to each of them. Sequence analysis showed that, the CBL proteins contain three calcium binding EF-hand domain that contains several conserved Asp and Glu amino acid residues. The third EF-hand of the CBL protein was found to posses the D/E-x-D calcium binding sensor motif. Phylogenetic analysis showed that, the CBL genes fall into six different groups. Additionally, except group B CBLs, all the CBL proteins were found to contain N-terminal palmitoylation and myristoylation sites. An evolutionary study showed that, CBL genes are evolved from a common ancestor and subsequently diverged during the course of evolution of land plants. Tajima’s neutrality test showed that, CBL genes are highly polymorphic and evolved via decreasing population size due to balanced selection. Differential expression analysis with cold and heat stress treatment led to differential modulation of OsCBL genes.
The basic architecture of plant CBL genes is conserved throughout the plant kingdom. Evolutionary analysis showed that, these genes are evolved from a common ancestor of lower eukaryotic plant lineage and led to broadening of the calcium signaling events in higher eukaryotic organisms.
In various biological processes, calcium signals play a vital role as intracellular secondary messengers because of their strong homeostatic mechanism, which maintains an intracellular free Ca2+ concentration . The concentration of calcium ions varies from 30 to 400 nM in resting cells and in millimolar range in organelles [2–4]. For cytosolic Ca2+ ion to be transported from cytosol to other parts of the cell, a low cellular level needs to be maintained. This can be achieved through the action of Ca2+-ATPase pump, which transports Ca2+ ions out of the cell across the plasma membrane, and sarco-endoplasmic reticulum Ca2+-ATPases that pump Ca2+ into the lumen of the endoplasmic reticulum . It has been reported that, once cells began to use high-efficiency phosphate compounds as metabolic currency, they faced great challenges in maintaining low levels of intracellular Ca2+  to prevent precipitation of calcium and phosphate salt in the cytosol, which ultimately forms a solid, bone-like structure. Since Ca2+ ion is a versatile signaling ion, it plays different roles across signaling cascades to regulate gene expression in plants . Indeed, Ca2+ signals are important regulator of growth, development, and biotic and abiotic stresses in plants . The signaling information encoded by Ca2+ ions is decoded and transmitted by calcium sensors of Ca2+-binding proteins [8, 9]. Such sensors binds Ca2+ ion and changes their conformation in a Ca2+ dependent manner in the presence of high levels of Mg2+ and monovalent cations [1, 10]. Some of the calcium sensor includes (i) calcium dependent protein kinases (CPKs), (ii) calmodulines (CaMs) and (iii) calcineurin B-like proteins (CBLs) [7, 11]. The CPKs are monomeric proteins with unique structures that contain five domains, the (i) N-terminal variable domain, (ii) kinase domain, (iii) an auto-inhibitory domain, (iv) a regulatory domain and (v) C-terminal domain. The regulatory domain of CPK is characterized by the presence of four Ca2+ binding EF (elongation factor)-hands. The EF-hands are calcium sensors characterized by the presence of a conserved Asp (D) or Glu (E) residue . The EF-hand motifs are highly conserved, with a helix-loop-helix structure of 36 amino acid residues in each EF-hand. Unlike CPKs, CaMs and CBLs are small proteins that lack effector kinase domain (Fig. 1). The CaMs contain four Ca2+ binding EF-hands, whereas CBL contains three (Fig. 1) . To transmit Ca2+ signals, CPKs, CBLs and CaMs interact with their target proteins, and regulate their gene expression . These target proteins are may be protein kinases, metabolic enzymes, or cyto-skeletal associated proteins. The CIPKs (CBL-interacting protein kinase) are important target proteins of CBLs .
Although a great deal of effort has been made to investigate of the role of CBL genes, there has been very little effort made to determine the exact characteristics of these genes. Therefore, in this study, we identified CBL gene family members from 38 different plant species and assigned a unique nomenclature system to them. Additionally, we investigated the gene expression, genomics, phylogenetics and evolutionary aspects of these CBL genes.
Results and discussion
Nomenclature of CBL genes
To date, different members of specific gene families have been named according to the serial number by which they were identified. If no CBL gene has been identified for a given plant species to date, the first one identified is named CBL1, the next one as CBL2 and so on, regardless of the orthologous sequence similarity with the known counterpart genes. The volume of genomic sequence data are increasing daily, providing an excellent platform for genomics study. However, lack of a systemic nomenclature system for specific genes or gene families has led to confusion and difficulty in understanding the ever increasing genomic information. For example, the AtCBL1 gene differentially regulates salt, drought, and cold responses in Arabidopsis , but it is not clear whether the OsCBL1 gene also confers the same functionality. In principle, sequence similarity confers the structural similarity and structural similarity confers the functional similarity of a gene [16, 17]. Accordingly, AtCBL1 and OsCBL1 may confer more or less similar function. However, lack of a proper nomenclature system makes it very difficult to understand its function properly. Orthology lends the legitimacy to transfer functional information from an experimentally characterized protein to an uncharacterized one [18, 19]. Accordingly, an orthology based nomenclature system was adopted to name all CBL genes identified during this study as proposed by different researchers [7, 20–23]. In this system, Arabidopsis thaliana and Oryza sativa CBL protein sequences were taken as orthologous query genes. In the naming system, the first letter of the genus was kept upper case and the first letter of the species was kept lower case followed by CBL and then A. thaliana or Oryza sativa based CBL gene number. The monocot plants were named according to O. sativa, while dicot and other plants were named according to A. thaliana. In the case of monocot plants, the CBL gene number was assigned according to the orthologous gene of Oryza sativa. If more than one ortholog was found in a particular species, additional numbers followed by a hyphen were used to distinguish between paralogs. When the first letter of the genus and species of an organism coincided with another organism, the first letter of the genus was kept constant and the first, second, third or fourth letter or including the first, second, third and fourth letter of the species were taken into consideration. For example, the CBL gene of Capsella rubella was named as CrCBL, while Chlamydomonas reinhardtii was named as CreinCBL. In this case, both the letter of the genus and species name coincided with each other; therefore, the CBL gene of C. reinhardtii was denoted as CreinCBL. This nomenclature system can also provide information about the related orthologous species. The unique orthologous gene of one species may resemble the orthologous counterpart gene of another species and have undergone similar cellular function. The same approaches are usually used to predict the potential function for a newly sequenced gene and its protein product. It is very difficult to investigate the roles of all CBL genes in all plant species with different functional aspects. Therefore, the orthology based nomenclature system of the CBL gene will help to provide the basic information required for the counterpart orthologous gene.
Genomics of CBL genes
The genome of a species is regarded as a bag of genes that contain all information’s necessary to bridge the gap between genotype and phenotype . In the next decade, the genome sequences of virtually all angiosperms as well as important green algae, bryophytes, pteridophytes and gymnosperms will be completed. These genome sequences will become valuable tools that can provide a powerful framework for relating genome-level events to decipher the morphological and physiological variations that have contributed to colonization from aquatic habitats to land habitats. Genome-wide analysis of CBL genes across 38 different plant species revealed the presence of 328 CBL genes (Table 1). Among these, G. raimondii was found to contain the highest number of CBL genes (13) among higher land plants. The lower algae like Chlamydomonas and Micromonas contain only 2 and 3 CBL genes, respectively, in their genome. The bryophyte plant, Physcomitrella patens, and the pteridophyte plant, Selagnella moellendorffii, only encodes four CBL genes. The numbers of CBL genes found in P. patens is in accordance with the study of Kleist et al. . The model gymnosperm plant, Picea abies, encodes 13 CBL genes. The genome size of an organism varies from species to species (Table 1). Among the monocot plant, Zea mays has the biggest genome (2500 Mbs) and encodes for 9 CBL genes where as among the dicot plants, Glycine max has the biggest genome (975 Mbs) and encodes for 9 CBL genes. The genome size of gymnosperm plant Picea abies is 1960 Mbs and encodes for 13 CBL genes. Similarly, the dicot plant Capsella rubella has the smallest genome (134.8 Mbs) and still contains 9 CBL genes in its genome. From this study, it is clear that, there is no correlation between the genome size and number of CBL genes in plants. In the case of blue green algae Micromonas pusila, its genome size is 22 Mbs and still contains 3 CBL genes whereas, the genome size of Chlamydomonas reinhardtii is 118.8 Mbs and only contains 2 CBL genes. The presence of specific numbers of CBL genes in its genome is independent of genome size and it might be correlated with functional evolutionary requirements of the plant. All CBL genes identified during this study contains only three calcium binding EF-hands. In our investigation, we did not find any CBL genes from green algae species Coccomyxa subellipsoidea, Ostreococcus lucimarinus or Volvox carteri. The CBL genes contain a maximum of six, seven, eight or nine introns in their gene; while only a few CBL genes are intronless (Additional file 1). The CBL genes of Picea abies are intronless. Other lower eukaryotic intronless CBL genes found during this study are from M. pusila (MpCBL2), P. patens (PpCBL3-3) and S. moellendorffii (SmCBL5), while higher eukaryotic intronless CBL genes were found from S. lycopersicum (SlCBL3-3) and S. tuberosum (StCBL3-3) (Additional file 1). The CBL gene of F. vesca FvCBL4 was found to be the largest CBL gene and posses an ORF (open reading frame) of 3048 nucleotides that encodes for 1015 amino acids. Similarly, the CBL gene of M. domestica MdCBL5 encodes the smallest CBL gene and that contains only 426 nucleotides ORF that encodes for 141 amino acids. The genome of Z. mays is the largest one, containing only nine CBL genes, whereas the genome of M. pusila is smallest one, with only two CBL genes. However, as shown in Table 1, larger genome size is not directly proportional to more CBL gene numbers. The molecular weights of CBL proteins are vary from 12.774 (PaCBL10) to 115.266 (FvCBL4) kDa, while the isoelectric point (pI) are ranges from 4.02 to 9.61. The majority of CBL proteins are acidic (Additional file 2). Based on the average amino acid composition of CBL proteins, the abundance of most important calcium sensing amino acids, Asp (D) and Glu (E) were found to be 8.07 and 8.94, respectively (Table 2). The average abundance of Trp and Cys amino acids in CBL proteins were 0.62 and 1.27, respectively.
The genome sizes of plants are remarkably diverse and vary from species to species, with sizes that range from 63 (Genlisea aurea) to 149,000 Mbs (Paris japonica), divided into n = 2 to approximately n = 600 chromosomes and remains constant within a species . In this study, we found that the dicot plant Arabidopsis thaliana and Carica papaya (135 Mbs) have the smallest genome size, whereas in the monocot plant Zea mays (2500 Mbs) have the largest genome size among the higher plants. The lower eukaryotic algae, Micromonas pusila (22 Mbs), contains the smallest genome among the investigated species. The gymnosperms are characterized by the presence of a very large genome (up to 35,000 Mb), and Picea abies contains 1960 Mbs genome . Despite their larger genome, gymnosperms do not have higher numbers of chromosomes, with the number ranging between 2n = 2x = 14-28. Arabidopsis genome sequencing was initiated based on the thinking that genes and gene sequences of Arabidopsis would be similar to those of other plants, which was later found to be true; however, the number of protein coding genes varied significantly. This also found to be true in this study as the numbers of protein coding genes vary in a specific gene family of a specific plant. The nuclear DNA of plant consists of a low copy number of coding sequences, introns, promoters and regulatory DNA sequences . In this study, the majority of CBL genes were found to have either six, seven or eight introns within it, suggesting, the presence of intron number within a specific gene family varies from species to species, as well as in their counterpart orthologous gene(s).
It is well known that individual genes and entire genome can vary significantly in nucleotide compositions [28, 29]. The mutational process and relationship between the primary structure and function of a protein is considered as the major determinants of amino acid composition and rate of protein evolution . The natural selection events usually enhances the protein specificity and stability by favouring codons that encodes particular amino acids in a specific genic region . However, metabolic constraints on protein structure and composition could include the energetic cost of amino acid biosynthesis. The biosynthesis of aromatic amino acids like Trp requires higher energy (74.3 unit) and hence the average abundance of Trp amino acid per CBL gene is only 0.62 amino acids . High energy is required to synthesize Trp amino acids, so plants have encoded only 0.6 amino acids per CBL protein to avoid extra energy expense. Similarly, 12.7 and 15.3 units of energy is require for biosynthesis of Asp and Glu amino acid, respectively. Biosynthesis of Asp and Glu amino acid is relatively less costly; hence, plants encoded 8.07 and 8.94 amino acids, respectively, per CBL protein. As plants use a substantial amount of energy for biosynthesis of amino acids, there is an advantage to encode less costly amino acid in their protein .
Multiple sequence alignment of the CBL proteins revealed the presence of several new conserved domains and motifs. The CBL proteins of the plant kingdom contain only three EF-hand domains and are conserved. Overall, each EF-hand is 36 amino acids in length and has a helix-loop-helix structure . Each helix loop contains 12 amino acids within it; hence, each EF-hand contains 36 amino acids. Multiple sequence alignment revealed that, Asp (D) amino acid is less significantly conserved at position 7 and 11 in the first EF-hand, but most significantly conserved at position 14 (Fig. 2, Additional file 3). Additionally, Asp (D)/Glu (E) amino acids are conserved at positions 22 and 25. Several other amino acids are also conserved in the first EF-hands. However, the major focus was given to calcium sensing Asp (D) and Glu (E) amino acid. If we consider the presence of conserved domains in CBL proteins, there is a conserved V-F-H-P-N domain at the end of the first EF-hand (Fig. 2). In the second EF-hand, Asp/Glu amino acids are slightly conserved at the 3, 4 and 7 position, but Asp is significantly conserved at the 14 position (Fig. 2). The Glu amino acid is most significantly conserved at position 22 and is less significantly conserved at position 25. The Glu amino acid is also significantly conserved at position 36. In the third EF-hand, Asp amino acid is conserved at position 7, 8 and 14; while Glu is conserved at position 11, 19, 20, 21 and 22 (Fig. 2). The Asp and Glu amino acids are present as a D/E-x-D motif at position 20, 21 and 22 of the third EF-hand. Another motif, D-x-E-E, is present at position 30, 31, 32 and 33 in the third EF-hand. Taken together, these findings indicate that, the third EF-hand contains the maximum Asp and Glu amino acids within it. In EF-hand loop, the calcium ion is coordinated in a pentagonal bi-pyramidal configuration. Earlier study in CPK EF-hand revealed that, six amino acid residues are involved in binding of calcium ion in each EF-hands and are present at position 1, 3, 5, 7, 9 and 12 . These residues are denoted by X, Y, Z, −Y, −X and –Z. The invariant Glu or Asp amino acid at position 12 provides two molecules of oxygen for liganding Ca2+ (bidentate ligand) ion . The position 1 (X), 3 (Y) and 12 (−Z) are the most conserved and plays critical role in calcium binding. In case of CBLs, the presence of Asp or Glu amino acids at position 7, 14 and 22 are very critical for binding calcium ion while other conserved Asp and Glu amino acids might provides the accessory affinity sites for strong calcium binding.
There is a presence of an upstream region immediately adjacent to the first EF-hand of the CBL protein (Fig. 3). This up-stream region is not significantly conserved, but contain several calcium binding Asp and Glu amino acids (Additional file 3). The Group D CBL protein was found to contain conserved Asp and Glu at position 16, 17 and 18 (E-E/D-P) in the N-terminal region (Fig. 3a). In the group A CBL protein, there is a D/E-x-E/D motif present at up-stream of the first EF-hand (N-terminal region) (Fig. 3b). A less conserved domain E/D-D-P-E-X4-E-X6-E is present at the N-terminal region of the CBL protein (Additional file 3). In the C-terminal region, there is a conserved P-S-F-V-F-x-S-E-V-D-E domain present downstream of the third EF-hand (Fig. 4).
The organisms are able to recognize sense and respond to their environment to survive. In plants, sensing mechanisms are evolved in response to hormonal and environmental signals . To elicit a cellular response, the perceived signal must be conveyed to its cellular machinery. One of the most important secondary messengers, Ca2+, perceives the stimulus and transduces it to the downstream protein to initiate Ca2+ mediated responses. The Ca2+ mediated stimuli causes plant to respond to hormone and external stimuli, which mediate and regulate diverse fundamental cellular processes such as cell division, cell elongation, cell differentiation, cell polarity, photo morphogenesis, plant defense and stress responses . The CBL protein is one of the several calcium sensing protein families, including calcium dependent protein kinase (CPK) and calmodulins. The CPK protein contains a kinase domain as well as a regulatory domain that has four calcium sensing EF-hands. The acidic amino acids Asp (D) and Glu (E) present in the EF-hands are important calcium sensors . The CBL proteins lack the kinase domain and contain only three calcium binding EF-hands. The CBL proteins of Arabidopsis thaliana and Oryza sativa were previously reported to contain four calcium binding EF-hands [35–37]. However, the scan prosite software study revealed that, CBL proteins of all plants contain only three calcium binding EF-hand domains (Figs. 1 and 5) . Investigations of the CBL proteins of Kudla et al. , Batistic and Kudla  and Gu et al.  using the scan prosite software revealed that, all CBL proteins reported to have four EF-hands actually contained only three EF-hands. They reported that, in some cases CBL protein contains four EF-hands while in other they contain incomplete four EF-hands. The prosite analysis of data provided by Weinl and Kudla  shows that, O. tauri protein contains clear four EF-hands where as S. moellendorfii protein shows only three EF-hands. One CBL protein contains four EF-hands whereas other contains three EF-hands. This is very contradicting. This proves that, the data provided by Weinl and Kudla are contradictory. Some other data provided in this manuscript belongs to genus Physcomitrella patens (FJ901251, FJ901252, FJ901253 and FJ901254). Here the P. patens FJ901254 protein contains four EF-hands while other contains only three EF-hands. The CBL genes are present from single celled Chlamydomonas to the modern land plants. The Chlamydomonas is considered as the basal evolutionary lineage of photosynthetic green plant that evolved since 3500 million years ago, which is far earlier than the evolution of land plants. So, it is highly unlikely that genome(s) will encode for incomplete functional EF-hands for more than 3500 million years. Genomes are very specific in nature. They would either encode for complete EF-hand or would remove the incomplete one. But nothing has happened; because there is not presence of such incomplete EF-hands in CBLs. Evolutionary pressure cannot allow transfer of incomplete and non-functional EF-hand for millions of years. This proves that CBLs protein contain only three calcium-binding EF-hands, not four or incomplete four.
Although there have been significant advances in our understanding of CBL proteins, no studies are available regarding their conserved domains and motifs. In this study, we found that the calcium binding EF-hands are highly conserved and contains the E/D-x-D motif in the third EF-hand (Fig. 2). In addition to this motif, CBLs also contain several C-terminal downstream conserved motifs, specifically conserved Asp and Glu amino acids (Fig. 3a and b). The high proportion of Asp and Glu amino acids in CBLs provides an opportunity for the accommodation of Ca2+ ions.
Myristoylation and palmitoylation sites
Protein myristoylation and palmitoylation are two important events necessary for protein trafficking, stability and aggregation . Addition of myristic acid to N-terminal Gly amino acid leads to protein myristoylation, while addition of palmitic acid to N-terminal Cys amino acid leads to protein palmitoylation . In most of the studied CBLs, N-terminal Gly amino acid is required for protein myristoylation and is conserved at the second position (Fig. 3). The N-terminal Gly amino acid in some other CBL proteins has been found to be conserved at the seventh position. Similarly, N-terminal Cys amino acid is required for protein palmitoylation and is conserved at the third position in group D CBL proteins (Fig. 3a) and at the fourth position in group A CBL proteins (Fig. 3b). The majority of group B CBLs don’t contain N-terminal Cys amino acids.
The protein palmitoylation is a widespread modification found in membrane bound protein that includes transmembrane-spanning protein synthesized in soluble ribosome . In general, protein palmitoylation increases the affinity of protein for membrane attachment and therefore affects protein localization and function. Proteins that undergo palmitoylation include RasGTPase , Rho GTPase  and CDPKs . The RasGTPase, Rho GTPase, and CDPKs contain N-terminal Cys residues at either the third, fourth or fifth position . All the 24 Arabidopsis CPKs are predicted to have a myristoylation consensus sequence and contain at least one Cys residue either at fourth, fifth or sixth position . This study revealed the presence of an N-terminal Cys residue at the third, fourth, fifth or sixth position in several CBLs (Figs. 3 and 5a and b). Except for group B CBLs (CBL10), all other group of CBL proteins (group A, C and D) contain the N-terminal Cys residue. These finding clearly demonstrates that, group B CBL protein does not undergo protein palmitoylation, and only selective CBL protein posse’s protein palmitoylation activity.
Co-translational addition of myristate to N-terminal glycine amino acid through amide bonds is known as myristoylation . Except in group B CBLs, all other CBLs contain N-terminal glycine residues at the second position (Fig. 5a). Additionally, all CBLs (except group B CBLs) that contain N-terminal cysteine amino acid concurrently possess N-terminal Gly amino acid at the second position (Figs. 3 and 5). The N-terminal myristoylation promotes protein-membrane attachment and protein-protein interactions. Mutation in the N-terminal Gly-abolishes lipid modification and thus prevents membrane association . Twenty-four of the Arabidopsis calcium sensing CDPK proteins were predicted to have the N-terminal myristoylation motif for membrane association. Among them, AtCPK2 has been experimentally confirmed to be myristoylated at the N-terminal Gly residue, and the first ten amino acids of the CPK protein are critical for localization to the ER (endoplasmic reticulum) membrane . In majority of cases, N-terminal myristoylation and palmitoylation events are complement to each other. Both N-terminal myristoylation in the Gly amino acid at position 2 and palmitoylation in the Cys amino acid at position 4 and 5 have been validated experimentally in membrane bound OsCPK2 . When N-terminal myristoylation was abolished by mutation at the Gly amino acid, the protein could no longer be palmitoylated, indicating that N-terminal myristoylation is the prerequisite for palmitoylation. Only protein myristoylation provides a weak affinity for membrane attachment, whereas palmitoylation and myristoylation provide very high affinity interactions .
Phylogeny and evolution
Protein families are defined as groups of protein with more than 50 % pairwise amino acid sequence similarity . Molecular evolution is generally studied at the level of individual gene or families of genes . However, there are still no models that can infer gene family evolution to enable the estimation of the ancestral state. Phylogenetic analysis can be a powerful tool to infer the relationships among genes and analyze their evolutionary events . Phylogenetic analyses of all CBL genes together revealed that they fall into six different groups (Fig. 6). Some lower eukaryotic specific CBL genes such as SmCBL9, PpCBL3-1, PpCBL9 and PpCBL3-2 are present as a cluster (group F) at the center of the phylogenetic tree, while group E CBLs are present at the distal end of the phylogenetic tree. The cluster of other CBL genes of higher eukaryotic plants (group A, B, C and D) was directly linked with the cluster of group F CBL genes (Fig. 6, Additional file 5). These findings indicate that CBL gene families of higher eukaryotic plants are derived from common ancestors of lower eukaryotic plants (Fig. 6). The lower eukaryotic plants are very simple, with unicellular to multi-cellular architecture. As complexity of an organism increases, it need to adapt from simpler aquatic habitats to complex terrestrial habitats, and hence the number of CBL genes per genome got increased . This indicates that these CBL genes might have been evolved for some unique and specific function responsible for adaptation to complex lifestyles. The CBL genes of lower eukaryotic plants such as algae, Physcomitrella, Selaginella and Pinus are fall in group E and F. These genes are probably evolved independently during evolution. Some of the CBL genes (SmCBL9, PpCBL3-1, PpCBL9 and PpCBL3-2) of lower eukaryotic plants fall in the middle of the phylogenetic tree, while CBL genes of higher angiosperm plants are phylogenetically linked with the cluster of CBL genes of lower eukaryotic plants. These findings indicates that, CBL genes of modern plants may have derived from a common ancestor of lower eukaryotic plant . The phylogenetic analysis revealed that, CBL2, CBL3, CBL6 and CBL7 fall in group A, CBL10 falls in group B, CBL1 and CBL9 fall in group C, and CBL4, CBL5 and CBL8 fall in group D. The lower eukaryotic CBL genes from Selaginella (SmCBL5), Micromonas (MpCBL2, MpCBL6), Chlamydomonas (CreinCBL8, CreinCBL9), and CBL genes of Picea abies fall in group E and F.
The significant similarities between the CBL gene sequences indicate that they arose relatively recently via gene duplication and might have similar or overlapping functions. The paralogous genes evolved due to the development of new function and provided the most probable role for adaptation. Gene duplication and diversification are considered to be the most important events in evolutionary biology. If a gene is duplicated from its original gene, the selective constraints become much lower for the extra copy, and it can evolve to have a slightly different function while the original function of the gene is kept in the other copy. Hence, gene duplication with subsequent diversification is one of the simplest ways to acquire new function. Because the role of the CBL gene is important for calcium sensing and there are several other calcium sensing gene families (CPK, CaM, etc.) present in the plant kingdom, duplicated genes are still being found for CBL genes. This may be due to the ploidy level, as well as some other aspects in different genomes. Some plant genomes that have undergone duplication during evolution contain few duplicated CBL genes including Brassica rapa, Eucalyptus grandis, Glycine max, Gossipium raimondii and Medicago truncatula.
Tajima’s molecular test hypothesis explains the significance and rate of evolution . Random analysis of CBL sequences was carried out in Tajima’s relative rate test and the p-value and X 2 -test was found to be significant (Table 3). Three random replicate analyses were carried out. In each analysis, three sequences were considered for the study by making them as group A, B and C. The first analysis contained sequences of MdCBL3 (group A), CsCBL3 (group B), and PerCBL3 (group C); the second analysis contained MeCBL3 (group A), BrCBL2-2 (group B), and PtCBL3 (group C); and the third analysis contained FvCBL10-1 (group A), BrCBL2-2 (group B), and MgCBL5 (group C). In the statistical analysis, the p-value was found to be 0.00666, 0.00284 and 0.00555 for the first, second and third analysis, respectively (Table 3). Similarly, the chi-square values for the first, second and third analysis was found to be 7.36, 8.91 and 7.69, respectively, with one degree of freedom (Table 3). These findings suggest that, the results presented herein are statistically significant. In Tajima’s test for neutrality, Tajima’s D value for CBLs was found to be 4.413697 (Table 4). In Tajima’s D-test, when D = 0, the average heterozygosity of a population becomes equal to the number of segregating sites. This occurred because the expected variation is similar to the observed variations [55, 56]. Hence, evolution of the population can be due to mutation-drift equilibrium, and there is no evidence of selection. When D < 0, the average heterozygosity is lower than the number of segregating sites [55, 56]. This indicates that, rare alleles are present at very low frequency and recent selective sweeps led to the expansion of the population size after recent bottleneck. When D > 0, the average heterozygosity is more than the segregating sites and can be considered as the presence of multiple alleles at high frequencies [55, 56]. This leads to balanced selection due to the sudden contraction in population size. Tajima’s negative D value signifies a very low frequency of polymorphism relative to expectation, indicating expansion in population by size via selective sweep or purifying selection. Tajima’s positive D value signifies a high frequency of polymorphism, indicating a decrease in population size by balancing selection. A Tajima’s D value greater than 2 or less than −2 is considered significant [55, 56]. In this study, Tajima’s D value is 4.413697 (Table 4), signifying that CBL genes have undergone high frequencies of polymorphism via decreasing population size due to balanced selection. Accordingly, the heterozygosity of CBLs is greater than the number of segregating sites and present as multiple alleles.
Differential expression of OsCBL genes
The plant have become the important target for genetic manipulation and provided an excellent platform for the investigation of different biological processes that control development. Analysis of these developmental processes at the molecular level requires isolation and characterization of important regulatory genes, including those are differentially expressed. Genes expressed in different developmental stages and specific tissues are of great interest. One of the major interests is whether the specific expression pattern of a gene in a specific cell or tissue type at a specific developmental stage can be used as a marker to study the development. Therefore, we investigated the relative expression patterns of the OsCBL gene at different developmental stages (Fig. 7). The relative expression of OsCBL genes in leaf tissue shows that OsCBL3-1, OsCBL3-2, OsCBL3-3, OsCBL4-3, OsCBL9 and OsCBL10-2 were upregulated at all four time points (Fig. 7). The expression of OsCBL4-1 undergone down regulation at the third and fourth week, while OsCBL4-2 undergone down regulation at weeks 1, 3 and 4. The major changes in the expression of the OsCBL genes were observed at 3 and 4 week. To better understand the role of CBL genes in stress responses, we conducted differential expression analysis of OsCBL genes by subjecting them to cold and heat stress at different time points (Fig. 8). The relative expression of OsCBL3-1, OsCBL3-2, and OsCBL9 was increased at all time points, whereas the expression of OsCBL4-1, OsCBL4-2, OsCBL4-3, OsCBL10-1 and OsCBL10-2 was down regulated at 24 h (Fig. 8). In heat treated plants, OsCBL3-1, OsCBL4-2, OsCBL4-3, and OsCBL10-2 had undergone up-regulation at all four time points (Fig. 9). The expression of OsCBL3-2 was down regulated at all the four time points (Fig. 9). The expression of OsCBL3-3 was down regulated at 3 and 6 h, and then gradually up-regulated at 12 and 24 h. Similarly, expression of OsCBL9 was down regulated at 3 h, but was gradually upregulated at 6, 12 and 24 h. Based on these findings, CBL genes are cold and heat stress responsive and differentially expressed upon exposure to different stresses.
This study revealed that the basic architecture of CBL genes are conserved among all plant species, including green algae, bryophytes, pteridophytes, gymnosperms and angiosperms. The CBL genes of lower eukaryotes such as green algae and pinus appear to have evolved independently. Based on these findings, the split between chlorophyta (green algae) and embryophyta (higher plants) played an important role in the evolution of CBL genes. During the course of evolution, CBL signaling events by land plants expanded significantly via gene duplication. Expression analysis shows that OsCBL3-1, OsBL3-2, OsCBL3-3 and OsCBL10-2 significantly modulated during different developmental stages in O. sativa. The differential expression of OsCBL3-1 was significantly modulated during cold and heat stress suggesting its important roles during these events.
The calcineurin B-like (CBL) gene family of Arabidopsis thaliana and Oryza sativa was downloaded from the Arabidopsis Information Resource (TAIR) database and the TIGR rice genome annotation project, respectively. Identified protein sequences of A. thaliana and O. sativa were then used to identify CBL gene family members in other plant species in the phytozome and spruce genome database . The BLASTP program (default) was used to identify the CBL gene family members of other plant species. The default statistical parameters used in the BLASTP analysis were as follows: BLASTP-protein query of the protein database; expected threshold (E): −1, comparison matrix: BLOSUM62; no. of alignments to show: 100. Sequences having an E-value up to 7.0e-200 were taken into consideration for further analysis to cover the maximum number of genes. Collected protein sequences were then subjected to the scan prosite software to analyze the presence of EF-hand domains. Sequences having only three calcium binding EF-hands were as considered calcineurin B-like genes because all CBL proteins of A. thaliana contain only three calcium binding EF-hands. These CBL proteins were again subjected to BLASTP analysis against the Arabidopsis genome database using the default parameters to reconfirm them as CBL proteins. Sequences of plant species that gave BLASTP hits against the A. thaliana CBL proteins were considered as CBL proteins. The CBL genes were numbered according to the Arabidopsis CBL gene they matched during BLAST search to ensure proper orthologous numbering.
Multiple sequence alignment and molecular modeling of CBL
The multiple sequence alignment of CBL genes of all species was carried out using the multalin software (http://multalin.toulouse.inra.fr/multalin/). The statistical parameters used to run the programs were as follows: protein weight matrix: BLOSSUM62, gap penalty at opening: default, gap penalty at extension: default, gap penalty at extremities: none, one iteration only: no, high consensus value: 90 % and low consensus value: 50 %.
Construction of phylogenetic tree
To construct a phylogenetic tree, a clustal file was generated in the CLUSTALW software using the protein sequences of all CBL genes . The parameters used to run the CLUSTALW program were as follows: protein weight matrix: BLOSSUM62, gap open: 10, gap extension: 0.2, iteration: none. The generated clustal file was downloaded and converted to MEGA file format using the MEGA5 software . The generated MEGA file was then run in the MEGA5 software to construct the phylogenetic tree. The statistical parameters used to construct the phylogenetic tree were as follows: analysis: phylogenetic reconstruction, statistical method: maximum likelihood, test of phylogeny: bootstrap method, no. of bootstrap replicates: 3000, substitution type: amino acids, model/methods: Jones-Taylor-Thornton (JTT) model, rates among sites: uniform rates, gaps/missing data treatment: partial deletion, site coverage cutoff: 95 % and branch swap filter: very strong.
Tajima’s relative rate test was carried out to investigate the significance and rate of evolution of plant CBL genes. The generated MEGA file used for the construction of the phylogenetic tree was subjected to MEGA5 to analyze Tajima’s relative rate test and Tajima’s test of neutrality. The statistical parameters used to run Tajima’s relative rate test were as follows: Tajima’s relative rate test scope: for three chosen sequences, substitution type: amino acid, and gaps/missing data treatment: complete deletion. The statistical parameters used to carry out Tajima’s test of neutrality were as follows: analysis: Tajima’s neutrality test, scope: all selected taxa, substitution type: amino acids, and gaps/missing data treatment: complete deletion.
Plant treatment and quantitative real time PCR
Wild types Oryza sativa japonica var. nipponbare were grown in soil in a greenhouse under a16 h light: 8 h dark cycle at 22–25 °C for 20 days. Cold treatment consisted of 4 °C, while drought/heat treatment consisted of 40 °C. The leaves were sampled at 0, 3, 6, 12 and 24 h and immediately transferred to liquid nitrogen for subsequent analysis. Untreated plants were used as control samples. Three biological replicates were employed during this study. Total RNA was isolated from the treated and control samples using Trizol reagent. The RNA was quantified using Nanodrop1000 and its integrity was checked by electrophoresis in 1.5 % (w/v) agarose gel. High quality RNA was subjected to the preparation of cDNA using a Fermentas RevertAid first strand cDNA synthesis kit. The reactions were prepared by adding 1 μg total RNA, 2 μl of 10× RT buffer, 2 μl 10 mMdNTPs mix, 2 μl of oligo (dT)18 primer, 1 μl of reverse transcriptase, 1 μl ribolock RNase inhibitor and nuclease free sterile water up to 20 μl. The reaction mixtures were then subjected to thermal incubation at 42 °C for 60 min followed by reaction termination at 70 °C for 5 min. The generated cDNA was diluted 10 times and kept for further use. The primers of O. sativa CBL genes were designed using primer3 software targeting either the extreme 5′ end, which is not conserved, or the 3′ UTR region, which generated an amplicon size between 120 and 200 bp (primer length between 20 and 24 bp) (Additional file 4). The specificity of primers was checked through regular PCR amplification followed by agarose gel electrophoresis, as well as by the primer test in a Mx3000P quantitative real time PCR machine by examining the melting curve. The melting curve analysis of the primers was conducted at 60–95 °C, with a temperature increasing step of 0.06 °C/s (five acquisitions per degree of Celsius) at the end of each run. The quantitative real-time PCR was carried out using a Mx3000P real-time PCR system with SYBR green master mix (2x) (Fermentas) and ROX as a passive reference standard to normalize the SYBR fluorescent signal. The PCR amplification was carried out in a 25 μl reaction mixture containing 1 μl cDNA as the template, 12.5 μl SYBR green master mix (2X), 1 μl of each forward and reverse primer and nuclease free water up to 25 μl. The thermal profile for quantitative real time PCR was: initial activation at 95 °C for 10 min, followed by 40 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 30 s. Analyses were conducted in triplicate using three biological replicates. The primers showing efficiency of 90–105 % were considered as significant. The relative expression of OsCBL genes was calculated using 2-ΔΔCt method .
Elongation factor hand
Calcium dependent protein kinase
CBL interacting protein kinase
Hamel L-P, Sheen J, Séguin A. Ancient signals: comparative genomics of green plant CDPKs. Trends Plant Sci. 2014;19:79–89.
Tuteja N, Mahajan S. Calcium Signaling Network in Plants. Plant Signal Behav. 2007;2(2):79–85.
Hogan PG, Lewis RS, Rao A. Molecular basis of calcium signaling in lymphocytes: STIM and ORAI. Annu Rev Immunol. 2010;28:491–533.
Mohanta TK, Sinha AK. Role of Calcium-Dependent Protein Kinases during Abiotic Stress Tolerance. In: Tuteja N, Gill S, editors. Abiotic Stress Response Plants. 2015th ed. Weinheim, Germany: Wiley-VCH Verlag GmbH & Co; 2015. p. 185–208.
Webb AR. Focus issue: calcium signaling. Plant Physiol. 2013;163:457–8.
Eprintsev AT, Fedorin DN, Igamberdiev AU. Ca2+ is involved in phytochrome A-dependent regulation of the succinate dehydrogenase gene sdh1-2 in Arabidopsis. J Plant Physiol. 2013;170:1349–52.
Kanchiswamy CN, Mohanta TK, Capuzzo A, Occhipinti A, Verrillo F, Maffei ME, et al. Differential expression of CPKs and cytosolic Ca2+ variation in resistant and susceptible apple cultivars (Malus x domestica) in response to the pathogen Erwinia amylovora and mechanical wounding. BMC Genomics. 2013;14:760.
Day IS, Reddy VS, Shad Ali G, Reddy AN. Analysis of EF-hand-containing proteins in Arabidopsis. Genome Biol. 2002;3:RESEARCH0056.
Luan S, Kudla J, Rodriguez-concepcion M, Yalovsky S, Gruissem W. Calmodulins and Calcineurin B – like proteins : calcium sensors for specific signal response coupling in plants. Plant Cell. 2002;14:389–400.
Sulmann S, Dell’orco D, Marino V, Behnen P, Koch K-W. Conformational changes in calcium-sensor proteins under molecular crowding conditions. Chemistry (Easton). 2014;20:6756–62.
Hrabak E, Chan C, Gribskov M. The arabidopsis CDPK-SnRK superfamily of protein kinases. Plant Physiol. 2003;132:666–80.
Lewit-Bentley A, Réty S. EF-hand calcium-binding proteins. Curr Opin Struct Biol. 2000;10:637–43.
Popescu SC, Popescu GV, Bachan S, Zhang Z, Seay M, Gerstein M, et al. Differential binding of calmodulin-related proteins to their targets revealed through high-density Arabidopsis protein microarrays. Proc Natl Acad Sci U S A. 2007;104:4730–5.
Nozawa A, Koizumi N, Sano H. An Arabidopsis SNF1-related protein kinase, AtSR1, interacts with a calcium-binding protein, AtCBL2, of which transcripts respond to light. Plant Cell Physiol. 2001;42:976–81.
Cheong YH, Kim K, Pandey GK, Gupta R, Grant JJ, Luan S. CBL1, a calcium sensor that differentially regulates salt, drought, and cold responses in arabidopsis. Plant Cell. 2003;15(8):1833–45.
Schlicker A, Domingues FS, Rahnenführer J, Lengauer T. A new measure for functional similarity of gene products based on Gene Ontology. BMC Bioinformatics. 2006;7:302.
Aravind L, Mazumder R, Vasudevan S, Koonin EV. Trends in protein evolution inferred from sequence and structure analysis. Curr Opin Struct Biol. 2002;12:392–9.
Illergård K, Ardell DH, Elofsson A. Structure is three to ten times more conserved than sequence–a study of structural response in protein cores. Proteins. 2009;77:499–508.
Choi I-G, Kim S-H. Evolution of protein structural classes and protein sequence families. Proc Natl Acad Sci U S A. 2006;103:14056–61.
Hamel L-P, Nicole M-C, Sritubtim S, Morency M-J, Ellis M, Ehlting J, et al. Ancient signals: comparative genomics of plant MAPK and MAPKK gene families. Trends Plant Sci. 2006;11:192–8.
Mohanta TK, Arora PK, Mohanta N, Parida P, Bae H. Identification of new members of the MAPK gene family in plants shows diverse conserved domains and novel activation loop variants. BMC Genomics. 2015;16:58.
Tapan KM, Mickael M, Nibedita M, Chidananda NK. In-silico identification and phylogenetic analysis of auxin efflux carrier gene family in Setaria italica L. African J Biotechnol. 2014;13:211–25.
Mohanta TK, Mohanta N. Genome wide identification of auxin efflux carrier gene family in physcomitrella patens. J Biotechnol Sci. 2013;1:54–64.
Huynen M, Bork P. Measuring genome evolution. Proc Natl Acad Sci U S A. 1998;95:5849–56.
Kleist TJ, Spencley AL, Luan S. Comparative phylogenomics of the CBL-CIPK calcium-decoding network in the moss Physcomitrella, Arabidopsis, and other green lineages. Front Plant Sci. 2014;5(May):187.
Heslop-Harrison JSP, Schwarzacher T. Organisation of the plant genome in chromosomes. Plant J. 2011;66:18–33.
Nystedt B, Street NR, Wetterbom A, Zuccolo A, Lin Y-C, Scofield DG, et al. The Norway spruce genome sequence and conifer genome evolution. Nature. 2013;497:579–84.
Singer G, Hickey D. Nucleotide bias causes a genomewide bias in the amino acid composition of proteins. Mol Biol Evol. 2000;17:1581–8.
Foster PG, Jermiin LS, Hickey D a. Nucleotide composition bias affects amino acid content in proteins coded by animal mitochondria. J Mol Evol. 1997;44:282–8.
Akashi H, Gojobori T. Metabolic efficiency and amino acid composition in the proteomes of Escherichia coli and Bacillus subtilis. Proc Natl Acad Sci U S A. 2002;99:3695–700.
Akashi H. Translational selection and yeast proteome evolution. Genetics. 2003;164:1291–303.
Revett SP, King G, Shabanowitz J, Hunt DF, Hartman KL, Laue TM, et al. Characterization of a helix-loop-helix (EF hand) motif of silver hake parvalbumin isoform B. Protein Sci. 1997;6:2397–408.
Xu T, Dai N, Chen J, Nagawa S, Cao M, Li H, et al. Cell surface ABP1-TMK auxin-sensing complex activates ROP GTPase signaling. Science (80-). 2014;343:1025–8.
Poovaiah BW, Du L, Wang H, Yang T. Recent advances in calcium/calmodulin-mediated signaling with an emphasis on plant-microbe interactions. Plant Physiol. 2013;163:531–42.
Kudla J, Xu Q, Harter K, Gruissem W, Luan S. Genes for calcineurin B-like proteins in Arabidopsis are differentially regulated by stress signals. Proc Natl Acad Sci U S A. 1999;96:4718–23.
Luan S, Kudla J, Rodriguez-Concepcio M, Yalovsky S, Gruissem W. Calmodulins and calcineurin B-like proteins: calcium sensors for specific signal response coupling in plants.Plant Cell. 2002; S389-S400.
Gu Z, Ma B, Jiang Y, Chen Z, Su X, Zhang H. Expression analysis of the calcineurin B-like gene family in rice (Oryza sativa L.) under environmental stresses. Gene. 2008;415:1–12.
De Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, et al. ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006;34(Web Server issue):W362–5.
Batistic O, Kudla J. Plant calcineurin B-like proteins and their interacting protein kinases. Biochim Biophys Acta. 2009;1793:985–92.
Weinl S, Kudla J. The CBL – CIPK Ca 2+−decoding signaling network: function and perspectives. New Phytol. 2009;184:517–28.
Linder ME, Deschenes RJ. Palmitoylation: policing protein stability and traffic. Nat Rev Mol Cell Biol. 2007;8:74–84.
Smotrys JE, Linder ME. Palmitoylation of intracellular signaling proteins: regulation and function. Annu Rev Biochem. 2004;73:559–87.
Blaskovic S, Blanc M, van der Goot FG. What does S-palmitoylation do to membrane proteins? FEBS J. 2013;280:2766–74.
Singaraja RR, Hadano S, Metzler M, Givan S, Wellington CL, Warby S, et al. HIP14, a novel ankyrin domain-containing protein, links huntingtin to intracellular trafficking and endocytosis. Hum Mol Genet. 2002;11:2815–28.
Michaelson D, Silletti J, Murphy G, D’Eustachio P, Rush M, Philips MR. Differential localization of Rho GTPases in live cells: regulation by hypervariable regions and RhoGDI binding. J Cell Biol. 2001;152:111–26.
Roberts PJ, Mitin N, Keller PJ, Chenette EJ, Madigan JP, Currin RO, et al. Rho Family GTPase modification and dependence on CAAX motif-signaled posttranslational modification. J Biol Chem. 2008;283:25150–63.
Cheng S, Willmann M, Chen H, Sheen J. Calcium signaling through protein kinases. The Arabidopsis calcium-dependent protein kinase gene family. Plant Physiol. 2002;129:469–85.
Martín ML, Busconi L. A rice membrane-bound calcium-dependent protein kinase is activated in response to low temperature. Plant Physiol. 2001;125:1442–9.
Lu SX, Hrabak EM. An Arabidopsis calcium-dependent protein kinase is associated with the endoplasmic reticulum 1. Plant Physiol. 2002;128:1008–21.
Huang IK, Pei J, Grishin NV. Defining and predicting structurally conserved regions in protein superfamilies. Bioinformatics. 2013;29:175–81.
Peña MI, Van Itallie E, Bennett MR, Shamoo Y. Evolution of a single gene highlights the complexity underlying molecular descriptions of fitness. Chaos. 2010;20:026107.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Prochnik SE, Umen J, Nedelcu AM, Hallmann A, Miller SM, Nishii I, et al. Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science (80-). 2010;329:223–6.
Galperin MY, Koonin EV. Divergence and convergence in enzyme evolution. J Biol Chem. 2012;287:21–8.
Tajima F. Simple methods for testing the molecular evolutionary clock hypothesis. Genetics. 1993;135:599–607.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;595:585–95.
Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, et al. Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012;40:D1178–86.
Larkin M, Blackshields G, Brown NP, Chenna R, McGettigan P, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.
This research was supported by the Yeungnam University Research Grant 214A367010.
The authors declare that there is no conflict of interests regarding the publication of this paper.
TKM: Conception and design of the experiment carried out the experiments, analyzed and interpreted the data, drafted the manuscript, NM: Analyzed the data, drafted manuscript, YM: carried out the experiments, PP: analyzed the data, HB: given approval for publication. All authors read and approved the final manuscript.
Table representing detailed genomic information of different CBL genes from 38 different plant species identified during this study.
Table representing molecular mass (kDa) and isoelectric point of different CBL genes from 38 plant species identified during this study.
Multiple sequence alignment of all CBL genes analyzed during this study. Sequence alignment was done using online available Multalin software (http://multalin.toulouse.inra.fr/multalin/ multalin.html) using default programme. Multiple alignments show presence of different conserved domains and motifs in CBL genes.
Table showing lists of OsCBL primers used during qRT-PCR analysis. Primers were designed using primer3 software (http://primer3.ut.ee/).
The details of phylogenetic datas are submitted to TreeBASE database and can be available in following link http://purl.org/phylo/treebase/phylows/study/TB2:S17414?x-access-code=1b88565e08ce238f8fc7928d2fa11a12&format=html .