Identification and comparative analysis of the CIPK gene family and characterization of the cold stress response in the woody plant Prunus mume

Understanding the salt overly sensitive pathway in Prunus: Identification and characterization of NHX, CIPK, and CBL genes

Salinity is a major abiotic stress factor that can significantly impact crop growth, and productivity. In response to salt stress, the plant Salt Overly Sensitive (SOS) signaling pathway regulates the homeostasis of intracellular sodium ion concentration. The SOS1, SOS2, and SOS3 genes play critical roles in the SOS pathway, which belongs to the members of Na+/H+ exchanger (NHX), CBL-interacting protein kinase (CIPK), and calcineurin B-like (CBL) gene families, respectively. In this study, we performed genome-wide identifications and phylogenetic analyses of NHX, CIPK, and CBL genes in six Rosaceae species: Prunus persica, Prunus dulcis, Prunus mume, Prunus armeniaca, Pyrus ussuriensis × Pyrus communis, and Rosa chinensis. NHX, CIPK, and CBL genes of Arabidopsis thaliana were used as controls for phylogenetic analyses. Our analysis revealed the lineage-specific and adaptive evolutions of Rosaceae genes. Our observations indicated the existence of two primary classes of CIPK genes: those that are intron-rich and those that are intron-less. Intron-rich CIPKs in Rosaceae and Arabidopsis can be traced back to algae CIPKs and CIPKs found in early plants, suggesting that intron-less CIPKs evolved from their intron-rich counterparts. This study identified one gene for each member of the SOS signaling pathway in P. persica: PpSOS1, PpSOS2, and PpSOS3. Gene expression analyses indicated that all three genes of P. persica were expressed in roots and leaves. Yeast two-hybrid-based protein-protein interaction analyses revealed a direct interaction between PpSOS3 and PpSOS2; and between PpSOS2 and PpSOS1C-terminus region. Our findings indicate that the SOS signaling pathway is highly conserved in P. persica.

AtCIPK16, a CBL-interacting protein kinase gene, confers salinity tolerance in transgenic wheat

Globally, wheat is the major source of staple food, protein, and basic calories for most of the human population. Strategies must be adopted for sustainable wheat crop production to fill the ever-increasing food demand. Salinity is one of the major abiotic stresses involved in plant growth retardation and grain yield reduction. In plants, calcineurin-B-like proteins form a complicated network with the target kinase CBL-interacting protein kinases (CIPKs) in response to intracellular calcium signaling as a consequence of abiotic stresses. The AtCIPK16 gene has been identified in Arabidopsis thaliana and found to be significantly upregulated under salinity stress. In this study, the AtCIPK16 gene was cloned in two different plant expression vectors, i.e., pTOOL37 having a UBI1 promoter and pMDC32 having a 2XCaMV35S constitutive promoter transformed through the Agrobacterium-mediated transformation protocol, in the local wheat cultivar Faisalabad-2008. Based on their ability to tolerate different levels of salt stress (0, 50, 100, and 200 mM), the transgenic wheat lines OE1, OE2, and OE3 expressing AtCIPK16 under the UBI1 promoter and OE5, OE6, and OE7 expressing the same gene under the 2XCaMV35S promoter performed better at 100 mM of salinity stress as compared with the wild type. The AtCIPK16 overexpressing transgenic wheat lines were further investigated for their K⁺ retention ability in root tissues by utilizing the microelectrode ion flux estimation technique. It has been demonstrated that after 10 min of 100 mM NaCl application, more K⁺ ions were retained in the AtCIPK16 overexpressing transgenic wheat lines than in the wild type. Moreover, it could be concluded that AtCIPK16 functions as a positive elicitor in sequestering Na⁺ ions into the cell vacuole and retaining more cellular K⁺ under salt stress to maintain ionic homeostasis.

Molecular Mechanism of Cold Tolerance of Centipedegrass Based on the Transcriptome

Low temperature is an important limiting factor in the environment that affects the distribution, growth and development of warm-season grasses. Transcriptome sequencing has been widely used to mine candidate genes under low-temperature stress and other abiotic stresses. However, the molecular mechanism of centipedegrass in response to low-temperature stress was rarely reported. To understand the molecular mechanism of centipedegrass in response to low-temperature stress, we measured physiological indicators and sequenced the transcriptome of centipedegrass under different stress durations. Under cold stress, the SS content and APX activity of centipedegrass increased while the SOD activity decreased; the CAT activity, POD activity and flavonoid content first increased and then decreased; and the GSH-Px activity first decreased and then increased. Using full-length transcriptome and second-generation sequencing, we obtained 38.76 G subreads. These reads were integrated into 177,178 isoforms, and 885 differentially expressed transcripts were obtained. The expression of AUX_IAA and WRKY transcription factors and HSF transcription-influencing factors increased during cold stress. Through KEGG enrichment analysis, we determined that arginine and proline metabolism, plant circadian rhythm, plant hormone signal transduction and the flavonoid biosynthesis pathways played important roles in the cold stress resistance of centipedegrass. In addition, by using weighted gene coexpression network analysis (WGCNA), we determined that the turquoise module was significantly correlated with SS content and APX activity, while the blue module was significantly negatively correlated with POD and CAT activity. This paper is the first to report the response of centipedegrass to cold stress at the transcriptome level. Our results help to clarify the molecular mechanisms underlying the cold tolerance of warm-season grasses.

Diverse roles of the CIPK gene family in transcription regulation and various biotic and abiotic stresses: A literature review and bibliometric study

CIPKs are a subclass of serine/threonine (Ser/Thr) protein kinases. CBLs are ubiquitous Ca2+ sensors that interact with CIPK with the aid of secondary Ca2+ messengers for regulation of growth and development and response to stresses faced by plants. The divergent roles of the CIPK-CBL interaction in plants include responding to environmental stresses (salt, cold, drought, pH, ABA signaling, and ion toxicity), ion homeostasis (K+, NH4 +, NO3 -, and microelement homeostasis), biotic stress, and plant development. Each member of this gene family produces distinct proteins that help plants adapt to diverse stresses or stimuli by interacting with calcium ion signals. CIPK consists of two structural domains-an N-terminal domain and a C-terminal domain-connected by a junction domain. The N-terminal domain, the site of phosphorylation, is also called the activation domain and kinase domain. The C-terminal, also known as the regulatory domain of CIPK, further comprises NAF/FISL and PPI. CBL comprises four EF domains and conserved PFPF motifs and is the site of binding with the NAF/FISL domain of CIPK to form a CBL-CIPK complex. In addition, we also performed a bibliometric analysis of the CIPK gene family of data extracted from the WoSCC. A total of 95 documents were retrieved, which had been published by 47 sources. The production over time was zigzagged. The top key terms were gene, CIPK, abiotic stress, and gene expression. Beijing Forestry University was the top affiliation, while The Plant Cell was the top source. The genomics and metabolomics of this gene family require more study.

Genome-wide transcriptome profiling revealed biological macromolecules respond to low temperature stress in Brassica napus L

Brassica napus L. (B. napus) is a vital oilseed crop cultivated worldwide; low temperature (LT) is one of the major stress factors that limit its growth, development, distribution, and production. Even though processes have been developed to characterize LT-responsive genes, only limited studies have exploited the molecular response mechanisms in B. napus. Here the transcriptome data of an elite B. napus variety with LT adaptability was acquired and applied to investigate the gene expression profiles of B. napus in response to LT stress. The bioinformatics study revealed a total of 79,061 unigenes, of which 3,703 genes were differentially expressed genes (DEGs), with 2,129 upregulated and 1,574 downregulated. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis pinpointed that the DEGs were enriched in LT-stress-responsive biological functions and metabolic pathways, which included sugar metabolism, antioxidant defense system, plant hormone signal transduction, and photosynthesis. Moreover, a group of LT-stress-responsive transcription factors with divergent expression patterns under LT was summarized. A combined protein interaction suggested that a complex interconnected regulatory network existed in all detected pathways. RNA-seq data was verified using real-time quantitative polymerase chain reaction analysis. Based on these findings, we presented a hypothesis model illustrating valuable information for understanding the LT response mechanisms in B. napus.

Sequence Characteristics and Expression Analysis of GhCIPK23 Gene in Upland Cotton (Gossypium hirsutum L.)

CIPK (calcineurin B-like-interacting protein kinase) is a kind of serine/threonine protein kinase widely existing in plants, and it plays an important role in plant growth and development and stress response. To better understand the biological functions of the GhCIPK23 gene in upland cotton, the coding sequence (CDS) of the GhCIPK23 gene was cloned in upland cotton, and its protein sequence, evolutionary relationship, subcellular localization, expression pattern and cis-acting elements in the promoter region were analyzed. Our results showed that the full-length CDS of GhCIPK23 was 1368 bp, encoding a protein with 455 amino acids. The molecular weight and isoelectric point of this protein were 50.83 KDa and 8.94, respectively. The GhCIPK23 protein contained a conserved N-terminal protein kinase domain and C-terminal regulatory domain of the CIPK gene family member. Phylogenetic tree analysis demonstrated that GhCIPK23 had a close relationship with AtCIPK23, followed by OsCIPK23, and belonged to Group A with AtCIPK23 and OsCIPK23. The subcellular localization experiment indicated that GhCIPK23 was located in the plasma membrane. Tissue expression analysis showed that GhCIPK23 had the highest expression in petals, followed by sepals, and the lowest in fibers. Stress expression analysis showed that the expression of the GhCIPK23 gene was in response to drought, salt, low-temperature and exogenous abscisic acid (ABA) treatment, and had different expression patterns under different stress conditions. Further cis-acting elements analysis showed that the GhCIPK23 promoter region had cis-acting elements in response to abiotic stress, phytohormones and light. These results established a foundation for understanding the function of GhCIPK23 and breeding varieties with high-stress tolerance in cotton.

Molecular and expression analysis indicate the role of CBL interacting protein kinases (CIPKs) in abiotic stress signaling and development in chickpea

Calcineurin B-like proteins (CBL)-interacting protein kinases (CIPKs) regulate the developmental processes, hormone signal transduction and stress responses in plants. Although the genome sequence of chickpea is available, information related to the CIPK gene family is missing in this important crop plant. Here, a total of 22 CIPK genes were identified and characterized in chickpea. We found a high degree of structural and evolutionary conservation in the chickpea CIPK family. Our analysis showed that chickpea CIPKs have evolved with dicots such as Arabidopsis and soybean, and extensive gene duplication events have played an important role in the evolution and expansion of the CIPK gene family in chickpea. The three-dimensional structure of chickpea CIPKs was described by protein homology modelling. Most CIPK proteins are localized in the cytoplasm and nucleus, as predicted by subcellular localization analysis. Promoter analysis revealed various cis-regulatory elements related to plant development, hormone signaling, and abiotic stresses. RNA-seq expression analysis indicated that CIPKs are significantly expressed through a spectrum of developmental stages, tissue/organs that hinted at their important role in plant development. The qRT-PCR analysis revealed that several CaCIPK genes had specific and overlapping expressions in different abiotic stresses like drought, salt, and ABA, suggesting the important role of this gene family in abiotic stress signaling in chickpea. Thus, this study provides an avenue for detailed functional characterization of the CIPK gene family in chickpea and other legume crops.

Genome-wide identification of calcineurin B-like protein-interacting protein kinase gene family reveals members participating in abiotic stress in the ornamental woody plant Lagerstroemia indica

Calcineurin B-like protein-interacting protein kinases (CIPKs) play important roles in plant responses to stress. However, their function in the ornamental woody plant Lagerstroemia indica is remains unclear. In this study, the LiCIPK gene family was analyzed at the whole genome level. A total of 37 LiCIPKs, distributed across 17 chromosomes, were identified. Conserved motif analysis indicated that all LiCIPKs possess a protein kinase motif (S_TKc) and C-terminal regulatory motif (NAF), while seven LiCIPKs lack a protein phosphatase interaction (PPI) motif. 3D structure analysis further revealed that the N-terminal and C-terminal 3D-structure of 27 members are situated near to each other, while 4 members have a looser structure, and 6 members lack intact structures. The intra- and interspecies collinearity analysis, synonymous substitution rate (K _s ) peaks of duplicated LiCIPKs, revealed that ∼80% of LiCIPKs were retained by the two whole genome duplication (WGD) events that occurred approximately 56.12-61.16 million year ago (MYA) and 16.24-26.34 MYA ago. The promoter of each LiCIPK contains a number of auxin, abscisic acid, gibberellic acid, salicylic acid, and drought, anaerobic, defense, stress, and wound responsive cis-elements. Of the 21 members that were successfully amplified by qPCR, 18 LiCIPKs exhibited different expression patterns under NaCl, mannitol, PEG8000, and ABA treatments. Given that LiCIPK30, the AtSOS2 ortholog, responded to all four types of stress it was selected for functional verification. LiCIPK30 complements the atsos2 phenotype in vivo. 35S:LiCIPK-overexpressing lines exhibit increased leaf area increment, chlorophyll a and b content, reactive oxygen species scavenging enzyme activity, and expression of ABF3 and RD22, while the degree of membrane lipid oxidation decreases under NaCl treatment compared to WT. The evolutionary history, and potential mechanism by which LiCIPK30 may regulate plant tolerance to salt stress were also discussed. In summary, we identified LiCIPK members involved in abiotic stress and found that LiCIPK30 transgenic Arabidopsis exhibits more salt and osmotic stress tolerance than WT. This research provides a theoretical foundation for further investigation into the function of LiCIPKs, and for mining gene resources to facilitate the cultivation and breeding of new L. indica varieties in coastal saline-alkali soil.

Identification of the CIPK-CBL family gene and functional characterization of CqCIPK14 gene under drought stress in quinoa

Background

Calcineurin-like Protein (CBL) and CBL interacting protein kinase (CIPK) play a key role in plant signal transduction and response to various environmental stimuli. Quinoa, as an important plant with high nutritional value, can meet the basic nutritional needs of human Cash crop, is also susceptible to abiotic stress. However, CBL-CIPK in quinoa have not been reported.

Results

In this study, 16 CBL and 41 CIPK genes were identified in quinoa. CBL-CIPK gene shows different intron-exon gene structure and motif, they participate in different biological processes, and form a complex regulatory network between CBL-CIPK proteins. Many cis-regulatory element associated with ABA and drought have been found. The expression patterns of CBL-CIPK showed different expression patterns in various abiotic stresses and tissues. RT-qPCR showed that most members of these two gene families were involved in drought regulation of quinoa, in particular, the expression levels of CqCIPK11, CqCIPK15, CqCIPK37 and CqCBL13 increased significantly under drought stress.

Conclusions

The structures and functions of the CBL-CIPK family in quinoa were systematically explored. Many CBL-CIPK may play vital roles in the regulation of organ development, growth, and responses to abiotic stresses. This research has great significance for the functional characterisation of the quinoa CBL-CIPK family and our understanding of the CBL-CIPK family in higher plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08683-6.

Inositol Improves Cold Tolerance Through Inhibiting CBL1 and Increasing Ca2+ Influx in Rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is an important oilseed crop worldwide. However, its productivity is significantly affected by various abiotic stresses, including cold stress. Among various stresses, cold stress is an important abiotic factor affecting plant growth, yield, and quality. The calcium channels are regarded as key pathways affecting cold tolerance in plants. Thus, improvement in cold tolerance is of great significance for crop improvement. The current study was designed to examine the beneficial role of exogenous inositol in improving cold stress tolerance in rapeseed. From the RNA-seq results, we identified 35 differently expressed genes encoding different inositol enzymes. The results show that inositol (a cyclic polyol) positively regulated cold tolerance by increasing calcium ion (Ca²⁺) influx in rapeseed. Furthermore, we found that the expression of calcineurin B-like (CBL1) gene was inhibited by inositol. On the other hand, overexpressed plant mediated the Ca²⁺ flux under cold stress suggesting the key role of inositol-Ca²⁺ pathway in cold tolerance. Moreover, the overexpression of BnCBL1-2 in Arabidopsis represented that transgenic plants mediated the Ca²⁺ flux highlighting the vital role of the inositol-Ca²⁺ pathway in conferring cold stress. Our study provides new insights into rapeseed cold tolerance mechanism and introduces a feasible method to improve the cold tolerance of rapeseed quickly.

Isolation and characterization of Salt Overly Sensitive family genes in spinach

The Salt Overly Sensitive (SOS) pathway regulates intracellular sodium ion homeostasis as a salt-stress response in plants. This pathway involves three main genes designated as SOS1, SOS2 and SOS3, which are members of the Na⁺ /H⁺ exchanger (NHX), CBL-interacting protein kinase (CIPK) and Calcineurin B-like (CBL) gene families, respectively. To identify and characterize SOS genes in spinach (Spinacia oleracea), a species of the Amaranthaceae family, we conducted genome-wide identification and phylogenetic analyses of NHX, CIPK and CBL genes from four Amaranthaceae species, Arabidopsis and rice. Most Amaranthaceae genes exhibited orthologous relationships with Arabidopsis and/or rice, except a clade of Vac-type Amaranthaceae NHX genes. Phylogenetic analyses also revealed gene gain/loss events in Amaranthaceae species and the intron-less to intron-rich evolution of CIPK genes. A bacterial protein-rooted CIPK tree allowed naming most of the phylogenetic clades based on their evolutionary history. Single S. oleracea (So) SOS1, SOS2 and SOS3 proteins were identified. Direct protein-protein interaction was observed between SoSOS2 and SoSOS3 but not between SoSOS2 and SoSOS1 based on yeast two-hybrid assay. This may suggest distinct modes of action of spinach SOS proteins compared to Arabidopsis SOS proteins. Unlike SoSOS1 and SoSOS2, which were expressed at similar or higher levels in leaves than roots, SoSOS3 expression was significantly higher in roots than leaves, suggesting its greater importance in roots. The expression of SoSOS3 was upregulated in both roots and leaves under salinity compared to the control; however, SoSOS1 was only upregulated in roots. Thus, this study demonstrated the conservation of SOS pathway genes in spinach and also highlighted the complexity of SOS signaling in Amaranthaceae species.

Physiological and molecular responses for long term salinity stress in common fig (Ficus carica L.)

Salinity stress in increasingly becoming a major challenge in current and expanding agricultural ecosystems. Unlike temporal abiotic stresses, plants are usually exposed to salinity stress for an entire lifespan. Therefore, a long term effect (10 weeks) of continuous salinity exposure was investigated for three common fig landraces (Zraki, Mwazi, and Khdari). Both relative water content and chlorophyll content decreased with elevated salinity stress, while stem length barely changed. The most prominent decline was observed in root biomass. The data would align common fig to moderately tolerant threshold slop with a C₅₀ range of 100 to 150 mM NaCl. A high and significant correlation was evident between root biomass and chlorophyll content (85%). Concurrently, differential expression of putative salinity responsive genes in common fig were determined; signal peptide peptidase-like 2B (FcSPPL2B), dehydration responsive element binding protein (FcDREB), calcineurin B-like protein (CBL)-CBL-interacting serine/threonine-protein kinase 11 (FcCIPK11), sorbitol dehydrogenase (FcSORD) and dehydrin (FcDHN). The data were discussed for each gene in respect of its potential role in salinity stress mitigation. The combined physiological and molecular data would conclude Zraki as the most salinity tolerant genotype. The major implication of the data emphasizes the tremendous genotype by environment (salinity stress) interaction in common fig.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at (10.1007/s12298-020-00921-z).

Genome-wide identification and expression analysis of the ERF transcription factor family in pineapple (Ananas comosus (L.) Merr.)

Pineapple (Ananas comosus (L.) Merr.) is an important tropical fruit with high economic value. The quality and yield of pineapple will be affected by various environmental conditions. Under adverse conditions, plants can produce a complex reaction mechanism to enhance their resistance. It has been reported that the member of ethylene responsive transcription factors (ERFs) plays a crucial role in plant developmental process and stress response. However, the function of these proteins in pineapple remains limited. In this study, a total of 74 ERF genes (AcoERFs) were identified in pineapple genome, named from AcoERF1 to AcoERF74, and divided into 13 groups based on phylogenetic analysis. We also analyzed gene structure, conserved motif and chromosomal location of AcoERFs, and the AcoERFs within the same group possess similar gene structures and motif compositions. Three genes (AcoERF71, AcoERF73 and AcoERF74) were present on unanchored scaffolds, so they could not be conclusively mapped on chromosome. Synteny and cis-elements analysis of ERF genes provided deep insight into the evolution and function of pineapple ERF genes. Furthermore, we analyzed the expression profiling of AcoERF in different tissues and developmental stages, and 22 AcoERF genes were expressed in all examined tissues, in which five genes (AcoERF13, AcoERF16, AcoERF31, AcoERF42, and AcoERF65) had high expression levels. Additionally, nine AcoERF genes were selected for functional verification by qRT-PCR. These results provide useful information for further investigating the evolution and functions of ERF family in pineapple.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.