Genome-Wide Identification and Analyses of Calmodulins and Calmodulin-like Proteins in Lotus japonicas

Evolution and diversification of CaM/CML gene family in green plants

Calmodulin (CaM) and calmodulin-like (CML) proteins are crucial Ca2+ sensors, which are widely involved in different biological processes of plants, including their growth and development, and stress responses. However, the origin and evolution of the CaM/CML gene family in plants remain elusive. In this study, 2133 CaM and 23094 CML genes were identified from the 1000 plants project (1 KP) species and the sequenced plants, covering algae, mosses, monilophytes, lycophytes, flowering plants, and all other green plant branches. Analysis showed that the size of the CML subfamily was correlated with the genome size of corresponding plant species, as well as the total gene number in the genome. Moreover, with the evolution from algae to angiosperms, the number of CML genes in plants increased gradually which could have been driven mainly by genome-wide segmental duplication events, while the number of CaMs remained basically stable at 2-3. Phylogenetic analysis demonstrated that CaM first appeared in green algae, while CML appeared earlier and has already been presented in dinoflagellates. Further analysis showed that the number and sequence of EF-hand domain in CaMs are highly conserved, while those of CMLs are diverse among different plant taxa. Expression analysis revealed that the expression level of CaMs was generally higher than that of CMLs, indicating that the high-expression genes have essential functions, while the low-expression genes are the main reasons for the functional diversity of the CaM/CML gene family in plants. The results might contribute to understanding the evolution of CaM/CML genes as well as their molecular functions.

Genome-Wide Characterization of Calmodulin and Calmodulin-like Protein Gene Families in Paulownia fortunei and Identification of Their Potential Involvement in Paulownia Witches' Broom

As significant Ca2+ sensors, calmodulin (CaM) and calmodulin-like proteins (CML), have been associated with a variety of environmental conditions in plants. However, whether CaMs/CMLs are related to the stress of phytoplasma infection has not been reported in Paulownia fortunei. In the current study, 5 PfCaMs and 58 PfCMLs were detected through a genome-wide investigation. The number of EF-hand motifs in all PfCaMs/CMLs varied. Bioinformatics analyses, including protein characteristics, conserved domain, gene structure, cis-elements, evolutionary relationship, collinearity, chromosomal location, post-translation modification site, subcellular localization and expression pattern analyses, represented the conservation and divergence of PfCaMs/CMLs. Furthermore, some PfCaMs/CMLs might be involved in plants' reaction to phytoplasma infection and exogenous calcium therapy, indicating these genes may play a role in abiotic as well as biotic stress responses. In addition, subcellular localization analysis showed that PfCML10 was located in the cell membrane and nucleus. In summary, these findings establish a stronger platform for their subsequent functional investigation in trees and further characterize their roles in Paulownia witches' broom (PaWB) occurrence.

Calcium decoders and their targets: The holy alliance that regulate cellular responses in stress signaling

Calcium (Ca²⁺) signaling is versatile communication network in the cell. Stimuli perceived by cells are transposed through Ca²⁺-signature, and are decoded by plethora of Ca²⁺ sensors present in the cell. Calmodulin, calmodulin-like proteins, Ca²⁺-dependent protein kinases and calcineurin B-like proteins are major classes of proteins that decode the Ca²⁺ signature and serve in the propagation of signals to different parts of cells by targeting downstream proteins. These decoders and their targets work together to elicit responses against diverse stress stimuli. Over a period of time, significant attempts have been made to characterize as well as summarize elements of this signaling machinery. We begin with a structural overview and amalgamate the newly identified Ca²⁺ sensor protein in plants. Their ability to bind Ca²⁺, undergo conformational changes, and how it facilitates binding to a wide variety of targets is further embedded. Subsequently, we summarize the recent progress made on the functional characterization of Ca²⁺ sensing machinery and in particular their target proteins in stress signaling. We have focused on the physiological role of Ca²⁺, the Ca²⁺ sensing machinery, and the mode of regulation on their target proteins during plant stress adaptation. Additionally, we also discuss the role of these decoders and their mode of regulation on the target proteins during abiotic, hormone signaling and biotic stress responses in plants. Finally, here, we have enumerated the limitations and challenges in the Ca²⁺ signaling. This article will greatly enable in understanding the current picture of plant response and adaptation during diverse stimuli through the lens of Ca²⁺ signaling.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in wheat (Triticum aestivum L.)

Calmodulin (CaM) and calmodulin-like (CML) genes are widely involved in plant growth and development and mediating plant stress tolerance. However, the whole genome scale studies about CaM and CML gene families have not been done in wheat, and the possible functions of most wheat CaM/CML gene members are still unknown. In this study, a total of 18 TaCaM and 230 TaCML gene members were identified in wheat genome. Among these genes, 28 TaCaM/CML gene members have 74 duplicated copies, while 21 genes have 48 transcript variants, resulting in 321 putative TaCaM/CML transcripts totally. Phylogenetic tree analysis showed that they can be classified into 7 subfamilies. Similar gene structures and protein domains can be found in members of the same gene cluster. The TaCaM/CML genes were spread among all 21 chromosomes with unbalanced distributions, while most of the gene clusters contained 3 homoeologous genes located in the same homoeologous chromosome group. Synteny analysis showed that most of TaCaM/CMLs gene members can be found with 1-4 paralogous genes in T. turgidum and Ae. Tauschii. High numbers of cis-acting elements related to plant hormones and stress responses can be observed in the promoters of TaCaM/CMLs. The spatiotemporal expression patterns showed that most of the TaCaM/TaCML genes can be detected in at least one tissue. The expression levels of TaCML17, 21, 30, 50, 59 and 75 in the root or shoot can be up-regulated by abiotic stresses, suggesting that TaCML17, 21, 30, 50, 59 and 75 may be related with responses to abiotic stresses in wheat. The spatiotemporal expression patterns of TaCaM/CML genes indicated they may be involved widely in wheat growth and development. Our results provide important clues for exploring functions of TaCaMs/CMLs in growth and development as well as responses to abiotic stresses in wheat in the future.

Identification of Salvia miltiorrhiza Bunge with high and low cadmium accumulation and insight into the mechanisms of cadmium accumulation

Salvia miltiorrhiza Bunge is used as a Chinese herbal medicine (CHM) particularly its roots. No relevant reports at home and abroad have been made on the mechanism of cadmium (Cd) accumulation in S. miltiorrhiza. The Cd accumulation characteristics of S. miltiorrhiza from main cultivation areas in China were evaluated for the first time to obtain high and low Cd accumulation in S. miltiorrhiza roots. Results showed obvious differences in the Cd enrichment capacity of S. miltiorrhiza from different cultivation areas. We took the lead in identifying the germplasm resources of S. miltiorrhiza with high and low Cd accumulation, that is, S. miltiorrhiza roots from Pingyi Shangdong (SDPY) belongs to the resource with high Cd accumulation (SDPY-H) and that from Zhongjiang Sichuan (SCZJ) is the resources with low Cd accumulation (SCZJ-L) based on relevant physiological and biochemical indexes. Although the Cd content of SDPY-H was apparently higher than that from SCZJ-L, its translocation factor from root to aboveground part is significantly lower than that in SCZJ-L. Therefore, planting SCZJ-L is not only an economic and effective way to use Cd in slightly and moderately polluted soil, but also its aboveground part can be used for phytoremediation. Changes in Cd content before and after the use of transpiration inhibitor indicate that SDPY-H enriched Cd through the symplastic pathway, whereas SCZJ-L mainly enriched Cd through the apoplastic pathways. In addition, the role of the symplastic pathway in SCZJ-L is weaker than that in SDPY-H, which were preliminarily revealed by fluorescent quantitative polymerase chain reaction. The significant reduction of the SmNramps transcription expression amount is one of the reasons for the low Cd accumulation of SCZJ-L.

A CCaMK/Cyclops response element in the promoter of Lotus japonicus calcium-binding protein 1 (CBP1) mediates transcriptional activation in root symbioses

Early gene expression in arbuscular mycorrhiza (AM) and the nitrogen-fixing root nodule symbiosis (RNS) is governed by a shared regulatory complex. Yet many symbiosis-induced genes are specifically activated in only one of the two symbioses. The Lotus japonicus T-DNA insertion line T90, carrying a promoterless uidA (GUS) gene in the promoter of Calcium Binding Protein 1 (CBP1) is exceptional as it exhibits GUS activity in both root endosymbioses. To identify the responsible cis- and trans-acting factors, we subjected deletion/modification series of CBP1 promoter : reporter fusions to transactivation and spatio-temporal expression analysis and screened ethyl methanesulphonate (EMS)-mutagenized T90 populations for aberrant GUS expression. We identified one cis-regulatory element required for GUS expression in the epidermis and a second element, necessary and sufficient for transactivation by the calcium and calmodulin-dependent protein kinase (CCaMK) in combination with the transcription factor Cyclops and conferring gene expression during both AM and RNS. Lack of GUS expression in T90 white mutants could be traced to DNA hypermethylation detected in and around this element. We concluded that the CCaMK/Cyclops complex can contribute to at least three distinct gene expression patterns on its direct target promoters NIN (RNS), RAM1 (AM), and CBP1 (AM and RNS), calling for yet-to-be identified specificity-conferring factors.

Genome-Wide Identification and Expression Analysis of CsCaM/CML Gene Family in Response to Low-Temperature and Salt Stresses in Chrysanthemum seticuspe

Calmodulin (CaM) and calmodulin-like proteins (CML) act as significant Ca²⁺ sensors binding Ca²⁺ with EF-hand motifs and have been reported to be involved in various environmental stresses in plants. In this study, calmodulin CsCaM/CML gene family members were identified based on the genome of Chrysanthemum seticuspe published recently; a phylogenetic tree was constructed; gene structures and chromosomal locations of CsCaM/CML were depicted; cis-acting regulatory elements were predicted; collinearity and duplicate events of CaM/CML were analyzed using MCScanX software; and the expression levels of CsCaM/CML in response to abiotic stress were analyzed, based on the published RNA-seq data. We identified 86 CsCaM/CML (4 CsCaMs and 82 CsCMLs) genes in total. Promoter sequences of CsCaM/CML contained elements related to abiotic stresses (including low-temperature and anaerobic stresses) and plant hormones (including abscisic acid (ABA), MeJA, and salicylic acid). CsCaM/CML genes were distributed on nine chromosomes unevenly. Collinearity analysis indicated that recent segmental duplications significantly enlarged the scale of the CML family in C. seticuspe. Four CsCMLs (CsCML14, CsCML50, CsCML65, and CsCML79) were statistically differentially regulated under low-temperature and salt stress compared with those in the normal condition. These results indicate diverse roles of CsCaM/CML in plant development and in response to environmental stimuli in C. seticuspe.

Overview of the roles of calcium sensors in plants’ response to osmotic stress signalling

Calcium signals serve an important function as secondary messengers between cells in various biological processes due to their robust homeostatic mechanism, maintaining an intracellular free Ca2+ concentration. Plant growth, development, and biotic and abiotic stress are all regulated by Ca2+ signals. Ca2+ binding proteins decode and convey the messages encoded by Ca2+ ions. In the presence of high quantities of Mg2+ and monovalent cations, such sensors bind to Ca2+ ions and modify their conformation in a Ca2+ -dependent manner. Calcium-dependent protein kinases (CPKs), calmodulins (CaMs), and calcineurin B-like proteins are all calcium sensors (CBLs). To transmit Ca2+ signals, CPKs, CBLs, and CaMs interact with target proteins and regulate the expression of their genes. These target proteins may be protein kinases, metabolic enzymes, or cytoskeletal-associated proteins. Beyond its role in plant nutrition as a macroelement and its involvement in the plant cell wall structure, calcium modulates many aspects of development, growth and adaptation to environmental constraints such as drought, salinity and osmotic stresses. This review summarises current knowledge on calcium sensors in plant responses to osmotic stress signalling.

A Comprehensive Analysis of Calmodulin-Like Proteins of Glycine max Indicates Their Role in Calcium Signaling and Plant Defense Against Insect Attack

The calcium (Ca²⁺) signaling is a crucial event during plant-herbivore interaction, which involves a transient change in cytosolic Ca²⁺ concentration, which is sensed by Ca²⁺-sensors, and the received message is transduced to downstream target proteins leading to appropriate defense response. Calmodulin-like proteins (CMLs) are calcium-sensing plant-specific proteins. Although CMLs have been identified in a few plants, they remained uncharacterized in leguminous crop plants. Therefore, a wide-range analysis of CMLs of soybean was performed, which identified 41 true CMLs with greater than 50% similarity with Arabidopsis CMLs. The phylogenetic study revealed their evolutionary relatedness with known CMLs. Further, the identification of conserved motifs, gene structure analysis, and identification of cis-acting elements strongly supported their identity as members of this family and their involvement in stress responses. Only a few Glycine max CMLs (GmCMLs) exhibited differential expression in different tissue types, and rest of them had minimal expression. Additionally, differential expression patterns of GmCMLs were observed during Spodoptera litura-feeding, wounding, and signaling compound treatments, indicating their role in plant defense. The three-dimensional structure prediction, identification of interactive domains, and docking with Ca²⁺ ions of S. litura-inducible GmCMLs, indicated their identity as calcium sensors. This study on the characterization of GmCMLs provided insights into their roles in calcium signaling and plant defense during herbivory.

The Same against Many: AtCML8, a Ca2+ Sensor Acting as a Positive Regulator of Defense Responses against Several Plant Pathogens

Calcium signals are crucial for the activation and coordination of signaling cascades leading to the establishment of plant defense mechanisms. Here, we studied the contribution of CML8, an Arabidopsis calmodulin-like protein in response to Ralstonia solanacearum and to pathogens with different lifestyles, such as Xanthomonas campestris pv. campestris and Phytophtora capsici. We used pathogenic infection assays, gene expression, RNA-seq approaches, and comparative analysis of public data on CML8 knockdown and overexpressing Arabidopsis lines to demonstrate that CML8 contributes to defense mechanisms against pathogenic bacteria and oomycetes. CML8 gene expression is finely regulated at the root level and manipulated during infection with Ralstonia, and CML8 overexpression confers better plant tolerance. To understand the processes controlled by CML8, genes differentially expressed at the root level in the first hours of infection have been identified. Overexpression of CML8 also confers better tolerance against Xanthomonas and Phytophtora, and most of the genes differentially expressed in response to Ralstonia are differentially expressed in these different pathosystems. Collectively, CML8 acts as a positive regulator against Ralstonia solanaceraum and against other vascular or root pathogens, suggesting that CML8 is a multifunctional protein that regulates common downstream processes involved in the defense response of plants to several pathogens.

Genome-wide identification and expression analysis of CaM/CML genes in Brassica napus under abiotic stress

Calmodulin (CaM) and calmodulin-like (CML) proteins are primary calcium (Ca²⁺) sensors and are involved in the regulation of plant development and stress responses by converting calcium signals into transcriptional responses, protein phosphorylation, or metabolic changes. However, the characterization and expression profiling of CaM/CML genes in Brassica napus remain limited. The present study reports that 25 BnaCaM and 168 BnaCML genes were identified in B. napus. The phylogenetics, gene structures, gene motifs, gene chromosomal locations, syntenic and Ka/Ks analysis, promoter cis-acting elements, and expression characteristics in various organs and under abiotic stress were evaluated. The phylogenetic results revealed a total of 11 subgroups, including one unique clade of CaMs distinct from CMLs. Most of group I (CaM), II, III, and X members are intron rich, while members from the other seven groups are intron-less. The majority of CaM/CML proteins have four EF-hands. Syntenic analysis showed that 91.3 % orthologous CaM/CML gene pairs between B. rapa and B. oleracea were retained as homologous gene pairs in B. napus. Ka/Ks analysis indicated that the majority of BnaCaM/CML experienced purifying selection. Expression analysis showed that BnaCaMs genes are highly and ubiquitously expressed in all of the organs and tissues examined, while distinct BnaCMLs are expressed specifically in particular organs and tissues. In total, 129 BnaCaM/CML were induced by abiotic stress and phytohormones. BnaCMLs from group IV, VI, VIII, and X were strongly induced by freezing treatment, but were not or just slightly induced by chilling treatment. The present study is the first to analyze the CaM/CML gene family in B. napus, which is useful for understanding the functions of the BnaCaM/CML in modulating plant responses to abiotic stress, especially freezing stress.

Revealing the full-length transcriptome of caucasian clover rhizome development

BACKGROUND

Caucasian clover (Trifolium ambiguum M. Bieb.) is a strongly rhizomatous, low-crowned perennial leguminous and ground-covering grass. The species may be used as an ornamental plant and is resistant to cold, arid temperatures and grazing due to a well-developed underground rhizome system and a strong clonal reproduction capacity. However, the posttranscriptional mechanism of the development of the rhizome system in caucasian clover has not been comprehensively studied. Additionally, a reference genome for this species has not yet been published, which limits further exploration of many important biological processes in this plant.

RESULT

We adopted PacBio sequencing and Illumina sequencing to identify differentially expressed genes (DEGs) in five tissues, including taproot (T1), horizontal rhizome (T2), swelling of taproot (T3), rhizome bud (T4) and rhizome bud tip (T5) tissues, in the caucasian clover rhizome. In total, we obtained 19.82 GB clean data and 80,654 nonredundant transcripts were analysed. Additionally, we identified 78,209 open reading frames (ORFs), 65,227 coding sequences (CDSs), 58,276 simple sequence repeats (SSRs), 6821 alternative splicing (AS) events, 2429 long noncoding RNAs (lncRNAs) and 4501 putative transcription factors (TFs) from 64 different families. Compared with other tissues, T5 exhibited more DEGs, and co-upregulated genes in T5 are mainly annotated as involved in phenylpropanoid biosynthesis. We also identified betaine aldehyde dehydrogenase (BADH) as a highly expressed gene-specific to T5. A weighted gene co-expression network analysis (WGCNA) of transcription factors and physiological indicators were combined to reveal 11 hub genes (MEgreen-GA3), three of which belong to the HB-KNOX family, that are up-regulated in T3. We analysed 276 DEGs involved in hormone signalling and transduction, and the largest number of genes are associated with the auxin (IAA) signalling pathway, with significant up-regulation in T2 and T5.

CONCLUSIONS

This study contributes to our understanding of gene expression across five different tissues and provides preliminary insight into rhizome growth and development in caucasian clover.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in apple (Malus × domestica)

Changes in intracellular calcium (Ca²⁺) levels in response to developmental processes or external stimuli serve as signals in eukaryotic cells. These Ca²⁺ signals are likely perceived through sensor proteins that bind Ca²⁺ by EF-hand (a helix-loop-helix structure) motif. Calmodulins (CaMs), a group of well-characterized Ca²⁺ sensors, and calmodulin-like (CMLs) are implicated in a large number of diverse cellular processes, including plant development and stress responses. In this study, apple (Malus × domestica) genes encoding CaM and CML proteins that only possess EF-hand motifs with no other functional domains were analyzed. A total of 4 MdCaM and 58 MdCML genes were identified, which are spread among 16 out of the 17 apple chromosomes. Bioinformatics analyses, including protein characteristics, conserved domain, evolutionary relationships and chromosomal locations, demonstrated the conservation and divergence of MdCaMs/CMLs. In addition, expression analysis showed that MdCaMs/CMLs are expressed in more than one tissue, including shoot tips, roots, mature leaves, flowers and fruit. Furthermore, the expression of some MdCaM/CML members responded to plant hormones (abscisic acid, jasmonic acid) and salt stress, suggesting a potential role of these genes in responses to biotic and abiotic stress. Overexpression of stress-induced MdCML3 gene significantly improved the tolerance of apple calli to salinity and ABA. The identification and characterization of MdCaMs/CMLs in apple lays a foundation for future functional studies of these genes.

Identification, characterization, and expression analysis of calmodulin and calmodulin-like genes in grapevine (Vitis vinifera) reveal likely roles in stress responses

Calcium (Ca²⁺) is an ubiquitous key second messenger in plants, where it modulates many developmental and adaptive processes in response to various stimuli. Several proteins containing Ca²⁺ binding domain have been identified in plants, including calmodulin (CaM) and calmodulin-like (CML) proteins, which play critical roles in translating Ca²⁺ signals into proper cellular responses. In this work, a genome-wide analysis conducted in Vitis vinifera identified three CaM- and 62 CML-encoding genes. We assigned gene family nomenclature, analyzed gene structure, chromosomal location and gene duplication, as well as protein motif organization. The phylogenetic clustering revealed a total of eight subgroups, including one unique clade of VviCaMs distinct from VviCMLs. VviCaMs were found to contain four EF-hand motifs whereas VviCML proteins have one to five. Most of grapevine CML genes were intronless, while VviCaMs were intron rich. All the genes were well spread among the 19 grapevine chromosomes and displayed a high level of duplication. The expression profiling of VviCaM/VviCML genes revealed a broad expression pattern across all grape organs and tissues at various developmental stages, and a significant modulation in biotic stress-related responses. Our results highlight the complexity of CaM/CML protein family also in grapevine, supporting the versatile role of its different members in modulating cellular responses to various stimuli, in particular to biotic stresses. This work lays the foundation for further functional and structural studies on specific grapevine CaMs/CMLs in order to better understand the role of Ca²⁺-binding proteins in grapevine and to explore their potential for further biotechnological applications.

A mongolian pine specific endoplasmic reticulum localized CALMODULIN-LIKE calcium binding protein enhances arabidopsis growth

Stress-adapted wild plants are natural sources of novel genes for molecular breeding. Here, we conducted a transcriptional analysis of Pinus sylvestris var. mongolica Litv, an evergreen pine in northeastern China, to identify a novel CALMODULIN-LIKE protein-encoding gene, PsCML1, no significant homologs found in other plant species. PsCML1 encodes a protein predicted to have a single trans-membrane domain at its N-terminal. Four EF-hand motifs (calcium [Ca]-binding structures) are located at its C-terminal and showed Ca²⁺-specific affinity in isothermal titration calorimetric analysis. Transient expression of PsCML1 in Nicotiana benthamiana showed that the PsCML1 localizes to the endoplasmic reticulum (ER). Heterologous expression of PsCML1 in Arabidopsis significantly promoted seedling growth, and increased resistance to stress from NaCl and Ca²⁺ deficiency. The roots of the transgenic seedlings had higher contents of cellulose and pectin, but less hemicellulose than those of the wild type (WT). The biosynthesis of cell wall components is linked with protein glycosylation in the ER and reactive oxygen species (ROS) homeostasis. No significant difference was found in the extent of protein glycosylation between the transgenic and WT plants. However, the transgenic roots had higher steady-state levels of ROS, NADPH oxidase activity, and endo-membrane dynamics than those of the WT. A working model was proposed to delineate the interaction among Ca²⁺, ROS homeostasis, and cell wall loosening-dependent cell division.

Towards Understanding Plant Calcium Signaling through Calmodulin-Like Proteins: A Biochemical and Structural Perspective

Ca2+ ions play a key role in a wide variety of environmental responses and developmental processes in plants, and several protein families with Ca2+-binding domains have evolved to meet these needs, including calmodulin (CaM) and calmodulin-like proteins (CMLs). These proteins have no catalytic activity, but rather act as sensor relays that regulate downstream targets. While CaM is well-studied, CMLs remain poorly characterized at both the structural and functional levels, even if they are the largest class of Ca2+ sensors in plants. The major structural theme in CMLs consists of EF-hands, and variations in these domains are predicted to significantly contribute to the functional versatility of CMLs. Herein, we focus on recent advances in understanding the features of CMLs from biochemical and structural points of view. The analysis of the metal binding and structural properties of CMLs can provide valuable insight into how such a vast array of CML proteins can coexist, with no apparent functional redundancy, and how these proteins contribute to cellular signaling while maintaining properties that are distinct from CaM and other Ca2+ sensors. An overview of the principal techniques used to study the biochemical properties of these interesting Ca2+ sensors is also presented.

Genome-Wide Analyses of Calcium Sensors Reveal Their Involvement in Drought Stress Response and Storage Roots Deterioration after Harvest in Cassava

Calcium (Ca²⁺) plays a crucial role in plant development and responses to environmental stimuli. Currently, calmodulins (CaMs), calmodulin-like proteins (CMLs), and calcineurin B-like proteins (CBLs), such as Ca²⁺ sensors, are not well understood in cassava (Manihotesculenta Crantz), an important tropical crop. In the present study, 8 CaMs, 48 CMLs, and 9 CBLs were genome-wide identified in cassava, which were divided into two, four, and four groups, respectively, based on evolutionary relationship, protein motif, and gene structure analyses. Transcriptomic analysis revealed the expression diversity of cassava CaMs-CMLs-CBLs in distinct tissues and in response to drought stress in different genotypes. Generally, cassava CaMs-CMLs-CBLs showed different expression profiles between cultivated varieties (Arg7 and SC124) and wild ancestor (W14) after drought treatment. In addition, numerous CaMs-CMLs-CBLs were significantly upregulated at 6 h, 12 h, and 48 h after harvest, suggesting their possible role during storage roots (SR) deterioration. Further interaction network and co-expression analyses suggested that a CBL-mediated interaction network was widely involved in SR deterioration. Taken together, this study provides new insights into CaMs-CMLs-CBLs-mediated drought adaption and SR deterioration at the transcription level in cassava, and identifies some candidates for the genetic improvement of cassava.