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

Deciphering the Atlas of Post-Translational Modification in Sugarcane

In plants, lysine acetylation (Kac), 2-hydroxyisobutyrylation (Khib), and lysine lactylation (Kla), the three new types of post-translational modification (PTM), play very important roles in growth, development, and resistance to adverse environmental stresses. Herein, we report the first global acetylome, 2-hydroxyisobutyrylome, and lactylome in sugarcane. A total of 8573 Kac, 4637 Khib, and 215 Kla sites across 3903, 1507, and 139 modified proteins were identified. Besides, homology analyses revealed the Kac, Khib, and Kla sites on histones were conserved between sugarcane and rice or poplar. Functional annotations demonstrated that the Kac, Khib, and Kla proteins were mainly involved in energy metabolism. In addition, a number of modified transcription factors and stress-related proteins, which were constitutively expressed in different tissues of sugarcane and induced by drought, cold or Sporisorium scitamineum stress, were identified. Finally, a proposed working mode on how PTM functions in sugarcane was depicted. We thus concluded that PTM should play a role in sugarcane growth, development, and response to biotic and abiotic stresses, but the mechanisms require further investigation. The present study provided the all-new comprehensive profile of proteins Kac, Khib, and Kla and a new perspective to understand the molecular mechanisms of protein PTMs in sugarcane.

Genome-wide survey of catalase genes in Brassica rapa, Brassica oleracea, and Brassica napus: identification, characterization, molecular evolution, and expression profiling of BnCATs in response to salt and cadmium stress

Catalase (CAT, EC 1.11.1.6), one of the most important antioxidant enzymes, can control excess levels of H₂O₂ produced under oxidative stress in plants. In this study, 16, 8, and 7 CAT genes in the genome of Brassica napus, B. rapa, and B. oleracea were identified, respectively. Phylogenetic studies showed that CATs could be divided into two main groups, each containing specific monocotyledon and dicotyledon subgroups. Motifs, gene structure, and intron phase of CATs in B. napus, Brassica rapa, and Brassica oleracea are highly conserved. Analysis of codon usage bias showed the mutation pressure and natural selection of the codon usage of CATs. Segmental duplication and polyploid were major factors in the expansion of this gene family in B. napus, and genes have experienced  negative selection during evolution. Existence of hormones and stress-responsive cis-elements and identifying miRNA molecules affecting CATs showed that these genes are complexly regulated at the transcriptional and posttranscriptional levels. Based on RNA-seq data, CATs are divided into two groups; the first group has moderate and specific expression in flowers, leaves, stems, and roots, while the second group shows expression in most tissues. qRT-PCR analysis showed that the expression of these genes is dynamic and has a specific expression consistent with other CAT genes in response to salinity and cadmium (Cd) stresses. These results provide information for further investigation of the function of CAT genes in response to stresses and the development of tolerant  plants.

Comprehensive Analysis of Calcium Sensor Families, CBL and CIPK, in Aeluropus littoralis and Their Expression Profile in Response to Salinity

Plants have acquired sets of highly regulated and complex signaling pathways to respond to unfavorable environmental conditions during evolution. Calcium signaling, as a vital mechanism, enables plants to respond to external stimuli, including abiotic and biotic stresses, and coordinate the basic processes of growth and development. In the present study, two calcium sensor families, CBL and CIPK, were investigated in a halophyte plant, Aeluropus littoralis, with a comprehensive analysis. Here, six AlCBL genes, and twenty AlCIPK genes were studied. The analysis of the gene structure and conserved motifs, as well as physicochemical properties, showed that these genes are highly conserved during evolution. The expression levels of AlCBL genes and AlCIPK genes were evaluated under salt stress in leaf and root tissue. Based on the real-time RT-PCR results, the AlCIPK gene family had a higher variation in mRNA abundance than the AlCBL gene family. AlCIPK genes were found to have a higher abundance in leaves than in roots. The results suggest that the correlation between AlCBL genes and AlCIPK is tissue-specific, and different correlations can be expected in leaves and roots. Based on these correlations, AlCIPK3.1-AlCBL4.1 and AlCIPK1.2-AlCBL4.4 can be co-expressed in the root tissue, while AlCBL10 has the potential to be co-expressed with AlCIPK5, AlCIPK26, and AlCIPK12.3 in the leaf tissue. Our findings reveal valuable information on the structure and function of calcium sensor families in A. littoralis, a halophyte plant, that can be used in future research on the biological function of CBLs and CIPKs on salt stress resistance.

Screening of Candidate Genes Associated with Brown Stripe Resistance in Sugarcane via BSR-seq Analysis

Sugarcane brown stripe (SBS), caused by the fungal pathogen Helminthosporium stenospilum, is one of the most serious threats to sugarcane production. However, its outbreaks and epidemics require suitable climatic conditions, resulting in the inefficient improvement of the SBS resistance by phenotype selection. The sugarcane F₁ population of SBS-resistant YT93-159 × SBS-susceptible ROC22 was used for constructing the bulks. Bulked segregant RNA-seq (BSR-seq) was then performed on the parents YT93-159 (T01) and ROC22 (T02), and the opposite bulks of 30 SBS-susceptible individuals mixed bulk (T03) and 30 SBS-resistant individuals mixed bulk (T04) collected from 287 F₁ individuals. A total of 170.00 Gb of clean data containing 297,921 SNPs and 70,426 genes were obtained. Differentially expressed genes (DEGs) analysis suggested that 7787 and 5911 DEGs were identified in the parents (T01 vs. T02) and two mixed bulks (T03 vs. T04), respectively. In addition, 25,363 high-quality and credible SNPs were obtained using the genome analysis toolkit GATK for SNP calling. Subsequently, six candidate regions with a total length of 8.72 Mb, which were located in the chromosomes 4B and 7C of sugarcane wild species Saccharum spontaneum, were identified, and 279 genes associated with SBS-resistance were annotated by ED algorithm and ΔSNP-index. Furthermore, the expression profiles of candidate genes were verified by quantitative real-time PCR (qRT-PCR) analysis, and the results showed that eight genes (LRR-RLK, DHAR1, WRKY7, RLK1, BLH4, AK3, CRK34, and NDA2) and seven genes (WRKY31, CIPK2, CKA1, CDPK6, PFK4, CBL2, and PR2) of the 20 tested genes were significantly up-regulated in YT93-159 and ROC22, respectively. Finally, a potential molecular mechanism of sugarcane response to H. stenospilum infection is illustrate that the activations of ROS signaling, MAPK cascade signaling, Ca²⁺ signaling, ABA signaling, and the ASA-GSH cycle jointly promote the SBS resistance in sugarcane. This study provides abundant gene resources for the SBS resistance breeding in sugarcane.

Genetic identification of SNP markers and candidate genes associated with sugarcane smut resistance using BSR-Seq

Sugarcane smut caused by Sporisorium scitamineum is one of the most severe fungal diseases worldwide. In this study, a cross was made between a smut-resistant variety YT93-159 and a smut-susceptible variety ROC22, and 312 progenies were obtained. Two bulks of progenies were then constructed, one consisted of 27 highly smut resistant progenies and the other 24 smut susceptible progenies. Total RNAs of the progenies of each bulk, were pooled and subject to bulked segregant RNA-sequence analysis (BSR-Seq). A total of 164.44 Gb clean data containing 2,341,449 SNPs and 64,999 genes were obtained, 7,295 of which were differentially expressed genes (DEGs). These DEGs were mainly enriched in stress-related metabolic pathways, including carbon metabolism, phenylalanine metabolism, plant hormone signal transduction, glutathione metabolism, and plant-pathogen interactions. Besides, 45,946 high-quality, credible SNPs, a 1.27 Mb region at Saccharum spontaneum chromosome Chr5B (68,904,827 to 70,172,982), and 129 candidate genes were identified to be associated with smut resistance. Among them, twenty-four genes, either encoding key enzymes involved in signaling pathways or being transcription factors, were found to be very closely associated with stress resistance. RT-qPCR analysis demonstrated that they played a positive role in smut resistance. Finally, a potential molecular mechanism of sugarcane and S. scitamineum interaction is depicted that activations of MAPK cascade signaling, ROS signaling, Ca²⁺ signaling, and PAL metabolic pathway and initiation of the glyoxalase system jointly promote the resistance to S. scitamineum in sugarcane. This study provides potential SNP markers and candidate gene resources for smut resistance breeding in sugarcane.

AcCIPK5, a pineapple CBL-interacting protein kinase, confers salt, osmotic and cold stress tolerance in transgenic Arabidopsis

Plant-specific calcineurin B-like proteins (CBLs) and their interacting kinases, CBL-interacting protein kinases (CIPKs) module, are essential for dealing with various biotic and abiotic stress. The kinases (CIPKs) of this module have been well studied in several plants; however, the information about pineapple CIPKs remains limited. To understand how CIPKs function against environmental cues in pineapple, the CIPK5 gene of pineapple was cloned and characterized. The phylogenetic analyses revealed that AcCIPK5 is homologous to the CIPK12 of Arabidopsis and other plant species. Quantitative real-time PCR (qRT-PCR) analysis revealed that AcCIPK5 responds to multiple stresses, including osmotic, salt stress, heat and cold. Under optimal conditions, AcCIPK5 gets localized to the cytoplasm and cell membrane. The ectopic expression of AcCIPK5 in Arabidopsis improved the germination under osmotic and salt stress. Furthermore, AcCIPK5 positively regulated osmotic, drought, salt and cold tolerance and negatively regulated heat and fungal stress in Arabidopsis. Besides, the expression of AcCIPK impacted ABA-related genes and ROS homeostasis. Overall, the present study demonstrates that AcCIPK5 contributes to multiple stress tolerance and has the potential to be utilized in the development of stress-tolerant crops.

Genome-Wide Identification of the Salvia miltiorrhiza SmCIPK Gene Family and Revealing the Salt Resistance Characteristic of SmCIPK13

Members of the CIPK (CBL-interacting protein kinases) gene family play important roles in calcium (Ca²⁺) signaling pathway-regulated plant resistance to abiotic stresses. Salvia miltiorrhiza, which is widely planted and grown in complex and diverse environments, is mainly focused on the transcriptional regulation of enzyme genes related to the biosynthesis of its bioactive components. However, the excavation of the genes related to the resistance of S.miltiorrhiza and the involved signaling pathways have not been deeply studied. In this study, 20 SmCIPK genes were identified and classified into two families and five subfamilies by biochemical means. Sequence characteristics and conserved motif analysis revealed the conservation and difference of SmCIPK protein in plants. Expression pattern analysis showed that SmCIPKs were mainly expressed in flowers and roots, and more than 90% of gene expression was induced by SA (salicylic acid), and MeJA (methyl jasmonate). Furthermore, the expression level of SmCIPK13 could be significantly increased after stress treatment with NaCl. SmCIPK13 expression in yeast reduces sensitivity to salt, while overexpression of it in Arabidopsis has the same effect and was localized in the cytoplasm, cell membrane and nucleus. In conclusion, the identification of the SmCIPK gene family and the functional characterization of the SmCIPK13 gene provides the basis for clarification of key genes in the Ca²⁺ signaling pathway and abiotic stress in S.miltiorrhiza.

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.

Characterization of silicon transporter gene family in Saccharum and functional analysis of the ShLsi6 gene in biotic stress

Silicon, one of the most prevalent elements in the soil, is beneficial for plant growth and defense against different stresses. The silicon transporter gene (Lsi) plays an important role in the uptake and transport of silicon in higher plants. In this study, a total of 32 Lsi genes, including 20 SsLsi in sugarcane wild species Saccharum spontaneum, 5 ShLsi in Saccharum hybrid cultivar R570 and 7 SbLsi in sugarcane related species Sorghum bicolor, were identified and classified into three groups. Bioinformatics analysis showed that instability, hydrophobicity, localization of cell membranes and vacuoles were the main features of the Lsi proteins. Whole genome and segmental duplication contributed to the main expansion of Lsi gene family. Collinearity analysis of the Lsi genes showed that S. spontanum and R570 had a collinear relationship with monocotyledonous plants S. bicolor and Oryza sativa, but not with dicotyledonous plants Arabidopsis thaliana and Vitis vinifera. The replicated Lsi genes were mainly subjected to strong selection pressure for purification. The diverse cis-regulatory elements in the promoter of SsLsi, ShLsi and SbLsi genes suggested that they were widely involved in the response of plants to various stresses and the regulation of the growth and development. Transcriptome data and real time quantitative PCR analysis showed that the Lsi genes exhibited different expression profiles in sugarcane tissues and under Sporisorium scitamineum, drought and cold stresses. In addition, the cDNA and genomic DNA sequences of ShLsi6 that was homologous to SsLsi1b gene was cloned from Saccharum hybrid cultivar ROC22. Transient expression analysis showed that, compared with the control, Nicotiana benthamiana leaves which overexpressed the ShLsi6 gene showed a high sensitivity after inoculation with tobacco pathogens Ralstonia solanacearum and Fusarium solani var. coeruleum. This study provides important information for further functional analysis of Lsi genes and resistant breeding in sugarcane.

New functions of CIPK gene family are continue to emerging

CIPK protein family is a key protein family in Ca²⁺ mediated plant signaling pathway, which plays an indispensable role in plant response to stress and development. Every gene in this family encodes specific proteins. They interact with calcium ion signals, make plants to deal with various stress or stimuli. This article mainly reviews the mechanism, positioning and physiological functions of the CIPK family in different species in recent years. According to our team's research, CIPK8 interacts with CBL5 to improve salt tolerance, and CIPK23 interacts with TGA1 to regulate nitrate uptake negatively in chrysanthemum. In addition, we discussed current limitations and future research directions. The article will enhance the understanding of the functional characteristics of the CIPK gene family under different stresses, provide insights for future breeding and the development of new crop varieties with enhanced stress tolerance.

Interspecific complementation-restoration of phenotype in Arabidopsis cuc2cuc3 mutant by sugarcane CUC2 gene

BACKGROUND

In plants, a critical balance between differentiation and proliferation of stem cells at the shoot apical meristem zone is essential for proper growth. The spatiotemporal regulation of some crucial genes dictates the formation of a boundary within and around budding organs. The boundary plays a pivotal role in distinguishing one tissue type from another and provides a defined shape to the organs at their developed stage. NAM/CUC subfamily of the NAC transcription factors control the boundary formation during meristematic development.

RESULTS

Here, we have identified the CUP-SHAPED COTYLEDON (CUC) genes in sugarcane and named SsCUC2 (for the orthologous gene of CUC1 and CUC2) and SsCUC3. The phylogenetic reconstruction showed that SsCUCs occupy the CUC2 and CUC3 clade together with monocots, whereas eudicot CUC2 and CUC3 settled separately in the different clade. The structural analysis of CUC genes showed that most of the CUC3 genes were accompanied by an intron gain during eudicot divergence. Besides, the study of SsCUCs expression in the RNA-seq obtained during different stages of ovule development revealed that SsCUCs express in developing young tissues, and the expression of SsCUC2 is regulated by miR164. We also demonstrate that SsCUC2 (a monocot) could complement the cuc2cuc3 mutant phenotype of Arabidopsis (eudicot).

CONCLUSIONS

This study further supports that CUC2 has diverged in CUC1 and CUC2 during the evolution of monocots and eudicots from ancestral plants. The functional analysis of CUC expression patterns during sugarcane ovule development and ectopic expression of SsCUC2 in Arabidopsis showed that SsCUC2 has a conserved role in boundary formation. Overall, these findings improve our understanding of the functions of sugarcane CUC genes. Our results reveal the crucial functional role of CUC genes in sugarcane.

Genome-wide identification and expression analysis of the coronatine-insensitive 1 (COI1) gene family in response to biotic and abiotic stresses in Saccharum

Background

The coronatine insensitive 1 (COI1) gene is the core member of jasmonate signaling pathway, which is closely related to plant biotic and abiotic resistance. However, there have been no reports on COI1 in sugarcane (Sacharum spp.). Hence, systematically investigating the characteristics of the COI1 multigene family in sugarcane can provide a means to study and manipulate the jasmonic acid signaling pathway.

Results

A total of 156 COI1 proteins were obtained from the genomes of 19 land plants, while none were obtained from five algae species. A phylogenetic tree demonstrated that these COI1 proteins were classified into four groups, while 31 proteins of SsCOI1 from Saccharum spontaneum, SbCOI1 from Sorghum bicolor, and ShCOI1 from Saccharum spp. hybrid cultivar R570 clustered into three groups. Synteny analysis and duplication patterns revealed that COI1 genes expanded through various genome replication events and could have experienced strong purifying selective pressure during evolution in S. spontaneum, S. bicolor, and R570. An investigation of cis-acting elements suggests that COI1 genes may be involved in plant growth and development and response to various stresses. Expression analysis implied that 21 SsCOI1 genes were constitutively expressed, and had positive responses to drought, cold, and Sporisorium scitamineum stresses with different expression patterns. Among them, seven SsCOI1 haplotype genes may play different roles in response to methyl jasmonate. Furthermore, the ShCOI1–4, ShCOI1–5, and ShCOI1–6 genes were cloned from Saccharum spp. hybrid cultivar ROC22. Real-time quantitative PCR (RT-qPCR) analysis demonstrated that these three ShCOI1 genes had divergent expression profiles in response to salicylic acid, abscisic acid, polyethylene glycol, cold, and S. scitamineum.

Conclusions

These results suggest that COI1 genes may act in sugarcane growth, development, and response to various stresses via different regulatory mechanisms, which laying a foundation for the functional identification of the sugarcane COI1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08255-0.

The CaCA superfamily genes in Saccharum: comparative analysis and their functional implications in response to biotic and abiotic stress

BACKGROUND

In plants, Calcium (Ca2+) acts as a universal messenger in various signal transduction pathways, including responses to biotic and abiotic stresses and regulation of cellular and developmental processes. The Ca2+/cation antiporter (CaCA) superfamily proteins play vital roles in the transport of Ca2+ and/or other cations. However, the characteristics of these superfamily members in Saccharum and their evolutionary and functional implications have remained unclear.

RESULTS

A total of 34 CaCA genes in Saccharum spontaneum, 5 CaCA genes in Saccharum spp. R570, and 14 CaCA genes in Sorghum bicolor were identified and characterized. These genes consisted of the H+/cation exchanger (CAX), cation/Ca2+ exchanger (CCX), EF-hand / CAX (EFCAX), and Mg2+/H+ exchanger (MHX) families, among which the CCX and EFCAX could be classified into three groups while the CAX could be divided into two groups. The exon/intron structures and motif compositions suggested that the members in the same group were highly conserved. Synteny analysis of CaCAs established their orthologous and paralogous relationships among the superfamily in S. spontaneum, R570, and S. bicolor. The results of protein-protein interactions indicated that these CaCA proteins had direct or indirect interactions. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis demonstrated that most members of Saccharum CaCA genes exhibited a similar expression pattern in response to hormonal (abscisic acid, ABA) treatment but played various roles in response to biotic (Sporisorium scitamineum) and abiotic (cold) stresses. Furthermore, ScCAX4, a gene encoding a cytoplasm, plasma membrane and nucleus positioning protein, was isolated from sugarcane. This gene was constitutively expressed in different sugarcane tissues and its expression was only induced at 3 and 6 h time points after ABA treatment, however was inhibited and indued in the whole process under cold and S. scitamineum stresses, respectively.

CONCLUSIONS

This study systematically conducted comparative analyses of CaCA superfamily genes among S. spontaneum, R570, and S. bicolor, delineating their sequence and structure characteristics, classification, evolutionary history, and putative functions. These results not only provided rich gene resources for exploring the molecular mechanism of the CaCA superfamily genes but also offered guidance and reference for research on other gene families in Saccharum.

Genome-wide identification, characterization and expression analysis of the carotenoid cleavage oxygenase (CCO) gene family in Saccharum

Carotenoid cleavage oxygenases (CCOs) play crucial roles in plant growth and development, as well as in the response to phytohormonal, biotic and abiotic stresses. However, comprehensive and systematic research on the CCO gene family has not yet been conducted in Saccharum. In this study, 47 SsCCO and 14 ShCCO genes were identified and characterized in Saccharum spontaneum and Saccharum spp. R570 cultivar, respectively. The SsCCOs consisted of 38 SsCCDs and 9 SsNCEDs, while ShCCOs contained 11 ShCCDs and 3 ShNCEDs. The SsCCO family could be divided into 7 groups, while ShCCO family into 5 groups. The genes/proteins contained similar compositions within the same group, and the evolutionary mechanisms differed between S. spontaneum and R570. Gene Ontology annotation implied that CCOs were involved in many physiological and biochemical processes. Additionally, 41 SsCCOs were regulated by 19 miRNA families, and 8 ShCCOs by 9 miRNA families. Cis-regulatory elements analysis suggested that CCO genes functioned in the process of growth and development or under the phytohormonal, biotic and abiotic stresses. qRT-PCR analysis indicated that nine CCO genes from different groups exhibited similar expression patterns under abscisic acid treatment, while more divergent profiles were observed in response to Sporisorium scitamineum and cold stresses. Herein, comparative genomics analysis of the CCO gene family between S. spontaneum and R570 was conducted to investigate its evolution and functions. This is the first report on the CCO gene family in S. spontaneum and R570, thus providing valuable information and facilitating further investigation into its function in the future.