Next Article in Journal
Erratum: Jimena Soledad Cadona, et al.; Pathogenicity Islands Distribution in Non-O157 Shiga Toxin-Producing Escherichia coli (STEC). Genes 2018, 9, 81
Previous Article in Journal
Full Mitogenomes in the Critically Endangered Kākāpō Reveal Major Post-Glacial and Anthropogenic Effects on Neutral Genetic Diversity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Xueyuan Road 4, Haikou 571101, Hainan, China
*
Authors to whom correspondence should be addressed.
Genes 2018, 9(4), 221; https://doi.org/10.3390/genes9040221
Submission received: 22 February 2018 / Revised: 1 April 2018 / Accepted: 12 April 2018 / Published: 19 April 2018
(This article belongs to the Section Plant Genetics and Genomics)

Abstract

:
Calcium (Ca2+) 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 Ca2+ sensors, are not well understood in cassava (Manihot esculenta 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.

1. Introduction

As a second messenger, calcium (Ca2+) plays a crucial role in various biological processes [1]. Cellular calcium concentrations are modulated in plants responding to signals, such as light, hormones, pathogens, and abiotic stress [2,3]. The transient cellular Ca2+ increase can be recognized and decoded by four classes of Ca2+ sensors, including calmodulins (CaMs), calmodulin-like proteins (CMLs), calcineurin B-like proteins (CBLs), and calcium-dependent protein kinases (CDPKs/CPKs), which results in the activation of downstream events, including protein phosphorylation and gene expression in various organisms [4]. All four types of Ca2+ sensors have the typical EF-hand domain that is responsible for Ca2+ binding.
Ca2+ sensors are widely involved in multiple processes related to plant growth and development, such as root hair elongation, guard cell regulation, pollen tube growth, and hormone signaling. Vacuolar CBL2/3-CIPK12 (CBL-interacting protein kinases) complexes are required for polarized pollen tube growth [5]. CBL2/3 also functions positively, affecting leaf, root, silique, and seed development and ion homeostasis through regulation of vacuolar-type H+-ATPase (V-ATPase) activity [6,7]. CML39 plays a positive role in transducing light signals and promoting early seedling establishment [8]. CBL1 and CBL9 regulate pollen germination and pollen tube growth through affecting K+ (potassium) homeostasis [9]. CBL1 also plays a role in plants responding to glucose and gibberellin signals during germination and seedling development [10]. Moreover, there is considerable evidence showing that Ca2+ sensors participate in the abiotic stress response by regulating abscisic acid (ABA) signaling, stomatal behavior, ion homeostasis, and osmotic balance. In Arabidopsis, many genes encoding Ca2+ sensors were confirmed to positively regulate plant tolerance of abiotic stress (osmotic, drought, cold, or salt), including CPK3, CPK4, CPK6, CPK10, CPK11, CPK21, CPK27, CPK32, CaM1, CaM3, CaM4, and CBL1 [11,12,13,14,15,16]. This is also validated in other plants [17,18]. In contrast, some Ca2+ sensors were found to negatively relate to abiotic stress response. cml9 knock-out mutants showed enhanced tolerance of drought and salt stresses in Arabidopsis [19]. Overexpression and mutation of CPK23 revealed its negative role in tolerance to drought and salt stresses [20]. Loss-of-function of CPK21 resulted in increased tolerance to hyperosmotic stress [13]. These evidences indicate that the CDPK-mediated abiotic stress response is complex.
Currently, based on genome-wide analysis, Ca2+ sensors have been identified in various species [21,22,23,24,25,26]. There are 10 CBLs in Arabidopsis and 10 in rice, seven CaMs in Arabidopsis and five in rice, and 50 CMLs in Arabidopsis and 32 in rice. Previously, we identified 27 CPKs from cassava (Manihot esculenta Crantz) [24]. However, less information is known about other Ca2+ sensors (CBLs, CaMs, and CMLs) in cassava. As a staple food crop, cassava is the sixth most important crop after wheat, rice, maize, potato, and barley [27]. Its edible storage roots (SR) supply a source of dietary carbohydrate for over 600 million people worldwide [28]. Cassava is also considered as a potential biofuel crop for production of ethanol and bioenergy, due to its high starch production [27]. Moreover, cassava has strong tolerance of drought and low-fertility environments, due to its effective use of light, heat, and water resources. Notably, the potential of cassava as a food and industrial crop is largely restricted by the rapid post-harvest physiological deterioration of the SR, which starts within 72 h after harvest [29]. The mechanism underlying cassava tolerance of drought stress and sensitivity to post-harvest physiological deterioration are largely unknown.
Considering the importance of Ca2+ sensors in plant growth, development, and response to abiotic stresses, an effort was made to identify CBLs, CaMs, and CMLs from cassava, and investigate their phylogenetic relationships, protein motifs, gene structure, and expression profiles in distinct tissues, in response to drought stress and during SR deterioration. This systematic study would increase our understanding of Ca2+ sensors associated with drought responses and SR deterioration, and lays a foundation for the genetic improvement of cassava.

2. Materials and Methods

2.1. Plant Materials and Treatments

The characteristics of SC124, Arg7, and W14 cassava genotypes were described in previous studies [30,31]. W14 is an accession of the wild ancestor, and had stronger drought resistance than the cultivated varieties SC124 and Arg7 [31]. Segments of cassava stems from mother plants were cultured in pots filled with soil and vermiculite (1:1) in a growth room with a 16 h/35 °C day and 8 h/20 °C night regime, and a relative humidity of 70%. Thereafter, stems (90-day-old), leaves (90-day-old), and SR (150-day-old) were cut from Arg7 and W14 under growth room conditions, to study the expression levels of CaM-CML-CBLs in distinct organs. To detect the transcriptional changes of CaM-CML-CBLs in response to drought, 90-day-old cassava plants of Arg7, SC124, and W14 were subjected to water withholding for 12 days, and the leaves and roots were collected for RNA sequencing. To study the transcriptional changes of CaM-CML-CBLs during SR deterioration, 10-month-old cassava SR with 5–6 cm diameter stems were cut into slices approximately 5 mm thick, and then moved into petri dishes containing a wet filter paper [29]. After incubation at 28 °C and 60% relative humidity in the dark for 0 h, 6 h, 12 h, and 48 h, the slices were sampled and frozen in liquid nitrogen until the extraction of total RNA.

2.2. Identification and Phylogenetic Analyses

The whole protein sequences of cassava were downloaded from the cassava genome database [32]. The CaM, CML, and CBL protein sequences from rice and Arabidopsis were obtained from the rice genome annotation project (RGAP) and UniPort databases, respectively [33,34]. The hidden markov models (HMM) profiles built from CaMs, CMLs, and CBLs in Arabidopsis and rice were used as queries to identify the putative cassava CaMs, CMLs, and CBLs with HMMER software [35]. In addition, protein sequences of CaMs, CMLs, and CBLs from Arabidopsis and rice were used in a basic local alignment search tool (BLAST [36]) search against cassava proteins, to detect other CaMs, CMLs, and CBLs that might be missed by the HMM profile. After removing the redundant sequences, the conserved domain (EF-hand) of predicted cassava CaMs, CMLs, and CBLs were further confirmed with Pfam and conserved domains database (CDD) databases [37,38]. The accession number of identified cassava CaMs, CMLs, and CBLs is displayed in Table S1. The phylogenetic tree was created with the CaMs, CMLs, and CBLs from Arabidopsis, rice, and cassava using MEGA 5.0 and Clustal X2.0 softwares (bootstrap values for 1000 replicates) [39,40].

2.3. Protein Properties and Sequence Analyses

The isoelectric points and molecular weights of cassava CaMs, CMLs, and CBLs were predicted using ExPASy database [41]. The conserved motifs of cassava CaMs, CMLs, and CBLs were identified by MEME database, and were further annotated with InterProScan database [42,43]. The optimum width of motifs ranged from 6 to 50 and the maximum number of motifs was 5. The gene structures of cassava CaMs, CMLs, and CBLs were analyzed by gene structure display server (GSDS) [44]. The interaction network of CaM, CML, and CBL members were built by STRING database [45], with interaction score > 0.7 and no more than 20 interactors.

2.4. Transcriptomic Analysis

The total RNA of each sample, extracted with a plant RNA extraction kit (DP432, TIANGEN Biotech, Beijing, China), was used to construct complementary DNA (cDNA) libraries. The sequencing was performed with an Illumina GAII (Illumina, San Diego, CA, USA), following the manufacturer’s instructions. Single-end libraries (non-strand-specific) were sequenced with a read length of 51 bp. Adapter sequences in the raw sequence reads were removed using the FASTX-toolkit [46]. After examining the sequence quality and removing low quality sequences by FastQC [47], clean reads were generated. Using TopHat v2.0.10, clean reads were mapped to the cassava reference genome (v4.1) [48]. The transcriptome assemblies were performed by Cufflinks [49]. Gene expression levels were calculated as fragments per kilobase of transcript per million fragments mapped (FPKM). Differentially expressed genes between treatments (after drought or after harvest) and control were identified based on absolute value of Log2 based fold changes > 1. The accession number of transcriptomic data was listed in Table S1.

3. Results

3.1. Genome-Wide Identification and Evolutionary Analyses of the Calmodulins (CaMs), Calmodulin-Like Proteins (CMLs), and Calcineurin B-Like Proteins (CBLs) in Cassava

Previously, we identified 27 calcium-dependent protein kinases and one type of calcium sensor from the cassava genome [24]. Here, a total of eight CaM, 48 CML, and nine CBL proteins were identified from the cassava genome, based on HMM and BLAST searches. The predicted CaM, CML, and CBL proteins varied from 84 to 318 amino acid residues, with relative molecular masses in the range of 9.3–34.9 kDa, and theoretical isoelectric points (PIs) ranging from 3.9 to 7.0 (Table S2).
To understand the phylogenetic relationship of CaM, CML, and CBL proteins, neighbor-joining (NJ) trees were created with CaM, CML, and CBL proteins from cassava, Arabidopsis and rice (Figure 1; Table S3). The cassava CaM, CML, and CBL families were classified into 2 (group A-B), 4 (group A-D), and 4 (group A-D) groups, respectively, based on their phylogenetic relationships. Generally, cassava CaMs, CMLs, and CBLs showed closer relationships with those in Arabidopsis than those in rice, supported by the orthologous genes between Arabidopsis and cassava.

3.2. Conserved Motifs and Gene Structure Analyses of calmodulins (CaMs), Calmodulin-Like Proteins (CMLs), and Calcineurin B-Like Proteins (CBLs) in Cassava

To investigate structural features of the cassava CaMs, CMLs, and CBLs, conserved motifs were analyzed based on the phylogenetic relationship. Four conserved motifs that were annotated as EF-hand domains for the CaM, CML, and CBL families were acquired with MEME and InterPro databases (Figure 2). For the CaM family, all the members contain four4 EF-hand domains. For the CML family, groups A, B, C, and D have 1–3, 2–4, 2–4 and 2–4 EF-hand domains, respectively. For the CBL family, all the members show four EF-hand domains, except for MeCBL7, which has three EF-hand domains. Based on the above results, all the identified CaMs, CMLs, and CBLs had at least one EF-hand domain, indicating their typical family features of calcium binding.
Additionally, the exon-intron organizations of these genes were tested using the GSDS database (Figure 3). For the CaM family, group A shows two to three exons and group B contains four exons. For the CML family, group A shows three to five exons, group B has four exons, group C has one to two exons, and group D displays one exon, except for MeCML24, with six exons. For the CBL family, groups A, B, C, and D contain seven to eight, eight, nine, and eight exons, respectively. Generally, the proteins classified into the same subgroup share similar motifs and exon-intron organizations, suggesting a link between evolution and conserved motifs/gene structure.

3.3. Expression Analyses of Calmodulins (CaMs), Calmodulin-Like Proteins (CMLs), and Calcineurin B-Like Proteins (CBLs) in Different Cassava Tissues

To study the expression profiles of CaMs, CMLs, and CBLs in different tissues, transcriptome analyses were carried out from samples of leaves, stems, and storage roots in wild ancestor (W14) and cultivated variety (Arg7) (Figure 4; Table S4). Generally, most of the CaMs, CMLs, and CBLs had similar expression profiles in different tissues between Arg7 and W14. For example, several genes (MeCaM-1, -4, MeCML-19 -24, -27, -47, -48 and MeCBL-8) showed high expression levels (Log2 FPKM value > 5) in various tissues of Arg7 and W14. In contrast, some genes (MeCaM-8, MeCML-1, -26, -41, -43 and MeCBL-7) had low transcript abundance (Log2 FPKM value < 1) in both Arg7 and W14.
Additionally, some CaM, CML, and CBL genes showed different expression profiles in different genotypes. MeCaM3 and MeCML-8, -22, -23, -40 showed high expression levels (Log2 FPKM > 5) in Arg7 stem, but low expression (Log2 FPKM < 3) in W14 stem. MeCaM-6, -7 had low transcript abundance (Log2 FPKM < 3) in Arg7 leaf, but high expression (Log2 FPKM > 5) in W14 leaf. MeCML-8, -11, -12, -13, -36 had high expression (Log2 FPKM > 5) in Arg7 storage root, but low expression (Log FPKM < 3) in W14 storage root. Together, the tissue expression profiles of CaMs, CMLs, and CBLs would be benefit for further study of tissue development and function.

3.4. Expression Analyses of Calmodulins (CaMs), Calmodulin-Like Proteins (CMLs), and Calcineurin B-Like Proteins (CBLs) Genes in Response to Drought Stress

Since previous reports revealed that the CaM, CML, and CBL families could participate in osmotic adjustment and drought stress response, the expression profiles of CaMs, CMLs, and CBLs in response to drought stress were further detected in three cassava genotypes by transcriptome analysis (Figure 4; Table S5). Totally, 14 and 10 MeCaMs-MeCMLs-MeCBLs showed induction in leaf and root of Arg7, respectively. The corresponding number is 12 and four in SC124, and five and 13 in W14. From the above results, we found that: (i) the total number of MeCaMs-MeCMLs-MeCBLs induced by drought is greater in Arg7 than in SC124 and W14; (ii) the number of MeCaMs-MeCMLs-MeCBLs induced by drought was greater in roots than in leaves in W14, whereas it was reduced in roots compared with leaves in Arg7 and SC124.

3.5. Expression Patterns of Cassava Calmodulins (CaMs), Calmodulin-Like Proteins (CMLs), and Calcineurin B-Like Proteins (CBLs) after Harvest

Transcriptomic analysis was performed to study the expression profiles of cassava CaMs, CMLs, and CBLs after harvest (Figure 4; Table S6). In comparison to cross-sessions of SR at 0 h after harvest, 38, 38, and 42 calcium sensor genes showed up-regulation at 6 h, 12 h, and 48 h after harvest, respectively, among which 30, 27, and 33 genes showed significant induction (Log2 based fold change > 1) at 6 h, 12 h, and 48 h after harvest, respectively. Moreover, 23 genes had significantly increased transcripts (Log2 based fold change > 1) during all the tested stages after harvest. Together, these results suggested the possible involvement of calcium signaling in the regulation of cassava SR deterioration.

3.6. Calmodulins (CaMs), Calmodulin-Like Proteins (CMLs), and Calcineurin B-Like Proteins (CBLs) Interaction Networks and Their Co-Expression after Harvest

To better understand the function of cassava CaM-CML-CBLs, the interaction networks and co-expression profiles of these proteins were detected according to experimentally validated interactions of CaM-CML-CBLs in Arabidopsis (Figure 5; Tables S7–S12). Firstly, CaM-, CML-, and CBL-mediated interaction networks in Arabidopsis were built using STRING database (Tables S7–S9). Then, the homologs of these proteins in the interaction network were identified in cassava using the reciprocal BLAST protein analysis (Table S10). Finally, the expression profiles of these cassava genes at 0 h and 48 h after harvest were searched from our previous RNA-seq. data sets. The gene pairs showing upregulation (Log2 based fold change > 1) or downregulation (Log2 based fold change < −1) at 48 h/0 h after harvest were displayed in Table S10.
For the CaM-mediated interaction network, no gene pair displayed co-expression (Figure S1A; Table S11). For the CML-mediated interaction network, only one gene pair (AT2G27480:Me014830-AT1G15130:Me001472) showed uniform upregulation (Figure S1B; Table S11). For the CBL-mediated interaction network, SOS3/MeCBL3-mediated seven gene pairs and CBL1/MeCBL5-mediated four gene pairs had uniform upregulation (Figure 5; Tables S11 and S12). These genes included eight CBL-interacting protein kinases (CIPKs), one phosphatidylinositol 4-kinase (PI-4KBETA1), and one SNF1-related protein kinase (AKINBETA1) (Table S10). These results indicated that CBL-mediated interaction and co-expression network may play a crucial role during cassava SR deterioration.

4. Discussion

Calcium, as an important second messenger, plays crucial roles in various aspects of biological processes in plants. Calcium signals are sensed and decoded by Ca2+ sensors, including CaMs, CMLs, CBLs, and CPKs [24,30]. Previously, we identified 27 CPKs from the cassava genome [24]. Here, eight MeCaMs, 48 MeCMLs, and nine MeCBLs were identified from the cassava genome, and their phylogenetic relationships, gene structures, and conserved motifs were investigated (Figure 1, Figure 2 and Figure 3). Phylogenetic analyses showed that cassava CaMs, CMLs, and CBLs could be divided into two, four and four groups, respectively. together with those in Arabidopsis and rice (Figure 1), which is consistent with previous phylogenetic classification of CaMs, CMLs, and CBLs in Arabidopsis, rice, strawberry, cotton, grapevine, and Lotus japonicus [22,23,25,50]. Moreover, this classification was also validated by conserved motif analysis, which indicated that all the identified CaMs, CMLs, and CBLs contained the typical EF-hand, and each subfamily shared similar motifs and exon-intron organizations (Figure 2 and Figure 3). The typical EF-hand organizations of CaMs, CMLs, and CBLs in cassava were the same as those in other plant species. The CaMs strictly contain four typical EF-hand domains in various plant species [22,25,51]; the CMLs generally had variable numbers of EF-hands from one to six [22,25,51]; and the CBLs generally had four EF-hands [21,23]. Thus, the conserved EF-hand organizations of CaMs, CMLs, and CBLs are conserved across plant species. Interestingly, CaMs in group B from cassava apparently separated from Arabidopsis and rice with strong bootstrap support, indicating an independent evolution of these CaMs in cassava relative to Arabidopsis and rice (Figure 1). This group of cassava CaMs showed similar conserved motifs and exon-intron organization, implying their consistent origin (Figure 2 and Figure 3). Notably, cassava CaMs in group A (MeCaM1 and MeCaM4) were clustered together with CaMs from Arabidopsis and rice. All of these genes are broadly expressed in various tissues of cassava, Arabidopsis, and rice (Figure 4B) [51,52]. Besides, some CMLs and CBLs with close phylogenetic relationship among cassava, Arabidopsis, and rice also showed similar expression profiles in different tissues [52,53]. This indicated a possible link between sequence conservation and gene expression for the CaMs, CMLs, and CBLs.
Accumulated evidences have documented the positive role of Ca2+ sensors in drought stress response [11,54,55,56,57]. In the present study, we found that some cassava MeCaMs-MeCMLs-MeCBLs were upregulated by drought stress in roots and leaves of different genotypes, suggesting the possible role of Ca2+ sensor genes in cassava response to drought stress (Figure 4C). Additionally, the number of MeCaMs-MeCMLs-MeCBLs induced by drought is greater in Arg7 than in SC124 and W14, suggesting the comprehensive activation of MeCaMs-MeCMLs-MeCBLs in the Arg7 variety (Figure 4C). Arg7 was confirmed to be a drought-sensitive variety, which showed lower tolerance to drought than W14 and SC124 [24]. It is concluded that Arg7 variety has strong response to drought stress, thus leading to the comprehensive transcriptional induction of MeCaMs-MeCMLs-MeCBLs. Interestingly, the number of MeCaMs-MeCMLs-MeCBLs induced by drought was greater in roots than in leaves in W14, whereas it was reduced in roots compared with leaves in Arg7 and SC124 (Figure 4C). Cassava maintains its robust drought tolerance through absorbing water stored in the deeper soil layers by the deep root systems [58]. Thus, it is speculated that significant induction of MeCaMs-MeCMLs-MeCBLs in roots of W14 may be involved in water uptake, thus benefit for maintaining its strong tolerance to drought stress.
During cassava harvest and storage, SR deterioration occurrence is difficult to prevent due to storage conditions, handling operations, and mechanical injury [29]. Physiological and genetic evidences revealed that oxidation caused by reactive oxygen species (ROS) lead to the SR deterioration symptoms, and reduction of ROS accumulation could attenuate SR deterioration [27,29,59,60]. Calcium, an important signaling molecular, plays a protective role in plants response to abiotic stress [61]. One of the important events underlying this protection is activating antioxidative system and repressing ROS accumulation [62,63,64,65]. These studies suggested that calcium signaling may play a role in regulating SR deterioration by influencing cellular ROS level. In the present study, we found that numerous calcium sensor genes were significantly induced in cross-sessions of SR at 6 h, 12 h, and 48 h compared with 0 h after harvest (Figure 4D). Further interaction network and co-expression analyses identified some gene pairs uniformly upregulated at 48 h/0 h in cross-sessions of SR after harvest (Figure 5). These results indicated the possible involvement of calcium sensors during SR deterioration of cassava. Based on the transcriptomic data, there were some Ca2+ sensor genes that showed induction after drought treatment and during SR deterioration, including MeCML-4, -10, -16, -17, -20, -22, -25, -26, and MeCBL2. These genes can serve as candidates for genetic improvement of cassava resistance to drought and SR deterioration.

5. Conclusions

This study identified 12 MeCaMs, 56 MeCMLs, and 9 MeCBLs from the cassava genome, and investigated their classification, protein motif, and gene structure. Transcriptome analysis showed the possible role of MeCaMs-MeCMLs-MeCBLs against drought stress and in cross-sessions of SR deterioration after harvest. Interaction network and co-expression analyses revealed the involvement of a CBL-mediated network in SR deterioration after harvest. Together, these results will advance the understanding of calcium signaling-mediated drought stress response and SR deterioration regulation in cassava.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/4/221/s1. Figure S1: Interaction network and co-expression analyses of CAMs (A) and CMLs (B) at 48 h/0 h after harvest, Table S1: The accession number of transcriptomic data in NCBI, Table S2: Characteristics of CaMs, CMLs, and CBLs in cassava, Table S3: The accession numbers of CaMs, CMLs, and CBLs in Arabidopsis and rice, Table S4: Expression of cassava CaMs, CMLs, and CBLs in different tissues of Arg7 and W14, Table S5: Expression of cassava CaMs, CMLs, and CBLs in response to drought stress in Arg7, SC124 and W14, Table S6: Expression of cassava CaMs, CMLs, and CBLs during SR deterioration, Table S7: The protein interaction relationship in the CaM-mediated interaction network in Arabidopsis, Table S8: The protein interaction relationship in the CML-mediated interaction network in Arabidopsis, Table S9: The protein interaction relationship in the CBL-mediated interaction network in Arabidopsis, Table S10: Expression data of the genes involved in CaM/CBL/CML-mediate interaction networks at CK48/CK0 after harvest, Table S11: The co-expression patterns of gene pairs involved in CAM/CBL/CML-mediate interaction networks at 48 h/0 h after harvest, Table S12: The accession number of the genes in the CBL-mediated interaction network.

Acknowledgments

This work was supported by the Key Research and Development Project of Hainan Province (ZDYF2018220), the National Natural Science Foundation of China (31771859), the Natural Science Foundation of Hainan Province (317255), the Central Public-interest Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences (1630052016005, 1630052016006, 1630052017021), the Central Public-interest Scientific Institution Basal Research Fund for Innovative Research Team Program of CATAS (17CXTD-28, 1630052017017), and the earmarked fund for Modern Agro-industry Technology Research System (CARS-11).

Author Contributions

J.G. and M.P. conceived and designed the experiments. W.H., Y.Y., Z.D., C.W., W.W. and Z.X. performed the experiments. W.H., W.T. and C.W. analyzed the data. W.T. and X.D. contributed reagents/materials/analysis tools. W.H. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Asano, T.; Tanaka, N.; Yang, G.; Hayashi, N.; Komatsu, S. Genome-wide identification of the rice calcium-dependent protein kinase and its closely related kinase gene families: Comprehensive analysis of the CDPKs gene family in rice. Plant Cell Physiol. 2005, 46, 356–366. [Google Scholar] [CrossRef] [PubMed]
  2. Harper, J.F. Dissecting calcium oscillators in plant cells. Trends Plant Sci. 2001, 6, 395–397. [Google Scholar] [CrossRef]
  3. Knight, H.; Knight, M.R. Abiotic stress signalling pathways: Specificity and cross-talk. Trends Plant Sci. 2001, 6, 262–267. [Google Scholar] [CrossRef]
  4. Sanders, D.; Brownlee, C.; Harper, J.F. Communicating with calcium. Plant Cell 1999, 11, 691–706. [Google Scholar] [CrossRef] [PubMed]
  5. Steinhorst, L.; Mahs, A.; Ischebeck, T.; Zhang, C.; Zhang, X.; Arendt, S.; Schultke, S.; Heilmann, I.; Kudla, J. Vacuolar CBL-CIPK12 Ca2+-sensor-kinase complexes are required for polarized pollen tube growth. Curr. Biol. 2015, 25, 1475–1482. [Google Scholar] [CrossRef] [PubMed]
  6. Eckert, C.; Offenborn, J.N.; Heinz, T.; Armarego-Marriott, T.; Schultke, S.; Zhang, C.; Hillmer, S.; Heilmann, M.; Schumacher, K.; Bock, R.; et al. The vacuolar calcium sensors CBL2 and CBL3 affect seed size and embryonic development in Arabidopsis thaliana. Plant J. 2014, 78, 146–156. [Google Scholar] [CrossRef] [PubMed]
  7. Tang, R.J.; Liu, H.; Yang, Y.; Yang, L.; Gao, X.S.; Garcia, V.J.; Luan, S.; Zhang, H.X. Tonoplast calcium sensors CBL2 and CBL3 control plant growth and ion homeostasis through regulating V-ATPase activity in Arabidopsis. Cell Res. 2012, 22, 1650–1665. [Google Scholar] [CrossRef] [PubMed]
  8. Bender, K.W.; Rosenbaum, D.M.; Vanderbeld, B.; Ubaid, M.; Snedden, W.A. The Arabidopsis calmodulin-like protein, CML39, functions during early seedling establishment. Plant J. 2013, 76, 634–647. [Google Scholar] [CrossRef] [PubMed]
  9. Mahs, A.; Steinhorst, L.; Han, J.P.; Shen, L.K.; Wang, Y.; Kudla, J. The calcineurin B-like Ca2+ sensors CBL1 and CBL9 function in pollen germination and pollen tube growth in Arabidopsis. Mol. Plant 2013, 6, 1149–1162. [Google Scholar] [CrossRef] [PubMed]
  10. Li, Z.Y.; Xu, Z.S.; Chen, Y.; He, G.Y.; Yang, G.X.; Chen, M.; Li, L.C.; Ma, Y.Z. A novel role for Arabidopsis CBL1 in affecting plant responses to glucose and gibberellin during germination and seedling development. PLoS ONE 2013, 8, e56412. [Google Scholar] [CrossRef] [PubMed]
  11. Wei, S.; Hu, W.; Deng, X.; Zhang, Y.; Liu, X.; Zhao, X.; Luo, Q.; Jin, Z.; Li, Y.; Zhou, S.; et al. A rice calcium-dependent protein kinase OsCPK9 positively regulates drought stress tolerance and spikelet fertility. BMC Plant Biol. 2014, 14, 133. [Google Scholar] [CrossRef] [PubMed]
  12. Zhou, S.; Jia, L.; Chu, H.; Wu, D.; Peng, X.; Liu, X.; Zhang, J.; Zhao, J.; Chen, K.; Zhao, L. Arabidopsis CaM1 and CaM4 promote nitric oxide production and salt resistance by inhibiting S-nitrosoglutathione reductase via direct binding. PLoS Genet. 2016, 12, e1006255. [Google Scholar] [CrossRef] [PubMed]
  13. Franz, S.; Ehlert, B.; Liese, A.; Kurth, J.; Cazale, A.C.; Romeis, T. Calcium-dependent protein kinase CPK21 functions in abiotic stress response in Arabidopsis thaliana. Mol. Plant 2011, 4, 83–96. [Google Scholar] [CrossRef] [PubMed]
  14. Zhao, R.; Sun, H.; Zhao, N.; Jing, X.; Shen, X.; Chen, S. The Arabidopsis Ca2+-dependent protein kinase CPK27 is required for plant response to salt-stress. Gene 2015, 563, 203–214. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, W.; Zhou, R.G.; Gao, Y.J.; Zheng, S.Z.; Xu, P.; Zhang, S.Q.; Sun, D.Y. Molecular and genetic evidence for the key role of AtCaM3 in heat-shock signal transduction in Arabidopsis. Plant Physiol. 2009, 149, 1773–1784. [Google Scholar] [CrossRef] [PubMed]
  16. Cheong, Y.H.; Kim, K.N.; Pandey, G.K.; Gupta, R.; Grant, J.J.; Luan, S. CBL1, a calcium sensor that differentially regulates salt, drought, and cold responses in Arabidopsis. Plant Cell 2003, 15, 1833–1845. [Google Scholar] [CrossRef] [PubMed]
  17. Asano, T.; Hayashi, N.; Kobayashi, M.; Aoki, N.; Miyao, A.; Mitsuhara, I.; Ichikawa, H.; Komatsu, S.; Hirochika, H.; Kikuchi, S.; et al. A rice calcium-dependent protein kinase OsCPK12 oppositely modulates salt-stress tolerance and blast disease resistance. Plant J. 2012, 69, 26–36. [Google Scholar] [CrossRef] [PubMed]
  18. Rao, S.S.; El-Habbak, M.H.; Havens, W.M.; Singh, A.; Zheng, D.; Vaughn, L.; Haudenshield, J.S.; Hartman, G.L.; Korban, S.S.; Ghabrial, S.A. Overexpression of GmCaM4 in soybean enhances resistance to pathogens and tolerance to salt stress. Mol. Plant Pathol. 2014, 15, 145–160. [Google Scholar] [CrossRef] [PubMed]
  19. Magnan, F.; Ranty, B.; Charpenteau, M.; Sotta, B.; Galaud, J.P.; Aldon, D. Mutations in AtCML9, a calmodulin-like protein from Arabidopsis thaliana, alter plant responses to abiotic stress and abscisic acid. Plant J. 2008, 56, 575–589. [Google Scholar] [CrossRef] [PubMed]
  20. Ma, S.Y.; Wu, W.H. AtCPK23 functions in Arabidopsis responses to drought and salt stresses. Plant Mol. Biol. 2007, 65, 511–518. [Google Scholar] [CrossRef] [PubMed]
  21. Kolukisaoglu, U.; Weinl, S.; Blazevic, D.; Batistic, O.; Kudla, J. Calcium sensors and their interacting protein kinases: Genomics of the Arabidopsis and rice CBL-CIPK signaling networks. Plant Physiol. 2004, 134, 43–58. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, K.; Yue, D.; Wei, W.; Hu, Y.; Feng, J.; Zou, Z. Characterization and functional analysis of calmodulin and calmodulin-like genes in Fragaria vesca. Front. Plant Sci. 2016, 7, 1820. [Google Scholar] [CrossRef] [PubMed]
  23. Xi, Y.; Liu, J.; Dong, C.; Cheng, Z.M. The CBL and CIPK gene family in grapevine (Vitis vinifera): Genome-wide analysis and expression profiles in response to various abiotic stresses. Front. Plant Sci. 2017, 8, 978. [Google Scholar] [CrossRef] [PubMed]
  24. Hu, W.; Hou, X.; Xia, Z.; Yan, Y.; Wei, Y.; Wang, L.; Zou, M.; Lu, C.; Wang, W.; Peng, M. Genome-wide survey and expression analysis of the calcium-dependent protein kinase gene family in cassava. Mol. Genet. Genomics 2016, 291, 241–253. [Google Scholar] [CrossRef] [PubMed]
  25. Liao, J.; Deng, J.; Qin, Z.; Tang, J.; Shu, M.; Ding, C.; Liu, J.; Hu, C.; Yuan, M.; Huang, Y.; et al. Genome-wide identification and analyses of calmodulins and calmodulin-like proteins in Lotus japonicas. Front. Plant Sci. 2017, 8, 482. [Google Scholar] [CrossRef] [PubMed]
  26. Yin, X.; Wang, Q.; Chen, Q.; Xiang, N.; Yang, Y. Genome-wide identification and functional analysis of the calcineurin B-like protein and calcineurin B-like protein-interacting protein kinase gene families in Turnip (Brassica rapa var. rapa). Front. Plant Sci. 2017, 8, 1191. [Google Scholar] [CrossRef] [PubMed]
  27. Zidenga, T.; Leyva-Guerrero, E.; Moon, H.; Siritunga, D.; Sayre, R. Extending cassava root shelf life via reduction of reactive oxygen species production. Plant Physiol. 2012, 159, 1396–1407. [Google Scholar] [CrossRef] [PubMed]
  28. Zhang, P.; Wang, W.Q.; Zhang, G.L.; Kaminek, M.; Dobrev, P.; Xu, J.; Gruissem, W. Senescence-inducible expression of isopentenyl transferase extends leaf life, increases drought stress resistance and alters cytokinin metabolism in cassava. J. Integr. Plant Biol. 2010, 52, 653–669. [Google Scholar] [CrossRef] [PubMed]
  29. Vanderschuren, H.; Nyaboga, E.; Poon, J.S.; Baerenfaller, K.; Grossmann, J.; Hirsch-Hoffmann, M.; Kirchgessner, N.; Nanni, P.; Gruissem, W. Large-Scale proteomics of the cassava storage root and identification of a target gene to reduce postharvest deterioration. Plant Cell 2014, 26, 1913–1924. [Google Scholar] [CrossRef] [PubMed]
  30. Hu, W.; Xia, Z.; Yan, Y.; Ding, Z.; Tie, W.; Wang, L.; Zou, M.; Wei, Y.; Lu, C.; Hou, X.; et al. Genome-wide gene phylogeny of CIPK family in cassava and expression analysis of partial drought-induced genes. Front. Plant Sci. 2015, 6, 914. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, W.; Feng, B.; Xiao, J.; Xia, Z.; Zhou, X.; Li, P.; Zhang, W.; Wang, Y.; Moller, B.L.; Zhang, P.; et al. Cassava genome from a wild ancestor to cultivated varieties. Nat. Commun. 2014, 5, 5110. [Google Scholar] [CrossRef] [PubMed]
  32. Prochnik, S.; Marri, P.R.; Desany, B.; Rabinowicz, P.D.; Kodira, C.; Mohiuddin, M.; Rodriguez, F.; Fauquet, C.; Tohme, J.; Harkins, T.; et al. The cassava genome: Current progress, future directions. Trop. Plant Biol. 2012, 5, 88–94. [Google Scholar] [CrossRef] [PubMed]
  33. The UniProt Consortium. UniProt: A hub for protein information. Nucleic Acids Res. 2015, 43, D204–D212. [Google Scholar] [CrossRef]
  34. Kawahara, Y.; de la Bastide, M.; Hamilton, J.P.; Kanamori, H.; McCombie, W.R.; Ouyang, S.; Schwartz, D.C.; Tanaka, T.; Wu, J.; Zhou, S.; et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice 2013, 6, 4. [Google Scholar] [CrossRef] [PubMed]
  35. Finn, R.D.; Clements, J.; Eddy, S.R. HMMER web server: Interactive sequence similarity searching. Nucleic Acids Res. 2011, 39, W29–W37. [Google Scholar] [CrossRef] [PubMed]
  36. Basic Local Alignment Search Tool. Available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 4 November 2017).
  37. Marchler-Bauer, A.; Derbyshire, M.K.; Gonzales, N.R.; Lu, S.; Chitsaz, F.; Geer, L.Y.; Geer, R.C.; He, J.; Gwadz, M.; Hurwitz, D.I.; et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015, 43, D222–D226. [Google Scholar] [CrossRef] [PubMed]
  38. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef] [PubMed]
  39. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X. version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed]
  40. 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–2739. [Google Scholar] [CrossRef] [PubMed]
  41. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef] [PubMed]
  42. Brown, P.; Baxter, L.; Hickman, R.; Beynon, J.; Moore, J.D.; Ott, S. MEME-LaB: Motif analysis in clusters. Bioinformatics 2013, 29, 1696–1697. [Google Scholar] [CrossRef] [PubMed]
  43. Mulder, N.; Apweiler, R. InterPro and InterProScan: Tools for protein sequence classification and comparison. Methods Mol Biol. 2007, 396, 59–70. [Google Scholar] [CrossRef] [PubMed]
  44. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [PubMed]
  45. Jensen, L.J.; Kuhn, M.; Stark, M.; Chaffron, S.; Creevey, C.; Muller, J.; Doerks, T.; Julien, P.; Roth, A.; Simonovic, M.; et al. STRING 8—A global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res. 2009, 37, D412–D416. [Google Scholar] [CrossRef] [PubMed]
  46. FASTX-Toolkit. FASTQ/A short-reads pre-processing tools. Available online: http://hannonlab.cshl.edu/fastx_toolkit/index.html (accessed on 5 December 2017).
  47. FastQC: A quality control tool for high throughput sequence data. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 15 December 2017).
  48. Trapnell, C.; Pachter, L.; Salzberg, S.L. TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  49. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed]
  50. Lu, T.; Zhang, G.; Sun, L.; Wang, J.; Hao, F. Genome-wide identification of CBL family and expression analysis of CBLs in response to potassium deficiency in cotton. PeerJ 2017, 5, e3653. [Google Scholar] [CrossRef] [PubMed]
  51. Wang, X.; Komatsu, S. Proteomic analysis of calcium effects on soybean root tip under flooding and drought stresses. Plant Cell Physiol. 2017, 58, 1405–1420. [Google Scholar] [CrossRef] [PubMed]
  52. Xu, C.; Li, X.; Zhang, L. The effect of calcium chloride on growth, photosynthesis, and antioxidant responses of Zoysia japonica under drought conditions. PLoS ONE 2013, 8, e68214. [Google Scholar] [CrossRef] [PubMed]
  53. Ranty, B.; Aldon, D.; Cotelle, V.; Galaud, J.P.; Thuleau, P.; Mazars, C. Calcium sensors as key hubs in plant responses to biotic and abiotic stresses. Front. Plant Sci. 2016, 7, 327. [Google Scholar] [CrossRef] [PubMed]
  54. Xu, J.; Tian, Y.S.; Peng, R.H.; Xiong, A.S.; Zhu, B.; Jin, X.F.; Gao, F.; Fu, X.Y.; Hou, X.L.; Yao, Q.H. AtCPK6, a functionally redundant and positive regulator involved in salt/drought stress tolerance in Arabidopsis. Planta 2010, 231, 1251–1260. [Google Scholar] [CrossRef] [PubMed]
  55. Zou, J.J.; Wei, F.J.; Wang, C.; Wu, J.J.; Ratnasekera, D.; Liu, W.X.; Wu, W.H. Arabidopsis calcium-dependent protein kinase CPK10 functions in abscisic acid- and Ca2+-mediated stomatal regulation in response to drought stress. Plant Physiol. 2010, 154, 1232–1243. [Google Scholar] [CrossRef] [PubMed]
  56. Li, Z.Y.; Xu, Z.S.; He, G.Y.; Yang, G.X.; Chen, M.; Li, L.C.; Ma, Y.Z. Overexpression of soybean GmCBL1 enhances abiotic stress tolerance and promotes hypocotyl elongation in Arabidopsis. Biochem. Biophys. Res. Commun. 2012, 427, 731–736. [Google Scholar] [CrossRef] [PubMed]
  57. Munir, S.; Liu, H.; Xing, Y.; Hussain, S.; Ouyang, B.; Zhang, Y.; Li, H.; Ye, Z. Overexpression of calmodulin-like (ShCML44) stress-responsive gene from Solanum habrochaites enhances tolerance to multiple abiotic stresses. Sci. Rep. 2016, 6, 31772. [Google Scholar] [CrossRef] [PubMed]
  58. Okogbenin, E.; Setter, T.L.; Ferguson, M.; Mutegi, R.; Ceballos, H.; Olasanmi, B.; Fregene, M. Phenotypic approaches to drought in cassava: Review. Front. Physiol. 2013, 4, 93. [Google Scholar] [CrossRef] [PubMed]
  59. Reilly, K.; Gomez-Vasquez, R.; Buschmann, H.; Tohme, J.; Beeching, J.R. Oxidative stress responses during cassava post-harvest physiological deterioration. Plant Mol. Biol. 2004, 56, 625–641. [Google Scholar] [CrossRef] [PubMed]
  60. Xu, J.; Duan, X.; Yang, J.; Beeching, J.R.; Zhang, P. Enhanced reactive oxygen species scavenging by overproduction of superoxide dismutase and catalase delays postharvest physiological deterioration of cassava storage roots. Plant Physiol. 2013, 161, 1517–1528. [Google Scholar] [CrossRef] [PubMed]
  61. Gilroy, S.; Suzuki, N.; Miller, G.; Choi, W.G.; Toyota, M.; Devireddy, A.R.; Mittler, R. A tidal wave of signals: Calcium and ROS at the forefront of rapid systemic signaling. Trends Plant Sci. 2014, 19, 623–630. [Google Scholar] [CrossRef] [PubMed]
  62. Deng, X.; Hu, W.; Wei, S.; Zhou, S.; Zhang, F.; Han, J.; Chen, L.; Li, Y.; Feng, J.; Fang, B.; et al. TaCIPK29, a CBL-interacting protein kinase gene from wheat, confers salt stress tolerance in transgenic tobacco. PLoS ONE 2013, 8, e69881. [Google Scholar] [CrossRef] [PubMed]
  63. Deng, X.; Zhou, S.; Hu, W.; Feng, J.; Zhang, F.; Chen, L.; Huang, C.; Luo, Q.; He, Y.; Yang, G.; et al. Ectopic expression of wheat TaCIPK14, encoding a calcineurin B-like protein-interacting protein kinase, confers salinity and cold tolerance in tobacco. Physiol. Plant 2013, 149, 367–377. [Google Scholar] [CrossRef] [PubMed]
  64. Verslues, P.E.; Batelli, G.; Grillo, S.; Agius, F.; Kim, Y.S.; Zhu, J.; Agarwal, M.; Katiyar-Agarwal, S.; Zhu, J.K. Interaction of SOS2 with nucleoside diphosphate kinase 2 and catalases reveals a point of connection between salt stress and H2O2 signaling in Arabidopsis thaliana. Mol. Cell Biol. 2007, 27, 7771–7780. [Google Scholar] [CrossRef] [PubMed]
  65. Zou, J.J.; Li, X.D.; Ratnasekera, D.; Wang, C.; Liu, W.X.; Song, L.F.; Zhang, W.Z.; Wu, W.H. Arabidopsis CALCIUM-DEPENDENT PROTEIN KINASE8 and CATALASE3 function in abscisic acid-mediated signaling and H2O2 homeostasis in stomatal guard cells under drought stress. Plant Cell 2015, 27, 1445–1460. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogenetic analysis of calmodulin-like proteins (CMLs) (A), calmodulins (CaMs) (B), and calcineurin B-like proteins (CBLs) (C) from cassava, Arabidopsis, and rice. The Neighbor-joining (NJ) tree was constructed using Clustal X 2.0 and MEGA 5.0 software with the pair-wise deletion option. One thousand bootstrap replicates were used to assess tree reliability. Blue triangle, calcium (Ca2+) sensors in Arabidopsis; Red square, Ca2+ sensors in cassava; Green circle, Ca2+ sensors in rice. Shade areas with different colors indicate distinct groups of CaM and CBL families.
Figure 1. Phylogenetic analysis of calmodulin-like proteins (CMLs) (A), calmodulins (CaMs) (B), and calcineurin B-like proteins (CBLs) (C) from cassava, Arabidopsis, and rice. The Neighbor-joining (NJ) tree was constructed using Clustal X 2.0 and MEGA 5.0 software with the pair-wise deletion option. One thousand bootstrap replicates were used to assess tree reliability. Blue triangle, calcium (Ca2+) sensors in Arabidopsis; Red square, Ca2+ sensors in cassava; Green circle, Ca2+ sensors in rice. Shade areas with different colors indicate distinct groups of CaM and CBL families.
Genes 09 00221 g001
Figure 2. The conserved motifs of cassava CaMs, CMLs, and CBLs according to phylogenetic relationship. All the grey boxes indicate EF-hand domains for each protein. All motifs were identified by MEME database with the complete amino acid sequences of cassava CaMs, CMLs, and CBLs. Each family showed similar motif organization. Shade areas with different colors indicate distinct groups of CaM, CML, and CBL families.
Figure 2. The conserved motifs of cassava CaMs, CMLs, and CBLs according to phylogenetic relationship. All the grey boxes indicate EF-hand domains for each protein. All motifs were identified by MEME database with the complete amino acid sequences of cassava CaMs, CMLs, and CBLs. Each family showed similar motif organization. Shade areas with different colors indicate distinct groups of CaM, CML, and CBL families.
Genes 09 00221 g002
Figure 3. Gene structure analyses of cassava CaMs, CMLs, and CBLs according to phylogenetic relationship. Exon-intron structure analyses were performed through the gene structure display server database. The blue boxes, yellow boxes, and black lines indicate upstream of initiation codon/downstream of stop codon, exons, and introns, respectively. Shade areas with different colors indicate distinct groups of CaM, CML, and CBL families.
Figure 3. Gene structure analyses of cassava CaMs, CMLs, and CBLs according to phylogenetic relationship. Exon-intron structure analyses were performed through the gene structure display server database. The blue boxes, yellow boxes, and black lines indicate upstream of initiation codon/downstream of stop codon, exons, and introns, respectively. Shade areas with different colors indicate distinct groups of CaM, CML, and CBL families.
Genes 09 00221 g003
Figure 4. Expression profiles of CaMs, CMLs, and CBLs in cassava, according to phylogenetic relationship. (A) Phylogenetic relationship of cassava CaMs, CMLs, and CBLs. (B) Expression profiles of cassava CaMs, CMLs, and CBLs in stems (90-day-old), leaves (90-day-old), and storage roots (150-day-old) of Arg7 and W14. Asterisks indicate the genes that show high expression levels (Log2 based fragments per kilobase of transcript per million fragments mapped (FPKM) > 5) in all the tissues examined. S, stem; L, leaf; SR, storage root. (C) Expression profiles of cassava CaMs, CMLs, and CBLs in response to drought stress in leaves and roots of Arg7, SC124, and W14. Ninety-day-old cassava plants were subjected to water withholding for 12 days, and leaves and roots were collected for RNA-seq. Asterisks indicate the genes that are induced (Log2 based fold changes > 1) after drought treatment. L, leaf; R, root; D, drought treatment; CK, control check showing leaves or roots of cassava varieties under normal conditions. (D) Expression profiles of cassava CaMs, CMLs, and CBLs in cross-sessions of SR at 6 h, 12 h, and 48 h compared with 0 h after harvest. Ten-month-old storage roots were cut into slices approximately 5 mm thick, and then were moved into petri dishes containing a wet filter paper. After incubation at 28 °C and 60% relative humidity in the dark for 0 h, 6 h, 12 h, and 48 h, the slices were sampled for RNA-seq. Asterisks indicate the genes that is induced (Log2 based fold changes > 1) after harvest. AH, after harvest. Log2 based FPKM (B) or fold change (C) and (D) was used to create the heat map. Changes in gene expression are shown in color as the scale.
Figure 4. Expression profiles of CaMs, CMLs, and CBLs in cassava, according to phylogenetic relationship. (A) Phylogenetic relationship of cassava CaMs, CMLs, and CBLs. (B) Expression profiles of cassava CaMs, CMLs, and CBLs in stems (90-day-old), leaves (90-day-old), and storage roots (150-day-old) of Arg7 and W14. Asterisks indicate the genes that show high expression levels (Log2 based fragments per kilobase of transcript per million fragments mapped (FPKM) > 5) in all the tissues examined. S, stem; L, leaf; SR, storage root. (C) Expression profiles of cassava CaMs, CMLs, and CBLs in response to drought stress in leaves and roots of Arg7, SC124, and W14. Ninety-day-old cassava plants were subjected to water withholding for 12 days, and leaves and roots were collected for RNA-seq. Asterisks indicate the genes that are induced (Log2 based fold changes > 1) after drought treatment. L, leaf; R, root; D, drought treatment; CK, control check showing leaves or roots of cassava varieties under normal conditions. (D) Expression profiles of cassava CaMs, CMLs, and CBLs in cross-sessions of SR at 6 h, 12 h, and 48 h compared with 0 h after harvest. Ten-month-old storage roots were cut into slices approximately 5 mm thick, and then were moved into petri dishes containing a wet filter paper. After incubation at 28 °C and 60% relative humidity in the dark for 0 h, 6 h, 12 h, and 48 h, the slices were sampled for RNA-seq. Asterisks indicate the genes that is induced (Log2 based fold changes > 1) after harvest. AH, after harvest. Log2 based FPKM (B) or fold change (C) and (D) was used to create the heat map. Changes in gene expression are shown in color as the scale.
Genes 09 00221 g004
Figure 5. Interaction network and co-expression analyses of CBLs at 48 h/0 h in cross-sessions of SR after harvest. The interacting relationship between CBLs and other proteins are indicated by the gray lines. The genes at upper and lower around each ball show CBL-mediated interaction network in Arabidopsis and cassava, respectively. The gene pairs that showed uniform upregulation (Log2 based fold change > 1) at 48 h/0 h after harvest in cassava were marked with yellow.
Figure 5. Interaction network and co-expression analyses of CBLs at 48 h/0 h in cross-sessions of SR after harvest. The interacting relationship between CBLs and other proteins are indicated by the gray lines. The genes at upper and lower around each ball show CBL-mediated interaction network in Arabidopsis and cassava, respectively. The gene pairs that showed uniform upregulation (Log2 based fold change > 1) at 48 h/0 h after harvest in cassava were marked with yellow.
Genes 09 00221 g005

Share and Cite

MDPI and ACS Style

Hu, W.; Yan, Y.; Tie, W.; Ding, Z.; Wu, C.; Ding, X.; Wang, W.; Xia, Z.; Guo, J.; Peng, M. Genome-Wide Analyses of Calcium Sensors Reveal Their Involvement in Drought Stress Response and Storage Roots Deterioration after Harvest in Cassava. Genes 2018, 9, 221. https://doi.org/10.3390/genes9040221

AMA Style

Hu W, Yan Y, Tie W, Ding Z, Wu C, Ding X, Wang W, Xia Z, Guo J, Peng M. Genome-Wide Analyses of Calcium Sensors Reveal Their Involvement in Drought Stress Response and Storage Roots Deterioration after Harvest in Cassava. Genes. 2018; 9(4):221. https://doi.org/10.3390/genes9040221

Chicago/Turabian Style

Hu, Wei, Yan Yan, Weiwei Tie, Zehong Ding, Chunlai Wu, Xupo Ding, Wenquan Wang, Zhiqiang Xia, Jianchun Guo, and Ming Peng. 2018. "Genome-Wide Analyses of Calcium Sensors Reveal Their Involvement in Drought Stress Response and Storage Roots Deterioration after Harvest in Cassava" Genes 9, no. 4: 221. https://doi.org/10.3390/genes9040221

APA Style

Hu, W., Yan, Y., Tie, W., Ding, Z., Wu, C., Ding, X., Wang, W., Xia, Z., Guo, J., & Peng, M. (2018). Genome-Wide Analyses of Calcium Sensors Reveal Their Involvement in Drought Stress Response and Storage Roots Deterioration after Harvest in Cassava. Genes, 9(4), 221. https://doi.org/10.3390/genes9040221

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop