Evolution of an intron-poor cluster of the CIPK gene family and expression in response to drought stress in soybean

Characterization of CBL-CIPK signaling networks and their response to abiotic stress in sugarcane

Calcineurin B-like proteins (CBLs) perceive calcium signals triggered by abiotic stress and interact with CBL-interacting protein kinases (CIPKs) to form a complex signal network. This study identified 21 SsCBL and 89 SsCIPK genes in Saccharum spontaneum, and 90 ScCBL and 367 ScCIPK genes in the sugarcane cultivar ZZ1. Phylogenetic analysis classified CBL genes into three groups and CIPK genes into twenty-five groups, with whole-genome duplication events promoting their expansion in sugarcane. RNA-seq analysis revealed their involvement in abiotic stress responses through ABA, JA, and SA pathways. Four ScCBLs and eight ScCIPKs were cloned from ZZ1. Three CBL-CIPK interactions were detected using a yeast two-hybrid system and Firefly luciferase complementation imaging, showing CBLs as membrane proteins and CIPKs as nuclear proteins. Spatial expression profiles indicate these genes are expressed in various tissues, with the highest expression in roots. Gene expression analyses suggested that CBL-CIPK signaling networks are involved in responses to drought, salt, and reactive oxygen species, possibly through Ca2+-induced hormone pathways. These findings establish three CBL-CIPK signaling networks responding to abiotic stress, providing a molecular basis for improving sugarcane stress resistance.

Identification and characterization of NHX gene family for their role under salt stress in Vigna mungo

In the current study, we have performed a comprehensive analysis of the Sodium Hydrogen Exchanger (NHX) gene family in Vigna mungo, and a total of 44 NHX genes were identified. A bimodal distribution based on domains, gene structure and phylogenetic analysis was evident. All intronpoor and intron-rich genes were clustered in clades I and II, respectively. Interestingly, all genes of subclade IIb were localized to vacuoles and possess only the NHX domain. The isoelectric point and trans-membrane domain analysis reflect the wide distribution of the NHX genes. Interestingly, Vm_NHX2 and Vm_NHX3 lacked trans-membrane domain but were found to interact with other NHX genes as well as vital salinity pathway genes, including calcium-mediated salt-responsive genes. The comparison of the mRNA sequences with that of V. marina, a halophytic species, reflects their independent evolution, majorly supporting the convergent evolution. The Ka/Ks ratio reflects the abundance of purifying selection supporting their conserved function during evolution. In our analysis, several abiotic stress and hormone-responsive elements and transcription factor binding sites were present in the promoter of the NHX genes. Further, the ion partitioning of a tolerant (K90) and a susceptible (K49) variety of V. mungo suggested that K90 managed the Na+/K+ ratio more affluently, which was also supported by profiling of superoxide radicals, hydrogen peroxide, phenol, peroxidase activity and superoxide dismutase activity. From the expression, we identified five candidate Vm_NHX genes, four of which, i.e. Vm_NHX16, Vm_NHX17, Vm_NHX29 and Vm_NHX33, were localized to the vacuolar and lysosomal membrane.

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

Genome-Wide Identification of Members of the Soybean CBL Gene Family and Characterization of the Functional Role of GmCBL1 in Responses to Saline and Alkaline Stress

Calcium ions function as key messengers in the context of intracellular signal transduction. The ability of plants to respond to biotic and abiotic stressors is highly dependent on the calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) signaling network. Here, a comprehensive effort was made to identify all members of the soybean CBL gene family, leading to the identification of 15 total genes distributed randomly across nine chromosomes, including 13 segmental duplicates. All the GmCBL gene subfamilies presented with similar gene structures and conserved motifs. Analyses of the expression of these genes in different tissues revealed that the majority of these GmCBLs were predominantly expressed in the roots. Significant GmCBL expression and activity increases were also observed in response to a range of stress-related treatments, including salt stress, alkaline stress, osmotic stress, or exposure to salicylic acid, brassinosteroids, or abscisic acid. Striking increases in GmCBL1 expression were observed in response to alkaline and salt stress. Subsequent analyses revealed that GmCBL1 was capable of enhancing soybean salt and alkali tolerance through the regulation of redox reactions. These results offer new insight into the complex mechanisms through which the soybean CBL gene family regulates the responses of these plants to environmental stressors, highlighting promising targets for efforts aimed at enhancing soybean stress tolerance.

Genome-Wide Identification and Expression Analysis of Sucrose Nonfermenting 1-Related Protein Kinase (SnRK) Genes in Salvia miltiorrhiza in Response to Hormone

The SnRK gene family is the chief component of plant stress resistance and metabolism through activating the phosphorylation of downstream proteins. S. miltiorrhiza is widely used for the treatment of cardiovascular diseases in Asian countries. However, information about the SnRK gene family of S. miltiorrhiza is not clear. The aim of this study is to comprehensively analyze the SnRK gene family of S. miltiorrhiza and its response to phytohormone. Here, 33 SmSnRK genes were identified and divided into three subfamilies (SmSnRK1, SmSnRK2 and SmSnRK3) according to phylogenetic analysis and domain. SmSnRK genes within same subgroup shared similar protein motif composition and were unevenly distributed on eight chromosomes of S. miltiorrhiza. Cis-acting element analysis showed that the promoter of SmSnRK genes was enriched with ABRE motifs. Expression pattern analysis revealed that SmSnRK genes were preferentially expressed in leaves and roots. Most SmSnRK genes were induced by ABA and MeJA treatment. Correlation analysis showed that SmSnRK3.15 and SmSnRK3.18 might positively regulate tanshinone biosynthesis; SmSnRK3.10 and SmSnRK3.12 might positively regulate salvianolic acid biosynthesis. RNAi-based silencing of SmSnRK2.6 down-regulated the biosynthesis of tanshinones and biosynthetic genes expression. An in vitro phosphorylation assay verified that SmSnRK2.2 interacted with and phosphorylated SmAREB1. These findings will provide a valuable basis for the functional characterization of SmSnRK genes and quality improvement of S. miltiorrhiza.

Identification of loci associated with water use efficiency and symbiotic nitrogen fixation in soybean

Soybean (Glycine max) production is greatly affected by persistent and/or intermittent droughts in rainfed soybean-growing regions worldwide. Symbiotic N2 fixation (SNF) in soybean can also be significantly hampered even under moderate drought stress. The objective of this study was to identify genomic regions associated with shoot carbon isotope ratio (δ13C) as a surrogate measure for water use efficiency (WUE), nitrogen isotope ratio (δ15N) to assess relative SNF, N concentration ([N]), and carbon/nitrogen ratio (C/N). Genome-wide association mapping was performed with 105 genotypes and approximately 4 million single-nucleotide polymorphism markers derived from whole-genome resequencing information. A total of 11, 21, 22, and 22 genomic loci associated with δ13C, δ15N, [N], and C/N, respectively, were identified in two environments. Nine of these 76 loci were stable across environments, as they were detected in both environments. In addition to the 62 novel loci identified, 14 loci aligned with previously reported quantitative trait loci for different C and N traits related to drought, WUE, and N2 fixation in soybean. A total of 58 Glyma gene models encoding for different genes related to the four traits were identified in the vicinity of the genomic loci.

Characterization of Dendrobium catenatum CBL-CIPK signaling networks and their response to abiotic stress

Dendrobium catenatum is a traditional Chinese medicine listing as rare and endangered due to environmental impacts. But little is known about its stress resistance mechanism. The CBL-CIPK signaling pathway played vital roles in various stress responses. In this study, we identified 9 calcineurin B-like (CBL) genes and 28 CBL-interacting protein kinase (CIPK) genes from D. catenatum. Phylogenetic analysis showed that DcCBL and DcCIPK families could be divided into four and six subgroups, respectively. Members in each subgroup had similar gene structures. Cis-acting element analyses showed that these genes were involved in stress responses and hormone signaling. Spatial expression profiles showed that they were tissue-specific, and expressed lower in vegetative organs than reproductive organs. Gene expression analyses revealed that these genes were involved in drought, heat, cold, and salt responses and depended on abscisic acid (ABA) and salicylic acid (SA) signaling pathways. Furthermore, we cloned 19 DcCIPK genes and 9 DcCBL genes and detected ten interacting CBL-CIPK combinations using yeast two-hybrid system. Finally, we constructed 20 CBL-CIPK signaling pathways based on their expression patterns and interaction relationships. These results established CBL-CIPK signaling pathway responding to abiotic stress and provided a molecular basis for improving D. catenatum stress resistance in the future.

CIPK-B is essential for salt stress signalling in Marchantia polymorpha

Calcium signalling is central to many plant processes, with families of calcium decoder proteins having expanded across the green lineage and redundancy existing between decoders. The liverwort Marchantia polymorpha has fast become a new model plant, but the calcium decoders that exist in this species remain unclear. We performed phylogenetic analyses to identify the calcineurin B-like (CBL) and CBL-interacting protein kinase (CIPK) network of M. polymorpha. We analysed CBL-CIPK expression during salt stress, and determined protein-protein interactions using yeast two-hybrid and bimolecular fluorescence complementation. We also created genetic knockouts using CRISPR/Cas9. We confirm that M. polymorpha has two CIPKs and three CBLs. Both CIPKs and one CBL show pronounced salt-responsive transcriptional changes. All M. polymorpha CBL-CIPKs interact with each other in planta. Knocking out CIPK-B causes increased sensitivity to salt, suggesting that this CIPK is involved in salt signalling. We have identified CBL-CIPKs that form part of a salt tolerance pathway in M. polymorpha. Phylogeny and interaction studies imply that these CBL-CIPKs form an evolutionarily conserved salt overly sensitive pathway. Hence, salt responses may be some of the early functions of CBL-CIPK networks and increased abiotic stress tolerance required for land plant emergence.

Identification and Expression Analysis of MPK and MKK Gene Families in Pecan (Carya illinoinensis)

Mitogen-activated protein kinases consist of three kinase modules composed of MPKs, MKKs, and MPKKKs. As members of the protein kinase (PK) superfamily, they are involved in various processes, such as developmental programs, cell division, hormonal progression, and signaling responses to biotic and abiotic stresses. In this study, a total of 18 MPKs and 10 MKKs were annotated on the pecan genome, all of which could be classified into four subgroups, respectively. The gene structures and conserved sequences of family members in the same branch were relatively similar. All MPK proteins had a conserved motif TxY, and D(L/I/V)K and VGTxxYMSPER existed in all MKK proteins. Duplication events contributed largely to the expansion of the pecan MPK and MKK gene families. Phylogenetic analysis of protein sequences from six plants indicated that species evolution occurred in pecan. Organ-specific expression profiles of MPK and MKK showed functional diversity. Ka/Ks values indicated that all genes with duplicated events underwent strong negative selection. Seven CiPawMPK and four CiPawMKK genes with high expression levels were screened by transcriptomic data from different organs, and these candidates were validated by qRT-PCR analysis of hormone-treated and stressed samples.

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.

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.

Genome-wide characterization of C2H2 zinc-finger gene family provides insight into the mechanisms and evolution of the dehydration-rehydration responses in Physcomitrium and Arabidopsis

Dehydration tolerance is a vital factor for land plant evolution and world agricultural production. Numerous studies enlightened that the plant-specific C2H2-type zinc-finger proteins (C2H2-ZFPs) as master regulators played pivotal roles in the abiotic stress responses of plants. However, a comprehensive understanding of the evolution of C2H2-ZFPs in terrestrial plants and its regulatory mechanism in dehydration and rehydration response remains a mystery. In this study, the genome-wide identification of C2H2-ZFP genes revealed 549 homologs in the representatives of terrestrial plant lineages from liverwort to angiosperms. Based on the characteristics of the conserved C2H2-ZF domains, four major C2H2-ZF types (M-, Z-, Q-, and D-type) were identified in the C2H2-ZFPs, with the dominants of M-type in all selected species and followed by Z-type in non-seed plants and Q-type in seed plants, respectively. Phylogenetic analyses of the identified C2H2-ZFPs supported four major groups in the land plant representatives, among which the members from the desiccation-tolerant Physcomitrium patens and the dehydration-sensitive Arabidopsis thaliana displayed different topological relationships in the phylogenies reconstructed for a single species. C2H2-ZFPs clustered in the same subclades shared similar features in their conserved domains and gene structures. Approximately, 81% of the C2H2-ZFP promoters of all 549 identified C2H2-ZFPs harbored the conserved ABA-responsive elements (ABREs) and/or dehydration-responsive elements (DREs). Comparative transcriptomic analyses showed that 50 PpZFPs and 56 AtZFPs significantly changed their transcripts abundance. Interestingly, most of the dehydration- and rehydration-responsive PpZPFs and AtZFPs had been predicted to contain the ABRE and DRE elements in their promoter regions and with over half of which phylogenetically belonging to group III. The differences in the expression patterns of C2H2-ZFPs in responses to dehydration and rehydration between P. patens and A. thaliana reflected their different strategies to adapt to dehydration. The identified candidate PpZFPs were specifically induced by moderate dehydration and reached the peak transcript abundance in severe dehydration. Our study lays the foundations for further functional investigation of C2H2-ZFPs in dehydration responses from an evolutionary perspective in land plants. The findings will provide us with genetic resources and potential targets for drought tolerance breeding in crops and beyond.

GmCIPK21, a CBL-interacting protein kinase confers salt tolerance in soybean (Glycine max. L)

Salt stress severely affects plant development and yield. Calcineurin B-like protein interacting protein kinases (CIPKs) play a crucial role in plant adaptation to environmental challenges. However, the biological functions of CIPKs in soybean remain poorly understood. Here, we identified GmCIPK21, a salt-responsive CIPK gene from soybean. Overexpression of GmCIPK21 in Arabidopsis and soybean hairy roots led to increased salt tolerance. The hairy roots with GmCIPK21 suppression by RNA interference exhibited salt-sensitive phenotypes. Further physiological analysis revealed that GmCIPK21 reduced the content of hydrogen peroxide (H2O2) and malondialdehyde (MDA) and increased the activity of the antioxidant enzymes under salt stress. Additionally, GmCIPK21 was found to enhance the ABA sensitivity of transgenic plants. GmCIPK21 was also implicated in increasing the activation of antioxidant-, salt-, and ABA-related genes upon salt stress. Interestingly, GmCIPK21 interacted with GmCBL4, promoting the scavenging salt-induced reactive oxygen species (ROS). These results collectively suggested that GmCIPK21 affects ROS homeostasis and ABA response to improve salt tolerance in soybean.

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.

Genome-Wide Identification and Functional Analysis of the Calcineurin B-like Protein and Calcineurin B-like Protein-Interacting Protein Kinase Gene Families in Chinese Cabbage (Brassica rapa ssp. pekinensis)

In plants, calcineurin B-like proteins (CBL) are a unique set of calcium sensors that decode calcium signals by activating a plant-specific protein kinase family called CBL-interacting protein kinases (CIPKs). The CBL–CIPK family and its interacting complexes regulate plant responses to various environmental stimuli. Chinese cabbage (Brassica rapa ssp. pekinensis) is an important vegetable crop in Asia; however, there are no reports on the role of the CBLs–CIPKs’ signaling system in response to abiotic stress during cabbage growth. In this study, 18 CBL genes and 47 CIPK genes were identified from the Chinese cabbage genome. Expansion of the gene families was mainly due to tandem repeats and segmental duplication. An analysis of gene expression patterns showed that different duplicate genes exhibited different expression patterns in response to treatment with Mg²⁺, K⁺, and low temperature. In addition, differences in the structural domain sequences of NAF/FISL and interaction profiles in yeast two-hybrid assays suggested a functional divergence of the duplicate genes during the long-term evolution of Chinese cabbage, a result further validated by potassium deficiency treatment using trans-BraCIPK23.1/23.2/23.3 Arabidopsis thaliana. Our results provide a basis for studies related to the functional divergence of duplicate genes and in-depth studies of BraCBL–BraCIPK functions in Chinese cabbage.

Insight into the CBL and CIPK gene families in pecan (Carya illinoinensis): identification, evolution and expression patterns in drought response

BACKGROUND

Calcium (Ca²⁺) serves as a ubiquitous second messenger and plays a pivotal role in signal transduction. Calcineurin B-like proteins (CBLs) are plant-specific Ca²⁺ sensors that interact with CBL-interacting protein kinases (CIPKs) to transmit Ca²⁺ signals. CBL-CIPK complexes have been reported to play pivotal roles in plant development and response to drought stress; however, limited information is available about the CBL and CIPK genes in pecan, an important nut crop.

RESULTS

In the present study, a total of 9 CBL and 30 CIPK genes were identified from the pecan genome and divided into four and five clades based on phylogeny, respectively. Gene structure and distribution of conserved sequence motif analysis suggested that family members in the same clade commonly exhibited similar exon-intron structures and motif compositions. The segmental duplication events contributed largely to the expansion of pecan CBL and CIPK gene families, and Ka/Ks values revealed that all of them experienced strong negative selection. Phylogenetic analysis of CIPK proteins from 14 plant species revealed that CIPKs in the intron-poor clade originated in seed plants. Tissue-specific expression profiles of CiCBLs and CiCIPKs were analysed, presenting functional diversity. Expression profiles derived from RNA-Seq revealed distinct expression patterns of CiCBLs and CiCIPKs under drought treatment in pecan. Moreover, coexpression network analysis helped to elucidate the relationships between these genes and identify potential candidates for the regulation of drought response, which were verified by qRT-PCR analysis.

CONCLUSIONS

The characterization and analysis of CBL and CIPK genes in pecan genome could provide a basis for further functional analysis of CiCBLs and CiCIPKs in the drought stress response of pecan.

Genome-wide identification, molecular characterization, and gene expression analyses of honeysuckle NHX antiporters suggest their involvement in salt stress adaptation

Background

Ion homeostasis is an essential process for the survival of plants under salt stress. Na⁺/H⁺ antiporters (NHXs) are secondary ion transporters that regulate Na⁺ compartmentalization or efflux reduce Na⁺ toxicity and play a critical role during plant development and stress responses.

Methods and Results

To gain insight into the functional divergence of NHX genes in honeysuckle, a total of seven LjNHX genes were identified on the whole genome level and were renamed according to their chromosomal positions. All LjNHXs possessed the Na⁺/H⁺ exchanger domain and the amiloride-binding site was presented in all NHX proteins except LjNHX4. The phylogenetic analysis divided the seven NHX genes into Vac-clade (LjNHX1/2/3/4/5/7) and PM-clade (LjNHX6) based on their subcellular localization and validated by the distribution of conserved protein motifs and exon/intron organization analysis. The protein-protein interaction network showed that LjNHX4/5/6/7 shared the same putatively interactive proteins, including SOS2, SOS3, HKT1, and AVP1. Cis-acting elements and gene ontology (GO) analysis suggested that most LjNHXs involve in the response to salt stress through ion transmembrane transport. The expression profile analysis revealed that the expression levels of LjNHX3/7 were remarkably affected by salinity. These results suggested that LjNHXs play significant roles in honeysuckle development and response to salt stresses.

Conclusions

The theoretical foundation was established in the present study for the further functional characterization of the NHX gene family in honeysuckle.

The Calcium Sensor Calcineurin B-Like Proteins -Calcineurin B-Like Interacting Protein Kinases Is Involved in Leaf Development and Stress Responses Related to Latex Flow in Hevea brasiliensis

Latex flow in Hevea brasiliensis (the Para rubber tree), the sole commercial source of natural rubber (cis-1,4-polyisoprene, NR), renders it uniquely suited for the study of plant stress responses. Calcineurin B-like interacting protein kinases (CIPK) serving as calcium-sensor protein kinases react with calcineurin B-like proteins (CBL) to play crucial roles in hormone signaling transduction and response to abiotic stress in plant developmental processes. However, little is known about their functions in Hevea. In this study, a total of twelve CBL (HbCBL) and thirty CIPK (HbCIPK) genes were identified from the Hevea genome. Structure and phylogenetic analysis assigned these CIPKs to five groups and CBLs to four groups, and mapped onto fourteen of the eighteen Hevea chromosomes. RNA-seq and qPCR analysis showed that the expressions of HbCBL and HbCIPK genes varied in the seven Hevea tissues examined, i.e., latex (cytoplasm of rubber-producing laticifers), bark, leaf, root, seed, female flower, and male flower. The expressions of two HbCBL and sixteen HbCIPK genes showed upward trends during leaf development. Following ethylene yield stimulation and the latex tapping treatment, both practices invoking stress, the expression levels of most latex-expressed genes were significantly altered. Yeast two-hybrid test revealed interactions for multiple combinations of HbCBLs and HbCIPKs with substantial gene expression in latex or other Hevea tissues. However, all the HbCBL-HbCIPK complexes examined did not recruit HbSOS1 or AtSOS1 to form functional salt tolerance SOS pathway in yeast cells. Taken together, the results suggested a role of the Hevea CBL-CIPK network as a point of convergence for several different signaling pathways in growth, development, and stress responses in relation to latex production.

Overexpression of CpWRKY75 from Chimonanthus praecox Promotes Flowering Time in Transgenic Arabidopsis

WRKY transcription factors play critical roles in the physiological processes of plants. Although the roles of WRKYs have been characterized in some model plants, their roles in woody plants, especially wintersweet (Chimonanthus praecox), are largely unclear. In this study, a wintersweet WRKY gene named CpWRKY75 belonging to group IIc was isolated and its characteristics were identified. CpWRKY75 is a nucleus-localized protein, and exhibited no transcriptional activation activity in yeast. CpWRKY75 was highly expressed in flowers at different bloom stages. Ectopic expression of CpWRKY75 significantly promoted the flowering time of transgenic Arabidopsis (Arabidopsis thaliana), as determined by the rosette leaf number and first flower open time. The expression levels of flowering-related genes were quantified by qRT-PCR, and the results suggested that CpWRKY75 had obvious influence on the expression level of MICRORNA156C (MIR156C), SQUAMOSA PROMOTER BINDING PROTEIN-LIKE3 (SPL3) and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9 (SPL9), FLOWERING LOCUS T (FT), LEAFY (LFY), SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), APETALA1 (AP1), CAULIFLOWER (CAL), and FRUITFULL (FUL). These results suggest that CpWRKY75 might have a flowering time regulation function, and additionally provide a new gene resource for the genetic engineering of woody flowering plants.

Genome-Wide Identification, Characterization and Expression Analysis of the CIPK Gene Family in Potato (Solanum tuberosum L.) and the Role of StCIPK10 in Response to Drought and Osmotic Stress

The potato (Solanum tuberosum L.), one of the most important food crops worldwide, is sensitive to environmental stresses. Sensor–responder complexes comprising calcineurin B-like (CBL) proteins and CBL-interacting protein kinases (CIPKs) not only modulate plant growth and development but also mediate numerous stress responses. Here, using a Hidden Markov Model and BLAST searches, 27 CIPK genes were identified in potato and divided into five groups by phylogenetic analysis and into two clades (intron-poor and intron-rich) by gene structure analysis. Quantitative reverse-transcription PCR (qRT-PCR) assays revealed that StCIPK genes play important roles in plant growth, development and abiotic stress tolerance. Up-regulated expression of StCIPK10 was significantly induced by drought, PEG6000 and ABA. StCIPK10 enhances both the ability of potato to scavenge reactive oxygen species and the content of corresponding osmoregulation substances, thereby strengthening tolerance to drought and osmotic stress. StCIPK10 is located at the intersection between the abscisic acid and abiotic stress signaling pathways, which control both root growth and stomatal closure in potato. In addition, StCIPK10 interacts with StCBL1, StCBL4, StCBL6, StCBL7, StCBL8, StCBL11 and StCBL12, and is specifically recruited to the plasma membrane by StCBL11.

Genome-Wide Identification of CBL-CIPK Gene Family in Honeysuckle (Lonicera japonica Thunb.) and Their Regulated Expression Under Salt Stress

In plants, calcineurin B-like proteins (CBLs) are a unique group of Ca²⁺ sensors that decode Ca²⁺ signals by activating a family of plant-specific protein kinases known as CBL-interacting protein kinases (CIPKs). CBL-CIPK gene families and their interacting complexes are involved in regulating plant responses to various environmental stimuli. To gain insight into the functional divergence of CBL-CIPK genes in honeysuckle, a total of six LjCBL and 17 LjCIPK genes were identified. The phylogenetic analysis along with the gene structure analysis divided both CBL and CBL-interacting protein kinase genes into four subgroups and validated by the distribution of conserved protein motifs. The 3-D structure prediction of proteins shown that most LjCBLs shared the same Protein Data Bank hit 1uhnA and most LjCIPKs shared the 6c9Da. Analysis of cis-acting elements and gene ontology implied that both LjCBL and LjCIPK genes could be involved in hormone signal responsiveness and stress adaptation. Protein-protein interaction prediction suggested that LjCBL4 is hypothesized to interact with LjCIPK7/9/15/16 and SOS1/NHX1. Gene expression analysis in response to salinity stress revealed that LjCBL2/4, LjCIPK1/15/17 under all treatments gradually increased over time until peak expression at 72 h. These results demonstrated the conservation of salt overly sensitive pathway genes in honeysuckle and a model of Ca²⁺-LjCBL4/LjSOS3-LjCIPK16/LjSOS2 module-mediated salt stress signaling in honeysuckle is proposed. This study provides insight into the characteristics of the CBL-CIPK gene families involved in honeysuckle salt stress responses, which could serve as a foundation for gene transformation technology, to obtain highly salt-tolerant medicinal plants in the context of the global reduction of cultivated land.

The soybean CBL-interacting protein kinase, GmCIPK2, positively regulates drought tolerance and ABA signaling

Calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) play important roles in plant environmental stress responses. However, the biological functions of the CBL-CIPK signaling pathway in the tolerance of soybean (Glycine max) to drought stress remain elusive. Here, we characterized the GmCIPK2 gene in soybean, and its expression was induced by drought stress and exogenous abscisic acid (ABA) treatments. The overexpression of GmCIPK2 enhanced drought tolerance in transgenic Arabidopsis and soybean hairy roots, whereas downregulation of GmCIPK2 expression in soybean hairy roots by RNA interference resulted in increased drought sensitivity. Further analysis showed that GmCIPK2 was involved in ABA-mediated stomatal closure in plants under drought stress conditions. GmCIPK2 increased the expression of ABA- and drought-responsive genes during drought stress. Additionally, yeast two-hybrid, pull-down, and bimolecular fluorescence complementation assays demonstrated that a positive regulator of drought stress, GmCBL1, physically interacted with GmCIPK2 on the plasma membrane. Collectively, our results demonstrated that GmCIPK2 positively regulates drought tolerance and ABA signaling in plants, providing new insights into the underlying mechanisms of how the CBL-CIPK signaling pathway contributes to drought tolerance in soybean.

Pecan kinome: classification and expression analysis of all protein kinases in Carya illinoinensis

Protein kinases (PKs) are involved in plant growth and stress responses, and constitute one of the largest superfamilies due to numerous gene duplications. However, limited PKs have been functionally described in pecan, an economically important nut tree. Here, the comprehensive identification, annotation and classification of the entire pecan kinome are reported. A total of 967 PK genes were identified from the pecan genome, and further classified into 20 different groups and 121 subfamilies using the kinase domain sequences, which were verified by phylogenetic analysis. The receptor-like kinase (RLK) group contained 565 members, which constituted the largest group. Gene duplication contributed to the expansion of pecan kinome, 169 segmental duplication events including 285 PK genes were found, and the Ka/Ks ratio revealed they experienced strong negative selection. The RNA-Seq data of PK genes in pecan were further analyzed at the subfamily level, and different PK subfamilies performed various expression patterns across pecan embryo development or drought treatment, suggesting PK genes in pecan are involved in embryo development and drought stress response. Taken together, this study provides insight into the classification, expansion, evolution, and expression of pecan PKs. Our findings regarding expansion, expression and co-expression analyses lay a good foundation for future research to understand the roles of pecan PKs, and more efficiently determine the key candidate genes.

Using Network-Based Machine Learning to Predict Transcription Factors Involved in Drought Resistance

Gene regulatory networks underpin stress response pathways in plants. However, parsing these networks to prioritize key genes underlying a particular trait is challenging. Here, we have built the Gene Regulation and Association Network (GRAiN) of rice (Oryza sativa). GRAiN is an interactive query-based web-platform that allows users to study functional relationships between transcription factors (TFs) and genetic modules underlying abiotic-stress responses. We built GRAiN by applying a combination of different network inference algorithms to publicly available gene expression data. We propose a supervised machine learning framework that complements GRAiN in prioritizing genes that regulate stress signal transduction and modulate gene expression under drought conditions. Our framework converts intricate network connectivity patterns of 2160 TFs into a single drought score. We observed that TFs with the highest drought scores define the functional, structural, and evolutionary characteristics of drought resistance in rice. Our approach accurately predicted the function of OsbHLH148 TF, which we validated using in vitro protein-DNA binding assays and mRNA sequencing loss-of-function mutants grown under control and drought stress conditions. Our network and the complementary machine learning strategy lends itself to predicting key regulatory genes underlying other agricultural traits and will assist in the genetic engineering of desirable rice varieties.

The pigeon pea CcCIPK14-CcCBL1 pair positively modulates drought tolerance by enhancing flavonoid biosynthesis

Calcineurin B-like (CBL)-interacting protein kinases (CIPKs) play a central role in Ca²⁺ signalling and promote drought tolerance in plants. The CIPK gene family in pigeon pea (Cajanus cajan L.), a major food crop affected by drought, has not previously been characterised. Here, we identified 28 CIPK genes in the pigeon pea genome. Five CcCIPK genes were strongly upregulated in roots upon drought treatment and were selected for further characterisation. Overexpression of CcCIPK13 and CcCIPK14 increased survival rates by two- to three-fold relative to controls after 14 days of drought. Furthermore, the three major flavonoids, genistin, genistein and apigenin, were significantly upregulated in the same transgenic plants. Using CcCIPK14 as bait, we performed a yeast two-hybrid screen and identified six interactors, including CcCBL1. CcCIPK14 exhibited autophosphorylation and phosphorylation of CcCBL1 in vitro. CcCBL1-overexpressed plants displayed higher survival rates upon drought stress as well as higher expression of flavonoid biosynthetic genes and flavonoid content. CcCIPK14-overexpressed plants in which CcCBL1 transcript levels were reduced by RNA interference had lower survival rates, which indicated CcCBL1 in the same pathway as CcCIPK14. Together, our results demonstrate a role for the CcCIPK14-CcCBL1 complex in drought stress tolerance through the regulation of flavonoid biosynthesis in pigeon pea.

Identification and Characterization of Abiotic Stress Responsive CBL-CIPK Family Genes in Medicago

The calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) play important roles in plant signal transduction and response to abiotic stress. Plants of Medicago genus contain many important forages, and their growth is often affected by a variety of abiotic stresses. However, studies on the CBL and CIPK family member and their function are rare in Medicago. In this study, a total of 23 CBL and 58 CIPK genes were identified from the genome of Medicago sativa as an important forage crop, and Medicaog truncatula as the model plant. Phylogenetic analysis suggested that these CBL and CIPK genes could be classified into five and seven groups, respectively. Moreover, these genes/proteins showed diverse exon-intron organizations, architectures of conserved protein motifs. Many stress-related cis-acting elements were found in their promoter region. In addition, transcriptional analyses showed that these CBL and CIPK genes exhibited distinct expression patterns in various tissues, and in response to drought, salt, and abscisic acid treatments. In particular, the expression levels of MtCIPK2 (MsCIPK3), MtCIPK17 (MsCIPK11), and MtCIPK18 (MsCIPK12) were significantly increased under PEG, NaCl, and ABA treatments. Collectively, our study suggested that CBL and CIPK genes play crucial roles in response to various abiotic stresses in Medicago.

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.

Genome-wide identification and expression pattern of SnRK gene family under several hormone treatments and its role in floral scent emission in Hedychium coronarium

The SnRK (Snf1-Related protein Kinase) gene family plays crucial roles in various plant signaling pathways and stress-adaptive responses including biotic and abiotic stresses via activating protein phosphorylation pathways. However, there is no information available on the role of the SnRK gene family in Hedychium coronarium. H. coronarium is an important crop widely cultivated as an ornamental plant, herb, spice, or condiment. In this study, 60 HcSnRK genes were identified from the H. coronarium genomic and transcriptome data. Phylogenetic and gene structure analysis showed that the HcSnRK genes were divided into three groups (HcSnRK1, HcSnRK2 and HcSnRK3) and among them HcSnRK3 subfamily was further subdivided into two clades according to the number of introns. Chromosome localization analysis showed that HcSnRK genes were unevenly mapped onto all chromosomes, and the Ka/Ks ratio of 24 paralogues includes four tandems and 20 segmental duplications indicated that the HcSnRK gene family underwent a purifying selection. Cis-regulatory elements analysis suggested that the HcSnRK genes respond to multiple hormones and other stresses. The responsiveness of HcSnRK genes to several hormones was analyzed by quantitative real-time PCR. Based on the different transcriptome data, two candidates HcSnRK genes (HcSnRK2.2 and HcSnRK2.9) were screened out for further characterization . The subcellular localization experiment revealed that both genes were located in the nucleus and cytoplasm. Moreover, virus-induced gene silencing (VIGS) of HcSnRK2.2 and HcSnRK2.9 significantly reduced the floral volatile contents by suppressing the expression of terpene synthase genes (HcTPS1, HcTPS3, and HcTPS5), indicating that HcSnRK2.2 and HcSnRK2.9 genes play an important role in the regulatory mechanism of floral aroma. These results will provide novel insights into the functional dissection of H. coronarium SnRK gene family.

Remodeling of the cell wall as a drought-tolerance mechanism of a soybean genotype revealed by global gene expression analysis

Drought stress is major abiotic stress that affects soybean production. Therefore, it is widely desirable that soybean becomes more tolerant to stress. To provide insights into regulatory mechanisms of the stress response, we compared the global gene expression profiles from leaves of two soybean genotypes that display different responses to water-deficit (BR 16 and Embrapa 48, drought-sensitive and drought-tolerant, respectively). After the RNA-seq analysis, a total of 5335 down-regulated and 3170 up-regulated genes were identified in the BR16. On the other hand, the number of genes differentially expressed was markedly lower in the Embrapa 48, 355 up-regulated and 471 down-regulated genes. However, induction and expression of protein kinases and transcription factors indicated signaling cascades involved in the drought tolerance. Overall, the results suggest that the metabolism of pectin is differently modulated in response to drought stress and may play a role in the soybean defense mechanism against drought. This occurs via an increase of the cell wall plasticity and crosslink, which contributed to a higher hydraulic conductance (K _f) and relative water content (RWC%). The drought-tolerance mechanism of the Embrapa 48 genotype involves remodeling of the cell wall and increase of the hydraulic conductance to the maintenance of cell turgor and metabolic processes, resulting in the highest leaf RWC, photosynthetic rate (A), transpiration (E) and carboxylation (A/C _i). Thus, we concluded that the cell wall adjustment under drought is important for a more efficient water use which promoted a more active photosynthetic metabolism, maintaining higher plant growth under drought stress.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00043-4.

The CBL-Interacting Protein Kinase NtCIPK23 Positively Regulates Seed Germination and Early Seedling Development in Tobacco (Nicotiana tabacum L.)

CBL-interacting protein kinase (CIPK) family is a unique group of serine/threonine protein kinase family identified in plants. Among this family, AtCIPK23 and its homologs in some plants are taken as a notable group for their importance in ions transport and stress responses. However, there are limited reports on their roles in seedling growth and development, especially in Solanaceae plants. In this study, NtCIPK23, a homolog of AtCIPK23 was cloned from Nicotiana tabacum. Expression analysis showed that NtCIPK23 is mainly expressed in the radicle, hypocotyl, and cotyledons of young tobacco seedlings. The transcriptional level of NtCIPK23 changes rapidly and spatiotemporally during seed germination and early seedling growth. To study the biological function of NtCIPK23 at these stages, the overexpressing and CRISPR/Cas9-mediated knock-out (ntcipk23) tobacco lines were generated. Phenotype analysis indicated that knock-out of NtCIPK23 significantly delays seed germination and the appearance of green cotyledon of young tobacco seedling. Overexpression of NtCIPK23 promotes cotyledon expansion and hypocotyl elongation of young tobacco seedlings. The expression of NtCIPK23 in hypocotyl is strongly upregulated by darkness and inhibited under light, suggesting that a regulatory mechanism of light might underlie. Consistently, a more obvious difference in hypocotyl length among different tobacco materials was observed in the dark, compared to that under the light, indicating that the upregulation of NtCIPK23 contributes greatly to the hypocotyl elongation. Taken together, NtCIPK23 not only enhances tobacco seed germination, but also accelerate early seedling growth by promoting cotyledon greening rate, cotyledon expansion and hypocotyl elongation of young tobacco seedlings.

The emergence and evolution of intron-poor and intronless genes in intron-rich plant gene families

Eukaryotic genes can be classified into intronless (no introns), intron-poor (three or fewer introns per gene) or intron-rich. Early eukaryotic genes were mostly intron-rich, and their alternative splicing into multiple transcripts, giving rise to different proteins, might have played pivotal roles in adaptation and evolution. Interestingly, extant plant genomes contain many gene families with one or sometimes few sub-families with genes that are intron-poor or intronless, and it remains unknown when and how these intron-poor or intronless genes have originated and evolved, and what their possible functions are. In this study, we identified 33 such gene families that contained intronless and intron-poor sub-families. Intronless genes seemed to have first emerged in early land plant evolution, while intron-poor sub-families seemed first to have appeared in green algae. In contrast to intron-rich genes, intronless genes in intron-poor sub-families occurred later, and were subject to stronger functional constraints. Based on RNA-seq analyses in Arabidopsis and rice, intronless or intron-poor genes in AP2, EF-hand_7, bZIP, FAD_binding_4, STE_STE11, CAMK_CAMKL-CHK1 and C2 gene families were more likely to play a role in response to drought and salt stress, compared with intron-rich genes in the same gene families, whereas intronless genes in the B_lectin and S_locus_glycop gene family were more likely to participate in epigenetic processes and plant development. Understanding the origin and evolutionary trajectory, as well as the potential functions, of intronless and intron-poor sub-families provides further insight into plant genome evolution and the functional divergence of genes.

A soybean calcineurin B-like protein-interacting protein kinase, GmPKS4, regulates plant responses to salt and alkali stresses

Calcineurin B-like protein-interacting protein kinases (CIPKs) are key elements of plant abiotic stress signaling pathways. CIPKs are SOS2 (Salt Overly Sensitive 2)-like proteins (protein kinase S [PKS] proteins) which all contain a putative FISL motif. It seems that the FISL motif is found only in the SOS2 subfamily of protein kinases. In this study, the full-length cDNA of a soybean CIPK gene (GmPKS4) was isolated and was revealed to have an important role in abiotic stress responses. A qRT-PCR analysis indicated that GmPKS4 expression is upregulated under saline conditions or when exposed to alkali, salt-alkali, drought, or abscisic acid (ABA). A subcellular localization assay revealed the presence of GmPKS4 in the nucleus and cytoplasm. Further studies on the GmPKS4 promoter suggested it affects soybean resistance to various stresses. Transgenic Arabidopsis thaliana and soybean hairy roots overexpressing GmPKS4 had increased proline content as well as high antioxidant enzyme activities but decreased malondialdehyde levels following salt and salt-alkali stress treatments. Additionally, GmPKS4 overexpression activated reactive oxygen species scavenging systems, thereby minimizing damages due to oxidative and osmotic stresses. Moreover, upregulated stress-related gene expression levels were detected in lines overexpressing GmPKS4 under stress conditions. In conclusion, GmPKS4 improves soybean tolerance to salt and salt-alkali stresses. The overexpression of GmPKS4 enhances the scavenging of reactive oxygen species, osmolyte synthesis, and the transcriptional regulation of stress-related genes.

Using Network-Based Machine Learning to Predict Transcription Factors Involved in Drought Resistance

Gene regulatory networks underpin stress response pathways in plants. However, parsing these networks to prioritize key genes underlying a particular trait is challenging. Here, we have built the Gene Regulation and Association Network (GRAiN) of rice (Oryza sativa). GRAiN is an interactive query-based web-platform that allows users to study functional relationships between transcription factors (TFs) and genetic modules underlying abiotic-stress responses. We built GRAiN by applying a combination of different network inference algorithms to publicly available gene expression data. We propose a supervised machine learning framework that complements GRAiN in prioritizing genes that regulate stress signal transduction and modulate gene expression under drought conditions. Our framework converts intricate network connectivity patterns of 2160 TFs into a single drought score. We observed that TFs with the highest drought scores define the functional, structural, and evolutionary characteristics of drought resistance in rice. Our approach accurately predicted the function of OsbHLH148 TF, which we validated using in vitro protein-DNA binding assays and mRNA sequencing loss-of-function mutants grown under control and drought stress conditions. Our network and the complementary machine learning strategy lends itself to predicting key regulatory genes underlying other agricultural traits and will assist in the genetic engineering of desirable rice varieties.

New insights into the evolution and functional divergence of the CIPK gene family in Saccharum

Background

Calcineurin B-like protein (CBL)-interacting protein kinases (CIPKs) are the primary components of calcium sensors, and play crucial roles in plant developmental processes, hormone signaling transduction, and in the response to exogenous stresses.

Results

In this study, 48 CIPK genes (SsCIPKs) were identified from the genome of Saccharum spontaneum. Phylogenetic reconstruction suggested that the SsCIPK gene family may have undergone six gene duplication events from the last common ancestor (LCA) of SsCIPKs. Whole-genome duplications (WGDs) served as the driving force for the amplification of SsCIPKs. The Nonsynonymous to synonymous substitution ratio (Ka/Ks) analysis showed that the duplicated genes were possibly under strong purifying selection pressure. The divergence time of these duplicated genes had an average duplication time of approximately 35.66 Mya, suggesting that these duplication events occurred after the divergence of the monocots and eudicots (165 Mya). The evolution of gene structure analysis showed that the SsCIPK family genes may involve intron losses. Ten ScCIPK genes were amplified from sugarcane (Saccharum spp. hybrids). The results of real-time quantitative polymerase chain reaction (qRT-PCR) demonstrated that these ten ScCIPK genes had different expression patterns under abscisic acid (ABA), polyethylene glycol (PEG), and sodium chloride (NaCl) stresses. Prokaryotic expression implied that the recombinant proteins of ScCIPK3, − 15 and − 17 could only slightly enhance growth under salinity stress conditions, but the ScCIPK21 did not. Transient N. benthamiana plants overexpressing ScCIPKs demonstrated that the ScCIPK genes were involved in responding to external stressors through the ethylene synthesis pathway as well as to bacterial infections.

Conclusions

In generally, a comprehensive genome-wide analysis of evolutionary relationship, gene structure, motif composition, and gene duplications of SsCIPK family genes were performed in S. spontaneum. The functional study of expression patterns in sugarcane and allogenic expressions in E. coli and N. benthamiana showed that ScCIPKs played various roles in response to different stresses. Thus, these results improve our understanding of the evolution of the CIPK gene family in sugarcane as well as provide a basis for in-depth functional studies of CIPK genes in sugarcane.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07264-9.

CaCBL1 Acts as a Positive Regulator in Pepper Response to Ralstonia solanacearum

Bacterial wilt caused by Ralstonia solanacearum is an important disease of pepper (Capsicum annuum), an economically important solanaceous vegetable worldwide, in particular, under high temperature (HT) conditions. However, the molecular mechanism underlying pepper immunity against bacterial wilt remains poorly understood. Herein, CaCBL1, a putative calcineurin B-like protein, was functionally characterized in the pepper response to R. solanacearum inoculation (RSI) under HT (RSI/HT). CaCBL1 was significantly upregulated by RSI at room temperature (RSI/RT), HT, or RSI/HT. CaCBL1-GFP fused protein targeted to whole epidermal cells of Nicotiana benthamiana when transiently overexpressed. CaCBL1 silencing by virus-induced gene silencing significantly enhanced pepper susceptibility to RSI under RT or HT, while its transient overexpression triggered hypersensitive response mimic cell death and upregulation of immunity-associated marker genes, including CabZIP63, CaWRKY40, and CaCDPK15, the positive regulators in the pepper response to RSI or HT found in our previous studies. In addition, by chromatin immunoprecipitation PCR and electrophoretic mobility shift assay, CaCBL1 was found to be directly targeted by CaWRKY40, although not by CaWRKY27 or CaWRKY58, via the W-box-2 within its promoter, and its transcription was found to be downregulated by silencing of CaWRKY40 while it was enhanced by its transient overexpression. These results suggest that CaCBL1 acts as a positive regulator in pepper immunity against R. solanacearum infection, constituting a positive feedback loop with CaWRKY40.

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.

Foxtail millet (Setaria italica (L.) P. Beauv) CIPKs are responsive to ABA and abiotic stresses

CBL-interacting protein kinases (CIPKs) have been shown to regulate a variety of environmental stress-related signalling pathways in plants. Foxtail millet (Setaria italica (L.) P. Beauv) is known worldwide as a relatively stress-tolerant C4 crop species. Although the foxtail millet genome sequence has been released, little is known about the functions of CIPKs in foxtail millet. Therefore, a systematic genome-wide analysis of CIPK genes in foxtail millet was performed. In total, 35 CIPK members were identified in foxtail millet and divided into four subgroups (I to IV) on the basis of their phylogenetic relationships. Phylogenetic and gene structure analyses clearly divided all SiCIPKs into intron-poor and intron-rich clades. Cis-element analysis subsequently indicated that these SiCIPKs may be involved in responses to abiotic stimuli, hormones, and light signalling during plant growth and development, and stress-induced expression profile analysis revealed that all the SiCIPKs are involved in various stress signalling pathways. These results suggest that the CIPK genes in foxtail millet exhibit the basic characteristics of CIPK family members and play important roles in response to abiotic stresses. The results of this study will contribute to future functional characterization of abiotic stress responses mediated by CIPKs in foxtail millet.

Adaptation of Plants to Salt Stress: Characterization of Na+ and K+ Transporters and Role of CBL Gene Family in Regulating Salt Stress Response

Salinity is one of the most serious factors limiting the productivity of agricultural crops, with adverse effects on germination, plant vigor, and crop yield. This salinity may be natural or induced by agricultural activities such as irrigation or the use of certain types of fertilizer. The most detrimental effect of salinity stress is the accumulation of Na+ and Cl− ions in tissues of plants exposed to soils with high NaCl concentrations. The entry of both Na+ and Cl− into the cells causes severe ion imbalance, and excess uptake might cause significant physiological disorder(s). High Na+ concentration inhibits the uptake of K+, which is an element for plant growth and development that results in lower productivity and may even lead to death. The genetic analyses revealed K+ and Na+ transport systems such as SOS1, which belong to the CBL gene family and play a key role in the transport of Na+ from the roots to the aerial parts in the Arabidopsis plant. In this review, we mainly discuss the roles of alkaline cations K+ and Na+, Ion homeostasis-transport determinants, and their regulation. Moreover, we tried to give a synthetic overview of soil salinity, its effects on plants, and tolerance mechanisms to withstand stress.

Genome-Wide Identification, Characterization, and Expression Analysis of the Grapevine Superoxide Dismutase (SOD) Family

Superoxide dismutase (SOD) is an essential enzyme of the plant antioxidant system that responds to oxidative damage caused by adverse conditions. However, little is known about the SOD gene family in Vitis vinifera (Vv). In the present study, ten SOD genes, including 6 copper/zinc SODs, 2 iron SODs, and 2 manganese SODs, were identified in the grapevine genome where they were unevenly distributed on 12 chromosomes. Ten VvSOD genes were divided into three main groups based on phylogenetic analysis, subcellular localization, and the distribution of conserved protein motifs. Additionally, many cis-elements related to different stresses were found in the promoters of the 10 VvSOD genes. Syntenic analysis revealed that VvMSD1 and VvMSD2 were derived from segmental duplication, and VvCSD4 and VvCSD5 belong to a pair of tandemly duplicated genes. Gene expression analysis based on microarray data showed that the 10 VvSOD genes were expressed in all the tested tissues. Interestingly, the segmentally duplicated gene pair (VvMSD1 and VvMSD2) exhibited differential expression patterns in various organs. In contrast, the tandemly duplicated gene pair (VvCSD4 and VvCSD5) displayed similar expression patterns in the tested organs. Our results provide a basis for further functional research on the SOD gene family in grapevine.

The Complex Fine-Tuning of K⁺ Fluxes in Plants in Relation to Osmotic and Ionic Abiotic Stresses

As the main cation in plant cells, potassium plays an essential role in adaptive responses, especially through its involvement in osmotic pressure and membrane potential adjustments. K⁺ homeostasis must, therefore, be finely controlled. As a result of different abiotic stresses, especially those resulting from global warming, K⁺ fluxes and plant distribution of this ion are disturbed. The hormone abscisic acid (ABA) is a key player in responses to these climate stresses. It triggers signaling cascades that ultimately lead to modulation of the activities of K⁺ channels and transporters. After a brief overview of transcriptional changes induced by abiotic stresses, this review deals with the post-translational molecular mechanisms in different plant organs, in Arabidopsis and species of agronomical interest, triggering changes in K⁺ uptake from the soil, K⁺ transport and accumulation throughout the plant, and stomatal regulation. These modifications involve phosphorylation/dephosphorylation mechanisms, modifications of targeting, and interactions with regulatory partner proteins. Interestingly, many signaling pathways are common to K⁺ and Cl^-/NO3^- counter-ion transport systems. These cross-talks are also addressed.

Circadian Regulation of Alternative Splicing of Drought-Associated CIPK Genes in Dendrobium catenatum (Orchidaceae)

Dendrobium catenatum, an epiphytic and lithophytic species, suffers frequently from perennial shortage of water in the wild. The molecular mechanisms of this orchid's tolerance to abiotic stress, especially drought, remain largely unknown. It is well-known that CBL-interacting protein kinase (CIPKs) proteins play important roles in plant developmental processes, signal transduction, and responses to abiotic stress. To study the CIPKs' functions for D. catenatum, we first identified 24 CIPK genes from it. We divided them into three subgroups, with varying intron numbers and protein motifs, based on phylogeny analysis. Expression patterns of CIPK family genes in different tissues and in response to either drought or cold stresses suggested DcaCIPK11 may be associated with signal transduction and energy metabolism. DcaCIPK9, -14, and -16 are predicted to play critical roles during drought treatment specifically. Furthermore, transcript expression abundances of DcaCIPK16 showed polar opposites during day and night. Whether under drought treatment or not, DcaCIPK16 tended to emphatically express transcript1 during the day and transcript3 at night. This implied that expression of the transcripts might be regulated by circadian rhythm. qRT-PCR analysis also indicated that DcaCIPK3, -8, and -20 were strongly influenced by circadian rhythmicity. In contrast with previous studies, for the first time to our knowledge, our study revealed that the major CIPK gene transcript expressed was not always the same and was affected by the biological clock, providing a different perspective on alternative splicing preference.

The kinome of pineapple: catalog and insights into functions in crassulacean acid metabolism plants

BACKGROUND

Crassulacean acid metabolism (CAM) plants use water 20-80% more efficiently by shifting stomata opening and primary CO₂ uptake and fixation to the nighttime. Protein kinases (PKs) play pivotal roles in this biological process. However, few PKs have been functionally analyzed precisely due to their abundance and potential functional redundancy (caused by numerous gene duplications).

RESULTS

In this study, we systematically identified a total of 758 predicted PK genes in the genome of a CAM plant, pineapple (Ananas comosus). The pineapple kinome was classified into 20 groups and 116 families based on the kinase domain sequences. The RLK was the largest group, containing 480 members, and over half of them were predicted to locate at the plasma membrane. Both segmental and tandem duplications make important contributions to the expansion of pineapple kinome based on the synteny analysis. Ka/Ks ratios showed all of the duplication events were under purifying selection. The global expression analysis revealed that pineapple PKs exhibit different tissue-specific and diurnal expression patterns. Forty PK genes in a cluster performed higher expression levels in green leaf tip than in white leaf base, and fourteen of them had strong differential expression patterns between the photosynthetic green leaf tip and the non-photosynthetic white leaf base tissues.

CONCLUSIONS

Our findings provide insights into the evolution and biological function of pineapple PKs and a foundation for further functional analysis of PKs in CAM plants. The gene duplication, expression, and coexpression analysis helped us to rapidly identify the key candidates in pineapple kinome, which may play roles in the carbon fixation process in pineapple and help engineering CAM pathway into C3 crops for improved drought tolerance.

Overexpression of the Tibetan Plateau annual wild barley (Hordeum spontaneum) HsCIPKs enhances rice tolerance to heavy metal toxicities and other abiotic stresses

BACKGROUND

The calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) signaling system plays a key regulatory role in plant stress signaling. The roles of plant-specific CIPKs, essential for CBL-CIPK functions, in the response to various abiotic stresses have been extensively studied so far. However, until now, the possible roles of the CIPKs in the plant response to heavy metal toxicities are largely unknown.

RESULTS

In this study, we used bioinformatic and molecular strategies to isolate 12 HsCIPK genes in Tibetan Plateau annual wild barley (Hordeum spontaneum C. Koch) and subsequently identified their functional roles in the response to heavy metal toxicities. The results showed that multiple HsCIPKs were transcriptionally regulated by heavy metal toxicities (e.g., Hg, Cd, Cr, Pb, and Cu) and other abiotic stresses (e.g., salt, drought, aluminum, low and high temperature, and abscisic acid). Furthermore, the ectopic overexpression of each HsCIPK in rice (Oryza sativa L. cv Nipponbare) showed that transgenic plants of multiple HsCIPKs displayed enhanced tolerance of root growth to heavy metal toxicities (Hg, Cd, Cr, and Cu), salt and drought stresses. These results suggest that HsCIPKs are involved in the response to heavy metal toxicities and other abiotic stresses.

CONCLUSIONS

Tibetan Plateau annual wild barley HsCIPKs possess broad applications in genetically engineering of rice with tolerance to heavy metal toxicities and other abiotic stresses.

Arabidopsis RNA processing factor SERRATE regulates the transcription of intronless genes

Intron splicing increases proteome complexity, promotes RNA stability, and enhances transcription. However, introns and the concomitant need for splicing extend the time required for gene expression and can cause an undesirable delay in the activation of genes. Here, we show that the plant microRNA processing factor SERRATE (SE) plays an unexpected and pivotal role in the regulation of intronless genes. Arabidopsis SE associated with more than 1000, mainly intronless, genes in a transcription-dependent manner. Chromatin-bound SE liaised with paused and elongating polymerase II complexes and promoted their association with intronless target genes. Our results indicate that stress-responsive genes contain no or few introns, which negatively affects their expression strength, but that some genes circumvent this limitation via a novel SE-dependent transcriptional activation mechanism. Transcriptome analysis of a Drosophila mutant defective in ARS2, the metazoan homologue of SE, suggests that SE/ARS2 function in regulating intronless genes might be conserved across kingdoms.

Identification of a novel fused gene family implicates convergent evolution in eukaryotic calcium signaling

Background

Both calcium signals and protein phosphorylation responses are universal signals in eukaryotic cell signaling. Currently three pathways have been characterized in different eukaryotes converting the Ca²⁺ signals to the protein phosphorylation responses. All these pathways have based mostly on studies in plants and animals.

Results

Based on the exploration of genomes and transcriptomes from all the six eukaryotic supergroups, we report here in Metakinetoplastina protists a novel gene family. This family, with a proposed name SCAMK, comprises SnRK3 fused calmodulin-like III kinase genes and was likely evolved through the insertion of a calmodulin-like3 gene into an SnRK3 gene by unequal crossover of homologous chromosomes in meiosis cell. Its origin dated back to the time intersection at least 450 million-year-ago when Excavata parasites, Vertebrata hosts, and Insecta vectors evolved. We also analyzed SCAMK’s unique expression pattern and structure, and proposed it as one of the leading calcium signal conversion pathways in Excavata parasite. These characters made SCAMK gene as a potential drug target for treating human African trypanosomiasis.

Conclusions

This report identified a novel gene fusion and dated its precise fusion time in Metakinetoplastina protists. This potential fourth eukaryotic calcium signal conversion pathway complements our current knowledge that convergent evolution occurs in eukaryotic calcium signaling.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4685-y) contains supplementary material, which is available to authorized users.

Expression Patterns and Identified Protein-Protein Interactions Suggest That Cassava CBL-CIPK Signal Networks Function in Responses to Abiotic Stresses

Cassava is an energy crop that is tolerant of multiple abiotic stresses. It has been reported that the interaction between Calcineurin B-like (CBL) protein and CBL-interacting protein kinase (CIPK) is implicated in plant development and responses to various stresses. However, little is known about their functions in cassava. Herein, 8 CBL (MeCBL) and 26 CIPK (MeCIPK) genes were isolated from cassava by genome searching and cloning of cDNA sequences of Arabidopsis CBLs and CIPKs. Reverse-transcriptase polymerase chain reaction (RT-PCR) analysis showed that the expression levels of MeCBL and MeCIPK genes were different in different tissues throughout the life cycle. The expression patterns of 7 CBL and 26 CIPK genes in response to NaCl, PEG, heat and cold stresses were analyzed by quantitative real-time PCR (qRT-PCR), and it was found that the expression of each was induced by multiple stimuli. Furthermore, we found that many pairs of CBLs and CIPKs could interact with each other via investigating the interactions between 8 CBL and 25 CIPK proteins using a yeast two-hybrid system. Yeast cells co-transformed with cassava MeCIPK24, MeCBL10, and Na⁺/H⁺ antiporter MeSOS1 genes exhibited higher salt tolerance compared to those with one or two genes. These results suggest that the cassava CBL-CIPK signal network might play key roles in response to abiotic stresses.

The grapevine kinome: annotation, classification and expression patterns in developmental processes and stress responses

Protein kinases (PKs) have evolved as the largest family of molecular switches that regulate protein activities associated with almost all essential cellular functions. Only a fraction of plant PKs, however, have been functionally characterized even in model plant species. In the present study, the entire grapevine kinome was identified and annotated using the most recent version of the grapevine genome. A total of 1168 PK-encoding genes were identified and classified into 20 groups and 121 families, with the RLK-Pelle group being the largest, with 872 members. The 1168 kinase genes were unevenly distributed over all 19 chromosomes, and both tandem and segmental duplications contributed to the expansion of the grapevine kinome, especially of the RLK-Pelle group. Ka/Ks values indicated that most of the tandem and segmental duplication events were under purifying selection. The grapevine kinome families exhibited different expression patterns during plant development and in response to various stress treatments, with many being coexpressed. The comprehensive annotation of grapevine kinase genes, their patterns of expression and coexpression, and the related information facilitate a more complete understanding of the roles of various grapevine kinases in growth and development, responses to abiotic stress, and evolutionary history.

Analysis of EF-Hand Proteins in Soybean Genome Suggests Their Potential Roles in Environmental and Nutritional Stress Signaling

Calcium ion (Ca2+) is a universal second messenger that plays a critical role in plant responses to diverse physiological and environmental stimuli. The stimulus-specific signals are perceived and decoded by a series of Ca2+ binding proteins serving as Ca2+ sensors. The majority of Ca2+ sensors possess the EF-hand motif, a helix-loop-helix structure which forms a turn-loop structure. Although EF-hand proteins in model plant such as Arabidopsis have been well described, the identification, classification, and the physiological functions of EF-hand-containing proteins from soybean are not systemically reported. In this study, a total of at least 262 genes possibly encoding proteins containing one to six EF-hand motifs were identified in soybean genome. These genes include 6 calmodulins (CaMs), 144 calmodulin-like proteins (CMLs), 15 calcineurin B-like proteins, 50 calcium-dependent protein kinases (CDPKs), 13 CDPK-related protein kinases, 2 Ca2+- and CaM-dependent protein kinases, 17 respiratory burst oxidase homologs, and 15 unclassified EF-hand proteins. Most of these genes (87.8%) contain at least one kind of hormonal signaling- and/or stress response-related cis-elements in their -1500 bp promoter regions. Expression analyses by exploring the published microarray and Illumina transcriptome sequencing data revealed that the expression of these EF-hand genes were widely detected in different organs of soybean, and nearly half of the total EF-hand genes were responsive to various environmental or nutritional stresses. Quantitative RT-PCR was used to confirm their responsiveness to several stress treatments. To confirm the Ca2+-binding ability of these EF-hand proteins, four CMLs (CML1, CML13, CML39, and CML95) were randomly selected for SDS-PAGE mobility-shift assay in the presence and absence of Ca2+. Results showed that all of them have the ability to bind Ca2+. This study provided the first comprehensive analyses of genes encoding for EF-hand proteins in soybean. Information on the classification, phylogenetic relationships and expression profiles of soybean EF-hand genes in different tissues and under various environmental and nutritional stresses will be helpful for identifying candidates with potential roles in Ca2+ signal-mediated physiological processes including growth and development, plant-microbe interactions and responses to biotic and abiotic stresses.