Skip to main content

Comprehensive identification and expression analyses of the SnRK gene family in Casuarina equisetifolia in response to salt stress



Sucrose nonfermenting-1 (SNF1)-related protein kinases (SnRKs) play crucial roles in plant signaling pathways and stress adaptive responses by activating protein phosphorylation pathways. However, there have been no comprehensive studies of the SnRK gene family in the widely planted salt-tolerant tree species Casuarina equisetifolia. Here, we comprehensively analyze this gene family in C. equisetifolia using genome-wide identification, characterization, and profiling of expression changes in response to salt stress.


A total of 26 CeqSnRK genes were identified, which were divided into three subfamilies (SnRK1, SnRK2, and SnRK3). The intron–exon structures and protein‑motif compositions were similar within each subgroup but differed among groups. Ka/Ks ratio analysis indicated that the CeqSnRK family has undergone purifying selection, and cis-regulatory element analysis suggested that these genes may be involved in plant development and responses to various environmental stresses. A heat map was generated using quantitative real‑time PCR (RT-qPCR) data from 26 CeqSnRK genes, suggesting that they were expressed in different tissues. We also examined the expression of all CeqSnRK genes under exposure to different salt concentrations using RT-qPCR, finding that most CeqSnRK genes were regulated by different salt treatments. Moreover, co-expression network analysis revealed synergistic effects among CeqSnRK genes.


Several CeqSnRK genes (CeqSnRK3.7, CeqSnRK3.16, CeqSnRK3.17) were up-regulated following salt treatment. Among them, CeqSnRK3.16 expression was significantly up-regulated under various salt treatments, identifying this as a candidate gene salt stress tolerance gene. In addition, CeqSnRK3.16 showed significant expression change correlations with multiple genes under salt stress, indicating that it might exhibit synergistic effects with other genes in response to salt stress. This comprehensive analysis will provide a theoretical reference for CeqSnRK gene functional verification and the role of these genes in salt tolerance.


Salinity, due mainly to sodium chloride (NaCl), seriously affects plant growth and development. The presence of redundant salt ions causes damage by inducing oxidation stress, osmotic stress, and ion toxicity [1, 2]. Casuarina equisetifolia is a sort of salt-resistant tree widely planted in southern China in coastal shelterbelts to stabilize moving sands, provide fuel wood, and reclaim coastal ecosystems due to its superior biological characteristics, such as rapid growth, wind and salt tolerance, and nitrogen fixation [3,4,5]. Therefore, it has become the most important shelterbelt tree species in this coastal area. According to a previous study, C. equisetifolia is highly salt tolerant and can survive in 500 mM NaCl solution [6]. It would therefore be useful to identify the determinants of salt tolerance in this species.

In general, plants respond to adverse environmental pressure, including salt stress, in two methods: regulation of gene expression and modification of proteins [7]. Among them, protein kinase-mediated phosphorylation and dephosphorylation is one way of protein modification [8]. SnRKs (sucrose nonfermenting 1 (SNF1)-related protein kinases) are Ser/Thr protein kinases involved in various physiological activities [9, 10]. The SnRK gene family was segmented into three subfamilies, SnRK1, SnRK2, and SnRK3 [11, 12]. The SnRK1 subfamily shares a strong identity due to an extremely conservative N-terminal catalytic domain [13, 14]. In contrast to the SnRK1 subfamily, the SnRK2 and SnRK3 subfamilies are plant-specific and more diverse [15, 16]. Apart from the identical kinase domain at the N-terminus, each SnRK2 subfamily member contains a C-terminal diverse regulatory domain and an adenosine triphosphate (ATP) binding domain [17]. SnRK3s, also referred to as calcineurin B-like protein-interacting protein kinases (CIPKs), interact with and regulate calcineurin B-like protein (CBL) for transient decoding of calcium signals [18,19,20]. They possess two conserved domains, NAF and PPI, at the C terminus [21].

Previous studies had shown that SnRK genes could make a response to manifold stresses, of which salt stress was an important one [14, 17, 22]. Over-expression of TaSnRK2.9 (Triticum aestivum) in tobacco (Nicotiana tabacum) plants increases scavenging of reactive oxygen species (ROS), thereby promoting the maintenance of ROS at a normal level to protect against abiotic stress [23]. The over-expression of TaSnRK2.3/-2.4/-2.7/-2.8 in Arabidopsis thaliana had been reported to heighten plant resistance to salt and other stresses [24,25,26]. In addition, studies had confirmed that PtSnRK2.5 and PtSnRK2.7 in Populus trichocarpa were closely related to the salt tolerance of transgenic A. thaliana [27]. Similarly, BdSnRK2.9 in Brachypodium distachyon was shown to increase the tolerance of transgenic tobacco to high NaCl concentration treatment [28], and ZmCIPK21 and OsCIPK15 were found could improve plant resistance to salt stress in maize (Zea mays) and rice (Oryza sativa), respectively [29, 30]. Moreover, SnRK3s participates in the SOS (salt overly sensitive) stress signal transduction pathway to regulate intracellular ion homeostasis [31]. For example, the calcium signal generated by salt stress could be intelligently recognized by SOS3 (AtCBL4), and then eliminated superfluous Na+ from root cells by co-phosphorylating SOS1 (Na+/H+ antitransporter) with SOS2 (AtCIPK24) [32]. AtCIPK24 could bind to the photoperiod and circadian clock regulator GI (GIGANTEA), further inhibiting the SOS1 signaling process to regulate the adaptability of plants to salinity [33]. In summary, SnRK genes are important in salt stress responses, and their genetic modification could potentially improve plant salt tolerance. However, a functional understanding of SnRK genes in C. equisetifolia has been lacking, leaving it unclear what role(s) SnRK genes might play in this species.

In our study, a total of 26 SnRK gene family members were identified in the genome of C. equisetifolia. We then evaluated their phylogenetic relationships, physicochemical property, gene structure, conserved domains, and promoter analysis using bioinformatics. Furthermore, differential expression patterns of CeqSnRK genes in different tissues and under various salt treatments were profiled using RT-qPCR. These results lay a substantial groundwork for further investigations of the molecular mechanism in C. equisetifolia resistance to salt stress.


Identification of SnRK genes in C. equisetifolia

As confirmed by SMART, Pfam, and manual screening, twenty-six SnRK proteins in C. equisetifolia were identified. According to the nomenclature in A. thaliana, we named these CeqSnRK1.1 to CeqSnRK1.2, CeqSnRK2.1 to CeqSnRK2.7, and CeqSnRK3.1 to CeqSnRK3.17. These 26 CeqSnRK genes were distributed across 21 scaffolds (the genome of C. equisetifolia is presently assembled to the scaffold level). Their coding sequences ranged from 247 (CeqSnRK2.4) to 541 amino acids (CeqSnRK3.7) in length. Moreover, the molecular weight and pI ranged from 27.71 kDa (CeqSnRK2.4) to 61.48 kDa (CeqSnRK3.7) and from 4.43 (CeqSnRK2.4) to 9.21 (CeqSnRK3.3), respectively. The detailed parameters were presented in Table S1. The results of subcellular localization prediction showed that the SnRK1 subfamily and the majority of SnRK2 subfamily proteins were situated in the cytoplasm and nucleus, but CeqSnRK2.1, CeqSnRK2.2, CeqSnRK2.5, and CeqSnRK2.6 were located in the cytoskeleton. The SnRK3 subfamily members were mostly located in the cytoplasm and chloroplast, while CeqSnRK3.2 was located in the endoplasmic reticulum. The details were provided in Table S2.

Phylogenetic tree and multiple alignment of CeqSnRK genes

To reveal the phylogenetic relation of SnRK proteins from A. thaliana, Eucalyptus grandis, P. trichocarpa, Oryza sativa, and C. equisetifolia, we then built a phylogenetic tree using MEGA 7.0 with the neighbor-joining (NJ) method (Fig. 1). Additionally, detailed information on the protein sequences of the above-mentioned plants can be found in Table S3. The 186 SnRK proteins were distinctly split into three subfamilies: SnRK1, SnRK2, and SnRK3. The phylogenetic analysis indicated that the SnRK1 subfamily had the smallest and the SnRK3 subfamily contained the largest number of genes. In addition, these SnRK genes in each subfamily was even-distributed, except that the rice SnRK genes showed clear clustering, which might be caused by differences between monocotyledons and dicotyledons.

Fig. 1
figure 1

Phylogenetic tree of SnRK genes from Casuarina equisetifolia, Arabidopsis thaliana, Populous trichocarpa, Eucalyptus grandis, and Oryza sativa. A total of 26 CeqSnRK genes, 34 EgrSnRK genes, 39 AtSnRK genes, 45 PtSnRK genes and 48 OsSnRK genes were clustered into three subgroups (SnRK1, SnRK2 and SnRK3). Details of the SnRK genes from all plant species are listed in Table S2. The tree was generated using ClustalX 2.0 with the neighbor-joining (N-J) method

The multiple sequence alignment of the CeqSnRKs was further probed using DNAMAN 8.0 software. Results showed that all of CeqSnRK2 genes encoded an ATP binding site and serine/threonine protein kinase active site at N-terminal regions and a domain known as domain I at the C terminus that is indispensable for enhancing osmotic stress–mediated endurance (Fig. 2). CeqSnRK3 genes not only contained a protein kinase domain at the N terminus, but specifically encoded a NAF domain and PPI domain at the C terminus. In summary, the multiple sequence alignment validated that the CeqSnRK genes have complete functional domains.

Fig. 2
figure 2

A Multiple sequence alignment of CeqSnRK2 genes. B Multiple sequence alignment of CeqSnRK3 genes. Sequences were aligned using DNAMAN 8.0. The scale indicates the number of amino acids. Blue and red indicate conserved amino acids. Red lines and boxes indicate conserved domains

Gene structure, motif composition, and protein structural analyses

The exon–intron structures of the 26 CeqSnRK genes were analyzed to facilitate understanding of gene evolution (Fig. 3). It turned out that the genetic structures of members of the same subfamily shared analogical characteristics. Members of the CeqSnRK1 subfamily had nine introns, consistent with those of E. grandis [34], while the number of introns in the CeqSnRK2 subfamily changed from five to eight. However, there was a noticeable difference in the intron numbers of CeqSnRK3 subfamily members, and they were segmented into two forms: intron-free and intron-rich. Nine CeqSnRK3 genes had no introns, 6 CeqSnRK3 genes had 11 to 15 introns, and the remaining two genes, CeqSnRK3.2 and CeqSnRK3.3, had either 3 or 1 intron (Fig. 3).

Fig. 3
figure 3

Phylogenetic relationships and gene structure of SnRK genes in Casuarina equisetifolia. Left panel shows the phylogenetic tree of 26 CeqSnRK genes built using the neighbor-joining method. Red, green, and blue indicate the SnRK1, SnRK2, and SnRK3 subgroups, respectively. Right panel shows the gene structure of CeqSnRK genes. Exons are indicated by yellow rectangles. Gray lines connecting two exons represent introns

To further reveal the structure of CeqSnRK proteins, we identified 20 conserved motifs and assessed their distribution (Fig. 4); the details were shown in Table S4. According to the Pfam annotation results, motifs 1, 2, and 3 encode a protein kinase domain; motifs 10 and 11 encode a NAF domain; and motif 20 encodes a KA1 domain. The other motifs had no functional annotations based on the Pfam database. We then used WebLogo and found that four motifs (motifs 1, 4, 7, and 14) were present in all of the CeqSnRK members (Fig. 5A); motifs 16, 18, and 20 were present in all members of the CeqSnRK1 subfamily (Fig. 5B); and motifs 17 and 10 were unique to the CeqSnRK2 and CeqSnRK3 subfamilies, respectively (Fig. 5C, D). In conclusion, the motifs in CeqSnRK genes of the same subfamily contained similar regularity, indicating that they might share the same genetic structure and functional domains.

Fig. 4
figure 4

Conserved motifs of SnRK genes in Casuarina equisetifolia. The distribution of 20 conserved motifs in CeqSnRK genes was analyzed using MEME. Rectangles with different colors represent different motifs

Fig. 5
figure 5

Sequence logos of repeated domains in the CeqSnRK gene family. A Motifs shared by SnRK1, SnRK2, and SnRK3 subfamilies. B Motifs unique to the SnRK1 subfamily. C Motif unique to the SnRK2 subfamily. D Motif unique to the SnRK3 subfamily. The overall height of each stack shows the conservation of the SnRK protein sequence at that position. Letters indicate different amino acid residue

The functions of a protein are closely related with its structure; therefore, the secondary and tertiary protein structures of each CeqSnRK gene were analyzed. We predicted the secondary structures of CeqSnRK genes using the Phyre 2 software, finding that alpha helix and random coil accounted for a major proportion (Table S5). Furthermore, 3D models of CeqSnRK genes were constructed using Swiss-Model online server (Fig. S1). As illustrated in the Fig. S1, the 3D structures of CeqSnRKs were variable among different subfamilies, indicating the existence of potential functional diversity.

Homology analysis in C. equisetifolia

We identified 43 SnRK orthologues and 8 paralogues on the basis of the topology of the phylogenetic tree and BLASTN results. To further explore the effect of selection pressure on the evolution of CeqSnRK genes, the synonymous substitutions (Ks), non-synonymous substitutions (Ka), and the Ka/Ks ratios of paralogues and orthologues were calculated using DnaSP 5.0 software. We built a sliding-window analysis for paralogous genes (Fig. 6); the Ka/Ks values of the orthologues are shown in Fig. S2. The ratio of Ka/Ks can be used to judge the selective pressure. Generally, when the ratio is equal to 1, it means neutral selection, intimating that the DNA mutation exerts no influence on the organism. Ka/Ks < 1 indicates purification selection (negative selection), while Ka/Ks > 1 indicates accelerated evolution (positive selection) [27]. We discovered that the Ka/Ks ratios of the paralogue pairs in C. equisetifolia ranged from 0.037 to 0.751 (Table S6), indicating that the CeqSnRK gene family has undergone purifying selection.

Fig. 6
figure 6

Sliding-window analysis of eight paralogous gene pairs. Window size, 150 bp; step size, 9 bp

Promoter analyses

In order to determine the cis-acting elements of the CeqSnRK genes, the promoter region of the CeqSnRK gene (the genomic DNA sequence 2 kb upstream of the translation start point) was submitted to the PlantCARE database for search, which will be conducive to further understanding of gene function and regulation (Table S7). The results identified three classes of cis-elements associated with stress responses, hormone responses, and plant development. In Fig. 7, a total of 98 ABA-responsive elements (ABRE) were found in the promoters of 24 CeqSnRK genes, of which CeqSnRK2 genes all contained ABRE elements, indicating that most of the CeqSnRK genes are involved in ABA signaling transduction pathway. The SA-responsive element (TCA-element) was found in 17 CeqSnRK genes, which was also a common cis-acting element.

Fig. 7
figure 7

Analysis of cis-acting elements in the promoter regions of CeqSnRK genes. Binding sites in the promoter region are represented by boxes of different colors, and the graph shows the number of binding sites

In addition to CeqSnRK3.6/-3.11/-3.16, cis-elements responding to methyl jasmonate (MeJA) were found in the promoters of the remaining 23 CeqSnRK genes. Moreover, elements responding to auxin (AuxRR-core and TGA element) and the gibberellin-responsive element (TATC-box, GARE-motif, and P-box) were also detected in the promoter sequences, with GARE-motif only present in a few CeqSnRK genes. Among the cis-acting elements associated with stress, drought (MBS), low temperature (LTR), anaerobic (ARE), and stress response (TC-rich Repeats) elements were found in the promoters of 16, 17, 22, and 9 CeqSnRK genes, respectively.

Furthermore, 54 elements partook in plant growth and development were also discovered, among them, the element related to seed-specific regulation (RY-element) was only present in CeqSnRK3.7 and CeqSnRK2.1. The differentiation of the palisade mesophyll cells (HD-Zip1), the meristem expression (CAT-box), and endosperm expression (GCN4-motif) were found in the CeqSnRK3.16 promoter. In this study, 278 hormone-responsive elements, 124 stress-related response elements, and 54 growth-related elements were discovered in the CeqSnRK gene promoters, indicating that the CeqSnRK genes may respond to various stresses, as well as be involved in important developmental processes.

Expression profile of SnRK genes in different tissues

A heat map was generated using RT-qPCR data from 24 CeqSnRK genes, finding that these patterns differed considerably between tissues (Fig. 8). Additionally, we did not detect the expression of CeqSnRK2.3 or CeqSnRK3.2. Most of the genes had low expression levels in roots, except for CeqSnRK2.5, CeqSnRK2.6, and CeqSnRK3.11 (Fig. 8). The expression levels of CeqSnRK3.5 and CeqSnRK3.11 in xylem, phloem, stem, inflorescence, spray, and mature branches were low, indicating an opposite expression pattern to that in root. Moreover, 7 of 24 CeqSnRK genes were highly expressed in xylem, especially CeqSnRK3.5, CeqSnRK3.10, and CeqSnRK3.15. Notably, the genes clustered together in the heat map presented similar expression patterns. For example, CeqSnRK2.2, CeqSnRK3.3, CeqSnRK2.4, and CeqSnRK3.1 were clustered together and showed lower expression levels in roots, shoots, xylem, and phloem, but higher in inflorescences. Similar patterns were also observed for CeqSnRK1.1, CeqSnRK1.2, CeqSnRK2.1, and CeqSnRK2.6.

Fig. 8
figure 8

Expression analysis of CeqSnRK genes in various tissues. The heatmap shows the hierarchical clustering of 24 CeqSnRK genes in different tissues. The color scale on the right-hand side represents log10 expression values; white represents low expression, and red indicates a high expression

CeqSnRK gene expression under salt treatment

It was obvious that the expression level of a few SnRK genes was up-regulated or down-regulated under diverse NaCl treatments in roots of C. equisetifolia (Fig. S3). In particular, the expression levels of 11/24 CeqSnRKs were increased following NaCl treatment (Fig. 9); for example, the expression of CeqSnRK3.4 and CeqSnRK3.6 peaked at 300 mM and 400 mM salt treatment, respectively (Fig. 9A, B). We obtained similar results for the paralogous CeqSnRK3.13/-3.16, which were markedly up-regulated in roots at high salinity treatment (300/400 mM) (Fig. S3). However, the expression levels of the paralogous CeqSnRK1.1/-1.2 were strongly down-regulated under high salt treatment (300 and 400 mM).

Fig. 9
figure 9

A The expression levels of 5 CeqSnRK genes were up-regulated and peaked under 300 mM NaCl. B The expression levels of 5 CeqSnRK genes were up-regulated and peaked under 400 mM NaCl. The Y-axis and X-axis indicate relative expression levels and salt concentration of stress treatment, respectively. Mean values and standard deviations (SDs) were obtained from three biological and three technical replicates. The error bars indicate standard deviation. **P < 0.01 and *P < 0.05

Subsequently, the CeqSnRK gene expression profiles in shoots after treatment with different concentrations of NaCl were analyzed (Fig. S4). Of the 24 CeqSnRK genes, compared with their expression in the control (non-treated seedlings), 11 were up-regulated under various salt concentrations, 12 were down-regulated, and only CeqSnRK1.2 did not show any significant change. Of note, the expression levels of seven CeqSnRKs were found to be significantly enhanced and reached a maximum at 100 mM after NaCl treatment, followed by a sharp downregulation (Fig. 10A). For example, the expression levels of CeqSnRK3.8 (more than 11 times) and CeqSnRK3.10 (about fourfold) were rapidly increased under 100 mM treatment compared with that of the control (Fig. 10A). Moreover, CeqSnRK1.1, CeqSnRK2.4, and CeqSnRK3.5 were greatly up-regulated under high concentrations of NaCl (Fig. S4). Interestingly, the expression pattern of CeqSnRK2.4 was distinctive from that of other SnRK genes, mainly showing a trend of first decreasing and then increasing. In other words, the expression of CeqSnRK2.4 was slightly depressed at 100 mM and then was increased about fivefold under 200 mM. With the further increase in salt concentration, the expression of CeqSnRK2.4 remained at a high level.

Fig. 10
figure 10

Relative expression of several CeqSnRK genes in shoots following NaCl treatment as determined by RT-qPCR. A The expression levels of seven CeqSnRK genes were up-regulated and peaked under 100 mM NaCl. B The expression levels of 11 CeqSnRK genes were up-regulated and peaked at 1 h following 200 mM NaCl treatment. C The expression levels of 11 CeqSnRK genes were up-regulated at different time points following 200 mM NaCl treatment. The Y-axis and X-axis indicate relative expression levels and the time courses of stress treatments, respectively. Mean values and standard deviations (SDs) were obtained from three biological and three technical replicates. The error bars indicate standard deviation. **P < 0.01 and *P < 0.05

We further analyzed the expression profile of CeqSnRK genes after different periods of salt treatments and found 17 genes to be distinctly up-regulated at several time points, while 8 genes were down-regulated or not changed at any time point, compared to the control values (Fig. S5). For example, the expression of 11 genes intensely increased and peaked at 1 h after NaCl treatment and then distinctly decreased at subsequent moments (Fig. 10B). Notably, the expression levels of CeqSnRK1.1, CeqSnRK2.1, CeqSnRK3.6, CeqSnRK3.7, CeqSnRK3.16, and CeqSnRK3.17 were higher at each time point under NaCl treatment than in the control (Fig. 10C). Specifically, CeqSnRK3.16 was drastically up-regulated (sevenfold) at 1 h, and the up-regulation level was maintained at subsequent time points. In contrast, the expression of four CeqSnRK genes was significantly descended at each time period. In addition, the effect of diurnal cycle on the experiment was not considered for the time being in the time course analysis of this paper.

Correlations and co-regulatory networks of CeqSnRK genes

To explore the relationship between the responses of these CeqSnRK genes to salt stress, the correlation and co-regulatory networks were established in accordance with the Pearson’s correlation coefficient data for relative expression levels (Fig. 11). For resisting salt pressure, the expression changes of most genes (60 and 57% of genes in roots and stems, respectively) were positively correlated. However, some showed a negative correlation trend. For example, there was a negative correlation between CeqSnRK3.7 and CeqSnRK3.17/-3.14/-3.13 expression in both roots and shoots (Fig. 11A and B). Moreover, several gene pairs exhibited positive correlations in the root system but negative correlations in shoots. For example, CeqSnRK2.7 expression was positively correlated with that of CeqSnRK2.1 and CeqSnRK2.2 in the root system, but negatively correlated in the shoots.

Fig. 11
figure 11

Correlations and co-regulatory networks of CeqSnRK genes under stress treatments. A, B Correlation analysis of CeqSnRK genes under NaCl treatment in roots and shoots, respectively. Each correlation is shown by a smooth curve. Red and gray indicate positive and negative correlations, respectively. C, D The coregulatory network of CeqSnRK genes under NaCl treatment in roots and shoots, respectively. Pearson’s correlation coefficients of co-regulatory gene pairs were considered significant at the 0.05 significance level (P-value), and the different correlation levels of the gene pairs are marked by edge lines with different colors

The correlation coefficient of more than 0.5 in the shoots was significantly higher than that in the roots (Fig. S6), according to the co-regulatory networks associated with salt stresses (Fig. 11C and D). Meanwhile, most of the positive correlations occurred between members belonging to the same subfamily, especially the SnRK3 subfamily. For example, CeqSnRK3.16, CeqSnRK3.15, and CeqSnRK3.17 showed significant positive correlations (P-value ≤ 0.01 and 0.7 < PCC) under NaCl stress. In addition, CeqSnRK3.16 expression was positively correlated with that of most CeqSnRK genes after salt treatments, suggesting it might respond to salt stress by interacting with other genes. This suggests that gene duplication not only leads to functional divergence but also enhances the synergistic interactions of homologues to help plants adapt to their complex habitats.


According to previous reports, 34, 48, 39, 44, and 60 SnRK genes have been identified from E. grandis [34], rice [35], A. thaliana [12], B. distachyon [28], and Hedychium coronarium [36], respectively. Various roles of SnRKs in plant growth, development, and resistance to biotic and abiotic stresses have been reported, but less is known about the functions of these proteins in C. equisetifolia. Hence, we investigated systematically CeqSnRK genes by combining bioinformatic analysis and RT-qPCR experiment.

The 26 CeqSnRK genes were identified and separated into three subfamilies based on a comprehensive phylogenetic tree. There were 2, 7, and 17 CeqSnRK genes in the CeqSnRK1, CeqSnRK2, and CeqSnRK3 subfamilies, respectively. Moreover, different SnRK gene subfamilies contained various conserved domains, but the N-terminal protein kinase domain was retained throughout the family. We also discovered that there was a NAF domain in the CeqSnRK3 subfamily at the C terminus (Fig. 2). Based on a previous study, SnRK3 genes can cooperate with CBLs in a calcium-dependent pattern because of the existence of NAF domain, implying that CeqSnRK3s might respond to stresses by interacting with CBLs [18, 19, 37, 38]. Furthermore, most CeqSnRK genes were clustered with the PtSnRK rather than the OsSnRK genes, implying that CeqSnRK and PtSnRK genes have a relatively close evolutionary relationship and that there are significant evolutionary differences between the SnRK genes in dicots and monocots.

Genetic structural diversity is the main source of multigene family evolution [29, 30]. Therefore, the conserved motifs and gene structure of 26 CeqSnRK genes were analyzed. The number and distribution of exon/intron structures differed, but most genes within the same subfamily shared similar gene structures. However, the CeqSnRK3 subfamily could be subdivided, based on the number of introns, into two kinds, the intron-rich and intron-deficient clades, which is similar to the differentiation of the HcSnRK3 genes of H. coronarium [36]. During evolution, intron-rich populations were inclined to lose introns, thus becoming intron-deficient [39]. Most (9/17) CeqSnRK3 genes had no intron, suggesting that this subfamily was prone to loss of introns and that the whole family was likely to be more conservative with evolution. We also observed that the most conserved motifs among members of the same subfamily shared certain similar characteristics (Fig. 4), which might give evidence of a closer evolutionary relation within subfamilies. Moreover, motifs 4, 1, 7, and 14 appeared in every CeqSnRK gene (Fig. 5A), implying their significance of the function in CeqSnRKs. The similarities in structure and motif composition among SnRK genes were coincident with our phylogenetic analysis. Meanwhile, the differences among subfamilies suggest that the functions of SnRK members are diverse.

Accumulating evidence suggests that gene activity is often relevant with differences in promoter regions, as cis-elements play an important role in controlling gene expression during development and environmental changes [32, 33]. In soybean (Glycine max), GmSnRK2.16 and GmSnRK2.18 can respond to salt stress due to the presence of two MeJA-responsive elements in their promoters, while GmSnRK2.6, which lacks such an element in its promoter region, cannot respond to salt stress [40]. Consistent with this, the expression level in roots of CeqSnRK3.11, which lacks a MeJA-responsive element in its promoter (Fig. 7), was not significantly changed under different concentrations of NaCl. Moreover, 24 out of 26 CeqSnRK genes contained ABRE cis-regulatory elements related to ABA responsiveness. Accumulating evidence suggests that ABRE-binding protein/ABRE-binding factor (AREB/ABF) can positively regulate plant responses and tolerance, and the SnRK-ABF pathway plays a crucial role in abiotic stress resistance [41,42,43]. In wheat, TaSnRK2.9 can interact with NtABF2 (N. tabacum) and up-regulate the expression of NtABF2 under mannitol (that is, osmotic stress) or NaCl treatment [23]. These results indicate that SnRK genes could heighten the resistance to salt though mediating plant hormone signaling.

Gene expression profiles can provide important clues to reveal gene functions [28, 44]. We analyzed the expression of all 26 CeqSnRK genes under diverse NaCl treatments. The expression of 11 genes was increased in both roots and shoots with the treatment of NaCl, and among these, CeqSnRK2.7, CeqSnRK3.8, and CeqSnRK3.16 were significantly up-regulated. In contrast, CeqSnRK3.9 was down-regulated in both roots and shoots under different concentrations of NaCl. Analogous expression patterns existed in some paralogous pairs in the same tissue. For example, under high NaCl concentrations (300 and 400 mM), the expression levels of both CeqSnRK3.13 and CeqSnRK3.16 were not only significantly up-regulated in roots but similar to control expression levels in shoots (Figs. 9 and 10). Moreover, the expression levels of some CeqSnRK genes, such as CeqSnRK3.5, CeqSnRK3.11, and CeqSnRK3.12, were increased in shoots but decreased in roots under salt stress. Similar results were observed in E. grandis [34], suggesting that some SnRK genes are specifically expressed in shoots. CeqSnRK3.12, EgrSnRK3.9, and AtCIPK21—the latter known to regulate osmotic and salt stress—clustered together [45]. We found that the expression of CeqSnRK3.12 was obviously increased in shoots under 100 mM NaCl treatment. Thus, CeqSnRK3.12 may respond to salt stress. Furthermore, most CeqSnRK genes were up-regulated at different time points under saline condition in shoots (Fig. S4). Specifically, the expression levels of CeqSnRK3.6, CeqSnRK3.7, CeqSnRK3.16, and CeqSnRK3.17 were significantly up-regulated at each time period under NaCl treatment compared with the non-treated control (Fig. 10). Previous studies have shown that AtCIPK24 transports excess intracellular Na+ to the extracellular compartment under salt stress [46], thereby increasing the salt tolerance of plants. CeqSnRK3.13 and CeqSnRK3.16 were clustered in the same group with AtCIPK24, indicating that CeqSnRK3.13 and CeqSnRK3.16 might respond to salt stress. Expression of AtCIPK3 increases the tolerance of A. thaliana to high salt concentrations, drought and other stress stimuli [47]. In this study, CeqSnRK3.7, which is a homolog of AtCIPK3, was strongly up-regulated at all time points under salt treatment, indicating that CeqSnRK3.7 could play important roles in responses to salt stress. Taken together, our results indicate that most CeqSnRK genes were up-regulated to a degree after treating by NaCl, suggesting that CeqSnRK genes might be conducive to improving salt tolerance in C. equisetifolia.


SnRK genes are involved in a variety of signaling pathways, including responses to biotic and abiotic stresses. In the present study, we identified 26 CeqSnRK genes and divided them into three subfamilies on the basis of motif composition and gene structural similarity. Promoter analysis revealed CeqSnRK genes involved in plant development and responses to stress and hormones. The expression levels of most SnRK genes were induced in shoots and roots under disparate salt treatments. Notably, CeqSnRK3.16 expression was up-regulated under salt treatment, suggesting that it might respond to salt stress. Moreover, CeqSnRK3.16 showed significant expression change correlations with multiple genes under salt stress, indicating that it might have synergistic effects with other genes in response to salt stress. Taken together, these results provided a rich theoretical basis for further validation of the functions of CeqSnRKs and their roles in salt tolerance.

Materials and methods

Identification of SnRK proteins in the C. equisetifolia

Download the whole genome protein sequence of Casuarina from the Casuarina database ( as a local database [36]. The SnRK protein sequences in A. thaliana downloaded from the Phytozome database ( was used as the target sequence, and then BLAST (E-value-0.5) homology alignment was performed with the local data to obtain the candidate SnRK protein sequences in C. equisetifolia [48]. Finally, all candidate SnRK genes were manually filtered based on the information from the Pfam database (, the NCBI Conserved Domain database (, and the SMART database ( [49]. Additionally, detailed resources on SnRK genomes in E. grandis, P. trichocarpa, and rice were obtained from a previous study [34]. Furthermore, the basic physicochemical property parameters of each SnRK gene, including information such as open reading frame (ORF) length, molecular weight (MW) and isoelectric point (pI) were explored using the online website Expasy ( The online website WoLP PSORT ( was used to predict the subcellular localization of the CeqSnRK genes.

Multiple sequence alignment and phylogenetic analyses

All SnRK protein sequences from A. thaliana, C. equisetifolia, E. grandis, P. trichocarpa, and rice were aligned using ClustalX 2.11, and a phylogenetic tree was constructed using the neighbor-joining (NJ) method [49, 50] in MEGA 7.0 with a bootstrap value of 1000 [51]. DNAMAN was used to show multiple sequence alignment of CeqSnRK genes [36].

Identification of conserved motifs and analyses of gene and protein structure

The gene structures of exon–intron distribution of CeqSnRKs were predicted using the Gene Structure Display Server (GSDS: based on the coding sequence of each SnRK gene and its corresponding genomic DNA sequence. To analyze the conserved motifs of CeqSnRK proteins, Multiple Expectation Maximization for Motif Elicitation (MEME) ( was used with the following parameters: number of repetition = any; maximum number of motifs = 20; and optimum motif length = 6–200 residues [14, 28]. The 3D structure of each CeqSnRK was determined using SWISS-MODEL ( [52].

Ka and Ks analyses of homologous pairs

The paralogs and orthologs were identified on the basis of the way described by Wang et al. [35]. That is, the coding sequences of SnRK genes of five species were aligned pairwise by local BLASTN, and the candidate gene pairs of paralogous genes of C. equisetifolia and orthologous pairs of C. equisetifolia and other four species were initially screened. Afterward, the definition in the previous study was used to further confirm. The synonymous (Ks) and non-synonymous (Ka) substitutions each locus between duplicated genes pairs were calculated using DnaSP v5.0 [53]. To further analyze Ka/Ks values, GraphPad Prism 5 was used to output sliding window graphs for analysis.

Cis-elements in the promoter regions of CeqSnRK genes

The promoter region, 2000 bp upstream of the translation start of each CeqSnRK gene, was searched using all Casuarina whole genome sequences and GFF files. The subsequent prediction and analysis of cis-acting elements were completed using the online site PlantCARE ( [54].

Plant materials and stress treatments

Clone A8 of C. equisetifolia is widely used for coastal shelterbelt construction in Guangdong province. Different tissues (root、stem、phloem、xylem、inflorescence、spray and mature branch) were taken from 10-year-old clone A8 in our nursery. Then, young branchlets of A8 were collected and vegetatively propagated to obtain experiment material for this study. After hydroponically rooting, clone A8 was transferred to a greenhouse and cultured for 3 months until the seedlings reached 60–70 cm in height. Subsequently, two salt treatments were included, one with 200 mM NaCl for 0, 1, 6, 24 and 168 h and the other with different salt concentrations (0, 100, 200, 300, and 400 mM) for 24 h. Plant samples (roots and shoots) were taken after watering with the various NaCl solution. All collected samples were immediately frozen with liquid nitrogen and transferred to − 80℃ for storage until RNA extraction.

RNA extraction and RT-qPCR analysis

Total RNA was isolated from each sample following the protocol as described previously and examined using 1% agarose gel electrophoresis and a NanoDrop™ One/OneC (Thermo Fisher Scientific, USA) [34]. According to the instructions, the first-strand cDNA synthesis and RT-qPCR were completed using PrimerScript RT MasterMix (Takara, Tokyo, Japan) and Swiss-made LightCycler480 II Real-Time PCR system, respectively. The loading system was: 10 μL Premix Ex Taq II, 2 μL cDNA template, 0.4 μL ROX Reference Dye (50X), 6 μL sterile water, and 0.8 μL forward and reverse primers each. qPCR reaction conditions were as follows: 95 ℃ for 30 s, followed by 40 cycles of denaturation at 95 ℃ for 5 s, and annealing at 55–60 ℃ for 34 s [54]. Based on literature reports, using EF1α as the reference gene, primers were designed using Primer 6 software [55] and then sent to the company (Ruibo, Guangzhou) for synthesis. The primer sequence was detailed in Table S8. Three biological and technical replicates were performed for each real-time PCR reaction [35]. The experimental data were processed by the 2−ΔΔct method [56], and the significant difference in the data was calculated by using IBM SPSS Statistics 25 software. Statistical analysis and graphing of gene expression patterns were completed using GraphPad 8 [57]. The obtained data were log2-transformed and visualized as a heatmap using TBtools [58].

Pearson’s correlation analyses

Pearson’s correlation coefficients of CeqSnRK gene expression levels were calculated and visualized using R ( based on the RT-qPCR results [59]. Pairs of genes that were positively correlated with each other were collected for a gene coregulatory network analysis. The co-expression networks were graphically visualized using Cytoscape based on the Pearson’s correlation coefficients of these gene pairs [9].

Availability of data and materials

The genome sequences of E. grandis were downloaded from Phytozome database (;

the genome sequences of P. trichocarpa were downloaded from Phytozome database (; and the genome sequences of A. thaliana were downloaded from Phytozome database ( Rice protein sequences and cDNA sequences was provided by the Rice Annotation Project (RAP) ( The Casuarina genome data were downloaded from online website ( The datasets supporting the results of this article are included in the article and Additional files.



Number of non-synonymous substitutions per non-synonymous site


Number of synonymous substitutions per synonymous site


Quantitative real‑time PCR


Sucrose nonfermenting 1


Calcineurin B-like protein-interacting protein kinases


Calcineurin B-like protein


Salt overly sensitive




  1. Hu HR, Liu H, Liu FH. Seed germination of hemp (Cannabis sativa L.) cultivars responds differently to the stress of salt type and concentration. Industrial Crops and Products. 2018;123:254–61.

    Article  CAS  Google Scholar 

  2. Liu JG, Han X, Yang T, Cui WH, Wu AM, Fu CX, et al. Genome-wide transcriptional adaptation to salt stress in Populus. BMC Plant Biol. 2019;19(1):14.

    Article  Google Scholar 

  3. Goel VL, Behl HM. Growth and productivity assessment of Casuarina glauca Sieb. ex. Spreng on sodic soil sites. Bioresour Technol. 2005;96(12):1399–404.

    Article  CAS  Google Scholar 

  4. Pinyopusarerk K, Williams ER. Variations in growth and morphological characteristics of Casuarina junghuhniana provenances grown in Thailand. J Trop Forest Sci. 2005;17(4):574–87.

    Google Scholar 

  5. Hu P, Zhong CL, Zhang Y, Jiang QB, Chen Y, Chen Z, et al. Geographic variation in seedling morphology of Casuarina equisetifolia subsp equisetifolia (Casuarinaceae). Aust J Bot. 2016;64(2):160–70.

    Article  Google Scholar 

  6. Tani C, Sasakawa H. Proline accumulates in Casuarina equisetifolia seedlings under salt stress. Soil Sci Plant Nutr. 2006;52(1):21–5.

    Article  CAS  Google Scholar 

  7. Bohnert HJ, Gong QQ, Li PH, Ma SS. Unraveling abiotic stress tolerance mechanisms - getting genomics going. Curr Opin Plant Biol. 2006;9(2):180–8.

    Article  CAS  Google Scholar 

  8. Hunter T. Protein kinases and phosphatases: the yin and yang of protein phosphorylation and signaling. Cell. 1995;80(2):225–36.

    Article  CAS  Google Scholar 

  9. Wu P, Wang WL, Duan WK, Li Y, Hou XL. Comprehensive analysis of the CDPK-SnRK superfamily genes in Chinese cabbage and its evolutionary implications in plants. Frontiers in Plant Sci. 2017;8:19.

    CAS  Google Scholar 

  10. Zhang Y, Wang S, Chen J, Wang W, Zhou T, Xiao B, et al. Cloning and expression analysis of CsSnRK2.1 and CsSnRK2.2 genes in tea plant (Camellia sinensis) under abiotic stress. Journal of Tea Sci. 2018;38(2):183–92.

    Google Scholar 

  11. Wrzaczek M, Hirt H. Plant MAP kinase pathways: how many and what for? Biol Cell. 2001;93(1–2):81–7.

    Article  CAS  Google Scholar 

  12. Hrabak EM, Chan CWM, Gribskov M, Harper JF, Choi JH, Halford N, et al. The Arabidopsis CDPK-SnRK superfamily of protein kinases. Plant Physiol. 2003;132(2):666–80.

    Article  CAS  Google Scholar 

  13. Monger WA, Thomas TH, Purcell PC, Halford NG. Identification of a sucrose nonfermenting-1-related protein kinase in sugar beet (Beta vulgaris L). Plant Growth Regulation. 1997;22(3):181–8.

    Article  CAS  Google Scholar 

  14. Wang CL, Zhang ZP, Miao MM. SNF1-related protein kinase (SnRK) 1 involved in the regulation of raffinose family oligosaccharide metabolism in cucumber (Cucumis sativus L.) Calli. J Plant Growth Regul. 2016;35(3):851–64.

    Article  CAS  Google Scholar 

  15. Halford NG, Hardie DG. SNF1-related protein kinases: global regulators of carbon metabolism in plants? Plant Mol Biol. 1998;37(5):735–48.

    Article  CAS  Google Scholar 

  16. Halford NG, Hey SJ. Snf1-related protein kinases (SnRKs) act within an intricate network that links metabolic and stress signalling in plants. Biochem J. 2009;419:247–59.

    Article  CAS  Google Scholar 

  17. Dong XF, Cui N, Wang L, Zhao XC, Qu B, Li TL, et al. The SnRK protein kinase family and the function of SnRK1 protein kinase. Int J Agric Biol. 2012;14(4):575–9.

    CAS  Google Scholar 

  18. Albrecht V, Ritz O, Linder S, Harter K, Kudla J. The NAF domain defines a novel protein-protein interaction module conserved in Ca2+-regulated kinases. EMBO J. 2001;20(5):1051–63.

    Article  CAS  Google Scholar 

  19. Ma X, Li QH, Yu YN, Qiao YM, Haq S, Gong ZH. The CBL-CIPK pathway in plant response to stress signals. Int J Mol Sci. 2020;21(16):27.

    Article  Google Scholar 

  20. Kim KN, Lee JS, Han H, Choi SA, Go SJ, Yoon IS. Isolation and characterization of a novel rice Ca2+-regulated protein kinase gene involved in responses to diverse signals including cold, light, cytokinins, sugars and salts. Plant Mol Biol. 2003;52(6):1191–202.

    Article  CAS  Google Scholar 

  21. Ohta M, Guo Y, Halfter U, Zhu JK. A novel domain in the protein kinase SOS2 mediates interaction with the protein phosphatase 2C AB12. Proc Natl Acad Sci USA. 2003;100(20):11771–6.

    Article  CAS  Google Scholar 

  22. Ma QJ, Sun MH, Kang H, Lu J, You CX, Hao YJ. A CIPK protein kinase targets sucrose transporter MdSUT2.2 at Ser(254) for phosphorylation to enhance salt tolerance. Plant Cell and Environment. 2019;42(3):918–30.

    Article  CAS  Google Scholar 

  23. Feng JL, Wang LZ, Wu YN, Luo QC, Zhang Y, Qiu D, et al. TaSnRK2.9, a sucrose non-fermenting 1-related protein kinase gene, positively regulates plant response to drought and salt stress in transgenic tobacco. Front Plant Sci. 2019;9:17.

    Article  Google Scholar 

  24. Zhang HY, Mao XG, Wang CS, Jing RL. Overexpression of a common wheat gene TaSnRK2.8 enhances tolerance to drought, salt and low temperature in Arabidopsis. Plos One. 2010;5(12):12.

    Article  Google Scholar 

  25. Mao XG, Zhang HY, Tian SJ, Chang XP, Jing RL. TaSnRK2.4, an SNF1-type serine/threonine protein kinase of wheat (Triticum aestivum L.), confers enhanced multistress tolerance in Arabidopsis. J Exp Bot. 2010;61(3):683–96.

    Article  CAS  Google Scholar 

  26. Zhang HY, Mao XG, Jing RL, Chang XP, Xie HM. Characterization of a common wheat (Triticum aestivum L.) TaSnRK2.7 gene involved in abiotic stress responses. J Exp Bot. 2011;62(3):975–88.

    Article  CAS  Google Scholar 

  27. Song XQ, Yu X, Hori C, Demura T, Ohtani M, Zhuge Q. Heterologous overexpression of poplar SnRK2 genes enhanced salt stress tolerance in Arabidopsis thaliana. Frontiers in Plant Sci. 2016;7:11.

    Article  Google Scholar 

  28. Wang LZ, Hu W, Sun JT, Liang XY, Yang XY, Wei SY, et al. Genome-wide analysis of SnRK gene family in Brachypodium distachyon and functional characterization of BdSnRK2.9. Plant Sci. 2015;237:33–45.

    Article  CAS  Google Scholar 

  29. Huai JL, Wang M, He JG, Zheng J, Dong ZG, Lv HK, et al. Cloning and characterization of the SnRK2 gene family from Zea mays. Plant Cell Rep. 2008;27(12):1861–8.

    Article  CAS  Google Scholar 

  30. Xiang Y, Huang YM, Xiong LZ. Characterization of stress-responsive CIPK genes in rice for stress tolerance improvement. Plant Physiol. 2007;144(3):1416–28.

    Article  CAS  Google Scholar 

  31. Guo Y, Xiong L, Song C-P, Gong D, Halfter U, Zhu J-K. A calcium sensor and its interacting protein kinase are global regulators of abscisic acid signaling in Arabidopsis. Dev Cell. 2002;3(2):233–44.

    Article  CAS  Google Scholar 

  32. Liu J, Ishitani M, Halfter U, Kim CS, Zhu JK. The Arabidopsis thaliana SOS2 gene encodes a protein kinase that is required for salt tolerance. Proc Natl Acad Sci USA. 2000;97(7):3730–4.

    Article  CAS  Google Scholar 

  33. Park HJ, Kim W-Y, Yun D-J. A role for GIGANTEA: keeping the balance between flowering and salinity stress tolerance. Plant Signal Behav. 2013;8(7): e24820.

    Article  Google Scholar 

  34. Wang YJ, Yan HF, Qiu ZF, Hu B, Zeng BS, Zhong CL, Fan CJ. Comprehensive analysis of SnRK gene family and their responses to salt stress in Eucalyptus grandis. Int J Mol Sci. 2019;20(11):21.

    Article  CAS  Google Scholar 

  35. Kobayashi Y, Yamamoto S, Minami H, Kagaya Y, Hattori T. Differential activation of the rice sucrose nonfermenting1-related protein kinase2 family by hyperosmotic stress and abscisic acid. Plant Cell. 2004;16(5):1163–77.

    Article  CAS  Google Scholar 

  36. Wang CT, Abbas F, Zhou YW, Ke YG, Li XY, Yue YC, et al. Genome-wide identification and expression pattern of SnRK gene family under several hormone treatments and its role in floral scent emission in Hedychium coronarium. PeerJ. 2021;9:27.

    Google Scholar 

  37. Drerup MM, Schlucking K, Hashimoto K, Manishankar P, Steinhorst L, Kuchitsu K, et al. The calcineurin B-like calcium sensors CBL1 and CBL9 together with their interacting protein kinase CIPK26 regulate the Arabidopsis NADPH oxidase RBOHF. Mol Plant. 2013;6(2):559–69.

    Article  CAS  Google Scholar 

  38. Tang RJ, Wang C, Li KL, Luan S. The CBL-CIPK calcium signaling network: unified paradigm from 20 years of discoveries. Trends Plant Sci. 2020;25(6):604–17.

    Article  CAS  Google Scholar 

  39. Colina F, Amaral J, Carbo M, Pinto G, Soares A, Canal MJ, et al. Genome-wide identification and characterization of CKIN/SnRK gene family in Chlamydomonas reinhardtii. Sci Rep. 2019;9:16.

    Article  Google Scholar 

  40. Zhao W, Cheng YH, Zhang C, Shen XJ, You QB, Guo W, et al. Genome-wide identification and characterization of the GmSnRK2 family in soybean. Int J Mol Sci. 2017;18(9):21.

    Article  Google Scholar 

  41. Fujita Y, Yoshida T, Yamaguchi-Shinozaki K. Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants. Physiol Plant. 2013;147(1):15–27.

    Article  CAS  Google Scholar 

  42. Furihata T, Maruyama K, Fujita Y, Umezawa T, Yoshida R, Shinozaki K, et al. Abscisic acid-dependent multisite phosphorylation regulates the activity of a transcription activator AREB1. Proc Natl Acad Sci USA. 2006;103(6):1988–93.

    Article  CAS  Google Scholar 

  43. Yoshida T, Fujita Y, Sayama H, Kidokoro S, Maruyama K, Mizoi J, et al. AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation. Plant J. 2010;61(4):672–85.

    Article  CAS  Google Scholar 

  44. Yu X, Takebayashi A, Demura T, Ohtani M. Differential expression of poplar sucrose nonfermenting1-related protein kinase 2 genes in response to abiotic stress and abscisic acid. J Plant Res. 2017;130(5):929–40.

    Article  CAS  Google Scholar 

  45. Pandey GK, Kanwar P, Singh A, Steinhorst L, Pandey A, Yadav AK, et al. Calcineurin B-like protein-interacting protein kinase CIPK21 regulates osmotic and salt stress responses in Arabidopsis. Plant Physiol. 2015;169(1):780.

    Article  CAS  Google Scholar 

  46. Chinnusamy V, Schumaker K, Zhu JK. Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. J Exp Bot. 2004;55(395):225–36.

    Article  CAS  Google Scholar 

  47. Kim KN, Cheong YH, Grant JJ, Pandey GK, Luan S. CIPK3, a calcium sensor-associated protein kinase that regulates abscisic acid and cold signal transduction in Arabidopsis. Plant Cell. 2003;15(2):411–23.

    Article  CAS  Google Scholar 

  48. Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44(D1):D279–85.

    Article  CAS  Google Scholar 

  49. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.

    CAS  Google Scholar 

  50. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.

    Article  CAS  Google Scholar 

  51. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

    Article  CAS  Google Scholar 

  52. Yang XY, Guo T, Li J, Chen Z, Guo B, An XM. Genome-wide analysis of the MYB-related transcription factor family and associated responses to abiotic stressors in populus. Int J Biol Macromol. 2021;191:359–76.

    Article  CAS  Google Scholar 

  53. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.

    Article  CAS  Google Scholar 

  54. Wang YJ, Zhang Y, Fan CJ, Wei YC, Meng JX, Li Z, et al. Genome-wide analysis of MYB transcription factors and their responses to salt stress in Casuarina equisetifolia. BMC Plant Biol. 2021;21(1):17.

    CAS  Google Scholar 

  55. Chen ZW, Zhou LH, Jiang PP, Lu RJ, Halford NG, Liu CH. Genome-wide identification of sucrose nonfermenting-1-related protein kinase (SnRK) genes in barley and RNA-seq analyses of their expression in response to abscisic acid treatment. BMC Genomics. 2021;22(1):16.

    Article  Google Scholar 

  56. Fan C, Qiu Z, Zeng B, Liu Y, Li X, Guo G. Selection of reference genes for quantitative real-time PCR in Casuarina equisetifolia under salt stress. Biol Plant. 2017;61(3):463–72.

    Article  CAS  Google Scholar 

  57. Bryfczynski SP, Pargas RP. GraphPad: a graph creation tool for CS2/CS7. 2009. p. 389.

    Google Scholar 

  58. Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    Article  CAS  Google Scholar 

  59. Zhang YH, Wan SQ, Wang WD, Chen JF, Huang LL, et al. Genome-wide identification and characterization of the CsSnRK2 family in Camellia sinensis. Plant Physiol Biochem. 2018;132:287–96.

    Article  CAS  Google Scholar 

Download references


Not applicable.


This study was financially supported by the Open Project of Key Laboratory of Ministry of Education on Tropical Island Ecology (HNSF-OP-202001), the National Natural Science Foundation of China (Grant No. 31901334) and the 7th Scientific and Technological Research project of Forest Seedlings in Fujian province (FJZM 2021–03). The funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



Conceived and designed the experiments: DA, YJW, JZ and YZ. Performed the experiments: DA, YCW. Analyzed the data: DA, YJW, JXM and JZ. Wrote the paper: DA. Participated in the design of this study and revised manuscript: DA, and YZ. The authors read and approved the final manuscript.

Authors’ information

Not applicable.

Corresponding author

Correspondence to Yong Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Figure S1. Tertiary structures of the CeqSnRKproteins. Protein models were obtained using the SWISS-MODEL online server.

Additional file 2:

 Figure S2.Ka/Ks ratios of orthologs.Left image represents orthologs of CeqSnRK with EgrSnRK and PtSnRK.Right image represents orthologs of CeqSnRKwith AtSnRK and OsSnRK.

Additional file 3:

 Figure S3.Relative expression of CeqSnRK genes in roots following different NaCltreatments as determined by RT-qPCR. The Y-axis and X-axisindicate relative expression levels and salt concentration of stress treatment,respectively. Mean values and standard deviations (SDs) were obtained from threebiological and three technical replicates. The error bars indicate standarddeviation. **P < 0.01 and *P < 0.05. 

Additional file 4:

 Figure S4.Relative expression of CeqSnRK genes in shoots following different NaCltreatments as determined by RT-qPCR. The Y-axis and X-axisindicate relative expression levels and salt concentration of stress treatment,respectively. Mean values and standard deviations (SDs) were obtained fromthree biological and three technical replicates. The error bars indicate standarddeviation. **P < 0.01 and *P < 0.05.

Additional file 5:

 Figure S5.Relative expression of 25 selected CeqSnRKgenes in shoots at different times following NaCl treatment as determined byRT-qPCR.The Y-axis and X-axis indicaterelative expression levels and the time course of stress treatment,respectively. Mean values and standard deviations (SDs) were obtained fromthree biological and three technical replicates. The error bars indicatestandard deviation. **P < 0.01 and *P < 0.05.

Additional file 6:

 Figure S6. Correlation matrixof expression among the CeqSnRK genesunder NaCl treatment. Correlation analysis of the expression of CeqSnRK genes under NaCl treatment inroots (left) and shoots (right). Correlations are indicated by the size andcolor of circles. * and ** represent correlations with P-value ≤ 0.05and P-value ≤ 0.01, respectively. 

Additional file 7:

 Table S1 Details of the identified CeqSnRK genes. Table S2. The prediction of subcellularlocalization. Table S3. Details of SnRKgenes from Arabidopsis, rice, grapevine, and Populous trichocarpa.Table S4. Detailed information for the 20 motifs in the SnRK proteins of Casuarina equisetifolia. Table S5. Ka/Ks ratios of paralogous andorthologous gene pairs. Table S6. Promoteranalysis of the CeqSnRK genes. TableS7. Promoter analysis of the CeqSnRK genes. Table S8. List of primer sequences used for RT-qPCRanalysis of the CeqSnRK genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ai, D., Wang, Y., Wei, Y. et al. Comprehensive identification and expression analyses of the SnRK gene family in Casuarina equisetifolia in response to salt stress. BMC Plant Biol 22, 572 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Phylogenetic analysis
  • Cis-regulatory element analysis
  • Expression pattern
  • Salt stress