Molecular Evolution of Calmodulin and Calmodulin-Like Genes in the Cephalochordate Branchiostoma

Calmodulin is a calcium-binding EF-hand protein that is an activator of many enzymes as well as ion pumps and channels. Due to its multiple targets and its central role in the cell, understanding the evolutionary history of calmodulin genes should provide insights into the origin of genetic complexity in eukaryotes. We have previously isolated and characterized a calmodulin gene from the early-diverging chordate Branchiostoma lanceolatum (CaM1). In this paper, we report the existence of a second calmodulin gene (CaM2) as well as two CaM-like genomic fragments (CaML-2, CaML-3) in B. lanceolatum and a CaM2 and three CaM-like genes (CaML-1, CaML-2, CaML-3) in B. floridae. The CaM-like genes were isolated using low-stringency PCR. Surprisingly, the nucleotide sequences of the B. lanceolatum CaM1 and CaM2 cDNAs differ by 19.3%. Moreover, the CaM2 protein differs at two positions from the amino acid sequence of CaM1; the latter is identical to calmodulins in Drosophila melanogaster, the mollusc Aplysia californica, and the tunicate Halocynthia roretzi. The two B. lanceolatum CaM-like genes are more closely related to the CaM2 than to the CaM1 gene. This relationship is supported by the phylogenetic analyses and the identical exon/intron organization of these three genes, a relationship unique among animal CaM sequences. These data demonstrate the existence of a CaM multigene family in the cephalochordate Branchiostoma, which may have evolved independently from the multigene family in vertebrates. Received: 2 November 1999 / Accepted: 25 April 2000     

Analysis of <Emphasis Type="Italic">Acropora muricata</Emphasis> Calmodulin (<Emphasis Type="Italic">CaM</Emphasis>) Indicates That Scleractinian Corals Possess the Ancestral Exon/Intron Organization of the Eumetazoan <Emphasis Type="Italic">CaM</Emphasis> Gene

Calmodulin (CaM), belonging to the tropinin C (TnC) superfamily, is one of the calcium-binding proteins that are highly conserved in their protein and gene structure. Based on the structure comparison among published vertebrate and invertebrate CaM, it is proposed that the ancestral form of eumetazoan CaM genes should have five exons and four introns (four-intron hypothesis). In this study, we determined the gene structure of CaM in the coral Acropora muricata, an anthozoan cnidarian representing the basal position in animal evolution. A CaM clone was isolated from a cDNA library constructed from the spawned eggs of A. muricata. This clone was composed of 908 nucleotides, including 162 base pairs (bp) of 5′-untranslated region (UTR), 296 bp of 3′-UTR, and an open reading frame 450 bp in length. The deduced amino acid indicated that the Acropora CaM protein is identical to that of the actiniarian, Metridinium senile, and has four putative calcium-binding domains highly similar to those of other vertebrate or invertebrate CaMs. Southern blot analysis revealed that Acropora CaM is a putative single-copy gene in the nuclear genome. Genomic sequencing showed that Acropora CaM was composed of five exons and four introns, with intron II not corresponding to any region in the actiniarian CaM gene, which possesses only four exons and three introns. Our results highlight that the coral CaM gene isolated from A. muricata has four introns at the predicted positions of the early metazoan CaM gene organization, providing the first evidence from the basal eumetazoan phylum to support the four-intron hypothesis.     

Characterization and Ca2+-induced expression of calmodulin (CaM) in marine dinoflagellates

Calmodulin (CaM) is one of the major Ca²⁺-binding proteins in the cells, and it plays multiple roles in several Ca²⁺ signaling pathways and regulating the activities of other proteins. In the present study, we characterized CaM genes from the marine dinoflagellates Amphidinium carterae, Cochlodinium polykrikoides, Prorocentrum micans, and P. minimum, and examined their expression patterns upon the addition and chelation of calcium. Their cDNAs had same ORF length (450 bp) and encoded the same protein, but with few nucleotide differences in the ORF and different 3′- and 5′ untranslated regions (UTRs). The four CaM proteins consist of four EF-hand Ca²⁺-binding motifs, two N-terminal domains and two C-terminal domains, and they were highly conserved within eukaryotes. The CaM gene expressions in the tested species increased by calcium treatments; however, they were significantly down-regulated by the calcium-chelator EGTA. The CaM genes of the test species were inducible and regulated by different calcium doses, suggesting their major role in calcium regulation in dinoflagellates.     

The cephalochordate Branchiostoma genome contains 26 intermediate filament (IF) genes: Implications for evolution of chordate IF proteins

We analyzed the draft genome of the cephalochordate Branchiostoma floridae (B. floridae) for genes encoding intermediate filament (IF) proteins. From 26 identified IF genes 13 were not reported before. Four of the new IF genes belong to the previously established Branchiostoma IF group A, four to the Branchiostoma IF group B, one is homologous to the type II keratin E2 while the remaining four new IF sequences N1 to N4 could not be readily classified in any of the previously established Branchiostoma IF groups. All eleven identified A and B2-type IF genes are located on the same genomic scaffold and arose due to multiple cephalochordate-specific duplications. Another IF gene cluster, identified in the B. floridae genome, contains three keratins (E1, Y1, D1), two keratin-like IF genes (C2, X1), one new IF gene (N1) and one IF unrelated gene, but does not show any similarities to the well defined vertebrate type I or type II keratin gene clusters. In addition, some type III sequence features were documented in the new IF protein N2, which, however, seems to share a common ancestry with the Branchiostoma keratins D1 and two keratin-related genes C. Thus, a few type I and type II keratin genes existed in a common ancestor of cephalochordates and vertebrates, which after separation of these two lineages gave rise to the known complexities of the vertebrate cytoplasmic type I–IV IF proteins, as well as to the multiple keratin and related IF genes in cephalochordates, due to multiple gene duplications, deletions and sequence divergences.     

Structure and expression of the chicken calmodulin I gene

The chicken calmodulin I (CaMI) gene has been isolated and characterized on the level of cDNA and genomic DNA. The deduced amino acid (aa) sequence is identical to the one of chicken CaMII which consists of 148 aa. The CaMI gene contains six exons. Its intron/exon organization is identical to that of the chicken CaMII and the CaMI and CaMIII genes of rat and human. Expression of the CaMI gene was detected in all chicken tissues examined, although at varying levels. The gene is transcribed into four mRNAs of 0.8, 1.4, 1.7 and 4.4 kb as determined by Northern blot analysis. Our results demonstrate that the “multigene-one-protein” principle of CaM synthesis is not only applicable to mammals whose CaM is encoded by three different genes, but also to chickens.     

Backbone resonance assignments of complexes of apo human calmodulin bound to IQ motif peptides of voltage-dependent sodium channels Na<Subscript>V</Subscript>1.1, Na<Subscript>V</Subscript>1.4 and Na<Subscript>V</Subscript>1.7

Human voltage-gated sodium (Na_V) channels are critical for initiating and propagating action potentials in excitable cells. Nine isoforms have different roles but similar topologies, with a pore-forming α-subunit and auxiliary transmembrane β-subunits. Na_V pathologies lead to debilitating conditions including epilepsy, chronic pain, cardiac arrhythmias, and skeletal muscle paralysis. The ubiquitous calcium sensor calmodulin (CaM) binds to an IQ motif in the C-terminal tail of the α-subunit of all Na_V isoforms, and contributes to calcium-dependent pore-gating in some channels. Previous structural studies of calcium-free (apo) CaM bound to the IQ motifs of Na_V1.2, Na_V1.5, and Na_V1.6 showed that CaM binding was mediated by the C-domain of CaM (CaM_C), while the N-domain (CaM_N) made no detectable contacts. To determine whether this domain-specific recognition mechanism is conserved in other Na_V isoforms, we used solution NMR spectroscopy to assign the backbone resonances of complexes of apo CaM bound to peptides of IQ motifs of Na_V1.1, Na_V1.4, and Na_V1.7. Analysis of chemical shift differences showed that peptide binding only perturbed resonances in CaM_C; resonances of CaM_N were identical to free CaM. Thus, CaM_C residues contribute to the interface with the IQ motif, while CaM_N is available to interact elsewhere on the channel.     

Identification and expression of an encoding steroid receptor coactivator (SRA) in amphioxus (Branchiostoma japonicum)

Steroid receptor coactivator (SRA), a class of genes encoding both functional RNA and protein, has been shown to be present in vertebrates but little is known in invertebrates. Here we isolated a cDNA encoding a SRA homolog from amphioxus Branchiostoma japonicum, named AmphiSRA. The cDNA contained a 525 bp open reading frame corresponding to a deduced protein of 174 amino acids with a predicted molecular mass of ~21 kDa. Phylogenetic analysis showed that AmphiSRA was located at the base of its vertebrate counterparts, suggesting that it represents the archetype of vertebrate SRA. The genomic DNA sequence of AmphiSRA contained four exons and three introns, which was similar to B. floridae SRA exon/intron organization. The recombinant SRAP expressed in vitro shows a band with a molecular mass of 21 kDa and western blot confirmed it, which proved it is an encoding isoform. AmphiSRA is found to display a tissue specific expression pattern, with a predominant expression in gill, intestine, testis, neural tube and notochord. The whole-mount in situ hybridization demonstrated the expression of AmphiSRA in all the stages of development assayed. These implicated that SRA maybe play an important role during embryonic development of cephalochordate amphioxus.     

Overexpression of calmodulin gene fragment from Antarctic notothenioid fish improves chilling tolerance in <Emphasis Type="Italic">Nicotiana benthamiana</Emphasis>

Calmodulin (CaM) is a highly conserved calcium sensor protein associated with chilling tolerance in living organisms. It has four EF-hand domains for binding of four Ca²⁺, two of them located in the N-terminus, and the other two in the C-terminus. A notothenioid CaM gene fragment (CaMm), which only codes for N-terminus of CaM (with two EF-hand domains), was introduced into Nicotiana benthamiana. Effects of its overexpression on chilling tolerance in plants were explored. During 4?C or 0?C chilling treatment, both CaMm and CaM transgenic plants showed higher PSII maximum quantum yield, actual quantum yield, and soluble protein content, lower electrolyte leakage and malondialdehyde content than that of the control. The changes in these physiological indices were comparable between the CaMm and CaM transgenic plants during the treatments. These results indicate that the N-terminus of calmodulin is likely the key functional domain involved in the adaptive response to cold stress.     

Comparing the Calcium Binding Abilities of Two Soybean Calmodulins: Towards Understanding the Divergent Nature of Plant Calmodulins

The discovery that plants contain multiple calmodulin (CaM) isoforms of variable sequence identity to animal CaM suggested an additional level of sophistication in the intracellular role of calcium regulation in plants. Past research has focused on the ability of conserved or divergent plant CaM isoforms to activate both mammalian and plant protein targets. At present, however, not much is known about how these isoforms respond to the signal of an increased cytosolic calcium concentration. Here, using isothermal titration calorimetry and NMR spectroscopy, we investigated the calcium binding properties of a conserved (CaM1) and a divergent (CaM4) CaM isoform from soybean (Glycine max). Both isoforms bind calcium with a semisequential pathway that favors the calcium binding EF-hands of the C-terminal lobe over those of the N-terminal lobe. From the measured dissociation constants, CaM4 binds calcium with a threefold greater affinity than CaM1 (K_d,Ca,mean of 5.0 versus 14.9 μM) but has a significantly reduced selectivity against the chemically similar magnesium cation that binds preferentially to EF-hand I of both isoforms. The implications of a potential magnesium/calcium competition on the activation of CaM1 and CaM4 are discussed in context with their ability to respond to stimulus-specific calcium signatures and their known physiological roles.     

Life-threatening arrhythmogenic CaM mutations disrupt CaM binding to a distinct RyR2 CaM-binding pocket

Calmodulin (CaM) modulates the activity of several proteins that play a key role in excitation-contraction coupling (ECC). In cardiac muscle, the major binding partner of CaM is the type-2 ryanodine receptor (RyR2) and altered CaM binding contributes to defects in sarcoplasmic reticulum (SR) calcium (Ca²⁺) release. Many genetic studies have reported a series of CaM missense mutations in patients with a history of severe arrhythmogenic cardiac disorders. In the present study, we generated four missense CaM mutants (CaM^N98I, CaM^D132E, CaM^D134H and CaM^Q136P) and we used a CaM-RyR2 co-immunoprecipitation and a [³H]ryanodine binding assay to directly compare the relative RyR2-binding of wild type and mutant CaM proteins and to investigate the functional effects of these CaM mutations on RyR2 activity. Furthermore, isothermal titration calorimetry (ITC) experiments were performed to investigate and compare the interactions of the wild-type and mutant CaM proteins with various synthetic peptides located in the well-established RyR2 CaM-binding region (3584-3602aa), as well as another CaM-binding region (4255-4271aa) of human RyR2. Our data revealed that all four CaM mutants displayed dramatically reduced RyR2 interaction and defective modulation of [³H]ryanodine binding to RyR2, regardless of LQTS or CPVT association. Moreover, our isothermal titration calorimetry ITC data suggest that RyR2 3584-3602aa and 4255-4271aa regions interact with significant affinity with wild-type CaM, in the presence and absence of Ca²⁺, two regions that might contribute to a putative intra-subunit CaM-binding pocket. In contrast, screening the interaction of the four arrhythmogenic CaM mutants with two synthetic peptides that correspond to these RyR2 regions, revealed disparate binding properties and signifying differential mechanisms that contribute to reduced RyR2 association.     

Calmodulin-Mediated Signal Transduction Pathways in Arabidopsis Are Fine-Tuned by Methylation

Calmodulin N-methyltransferase (CaM KMT) is an evolutionarily conserved enzyme in eukaryotes that transfers three methyl groups to a highly conserved lysyl residue at position 115 in calmodulin (CaM). We sought to elucidate whether the methylation status of CaM plays a role in CaM-mediated signaling pathways by gene expression analyses of CaM KMT and phenotypic characterization of Arabidopsis thaliana lines wherein CaM KMT was overexpressed (OX), partially silenced, or knocked out. CaM KMT was expressed in discreet spatial and tissue-specific patterns, most notably in root tips, floral buds, stamens, apical meristems, and germinating seeds. Analysis of transgenic plants with genetic dysfunction in CaM KMT revealed a link between the methylation status of CaM and root length. Plants with suppressed CaM methylation had longer roots and CaM KMT OX lines had shorter roots than wild type (Columbia-0). CaM KMT was also found to influence the root radial developmental program. Protein microarray analyses revealed a number of proteins with specificity for methylated forms of CaM, providing candidate functional intermediates between the observed phenotypes and the target pathways. This work demonstrates that the functionality of the large CaM family in plants is fine-tuned by an overarching methylation mechanism.     

Both N- and C-lobes of calmodulin are required for Ca-dependent regulations of CaV1.2 Ca channels

We investigated the concentration- and Ca²⁺-dependent effects of CaM mutants, CaM₁₂ and CaM₃₄, in which Ca²⁺-binding to its N- and C-lobes was eliminated, respectively, on the Ca_V1.2 Ca²⁺ channel by inside-out patch clamp in guinea-pig cardiomyocytes. Both CaM₁₂ and CaM₃₄ (0.7-10 μM) applied with 3 mM ATP produced channel activity after “rundown”. Concentration-response curves were bell-shaped, similar to that for wild-type CaM. However, there was no obvious leftward shift of the curves by increasing [Ca²⁺], suggesting that both functional lobes of CaM were necessary for the Ca²⁺-dependent shift. However, channel activity induced by the CaM mutants showed Ca²⁺-dependent decrease, implying a Ca²⁺ sensor existing besides CaM. These results suggest that both N- and C-lobes of CaM are required for the Ca²⁺-dependent regulations of Ca_V1.2 Ca²⁺ channels.     

Calmodulin isoforms in Arabidopsis encoded by multiple divergent mRNAs

Three new, unique cDNA sequences encoding isoforms of calmodulin (CaM) were isolated from an Arabidopsis cDNA library cloned in gt10. These sequences (ACaM-4, -5, and -6) represent members of the Arabidopsis CaM gene family distinct from the three DNA sequences previously reported. ACaM-4 and -6 encode full-length copies of CaM mRNAs of ca. 0.75 kb. The ACaM-5 sequence encodes a partial length copy of CaM mRNA that is lacking sequences encoding the amino-terminal 10 amino acids of mature CaM and the initiator methionine. The derived amino acid sequence of ACaM-5 is identical to the sequences encoded by two of the previously characterized ACaM cDNAs, and is identical to TCH-1 mRNA, whose accumulation was increased by touch stimulation. The polypeptides encoded by ACaM-4 and -6 differ from that encoded by ACaM-5 by six and two amino acid substititions, respectively. Most of the deduced amino acid sequence substitutions in the Arabidopsis CaM isoforms occurred in the fourth Ca²⁺-binding domain. Polymerase chain reaction amplification assays of ACaM-4, -5 and -6 mRNA sequences indicated that each accumulated in Arabidopsis leaf RNA fractions, but only ACaM-4 and -5 mRNAs were detected in silique total RNA. The six different CaM cDNA sequences each hybridize with unique Eco RI restriction fragments in genomic Southern blots of Arabidopsis DNA, indicating that these sequences were derived from distinct structural genes. Our results suggest that CaM isoforms in Arabidopsis may have evolved to optimize the interaction of this Ca²⁺-receptor protein with specific subsets of response elements.     

Characterization of GmCaMK1, a member of a soybean calmodulin-binding receptor-like kinase family

Calmodulin(CaM)-regulated protein phosphorylation forms an important component of Ca²⁺ signaling in animals but is less understood in plants. We have identified a CaM-binding receptor-like kinase from soybean nodules, GmCaMK1, a homolog of Arabidopsis CRLK1. We delineated the CaM-binding domain (CaMBD) of GmCaMK1 to a 24-residue region near the C-terminus, which overlaps with the kinase domain. We have demonstrated that GmCaMK1 binds CaM with high affinity in a Ca²⁺-dependent manner. We showed that GmCaMK1 is expressed broadly across tissues and is enriched in roots and developing nodules. Finally, we examined the CaMBDs of the five-member GmCaMK family in soybean, and orthologs present across taxa.

Structured summary

MINT-8051564: AtCRLK2 (uniprotkb:Q9LFV3) binds (MI:0407) to CaM (uniprotkb:P62199) by filter binding (MI:0049)MINT-8051416: GmCaMK3 (uniprotkb:C6ZRS6) binds (MI:0407) to CaM (uniprotkb:P62199) by filter binding (MI:0049)MINT-8051258: CaM (uniprotkb:P62199) and GmCaMK1 (genbank_protein_gi:223452504) bind (MI:0407) by isothermal titration calorimetry (MI:0065)MINT-8051400: GmCaMK2 (uniprotkb: C6ZRY5) binds (MI:0407) to CaM (uniprotkb:P62199) by filter binding (MI:0049)MINT-8051242, MINT-8051295, MINT-8051313, MINT-8051327, MINT-8051341, MINT-8051355: GmCaMK1 (genbank_protein_gi:223452504) binds (MI:0407) to CaM (uniprotkb:P62199) by filter binding (MI:0049)MINT-8051467: GmCaMK4 (uniprotkb: C6TIQ0) binds (MI:0407) to CaM (uniprotkb:P62199) by filter binding (MI:0049)MINT-8051276: CaM (uniprotkb:P62199) and GmCaMK1 (genbank_protein_gi:223452504) bind (MI:0407) by comigration in non denaturing gel electrophoresis (MI:0404)MINT-8051374: CaM (uniprotkb:P62199) and GmCaMK1 (genbank_protein_gi:223452504) bind (MI:0407) by mass spectrometry studies of complexes (MI:0069)     

Focus Issue on Calcium Signaling: Recent Advances in Calcium/Calmodulin-Mediated Signaling with an Emphasis on Plant-Microbe Interactions

Calcium/calmodulin-mediated signaling contributes in diverse roles in plant growth, development, and response to environmental stimuli.During calcium (Ca2+) signaling, decoding the stimulus-response coupling involves a set of Ca2+ sensor proteins or Ca2+-binding proteins (DeFalco et al., 2010a; Kudla et al., 2010). These proteins usually possess one or more classical helix-loop-helix elongation factor (EF) hand motifs. Three major types of Ca2+-sensor proteins in plants are calmodulin (CaM)/CaM-like proteins, calcium-dependent protein kinases (CDPKs), and calcineurin B-like proteins. As compared with animals, plant genomes encode more diversified Ca2+ sensors; with the exception of canonic CaM, all other types of Ca2+ sensors (CaM-like proteins, CDPKs, and calcineurin B-like proteins) are plant specific. The large population and unique structural composition of Ca2+-binding proteins and the diversity of the target proteins regulated by the Ca2+ sensors reflect the complexity of Ca2+ signaling, which helps plants adapt to the changing environment. This update will be limited primarily to discussions on CaM and CaM-binding proteins and the recent advances in Ca2+/CaM-mediated signaling.CaM is a conserved Ca2+-binding protein found in all eukaryotes. The discovery of CaM can be traced back to the 1970s. An activator of cyclic nucleotide phosphodiesterase was shown to be involved in the regulation of cAMP concentration, which was stimulated by Ca2+ (Kakiuchi and Yamazaki, 1970; Cheung, 1971). The activator was found to bind Ca2+ and was eventually named “calmodulin,” an abbreviation of Ca2+-modulated protein. Since its discovery over 40 years ago, CaM has been regarded as a model Ca2+-binding protein and has been subjected to intensive studies in biochemistry, cell biology, and molecular biology because of its importance in almost all aspects of cellular regulation (Poovaiah and Reddy, 1987, 1993; Bouche et al., 2005; DeFalco et al., 2010a; Du et al., 2011; Reddy et al., 2011b). Disruption or depletion of the single copy of the CaM gene in yeast (Saccharomyces cerevisiae) results in a recessive lethal mutation (Davis et al., 1986), suggesting that CaM has a critical role in eukaryotic cells.The structure of CaM has been well studied, and the prototype of CaM found in all eukaryotes has 149 amino acids with two globular domains, each containing two EF hands connected by a long flexible helix (Meador et al., 1993; Zhang et al., 1995; Yun et al., 2004; Ishida et al., 2009). As more and more genomes are sequenced, it is becoming clear that CaM belongs to a small gene family in plants. In the model plant Arabidopsis (Arabidopsis thaliana), seven CaM genes encode for four highly conserved isoforms (CaM1/4, CaM2/3/5, CaM6, and CaM7) that differ in only one to five amino acid residues. Loss-of-function mutations of individual CaMs indicate that the different CaMs may have overlapping yet different functions. For example, a loss of function in Arabidopsis AtCaM2 affects pollen germination (Landoni et al., 2010). Phenotypic analysis showed that in normal growth conditions, atcam2-2 plants were indistinguishable from the wild type, while genetic analysis showed a reduced transmission of the atcam2-2 allele through the male gametophyte, and in vitro pollen germination revealed a reduced level of germination in comparison with the wild type. However, the atcam3 knockout mutant showed a clear reduction in thermotolerance after heat treatment at 45°C for 50 min (Zhang et al., 2009). Overexpression of AtCaM3 in either the atcam3 knockout or wild-type background significantly rescued or increased the thermotolerance, respectively. Further analysis of individual CaM mutants under different stress conditions should reveal more on the functional significance of individual CaM genes.     

Bacterial expression and characterization of proteins derived from the chicken calmodulin cDNA and a calmodulin processed gene

Both normal chicken calmodulin (CaM) and a CaM-like mutant protein have been expressed in bacteria, isolated and evaluated with respect to several physical and biological properties. The mutant CaM is derived from a CaM-like gene that lacks intervening sequences and probably evolved from a CaM-processed gene (Stein, J. P., Munjaal, R. P., Lagacé, L., Lai, E. C., O'Malley, B. W., and Means, A. R. (1983) Proc. Natl. Acad. Sci. U. S. A. 80, 6485-6489). The mutant CaM protein contains 16 of the 19 amino acids encoded by the CaM-like gene. Normal chicken CaM produced in bacteria is identical to rat CaM by all criteria tested except that it is not trimethylated. The protein product of the CaM-like gene has been termed CaML and exhibits properties which are very similar to CaM despite the presence of 16 amino acid substitutions. CaML binds Ca2+ as evidenced by Ca2+-dependent binding to phenothiazine- and phenyl-Sepharose affinity resins and a Ca2+-dependent electrophoretic mobility shift which is similar to but distinct from CaM. CaML cross-reacts with a monospecific CaM antibody and has an immunodilution curve which is identical to bacterially synthesized CaM. Finally, CaML can maximally activate rat brain phosphodiesterase but with altered kinetic parameters as compared to CaM. These data suggest that the nucleotide substitutions in the putative CaM processed gene are not random but are selected to retain CaM-like functions in the encoded protein. Such a mechanism may exist for other processed genes.