Next Article in Journal
A Simple and Efficient Molecularly Imprinted Electrochemical Sensor for the Selective Determination of Tryptophan
Previous Article in Journal
Antihypertensive Indigenous Lebanese Plants: Ethnopharmacology and a Clinical Trial
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response

1
Key Laboratory of Genetics, Breeding and Multiple Utilization of Crops, Ministry of Education, Fujian Provincial Key Laboratory of Haixia Applied Plant Systems Biology, College of Crop Science, Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China
2
State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou 350002, China
3
Life Science College, Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China
4
College of Forestry, Fujian Agriculture and Forestry University, Fuzhou 350002, Fujian, China
5
State Key Laboratory for Conservation and Utilization of Subtropical Agro-Bioresources, Guangxi Key Lab of Sugarcane Biology, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, China
*
Authors to whom correspondence should be addressed.
Biomolecules 2019, 9(7), 293; https://doi.org/10.3390/biom9070293
Submission received: 3 July 2019 / Revised: 19 July 2019 / Accepted: 19 July 2019 / Published: 20 July 2019

Abstract

:
Ca2+ serves as a ubiquitous second messenger regulating several aspects of plant growth and development. A group of unique calcium sensor proteins, calcineurin B-like (CBL), interact with CBL-interacting protein kinases (CIPKs) to decode the Ca2+ signature inside the cell. Although CBL-CIPK signaling toolkit has been shown to play significant roles in the responses to numerous stresses in different plants, the information about pineapple CBL-CIPK remains obscure. In the present study, a total of eight AcCBL and 21 AcCIPK genes were identified genome-wide in pineapple. The identified genes were renamed on the basis of gene ID in ascending order and phylogenetic analysis divided into five groups. Transcriptomic data analysis showed that AcCBL and AcCIPK genes were expressed differentially in different tissues. Further, the expression analysis of AcCBL1 in different tissues showed significant changes under various abiotic stimuli. Additionally, the ectopic expression of AcCBL1 in Arabidopsis resulted in enhanced tolerance to salinity, osmotic, and fungal stress. The present study revealed the crucial contribution of the CBL-CIPK gene in various biological and physiological processes in pineapple.

1. Introduction

In living organisms, the precise perception and timely decoding of environmental and developmental signals are essential for survival. The signal transduction pathway involved in this process includes several indispensable components including calcium (Ca2+), which serves as a ubiquitous second messenger in all eukaryotes [1,2,3]. Various organelles such as endoplasmic reticulum (ER), mitochondria, and vacuoles work as Ca2+ store and help in the maintenance of a critical balance of Ca2+ inside the cell [4]. Acknowledging the stress or developmental signals, the cytosolic calcium concentration rapidly shoots up as an early response [5]. Consistently, in plants, the rapid increase of cytosolic Ca2+ concentration is well documented against several stress factors such as salinity, drought, and cold [3,6,7]. In addition, pH dynamics, pathogens, phytohormones, light signal, and drugs, as well as the development of a pollen tube and regulation of guard cells are also linked with changes in the concentration of Ca2+ [8,9,10,11,12,13]. The elevated Ca2+ level inside the cell is then recognized by several calcium-binding proteins or calcium sensors [5]. After sensing the elevation in Ca2+ level these calcium sensors bind to calcium in their specific motifs (such as EF-hands), which changes their phosphorylation status. The change in the phosphorylation status of calcium sensors activates several protein kinases which sometimes lead to a protein phosphorylation cascade [14]. A number of calcium sensors have been reported in plants including calmodulin (CaM), calcineurin B-like (CBLs), and calcium-dependent protein kinases (CDPKs) [15].
CBLs are plant-specific genes that encode a protein similar to the calcineurin B subunit of protein phosphatase of yeast and animal cells [16]. CBLs were initially identified in Arabidopsis as an essential player of the salt stress response [17]. CBLs specifically bind with the conserved C-terminal NAF/FISL domain of CBL-interacting protein kinases (CIPKs) [18]. Both, the CBLs and CIPKs are encoded by the multigene family in plants. The Arabidopsis genome encodes 10 members of CBLs and 26 members of CIPKs [19]. In addition, individual CBLs can physically interact with multiple CIPKs after sensing the specific calcium signature of a specific signal in the cellular environment [20]. This interaction activates the downstream components of the CBL-CIPK signaling pathway, including ion channels and transporters involved in plant physiological and developmental response [18,21]. The complex signaling network of CBL-CIPK in plants participates in several essential functions of plant growth and development. Calcineurin B-like proteins have been shown to play key roles in biotic and abiotic stress responses [3]. For example, the CBL7 mediates the alkaline stress response in Arabidopsis by promoting the self-inhibition of plasma membrane H+-ATPase [22]. Consistently, CBL1 plays key roles in cold stress [23,24], aluminum stress [25], activation of AKT1 channel [26], glucose, and gibberellic acid signal [27]. In addition, CBL1 along with CBL9, is also involved in pollen germination and pollen tube growth [28], and in the regulation of potassium transport [29]. Similarly, CBL2 and CBL3 are involved in seed development and seed morphology [30], CBL3 in potassium homeostasis [31], and CBL4-CIPK6 complex modulates the activity and plasma membrane targeting of the K+ channel AKT2 [32].
Abiotic stresses greatly limit plant growth and development. Among them, salt stress is a prime abiotic stress which severely affects the plant productivity and significantly decreases the agricultural yield. Roughly, 20% of the irrigated agricultural land is under salt stress [33]. In plants, salt stress initiates the highly conserved pathways which export Na+ ions outside the cell, known as salt overlay sensitive (SOS) pathway [17,33]. SOS pathway key players SOS2 kinase and SOS1 Na+ antiporter are activated within 2 h of salt stress [34]. In Arabidopsis, CBL10 reads and decodes the salt-induced cytosolic changes in calcium concentration and further activates the CIPK24 (SOS2 kinase) [3,17,35]. Salinity induced SOS2 kinase then phosphorylates plasma membrane Na+/H+-ATPase (SOS1), which transports the Na+ to apoplast from the cytoplasm [34,35].
Pineapple, Ananas comosus (L.) Merr., belonging to the family Bromeliaceae is grown commercially in tropical and subtropical regions of the world. Pineapple is highly appreciated for its delicious fruit and aroma, ranking it among the top three economically important tropical fruits in the world, with banana and citrus [36]. The growth and development of the pineapple plant is seriously affected by the biotic and abiotic stresses. However, it can withstand up to some level of drought conditions owing to its morphological features and plants crassulacean acid metabolism [37]. Due to its worldwide cultivation and immense economic and research significance, researchers have shown considerable interest in fishing out crucial functional genes from pineapple.
CBL-CIPK gene family has been studied in several plants including Arabidopsis [19], Oryza [38], eggplant [39], grapevine [40], populus [41], brassica [42], and Zea mays [38]. However, the information about this crucial gene family in pineapple is still obscure. Due to the great significance of CBL-CIPK genes in various physiological and developmental processes, it is of great interest to investigate the CBL-CIPK genes in pineapple. Genome sequence of pineapple [43] provides an opportunity for genome-wide investigation of the CBL-CIPK genes. In this study, we identified eight CBL and 21 CIPK genes and performed a comprehensive analysis including exon-intron organization, chromosome distribution, and phylogenetic analysis. Additionally, we have performed expression analysis to identify the involvement of these genes in different biological and developmental stages in pineapple. We further characterized AcCBL1 by overexpressing it in Arabidopsis. Our study provides a basis for further utilization of CBL-CIPK genes to develop next-generation crops especially in the context of a globally changing environment.

2. Materials and Methods

2.1. Identification of CBL-CIPK Genes in Pineapple Genome

The CBL and CIPK gene sequences from Arabidopsis and Oryza sativa were downloaded from TAIR (http://www.arabidopsis.org) and China Rice Data Center (http://www.ricedata.cn/gene/index.htm) respectively. Further, the CBL and CIPK genes were used as keywords against the pineapple genome database of phytozome (https://phytozome.jgi.doe.gov/pz/portal.html) to identify CBL-CIPK genes from pineapple. Additionally, Hidden Markov Model (HMM) profile of the NAF domain (PF03822) from the Pfam database (http://pfam.xfam.org/) and EF-hand calcium-binding domain (PS50222) from prosite (https://prosite.expasy.org/) were downloaded. We then performed searches from the pineapple genome using the HMM profiles by BLAST-P with the e-value set at 0.01. The identified sequences were further verified and redundant sequences were removed from the list. Finally, ExPASy (http://web.expasy.org/compute_pi/) was used determine the isoelectric point (pI) and molecular weight (MW) of CBL and CIPK proteins of pineapple.

2.2. Phylogenetic Analysis

The phylogenetic relationship of CBL-CIPK proteins between pineapple, Arabidopsis, and rice were studied using identified CBL-CIPK amino acid sequences from pineapple, Arabidopsis, and rice. The multiple sequence alignments were performed by using MUSCLE 3.7 and the phylogenetic trees were constructed by MEGA7 using the Maximum likelihood (ML) method with default parameters with the bootstrap option n = 1000.

2.3. Gene Structure Analysis and Conserved Motif Identification

The schematic AcCBL-AcCIPK genes structures were drawn by the Gene Structure Display Server 130 (http://gsds.cbi.pku.edu.cn/). The conserved motifs in the pineapple AcCBL-AcCIPK proteins were identified by MEME program (http://meme.nbcr. net/meme/cgi-bin/meme.cgi).

2.4. Chromosome Location of CBL and CIPK Genes in Pineapple

The information about the location of AcCBL and AcCIPK genes on chromosomes were collected from phytozome. MapChart software was used to visualize the AcCBL and AcCIPK genes mapped on chromosomes based on the gene start position and the length of related chromosomes.

2.5. Plant Materials and Growth Condition and Treatments

2.5.1. Pineapple Growth and Treatments

Two-month-old tissue culture raised pineapple (Ananas comosus) variety MD2 plants were obtained from Qin lab (http://www.qinlab.net) [44] and grown in plastic pots containing soil mix (peat moss:perlite = 2:1 (v/v)) placed at 30 °C with 16 h light/8 h dark photoperiod under an intensity of 70 μmol m−2 s−1 70% humidity as described earlier by Chen et al. [45]. For the various stress treatment experiments, the pineapple plants were exposed to different treatments for 24, 48, and 72 h. The treatments included salinity stress (150 mM NaCl), osmotic stress (350 mM mannitol), heat (45 °C), and cold stress (4 °C) respectively. The leaf and root tissues from whole plants were harvested and frozen immediately in liquid nitrogen and stored at −80 °C until used.

2.5.2. Arabidopsis Treatments and Root Growth Assay

The wild-type Arabidopsis thaliana (L.) Heynh (Col-0; CS60000) was obtained from the Arabidopsis Biological Resource Center (Columbus, OH, USA). All plants used were Columbian in background. Surface sterilized Arabidopsis seeds were placed in a round, 9 cm Petri plates on modified Hoagland’s medium [46,47] containing 1% (w/v) sucrose and 1% (w/v) agar. The plates were kept at 4 °C in the dark for 2 day for seed stratification. After stratification, the plates were transferred to a growth room at 22 °C with 16 h light/8 h dark photoperiod under an intensity of 100 μmol m−2 sec−1 and seedlings were grown vertically.
Salinity stress treatment and growth recovery were performed from T3 generation of 3 different AcCBL1 overexpressing lines. 5-day-old seedlings were transferred to new plates containing the modified Hoagland medium supplemented with 150 mM NaCl and kept 10 days in a growth room. Primary roots were analyzed after 2 days, 5 days, and 10 days of transfer. The experiments were repeated at least three times with 8–10 seedlings per treatment. To measure salinity induced root growth inhibition seedlings were photographed by digital camera (NIKON D750, Japan) and root growth was analyzed by ImageJ software (http://rsbweb.nih.gov/ij/).
For in-situ detection of hydrogen peroxide induced against biotic stress after exposure to fungus (Sclerotinia sclerotiorum) mature Arabidopsis rosette leaves were exposed to fungus for 24 h followed by staining with 3,3-diaminobenzidine (DAB) as reported earlier by Daudi et al. [48].

2.6. RNA Isolation and Quantitative Real-Time PCR (qRT-PCR) Analysis

Total RNA from desired tissues was isolated by using the RNeasy kit (Qiagen, MD, USA), followed by treatment with DNase I (Thermo Fisher Scientific, CA, USA) and first-strand cDNA synthesis by reverse transcription using the ThermoScript RT-PCR kit (Thermo Fisher Scientific, CA, USA). Quantitative PCR was performed by using FastStart DNA Master SYBR Green I master mix (Takara, Shiga, Japan) in a CFX96 qPCR system (Bio-Rad, Singapore). Primers used in qRT-PCR are listed in supplementary Table S3. The PCR cycle conditions for quantitative PCR were as follows: 95 °C for 2 min followed by 40 cycles of 95 °C for 10 s, 60 °C for 15 s, and 72 °C for 15 s. For each analysis, two technical replicates from three biological replicates were performed and pineapple EF1a was used as an internal control to normalize the mRNA levels, and finally, the fold change of genes was calculated using the 2−ΔΔCT method.

2.7. RNA-Seq and Data Analysis

RNA was extracted from different developmental stages as described earlier by Chen et al., [45], pineapple of MD2 variety through Plant RNeasy Mini Kit (Qiagen MD, USA). Illumina sequencing was performed as previously described [45], with 1 μg RNA per sample and two independent biological replicates per genotype. cDNA libraries were constructed using the NEBNext Ultra™ RNA Library Prep Kit for Illumina (NEB, MA, USA) following standard protocols. The transcript abundance of pineapple CBL-CIPK genes were calculated as fragments per kilobase of exon model per million mapped reads (FPKM). The heatmap was created by pheatmap packages in R based the log2 (FPKM + 1).

2.8. Vector Constructs

AcCBL1:GFP was generated by amplifying 678 bp of coding sequence without the stop codon from pineapple cDNA using the primers listed in Table S3. The amplified PCR fragment was then cloned into the pENTR/D-TOPO vector (Invitrogen, Carlsbad, CA, USA). pENTR/D-TOPO clones were then recombined into the destination vector pGWB505 by using LR Clonase II (Invitrogen, Carlsbad, CA, USA). Before performing the bacterial and plant transformation, the construct was confirmed by sequencing.

2.9. Live Cell Imaging

For live cell microscopy, five-day-old AcCBL1 overexpressing plants were treated either by 150 mM salt or wounded manually. After the treatment, the seedlings were mounted in the water on a cover glass for observation on a TCS SP8 microscope (Leica, Wetzlar, Germany) and imaged.

2.10. Statistical Analysis

Results are expressed as the means ± standard error (SE) from at least 3 experiments. To analyze statistical significance, a two-tailed Student’s t-test was used.

3. Results

3.1. Identification of the CBL-CIPK Genes in Pineapple

A total of 8 CBL and 26 CIPK genes were initially obtained from the pineapple genome using BLAST, HMM search, and available pineapple annotation. Five redundant CIPK genes were removed after the Pfam NAF domain (PF03822) investigation identified putative CIPK genes. Finally, a total of 8 CBLs and 21 CIPK gene models were selected and annotated based on their gene ID in ascending order (Table 1a,b). All the CBLs contained the EF-hand and CIPKs possess the conserved NAF/FISL motif. The detailed information of the identified CBL-CIPK genes such as gene name, gene ID, protein length, isoelectric points, and molecular weights is listed in Table 1. Among 8 identified CBLs, AcCBL6 was the smallest having 193 amino acid and AcCBL4 was largest with 359 amino acid, whereas AcCIPK8 was smallest with 382 amino acids and AcCIPK20 with 506 amino acid was largest in the group. The molecular weight of the AcCBL proteins ranged from 21.9 kDa to 40.8 kDa and 43.7 kDa to 55.8 kDa for CIPK proteins.

3.2. Gene Structure and Phylogenetic Analysis of Pineapple CBL-CIPK Genes

The evolutionary aspect and structural diversity of the CBL-CIPK genes in pineapple were explored by studying the exon-intron organization. The exon-intron structure for AcCBLs and AcCIPKs were constructed using the online Gene Structure Display Server 2.0 [49]. All the AcCBLs possessed introns which ranged from six to ten in numbers (AcCBL6 had six introns, AcCBL1, AcCBL2, AcCBL3, and AcCBL7 possessed seven introns, AcCBL5 had eight introns, AcCBL8 had nine introns, and AcCBL4 possessed maximum ten introns) (Figure 1a). The result also showed that all the AcCBLs contained the UTR, though AcCBL8 had only 3’ UTR. The distribution of introns in AcCIPK genes ranged from zero to thirteen where AcCIPK2, AcCIPK4, AcCIPK11, AcCIPK12, AcCIPK13, and AcCIPK18 did not possess introns. AcCIPK5, AcCIPK9, AcCIPK10, AcCIPK17, and AcCIPK19 possessed single intron, AcCIPK1 and AcCIPK6 had two introns, AcCIPK21 had eleven introns, AcCIPK16, and AcCIPK7 had 12 introns, AcCIPK3, AcCIPK8, AcCIPK14, AcCIPK15, and AcCIPK19 possessed maximum thirteen introns (Figure 1b). Moreover, no UTRs were observed in ten AcCIPKs (AcCIPK2, AcCIPK4, AcCIPK5, AcCIPK6, AcCIPK8, AcCIPK9, AcCIPK11, AcCIPK12, AcCIPK13, and AcCIPK18) whereas two AcCIPKs (AcCIPK1 and AcCIPK20) possessed only one UTR (5’ UTR and 3’ UTR respectively). Additionally, pineapple CBL-CIPK protein structures were further analyzed by performing MEME analysis. Overall, ten conserved motifs were identified in both AcCBL and AcCIPK protein (Supplementary Figures S1 and S2).
The evolutionary relationships and functional associations of CBL-CIPK genes in pineapple were investigated by constructing a multi-species phylogenetic tree using full-length protein sequences of CBL (Figure 2) and CIPK (Figure 3) from pineapple, Arabidopsis, and rice using the Maximum Likelihood method. AcCBLs in phylogenetic trees formed four subfamilies more closely related to rice compared to Arabidopsis. The tree also divided the AcCBLs into 5 subfamilies viz. AcCBL1 and AcCBL7 formed a group with OsaCBL2 and OsaCBL3, AcCBL5 and AcCBL6 formed a group with OsaCBL1, and AcCBL6 paired with OsaCBL6 (Figure 2). Interestingly, some groups (depicted without any color in Figure 2 and Figure 3) both among CBL and CIPK did not have any members from pineapple, indicating the evolutionary divergence between monocotyledonous and dicotyledonous plants species.

3.3. Chromosomal Distribution of Pineapple CBL-CIPK Genes

The pineapple CBL-CIPK genes were unevenly distributed on 15 linkage groups (LG) out of 25 linkage groups present in the genome. The linkage groups LG 06, 09, 14, 15, 17, 21, 22, 23, 24, and 25 did not contain any CBL or CIPK genes (Figure 4). All 8 CBLs were located on different linkage groups, whereas CIPKs were located on different linkage groups ranging from 1-4 in number. Linkage groups LG01, 03, 07, 12, 16, and 19 had a single CIPK gene, LG 02, 08, 13, and 18 had two CIPK genes, LG05 had 3 CIPK genes and maximum 4 CIPK genes were located on LG 04 (Figure 4). The result showed no correlation between the length of linkage groups and number of CBL-CIPK gene distribution.

3.4. Expression Profiling of Pineapple CBL-CIPK Genes in Different Developmental Stages

The expression profiling of all the CBL-CIPK genes was carried out from RNA-Seq data generated from different organs of the pineapple plant. The expression profiles of all the pineapple CBL-CIPK genes were generated and hierarchically clustered using average log values of each gene (Figure 5). The hierarchical clustering of RNA-Seq data divided them into three major groups. In group I, low expressed genes (AcCIPK4, AcCIPK5, AcCIPK9, and AcCIPK11) were clustered together. In contrast, group II contained genes (6 AcCBLs and 13 AcCIPK genes) with high level of expression in all the selected stages (Figure 5). Whereas, group III had genes (2 AcCBLs and 4 AcCIPK genes) with relatively low levels of expression. The expression pattern of CBL-CIPK genes suggests that they are expressed differentially during different developmental stages and contribute essentially to plant growth and development.

3.5. Expression Profiling of Pineapple CBL1 Gene in Response to Different Treatments

In the RNA-Seq experiment AcCBL1 gene showed the highest expression, so we selected it for further study and investigated the role of AcCBL1 gene during abiotic stress response. For this, we used quantitative real-time PCR (qRT-PCR) and determined the spatial and temporal expression profile of AcCBL1 in different tissues (root and leaf) of pineapple at different time intervals (0, 24, 48, and 72 h) after salinity stress (150 mM NaCl), osmotic stress (350 mM mannitol), heat (45 °C), and cold stress (4 °C) treatments (Figure 6). The qPCR data showed that compared to untreated plants the transcript level of AcCBL1 significantly increased during salinity, osmotic, and heat treatments whereas, the cold stress did not show significant response. The results also showed the differential expression of AcCBL1 when it was compared between root and leaf (Figure 6). Collectively, these results strongly suggest that the AcCBL1 gene contributes to various abiotic stresses in pineapple.

3.6. Salt Stress Increases the Localization of AcCBL1 in Arabidopsis Roots

Having established the role of AcCBL1 in salinity, osmotic, and heat stress in pineapple. We further investigated the biological function of AcCBL1 by ectopically expressing AcCBL1 gene in Arabidopsis under CaMV 35S promoter and GFP tag. Using a confocal microscope we monitored the localization of AcCBL1:GFP after transfer experiment. Transfer experiment was carried out with five-day-old seedlings by transferring them to a fresh plate consisting of Hoagland medium supplemented with 150 mM NaCl for 48 h. The results showed that the AcCBL1:GFP mainly localized in the cytosol and its expression increased after salt stress treatment (Figure 7 and Figure 8b). The increased localization of AcCBL1 in Arabidopsis roots clearly indicates its involvement in salt stress. Taken together, our results indicate that the AcCBL1 helps the plant to thrive well under salt stress.

3.7. Injury Induces AcCBL1 Translocation from Cytosol to the Nucleus

Previous studies have established that Arabidopsis CBLs are induced during wounding and injury [50]. In order to investigate the possible role of AcCBL1 during injury, we monitored the expression of AcCBL1:GFP after manually injuring the AcCBL1 overexpressing plants. Interestingly, we found that the expression of AcCBL1 was regulated by wounding and AcCBL1 was translocated to the nucleus after the injury (Figure 8a). The results clearly indicate that the AcCBL1 is regulated by injury and may be involved in the protection against injury stress.

3.8. Overexpression of AcCBL1 Results in Resistance to Salinity, Osmotic, and Biotic Stress

As shown in Figure 6 and Figure 7, the expression of AcCBL1 is significantly increased during osmotic and salinity stress, suggesting that this gene may be participating in abiotic stress responses including salinity stress. Previously, the involvement of AtCBL1 and AtCBL10 in response to salt stress has been well documented [17,35,51]. To test whether AcCBL1 contributes to salinity stress we took different approaches: (1) We performed seedling transfer experiment of AcCBL1 overexpressing lines and observed the salinity effect on primary root growth and percent root growth compared to wild-type plants; (2) tested the germination and growth phenotype under continuous salt and osmotic stress; (3) checked the performance of AcCBL1 overexpressing lines under saline irrigation and (4) performed the fresh weight analysis of seedling after 10 days (d) of transfer experiment. In agreement to our hypothesis, AcCBL1 overexpressing Arabidopsis lines showed better growth performance on 150 mM NaCl containing plates when compared to wild-type plants in terms of percent root growth (Figure 9a), and primary root growth (Figure 9b,c), after 2 d and 5 d of the transfer. The AcCBL1 overexpressing lines also showed better germination rates under salt and osmotic (300 mM mannitol) stress (Figure 9d). Consistently, fresh weights of the seedling grown on medium containing NaCl were significantly more compared to wild type plants (Figure 9e,g). Moreover, the saline irrigation of 15 d old plants with 250 mM NaCl for 10 d also showed better growth performance compared to wild-type (Figure 9f). Additionally, to investigate the contribution of AcCBL1 overexpressing lines to resistance to biotic stresses, we infected the mature rosette leaves with plant pathogenic fungus (Sclerotinia sclerotiorum) for 24 h and performed DAB staining. The DAB staining results indicated that the AcCBL1 lines had less accumulation and distribution of hydrogen peroxide (Figure 10b), indicating the increased tolerance of these plants under biotic stress. Collectively, the above results convincingly show that the AcCBL1 contributes to abiotic and biotic stress and AcCBL1 overexpressing plants display better performance under unfavorable stress conditions.

4. Discussion

Due to their sessile nature, plants face several biotic and abiotic stresses during their life cycle. Several genes contribute in defense and protective measure for plants to thrive well under hostile conditions [37,52,53]. CBL-CIPK gene family consists of plant-specific and calcium involving signaling modules indispensable to various stress signaling pathways. Due to availability of sequenced genome, CBL-CIPK gene families have been extensively studied in several plant species [19,38,39,40,41,42,54] however, the information about pineapple CBL-CIPK genes remain elusive. In the present study, a comprehensive search of CBL-CIPK genes from pineapple genome resulted in the identification of eight CBL and 21 CIPK genes, which were designated as AcCBL1 to AcCBL8 and AcCIPK1 to AcCIPK21 on the basis of their gene ID in the order of increasing number.
In plants, the CBL-CIPK network represents an excellent module for decoding Ca2+ signals during stress. This network is also a classic example of a diverged Ca2+ signal decoding system for sensing differential Ca2+ signals in numerous stress signaling pathways [42]. CBL proteins interact with different CIPK proteins depending on the cellular Ca2+ signals (concentration) that determine the precise outcome due to specific complex formation [55]. A 24 amino acid domain, referred to as NAF domain (due to the presence of conserved Asn-Ala-Phe), is specific and sufficient to mediate the interaction between CBL and CIPK [18]. We also utilized NAF domain as a query to identify CIPK proteins in pineapple genome.
Gene structure (intron-exon structures) in the gene families are often considered as imprints of evolution [56]. During the genome duplication, the deletion and insertion or both could happen to the genes. The intron–exon structure of the identified AcCBL genes showed that the CDS of all the AcCBLs were discontinuous by the presence of introns (Figure 1a). Whereas, AcCIPK genes can be grouped in intron less or intron poor and intron rich genes (Figure 1b). The intron less or intron-poor and intron-rich genes of CIPKs are conserved among other genomes such as Arabidopsis, maize, grapevine, and Brassica [19,38,40,42]. The intron numbers of AcCBL and AcCIPK genes ranged from 0 to 13 showing the large differences between the structures of the paralogs. This difference could be one of the attributes of functionality difference among the paralogs of CBL-CIPK genes. The phylogenetic analysis serves as an excellent method to analyze evolutionary relationships among genes [57]. In the present study, phylogenetic analysis of AcCBL and AcCIPK genes, together with Arabidopsis and rice, categorized AcCBLs mainly into five different groups (Figure 2) and CIPKs also into five different groups (Figure 3). Consistent with the current information of plant evolution, in the phylogenetic tree AcCBLs were more closely related to rice CBLs (monocot) compared to Arabidopsis (dicot). These findings indicate that they may have arisen together by gene duplication via a common ancestor and could have similar functions. It is interesting to note that the lower plants contain a lesser number of CBL-CIPK genes; whereas their number remarkably increases in higher land plants. The abundance of CBL-CIPK genes could be due to the increase in the complexity of the land plants during their evolution [58]. We found that the pineapple genome contains eight CBL genes (Table 1a) compared to 10 CBLs of Arabidopsis and rice whereas, 21 CIPK genes (Table 1b) compared to 26 CIPKs in Arabidopsis and 34 CIPKs in rice. The multiple CBL-CIPK genes present in the genome of the pineapple indicates the significance of the CBL-CIPK gene family in growth, development, and response to various stimuli [6].
Considering the importance of the CBL-CIPK genes in overall plant growth and response, we studied their expression profile in different tissues and development stages of pineapple. The expression profile suggests that AcCBL and AcCIPK genes express differentially in different tissues and stages (Figure 5). The hierarchical clustering of the gene expression data grouped them into three different categories (Figure 5). RNA-Seq data from different tissues and stages clearly indicated that the CBL-CIPK genes are involved in the regulation of different biological functions correlated to tissue or stage. For example, AcCIPK18 was highly expressed during the different stages of pineapple fruit development, suggesting that it may be involved in the regulation of pineapple fruit development (Figure 5).
Increasing evidence suggests that the CBL-CIPK toolkit performs a plethora of functions in plant growth and development under optimal and stress conditions [6,20,21,22,59]. In our study, we found that both the AcCBL1 and AcCBL7 were highly expressed in RNA-seq (Figure 5). Therefore, we selected AcCBL1 to further investigate its expression under different stress conditions. Real-time quantitative RT-PCR analysis showed that the transcript level of AcCBL1 significantly changed during all the stress treatments applied. The AcCBL1 exhibited significant changes in root and leaf at all the time points under salt stress, osmotic (mannitol), cold, and heat stress (Figure 6). The real-time result indicates that AcCBL1 could be playing a key role in adaption to extreme conditions as Arabidopsis CBL1 contributes to salt, cold, and osmotic stress [29]. Several studies indicate that CBLs contribute to salt stress [6,29,35]. In Arabidopsis, AtCBL4 with CIPK24, along with Na+/H+ antiporter SOS1, function in exclusion of Na+ from the cytoplasm [3,6]. To further validate the hypothesis, we generated transgenic Arabidopsis plants overexpressing AcCBL1 tagged with GFP. Consistently, the overexpression lines also show the increased localization of AcCBL1 in roots under salt stress (Figure 7 and Figure 8b). Previously, AtCBL1 was reported to be induced after wounding [16,54]. Consistent with these reports, we also found that AcCBL1 is associated with wounding. AcCBL1 protein was induced and translocated to the nucleus after wounding (Figure 8a). Moreover, the independent homozygous T3 transgenic lines (designated as Ox1, Ox2, Ox3) showed better performance under salt stress in terms of present root growth (Figure 9a), primary root growth (Figure 9b,c), germination (Figure 9d), fresh weight (Figure 9e), and overall growth phenotype (Figure 9f,g).
Additionally, fungal infection followed by DAB (3,3-diaminobenzidine) staining resulted in the better performance of AcCBL1 overexpressing transgenic plants compared to wild-type plants. In the presence of peroxidases, DAB is oxidized to generate a dark brown precipitate which is easily monitored to detect the existence and dispersal of hydrogen peroxide in cells [48]. The biotic stress would result in the generation of reactive oxygen species (ROS), resulting in oxidative stress to plants. After 24 h of fungal stress the transgenic plants displayed less damage to leaves (Figure 10b) and less accumulation of hydrogen peroxide (Figure 10c) than wild-type. In plants, the elevation in the cytosolic calcium (Ca2+) has been previously documented [60]. The AcCBL1 could be sensing the change in the cytosol after salt and wound stress, and relaying the signal downstream would be helping the AcCBL1 overexpressing plant to perform better under stress conditions.
Taken together, our results clearly indicate that ectopic expression of AcCBL1 enhances salt and biotic stress tolerance by facilitating scavenging of the Na+ and ROS from the cell.

5. Conclusions

In the present study, a comprehensive analysis of CBL-CIPK gene family was carried out in the economically important pineapple plant. A total of 8 CBL and 21 CIPK proteins were identified and characterized. The expression profile of pineapple CBL-CIPK genes suggest that AcCBL and AcCIPK proteins play key roles in the pineapple growth, development, and response to abiotic stresses. Further, the functional characterization of AcCBL1 shed light on its biological role and contribution during salt stress and wounding. The present study could be a foundation for further research and serve as a basis to utilize CBL-CIPK genes for development of next-generation crops.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-273X/9/7/293/s1, Figure S1: The distribution of conserved motifs in pineapple CBL proteins, Figure S2: The distribution of conserved motifs in pineapple CIPK proteins, Table S3: List of primers used in present study.

Author Contributions

M.A. and Y.Q. conceived and designed the research. M.A., B.F., B.H.J., performed the experiments. M.A., B.F., B.H.J., L.Z., S.C. and Y.C. analyzed data. M.A., B.F. and Y.Q. wrote and revised the paper. All of the authors read and approved the final manuscript.

Funding

This work was supported by National Natural Science Foundation of China (U1605212; 31761130074) and a Newton Advanced Fellowship (NA160391) to Y.Q.

Acknowledgments

We thank to Qin lab members for assistance in the experiments.

Conflicts of Interest

The authors declare that they have no competing interests.

Abbreviations

CBLCalcineurin B-like
CIPKCBL-interacting protein kinases
Ca2+Calcium ion
EREndoplasmic reticulum
CaMCalmodulin
CDPKsCalcium-dependent protein kinases
SOSSalt overlay sensitive
UTRUntranslated region
CDSCoding sequence
MEMEMultiple Em for Motif Elicitation
LGLinkage groups
DAB3,3-diaminobenzidine

References

  1. Dodd, A.N.; Kudla, J.; Sanders, D. The language of calcium signaling. Annu. Rev. Plant Biol. 2010, 61, 593–620. [Google Scholar] [CrossRef] [PubMed]
  2. Tang, R.J.; Luan, S. Regulation of calcium and magnesium homeostasis in plants: From transporters to signaling network. Curr. Opin. Plant Biol. 2017, 39, 97–105. [Google Scholar] [CrossRef] [PubMed]
  3. Zhu, J.K. Abiotic Stress Signaling and Responses in Plants. Cell 2016, 167, 313–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Costa, A.; Navazio, L.; Szabo, I. The contribution of organelles to plant intracellular Calcium signalling. J. Exp. Bot. 2018. [Google Scholar] [CrossRef] [PubMed]
  5. Tuteja, N.; Mahajan, S. Calcium signaling network in plants: An overview. Plant Signal. Behav. 2007, 2, 79–85. [Google Scholar] [CrossRef] [PubMed]
  6. Kudla, J.; Becker, D.; Grill, E.; Hedrich, R.; Hippler, M.; Kummer, U.; Parniske, M.; Romeis, T.; Schumacher, K. Advances and current challenges in calcium signaling. New Phytol. 2018, 218, 414–431. [Google Scholar] [CrossRef] [PubMed]
  7. Ding, Y.; Shi, Y.; Yang, S. Advances and challenges in uncovering cold tolerance regulatory mechanisms in plants. New Phytol. 2019, 222, 1690–1704. [Google Scholar] [CrossRef] [Green Version]
  8. Behera, S.; Zhaolong, X.; Luoni, L.; Bonza, M.C.; Doccula, F.G.; De Michelis, M.I.; Morris, R.J.; Schwarzlander, M.; Costa, A. Cellular Ca(2+) Signals Generate Defined pH Signatures in Plants. Plant Cell. 2018, 30, 2704–2719. [Google Scholar] [CrossRef]
  9. Meena, M.K.; Prajapati, R.; Krishna, D.; Divakaran, K.; Pandey, Y.; Reichelt, M.; Mathew, M.K.; Boland, W.; Mithofer, A.; Vadassery, J. The Ca2+ Channel CNGC19 Regulates Arabidopsis defense against Spodoptera Herbivory. Plant Cell. 2019. [Google Scholar] [CrossRef]
  10. Van Loon, L.C. The Intelligent Behavior of Plants. Trend. Plant Sci. 2016, 21, 286–294. [Google Scholar] [CrossRef]
  11. De Vriese, K.; Himschoot, E.; Dünser, K.; Nguyen, L.; Drozdzecki, A.; Costa, A.; Nowack, M.K.; Kleine-Vehn, J.; Audenaert, D.; Beeckman, T.; et al. Identification of Novel Inhibitors of Auxin-Induced Ca2+ Signaling via a Plant-Based Chemical Screen. Plant Physiol. 2019, 180, 480–496. [Google Scholar] [CrossRef] [PubMed]
  12. Blatt, M.R.; Grabov, A. Signal redundancy, gates and integration in the control of ion channels for stomatal movement. J. Exp. Bot. 1997, 48, 529–537. [Google Scholar] [CrossRef] [PubMed]
  13. Michard, E.; Simon, A.A.; Tavares, B.; Wudick, M.M.; Feijo, J.A. Signaling with Ions: The Keystone for Apical Cell Growth and Morphogenesis in Pollen Tubes. Plant Physiol. 2017, 173, 91–111. [Google Scholar] [CrossRef] [PubMed]
  14. Bender, K.W.; Zielinski, R.E.; Huber, S.C. Revisiting paradigms of Ca(2+) signaling protein kinase regulation in plants. Biochem. J. 2018, 475, 207–223. [Google Scholar] [CrossRef] [PubMed]
  15. 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]
  16. Kudla, J.; Xu, Q.; Harter, K.; Gruissem, W.; Luan, S. Genes for calcineurin B-like proteins in Arabidopsis are differentially regulated by stress signals. Proc. Natl. Acad. Sci. USA 1999, 96, 4718–4723. [Google Scholar] [CrossRef]
  17. Liu, J.; Zhu, J.K. A calcium sensor homolog required for plant salt tolerance. Science 1998, 280, 1943–1945. [Google Scholar] [CrossRef]
  18. Albrecht, V.; Ritz, O.; Linder, S.; Harter, K.; Kudla, J. The NAF domain defines a novel protein-protein interaction module conserved in Ca2+-regulated kinases. EMBO J. 2001, 20, 1051–1063. [Google Scholar] [CrossRef] [Green Version]
  19. 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]
  20. Zhang, Y.; Lv, Y.; Jahan, N.; Chen, G.; Ren, D.; Guo, L. Sensing of Abiotic Stress and Ionic Stress Responses in Plants. Int. J. Mol. Sci. 2018, 19. [Google Scholar] [CrossRef]
  21. Xuan, Y.H.; Kumar, V.; Han, X.; Kim, S.H.; Jeong, J.H.; Kim, C.M.; Gao, Y.; Han, C.D. CBL-INTERACTING PROTEIN KINASE 9 regulates ammonium-dependent root growth downstream of IDD10 in rice (Oryza sativa). Ann. Bot. 2019. [Google Scholar] [CrossRef]
  22. Yang, Y.; Wu, Y.; Ma, L.; Yang, Z.; Dong, Q.; Li, Q.; Ni, X.; Kudla, J.; Song, C.-P.; Guo, Y. The Ca2+ sensor SCaBP3/CBL7 Fine Tunes Arabidopsis Alkali Tolerance and Modulats Plasma Membrane H+-ATPase Activity. Plant Cell 2019. [Google Scholar] [CrossRef] [PubMed]
  23. Li, Z.Y.; Xu, Z.S.; He, G.Y.; Yang, G.X.; Chen, M.; Li, L.C.; Ma, Y. The voltage-dependent anion channel 1 (AtVDAC1) negatively regulates plant cold responses during germination and seedling development in Arabidopsis and interacts with calcium sensor CBL1. Int. J. Mol. Sci. 2013, 14, 701–713. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, C.; Ding, S.; Zhang, H.; Du, H.; An, L. CIPK7 is involved in cold response by interacting with CBL1 in Arabidopsis thaliana. Plant Sci. 2011, 181, 57–64. [Google Scholar] [CrossRef] [PubMed]
  25. Ligaba-Osena, A.; Fei, Z.; Liu, J.; Xu, Y.; Shaff, J.; Lee, S.C.; Luan, S.; Kudla, J.; Kochian, L.; Pineros, M. Loss-of-function mutation of the calcium sensor CBL1 increases aluminum sensitivity in Arabidopsis. New Phytol. 2017, 214, 830–841. [Google Scholar] [CrossRef]
  26. Lan, W.Z.; Lee, S.C.; Che, Y.F.; Jiang, Y.Q.; Luan, S. Mechanistic analysis of AKT1 regulation by the CBL-CIPK-PP2CA interactions. Mol. Plant 2011, 4, 527–536. [Google Scholar] [CrossRef] [PubMed]
  27. 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]
  28. 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]
  29. Cheong, Y.H.; Pandey, G.K.; Grant, J.J.; Batistic, O.; Li, L.; Kim, B.G.; Lee, S.C.; Kudla, J.; Luan, S. Two calcineurin B-like calcium sensors, interacting with protein kinase CIPK23, regulate leaf transpiration and root potassium uptake in Arabidopsis. Plant J. 2007, 52, 223–239. [Google Scholar] [CrossRef]
  30. 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]
  31. Liu, L.L.; Ren, H.M.; Chen, L.Q.; Wang, Y.; Wu, W.H. A protein kinase, calcineurin B-like protein-interacting protein Kinase9, interacts with calcium sensor calcineurin B-like Protein3 and regulates potassium homeostasis under low-potassium stress in Arabidopsis. Plant Physiol. 2013, 161, 266–277. [Google Scholar] [CrossRef] [PubMed]
  32. Held, K.; Pascaud, F.; Eckert, C.; Gajdanowicz, P.; Hashimoto, K.; Corratge-Faillie, C.; Offenborn, J.N.; Lacombe, B.; Dreyer, I.; Thibaud, J.B.; et al. Calcium-dependent modulation and plasma membrane targeting of the AKT2 potassium channel by the CBL4/CIPK6 calcium sensor/protein kinase complex. Cell Res. 2011, 21, 1116–1130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yang, Y.; Guo, Y. Elucidating the molecular mechanisms mediating plant salt-stress responses. New Phytol. 2018, 217, 523–539. [Google Scholar] [CrossRef] [PubMed]
  34. Lin, H.; Yang, Y.; Quan, R.; Mendoza, I.; Wu, Y.; Du, W.; Zhao, S.; Schumaker, K.S.; Pardo, J.M.; Guo, Y. Phosphorylation of SOS3-LIKE CALCIUM BINDING PROTEIN8 by SOS2 protein kinase stabilizes their protein complex and regulates salt tolerance in Arabidopsis. Plant Cell. 2009, 21, 1607–1619. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, B.G.; Waadt, R.; Cheong, Y.H.; Pandey, G.K.; Dominguez-Solis, J.R.; Schultke, S.; Lee, S.C.; Kudla, J.; Luan, S. The calcium sensor CBL10 mediates salt tolerance by regulating ion homeostasis in Arabidopsis. Plant J. 2007, 52, 473–484. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, J.; Liu, J.; Ming, R. Genomic analyses of the CAM plant pineapple. J. Exp. Bot. 2014, 65, 3395–3404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Xie, T.; Chen, C.; Li, C.; Liu, J.; Liu, C.; He, Y. Genome-wide investigation of WRKY gene family in pineapple: Evolution and expression profiles during development and stress. BMC Genom. 2018, 19, 490. [Google Scholar] [CrossRef]
  38. Chen, X.; Gu, Z.; Xin, D.; Hao, L.; Liu, C.; Huang, J.; Ma, B.; Zhang, H. Identification and characterization of putative CIPK genes in maize. J. Genet. Genom. 2011, 38, 77–87. [Google Scholar] [CrossRef]
  39. Li, J.; Jiang, M.M.; Ren, L.; Liu, Y.; Chen, H.Y. Identification and characterization of CBL and CIPK gene families in eggplant (Solanum melongena L.). Mol. Genet. Genom. 2016, 291, 1769–1781. [Google Scholar] [CrossRef]
  40. 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]
  41. Zhang, H.; Yin, W.; Xia, X. Calcineurin B-Like family in Populus: Comparative genome analysis and expression pattern under cold, drought and salt stress treatment. Plant Growth Regul. 2008, 56, 129–140. [Google Scholar] [CrossRef]
  42. Yin, X.; Wang, Q.; Chen, Q.; Xiang, N.; Yang, Y.; 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]
  43. Ming, R.; Wai, C.M.; Guyot, R. Pineapple Genome: A Reference for Monocots and CAM Photosynthesis. Trends Genet. 2016, 32, 690–696. [Google Scholar] [CrossRef] [PubMed]
  44. Priyadarshani, S.; Hu, B.; Li, W.; Ali, H.; Jia, H.; Zhao, L.; Ojolo, S.P.; Azam, S.M.; Xiong, J.; Yan, M.; et al. Simple protoplast isolation system for gene expression and protein interaction studies in pineapple (Ananas comosus L.). Plant Method. 2018, 14, 95. [Google Scholar] [CrossRef] [PubMed]
  45. Chen, P.; Li, Y.; Zhao, L.; Hou, Z.; Yan, M.; Hu, B.; Liu, Y.; Azam, S.M.; Zhang, Z.; Rahman, Z.U.; et al. Genome-Wide Identification and Expression Profiling of ATP-Binding Cassette (ABC) Transporter Gene Family in Pineapple (Ananas comosus (L.) Merr.) Reveal the Role of AcABCG38 in Pollen Development. Front. Plant Sci. 2017, 8, 2150. [Google Scholar] [CrossRef] [PubMed]
  46. Baskin, T.I.; Wilson, J.E. Inhibitors of protein kinases and phosphatases alter root morphology and disorganize cortical microtubules. Plant Physiol. 1997, 113, 493–502. [Google Scholar] [CrossRef] [PubMed]
  47. Rahman, A. Auxin: A regulator of cold stress response. Physiol. Plant 2013, 147, 28–35. [Google Scholar] [CrossRef]
  48. Daudi, A.; O’Brien, J.A. Detection of Hydrogen Peroxide by DAB Staining in Arabidopsis Leaves. Bio Protoc. 2012, 2. [Google Scholar] [CrossRef]
  49. 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]
  50. Yadav, A.K.; Jha, S.K.; Sanyal, S.K.; Luan, S.; Pandey, G.K. Arabidopsis calcineurin B-like proteins differentially regulate phosphorylation activity of CBL-interacting protein kinase 9. Biochem. J. 2018, 475, 2621–2636. [Google Scholar] [CrossRef]
  51. Manik, S.M.; Shi, S.; Mao, J.; Dong, L.; Su, Y.; Wang, Q.; Liu, H. The Calcium Sensor CBL-CIPK Is Involved in Plant’s Response to Abiotic Stresses. Int. J. Genom. 2015, 2015, 493191. [Google Scholar] [CrossRef] [PubMed]
  52. Liu, C.; Xie, T.; Chen, C.; Luan, A.; Long, J.; Li, C.; Ding, Y.; He, Y. Genome-wide organization and expression profiling of the R2R3-MYB transcription factor family in pineapple (Ananas comosus). BMC Genom. 2017, 18, 503. [Google Scholar] [CrossRef] [PubMed]
  53. Tak, H.; Negi, S.; Ganapathi, T.R. Banana NAC transcription factor MusaNAC042 is positively associated with drought and salinity tolerance. Protoplasma 2017, 254, 803–816. [Google Scholar] [CrossRef] [PubMed]
  54. Niu, L.; Dong, B.; Song, Z.; Meng, D.; Fu, Y. Genome-Wide Identification and Characterization of CIPK Family and Analysis Responses to Various Stresses in Apple (Malus domestica). Int. J. Mol. Sci. 2018, 19. [Google Scholar] [CrossRef] [PubMed]
  55. Zhu, J.K. Regulation of ion homeostasis under salt stress. Curr. Opin. Plant Biol. 2003, 6, 441–445. [Google Scholar] [CrossRef]
  56. Boudet, N.; Aubourg, S.; Toffano-Nioche, C.; Kreis, M.; Lecharny, A. Evolution of intron/exon structure of DEAD helicase family genes in Arabidopsis, Caenorhabditis, and Drosophila. Genome Res. 2001, 11, 2101–2114. [Google Scholar] [CrossRef] [PubMed]
  57. 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]
  58. Edel, K.H.; Marchadier, E.; Brownlee, C.; Kudla, J.; Hetherington, A.M. The Evolution of Calcium-Based Signalling in Plants. Curr. Biol. 2017, 27, R667–R679. [Google Scholar] [CrossRef]
  59. Lyzenga, W.J.; Sullivan, V.; Liu, H.; Stone, S.L. The Kinase Activity of Calcineurin B-like Interacting Protein Kinase 26 (CIPK26) Influences Its Own Stability and that of the ABA-regulated Ubiquitin Ligase, Keep on Going (KEG). Front. Plant Sci. 2017, 8, 502. [Google Scholar] [CrossRef] [Green Version]
  60. Mousavi, S.A.; Chauvin, A.; Pascaud, F.; Kellenberger, S.; Farmer, E.E. GLUTAMATE RECEPTOR-LIKE genes mediate leaf-to-leaf wound signalling. Nature 2013, 500, 422–426. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Exon-intron structure of pineapple CBL-CIPK genes. (a) Exon-intron structure of pineapple CBL genes (b) exon-intron structure of pineapple CIPK genes; green boxes indicate untranslated 5′- and 3′-regions; red boxes indicate exons; black lines indicate introns. Prefix ‘Ac’ indicates Ananas comosus.
Figure 1. Exon-intron structure of pineapple CBL-CIPK genes. (a) Exon-intron structure of pineapple CBL genes (b) exon-intron structure of pineapple CIPK genes; green boxes indicate untranslated 5′- and 3′-regions; red boxes indicate exons; black lines indicate introns. Prefix ‘Ac’ indicates Ananas comosus.
Biomolecules 09 00293 g001
Figure 2. Phylogenetic tree depicting the relationships among CBL proteins from pineapple, rice, and Arabidopsis. The different colored arcs indicate different subgroups. Prefix ‘Ath’, Osa, and ‘Ac’ indicate CBL proteins from Arabidopsis, Oryza sativa, and Ananas comosus respectively.
Figure 2. Phylogenetic tree depicting the relationships among CBL proteins from pineapple, rice, and Arabidopsis. The different colored arcs indicate different subgroups. Prefix ‘Ath’, Osa, and ‘Ac’ indicate CBL proteins from Arabidopsis, Oryza sativa, and Ananas comosus respectively.
Biomolecules 09 00293 g002
Figure 3. Phylogenetic tree depicting the relationships among CIPK proteins from pineapple, rice, and Arabidopsis. The different colored arcs indicate different subgroups. Prefix ‘Ath’, Osa, and ‘Ac’ indicate CBL proteins from Arabidopsis, Oryza sativa, and Ananas comosus respectively.
Figure 3. Phylogenetic tree depicting the relationships among CIPK proteins from pineapple, rice, and Arabidopsis. The different colored arcs indicate different subgroups. Prefix ‘Ath’, Osa, and ‘Ac’ indicate CBL proteins from Arabidopsis, Oryza sativa, and Ananas comosus respectively.
Biomolecules 09 00293 g003
Figure 4. Chromosomal locations of pineapple CBL and CIPK genes. The 8 CBL and 21CIPK genes of pineapple were mapped to different chromosomes using MapChart.
Figure 4. Chromosomal locations of pineapple CBL and CIPK genes. The 8 CBL and 21CIPK genes of pineapple were mapped to different chromosomes using MapChart.
Biomolecules 09 00293 g004
Figure 5. Expression profiles of the pineapple CBL and CIPK genes. Hierarchical clustering of expression profiles of pineapple CBL-CIPK genes in different organs and developmental stages. Red color indicates high levels of transcript abundance, and blue color indicates low transcript abundance. The color scale is shown at right side of the figure. Sample details are mentioned at the bottom of each lane: sepal Se1–Se4, gynoecium Gy1- Gy7, ovule Ov1–Ov7, petal Pe1–Pe3, stamen St1–St6, fruit ‘Fr_S1–Fr_S7’, root ‘Ro’, leaf ‘Le’, and flower ‘Fl’. Where ‘S’ is the abbreviation for ‘stage’. Pineapple CBL gene names are represented in green color.
Figure 5. Expression profiles of the pineapple CBL and CIPK genes. Hierarchical clustering of expression profiles of pineapple CBL-CIPK genes in different organs and developmental stages. Red color indicates high levels of transcript abundance, and blue color indicates low transcript abundance. The color scale is shown at right side of the figure. Sample details are mentioned at the bottom of each lane: sepal Se1–Se4, gynoecium Gy1- Gy7, ovule Ov1–Ov7, petal Pe1–Pe3, stamen St1–St6, fruit ‘Fr_S1–Fr_S7’, root ‘Ro’, leaf ‘Le’, and flower ‘Fl’. Where ‘S’ is the abbreviation for ‘stage’. Pineapple CBL gene names are represented in green color.
Biomolecules 09 00293 g005
Figure 6. Expression profiles of the pineapple CBL1 gene in response to different stress treatments. qRT-PCR of CBL1 gene in at different time point (0 h, 24 h, 48 h, and 72 h) and samples (leaf and root) after (a) salt (NaCl 150 mM), (b) cold (4 °C) stress, (c) osmotic (Mnnitol 350 mM) stress, and (d) heat (45 °C) stress. Data were normalized to EF1a gene. Vertical bars indicate SE and *** indicates significantly different values compared to control treatment (p < 0.001). All experiments were performed with three technical and three biological repeats.
Figure 6. Expression profiles of the pineapple CBL1 gene in response to different stress treatments. qRT-PCR of CBL1 gene in at different time point (0 h, 24 h, 48 h, and 72 h) and samples (leaf and root) after (a) salt (NaCl 150 mM), (b) cold (4 °C) stress, (c) osmotic (Mnnitol 350 mM) stress, and (d) heat (45 °C) stress. Data were normalized to EF1a gene. Vertical bars indicate SE and *** indicates significantly different values compared to control treatment (p < 0.001). All experiments were performed with three technical and three biological repeats.
Biomolecules 09 00293 g006
Figure 7. AcCBL1-GFP localizes in cytosol under salt stress. Confocal analyses show AcCBL1-GFP localization in the cytosol when of 5-day-old seedlings grown in the standard conditions were transferred in a medium containing 150 mM NaCl. GFP fluorescence is represented in green and propidium iodide (PI)fluorescence channel is represented in red. Scale bars = 50 μm.
Figure 7. AcCBL1-GFP localizes in cytosol under salt stress. Confocal analyses show AcCBL1-GFP localization in the cytosol when of 5-day-old seedlings grown in the standard conditions were transferred in a medium containing 150 mM NaCl. GFP fluorescence is represented in green and propidium iodide (PI)fluorescence channel is represented in red. Scale bars = 50 μm.
Biomolecules 09 00293 g007
Figure 8. AcCBL1-GFP translocalizes to nucleus from cytosol after injury. (a) Localization AcCBL1-GFP in the nucleus when the 7-day-old root of transgenic seedlings is injured. (b) AcCBL1-GFP localization in cytosol after salt stress. GFP fluorescence is represented in green and PI fluorescence channel is represented in red. Scale bars = 50 μm.
Figure 8. AcCBL1-GFP translocalizes to nucleus from cytosol after injury. (a) Localization AcCBL1-GFP in the nucleus when the 7-day-old root of transgenic seedlings is injured. (b) AcCBL1-GFP localization in cytosol after salt stress. GFP fluorescence is represented in green and PI fluorescence channel is represented in red. Scale bars = 50 μm.
Biomolecules 09 00293 g008
Figure 9. Overexpression of AcCBL1 results in resistance to salinity and osmotic stress. (a) Percent root growth in 5-day-old wild-type and AcCBL1 overexpressing seedlings transferred to new plates supplemented with 150 mm NaCl and analyzed after 2 days and 5 days of transfer. (b) Primary root growth after 2 days and (c) primary root growth after 5 days and (d) phenotype of 7-day-old AcCBL1 overexpressing seedlings germinated under salinity (150 mM NaCl) and osmotic stress (300 mM mannitol). (e) Fresh weight of 5-day-old wild-type and AcCBL1 overexpressing lines grown for an additional 10 days on 150 mM NaCl after transfer to new plate and analysis. (f) 15-day-old seedlings were subjected to saline irrigation with 250 mM NaCl for an additional 10 days and then photographed. (g) Seedlings grown for 5 days and transferred to new plate with (lower panel), or without (upper panel), 150 mM NaCl for 10 additional days and photographed. Vertical bars represent ± S.E., and *** indicates significantly different values between treatments (p < 0.001).
Figure 9. Overexpression of AcCBL1 results in resistance to salinity and osmotic stress. (a) Percent root growth in 5-day-old wild-type and AcCBL1 overexpressing seedlings transferred to new plates supplemented with 150 mm NaCl and analyzed after 2 days and 5 days of transfer. (b) Primary root growth after 2 days and (c) primary root growth after 5 days and (d) phenotype of 7-day-old AcCBL1 overexpressing seedlings germinated under salinity (150 mM NaCl) and osmotic stress (300 mM mannitol). (e) Fresh weight of 5-day-old wild-type and AcCBL1 overexpressing lines grown for an additional 10 days on 150 mM NaCl after transfer to new plate and analysis. (f) 15-day-old seedlings were subjected to saline irrigation with 250 mM NaCl for an additional 10 days and then photographed. (g) Seedlings grown for 5 days and transferred to new plate with (lower panel), or without (upper panel), 150 mM NaCl for 10 additional days and photographed. Vertical bars represent ± S.E., and *** indicates significantly different values between treatments (p < 0.001).
Biomolecules 09 00293 g009
Figure 10. Overexpression of AcCBL1 results in resistance to biotic stress. (a) Morphology of four-week-old wild-type (Col-0) and AcCBL1 overexpressing plant. (b) Damage in leaves by exposure to fungus (Sclerotinia sclerotiorum) to mature rosette leaves of wild-type and 3 different lines (Ox1, Ox2, and Ox3) of AcCBL1 overexpressing plants for 24 h. (c) Detection of H2O2 by DAB staining in leaves of wild-type (Col-0) and AcCBL1 overexpressing lines after 24 h infection of fungus (Sclerotinia sclerotiorum). The brown color reflects the reactive oxygen species (ROS) accumulation after the fungus treatment. Scale bar = 1 cm.
Figure 10. Overexpression of AcCBL1 results in resistance to biotic stress. (a) Morphology of four-week-old wild-type (Col-0) and AcCBL1 overexpressing plant. (b) Damage in leaves by exposure to fungus (Sclerotinia sclerotiorum) to mature rosette leaves of wild-type and 3 different lines (Ox1, Ox2, and Ox3) of AcCBL1 overexpressing plants for 24 h. (c) Detection of H2O2 by DAB staining in leaves of wild-type (Col-0) and AcCBL1 overexpressing lines after 24 h infection of fungus (Sclerotinia sclerotiorum). The brown color reflects the reactive oxygen species (ROS) accumulation after the fungus treatment. Scale bar = 1 cm.
Biomolecules 09 00293 g010
Table 1. Characteristics of calcineurin B-like (CBL) and CBL-interacting protein kinases (CIPK) genes identified from pineapple. Prefix ‘Ac’ Indicates Ananas comosus.
Table 1. Characteristics of calcineurin B-like (CBL) and CBL-interacting protein kinases (CIPK) genes identified from pineapple. Prefix ‘Ac’ Indicates Ananas comosus.
(a) Characteristics of Pineapple CBL Gene
NameTranscript IDChromosome LocationLengthMWpI
AcCBL1Aco000875LG02:16012605-1601840322625,958.744.81
AcCBL2Aco004447LG05:3255687-325866821324,478.924.81
AcCBL3Aco005668LG11:12099528-1210304721624,888.394.86
AcCBL4Aco005877LG16:10709202-1071792735940,867.514.68
AcCBL5Aco011645LG08:12541062-1254569323426,797.474.72
AcCBL6Aco012888LG03:15510233-1551637919321,928.834.52
AcCBL7Aco015301LG20:10385809-1039110922425,763.524.81
AcCBL8Aco024403LG12:9915463-992754625228,939.196.16
(b) Characteristics of Pineapple CIPK Genes
AcCIPK1Aco001621LG18:9219002-922061240943,926.739.35
AcCIPK2Aco001625LG18:9192567-919389544250,063.849.26
AcCIPK3Aco002215LG04:3262495-326816345050,565.198.78
AcCIPK4Aco004322LG05:2160825-216224347252,996.596.69
AcCIPK5Aco004323LG05:2169056-217082348754,823.098.75
AcCIPK6Aco005225LG07:3071294-307266740542,654.619.99
AcCIPK7Aco006700LG01:22600010-2261024538544,053.777.12
AcCIPK8Aco007507LG12:11547144-1155834738243,701.046.87
AcCIPK9Aco008201LG19:8409312-841120042445,918.028.61
AcCIPK10Aco010015LG10:2187779-218958147854,071.926.51
AcCIPK11Aco010017LG10:2198749-220014946652,913.089.06
AcCIPK12Aco011132LG04:13486171-1348753845549,351.156.94
AcCIPK13Aco012533LG13:1862440-186374443449,113.689.01
AcCIPK14Aco014260LG05:45295-6006442247,747.217.71
AcCIPK15Aco015008LG04:11703020-1171039644950,875.828.01
AcCIPK16Aco015525LG03:12649272-1265417239645,072.728.99
AcCIPK17Aco016931LG02:2099383-210147345150,452.218.09
AcCIPK18Aco016932LG02:2036912-203829145950,213.689.08
AcCIPK19Aco019253LG13:3471010-347847744150,378.767.2
AcCIPK20Aco021653LG16:7092847-709589250655,836.966.63
AcCIPK21Aco022003LG04:312075-31654544449,867.096.87

Share and Cite

MDPI and ACS Style

Aslam, M.; Fakher, B.; Jakada, B.H.; Zhao, L.; Cao, S.; Cheng, Y.; Qin, Y. Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response. Biomolecules 2019, 9, 293. https://doi.org/10.3390/biom9070293

AMA Style

Aslam M, Fakher B, Jakada BH, Zhao L, Cao S, Cheng Y, Qin Y. Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response. Biomolecules. 2019; 9(7):293. https://doi.org/10.3390/biom9070293

Chicago/Turabian Style

Aslam, Mohammad, Beenish Fakher, Bello Hassan Jakada, Lihua Zhao, Shijiang Cao, Yan Cheng, and Yuan Qin. 2019. "Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response" Biomolecules 9, no. 7: 293. https://doi.org/10.3390/biom9070293

APA Style

Aslam, M., Fakher, B., Jakada, B. H., Zhao, L., Cao, S., Cheng, Y., & Qin, Y. (2019). Genome-Wide Identification and Expression Profiling of CBL-CIPK Gene Family in Pineapple (Ananas comosus) and the Role of AcCBL1 in Abiotic and Biotic Stress Response. Biomolecules, 9(7), 293. https://doi.org/10.3390/biom9070293

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