Identification and biochemical characteristics of CDPK and CRK genes in Solanum Pennellii
In this study, 39 non-redundant sequences were collected in the whole genome of Solanum pennellii, including 32 SpCDPK and 7 SpCRK. All of the SpCRK and SpCDPK proteins were designated separately as SpCRK1-SpCRK7 and SpCDPK1-SpCDPK32 based on their accession locus. By analyzing the physicochemical properties of SpCDPK and SpCRK proteins, it was found that the length of open reading frame of the SpCDPK proteins ranged from 511–607 aa, the molecular weight ranged from 55.9–68.0 kDa and the isoelectric point ranged from 5.0–9.0. The length of open reading frame, molecular weight, isoelectric point of SpCRK protein ranged from 501–607 aa, 56.4–68.2 kDa, 6.0-9.3, respectively (Table 1).
The feature domain analysis showed that the SpCRK carried ATP-binding region and Ser/Thr protein kinase region but lacked the EF-hand calcium-binding domain. However, except for SpCDPK3, SpCDPK5, SpCDPK27 and SpCDPK31, the ATP-binding region, Ser/Thr protein kinase region and four EF-hand calcium-binding domains existed in the remaining SpCDPK. Subcellular localization analysis showed that these genes acted in different organelles (Table 1).
Table 1
Physicochemical characterization of SpCRK and SpCDPK gene family proteins in Solanum Pennellii
Gene name
|
Accession locus
|
ORF length
(aa)
|
Protein kinase domain
|
The number of EF-Hand calcium-binding region
|
Mol Wt (kDa)
|
pi
|
N-myristoylation site
|
N-Palmitoylation site
|
Sub-cell localizationa
|
Protein kinases ATP-binding region
|
Ser/Thr protein kinases region
|
SpCRK1
|
XP_015055812
|
574
|
+
|
+
|
|
64.3
|
8.7
|
+
|
+
|
cyto
|
SpCRK2
|
XP_015062453
|
588
|
+
|
+
|
|
65.9
|
9.0
|
+
|
+
|
chlo
|
SpCRK3
|
XP_015063635
|
589
|
+
|
+
|
|
66.1
|
7.7
|
+
|
+
|
E.R.
|
SpCRK4
|
XP_015067044
|
607
|
+
|
+
|
|
68.0
|
9.0
|
+
|
+
|
chlo
|
SpCRK5
|
XP_015070667
|
511
|
+
|
+
|
|
55.9
|
6.0
|
|
+
|
cyto
|
SpCRK6
|
XP_015071034
|
598
|
+
|
+
|
|
66.9
|
9.0
|
+
|
+
|
nucl
|
SpCRK7
|
XP_015076220
|
587
|
+
|
+
|
|
65.8
|
8.5
|
+
|
+
|
cyto
|
SpCDPK1
|
XP_015055203
|
579
|
+
|
+
|
4
|
63.5
|
5.0
|
|
+
|
chlo
|
SpCDPK2
|
XP_015055554
|
544
|
+
|
+
|
4
|
60.4
|
5.5
|
+
|
+
|
cyto
|
SpCDPK3
|
XP_015055614
|
525
|
+
|
+
|
3
|
59.6
|
5.9
|
+
|
+
|
chlo
|
SpCDPK4
|
XP_015055658
|
557
|
+
|
+
|
4
|
62.2
|
5.6
|
|
+
|
nucl
|
SpCDPK5
|
XP_015056683
|
533
|
+
|
+
|
3
|
59.7
|
6.0
|
+
|
+
|
cyto
|
SpCDPK6
|
XP_015057806
|
528
|
+
|
+
|
4
|
59.4
|
5.7
|
+
|
+
|
cyto
|
SpCDPK7
|
XP_015058701
|
578
|
+
|
+
|
4
|
64.8
|
5.5
|
|
+
|
chlo
|
SpCDPK8
|
XP_015058992
|
505
|
+
|
+
|
4
|
56.8
|
5.4
|
|
+
|
nucl
|
SpCDPK9
|
XP_015059608
|
535
|
+
|
+
|
4
|
60.0
|
5.6
|
|
+
|
cyto
|
SpCDPK10
|
XP_015061184
|
525
|
+
|
+
|
4
|
58.4
|
6.5
|
+
|
+
|
mito
|
SpCDPK11
|
XP_015061185
|
518
|
+
|
+
|
4
|
57.5
|
6.5
|
+
|
+
|
mito
|
SpCDPK12
|
XP_015061191
|
535
|
+
|
+
|
4
|
59.6
|
5.4
|
+
|
+
|
cyto
|
SpCDPK13
|
XP_015063064
|
539
|
+
|
+
|
4
|
60.9
|
5.5
|
+
|
+
|
cyto
|
SpCDPK14
|
XP_015064109
|
570
|
+
|
+
|
4
|
64.3
|
9.3
|
|
+
|
chlo
|
SpCDPK15
|
XP_015064388
|
526
|
+
|
+
|
4
|
60.1
|
6.1
|
+
|
+
|
cyto
|
SpCDPK16
|
XP_015070509
|
573
|
+
|
+
|
4
|
64.9
|
9.0
|
+
|
+
|
chlo
|
SpCDPK17
|
XP_015070510
|
565
|
+
|
+
|
4
|
64.0
|
9.0
|
+
|
+
|
chlo
|
SpCDPK18
|
XP_015070549
|
582
|
+
|
+
|
4
|
64.6
|
5,6
|
|
+
|
chlo
|
SpCDPK19
|
XP_015070664
|
553
|
+
|
+
|
4
|
62.8
|
6.2
|
+
|
+
|
cyto
|
SpCDPK20
|
XP_015071018
|
538
|
+
|
+
|
4
|
60.9
|
6.4
|
|
+
|
cyto
|
SpCDPK21
|
XP_015071199
|
598
|
+
|
+
|
4
|
67.5
|
5.3
|
|
+
|
cyto
|
SpCDPK22
|
XP_015073898
|
527
|
+
|
+
|
4
|
59.5
|
5.3
|
+
|
+
|
cyto
|
SpCDPK23
|
XP_015074249
|
581
|
+
|
+
|
4
|
64.6
|
5.5
|
|
+
|
chlo
|
SpCDPK24
|
XP_015074284
|
508
|
+
|
+
|
4
|
57.2
|
5.1
|
|
+
|
cyto
|
SpCDPK25
|
XP_015076006
|
503
|
+
|
+
|
4
|
56.5
|
5.0
|
|
+
|
cyto
|
SpCDPK26
|
XP_015077134
|
501
|
+
|
+
|
4
|
56.4
|
5.6
|
|
+
|
chlo
|
SpCDPK27
|
XP_015078365
|
533
|
+
|
+
|
3
|
60.0
|
7.0
|
+
|
+
|
pero
|
SpCDPK28
|
XP_015079858
|
536
|
+
|
+
|
4
|
61.0
|
5.8
|
+
|
+
|
mito
|
SpCDPK29
|
XP_015082059
|
521
|
+
|
+
|
4
|
57.8
|
6.6
|
+
|
+
|
pero
|
SpCDPK30
|
XP_015085844
|
516
|
|
+
|
4
|
57.7
|
5.8
|
+
|
+
|
chlo
|
SpCDPK31
|
XP_015088080
|
529
|
+
|
+
|
3
|
59.6
|
6.0
|
|
+
|
cyto
|
SpCDPK32
|
XP_015089561
|
607
|
+
|
+
|
4
|
68.2
|
5.4
|
|
+
|
chlo
|
a Nucl, nuclear; ER, endoplasmic reticulum; Mito, mitochondria; Cyto, cytosol; Chlo, chloroplast; Pero, Pero, peroxisomesa Nucl, nuclear; ER, endoplasmic reticulum; Mito, mitochondria; Cyto, cytosol; Chlo, chloroplast; Pero, Pero, peroxisomes |
Chromosomal localization of SpCRK and SpCDPK genes
The chromosomal localization of SpCRK and SpCDPK genes were shown in Figure 2. The number of genes on chromosome 1 to chromosome 12 were 7, 3, 6, 3, 2, 2, 1, 1, 1, 6, 4 and 3, respectively. Chromosome 7, 8 and 9 involved only one SpCDPK gene. Chromosome 1 possessed the largest number of genes, with 7 genes (5 SpCDPK and 2 SpCRK). Chromosomes 2 and 3 each had 2 SpCRK genes, followed by chromosomes 2, 5 and 10 each had one SpCRK gene.
Protein sequence comparison of SpCRKs and SpCDPKs
In order to investigate the sequence characteristics of SpCRK and SpCDPK proteins, we performed multiple sequence alignments on the protein sequences of SpCRK and SpCDPK (Figure 2). As shown in Figure 2a, SpCRK and SpCDPK proteins contained protein kinase domains. Two highly conserved domains, the protein kinases ATP-binding region (LGxGxFGxTxCGxACKxIxK) and the Ser/Thr protein kinases region (VxHDRLKPENFLx), existed in the protein kinase domain.
The amino acid sequence analysis of the ATP-binding region of SpCDPK and SpCRK proteins was found that, compared with the SpCDPK proteins, the first amino acid in the ATP-binding region of SpCRK genes were mutated from L (leucine) to I (isoleucine)/V (Valine). And the F (Phenylalanine) mutation occurred in SpCDPK15, T (Threonine) mutation in SpCRK3 and SpCDPK15, two C (cysteine) mutations in SpCRK1-4, SpCRK6, SpCRK7, SpCDPK14, SpCDPK16, and SpCDPK17, only one C mutation in SpCDPK3 and SpCDPK30. Obviously, only the fourth G (glycines) in the ATP-binding region was mutated, and this phenomenon occurred in 10 genes. The last amino acid in the ATP-binding region, K (lysine), was mutated in SpCDPK22 and SpCDPK30. There were also a few proteins with amino acid mutations in the Ser/Thr protein kinases region. Two amino acid mutations existed in the Ser/Thr protein kinases region of SpCRK5, SpCDPK14, SpCDPK16 SpCDPK17, and SpCDPK32, while the Ser/Thr protein kinases region of the other genes were intact. Therefore, in the SpCRK and SpCDPK families, the Ser/Thr protein kinases region was more conserved than the ATP-binding region.
As shown in Figure 2b, all SpCDPKs contained the EF-hand calcium-binding domain (DxD/NxGxE), while SpCRKs lacked it. Significantly, the second D (Asparticacid) was mutated to G in the first EF-hand calcium-binding region of SpCDPK3, SpCDPK5, and SpCDPK27, which resulted in the deletion of an EF-hand calcium-binding region. It was also found that the second EF-hand calcium-binding region of SpDPK31 was incorrect, although it was consistent with the characteristic structure of the EF-hand region.
Phylogentic analysis of SpCRKs and SpCDPKs
The phylogenetic relationship about CDPK and CRK family members of Solanum Pennellii (32 SpCDPK and 7 SpCRK), Arabidopsis (34 CPK and 8 CRK) and rice (29 CPK and 5 CRK) were constructed via using the neighbor-joining method of MEGA6.0 (Figure 3). As Figure 3 showed, CDPK and CRK gene families were divided into four groups (I, II, III, and IV). The size of four groups were similar. The 10 CDPK in Solanum Pennellii, 13 CPK in Arabidopsis and 8 CPK in rice were put into group I. The group II included 13 SpCDPK, 10 AtCPK and 11 OsCPK. The group III comprised 6 Solanum Pennellii CDPK, 8 Arabidopsis CPK, and 8 rice CPK. The 20 CRK proteins from the three species were put into group IV. In addition to CRK, 3 CDPK (SpCDPK14, SpCDPK16 and SpCDPK17) of Solanum Pennellii, 2 CPK (AtCPK16, AtCPK18 and AtCPK28) of Arabidopsis and 3 CPK (OsCPK4, OsCPK5 and OsCPK18) of rice existed in group IV. This phenomenon was not discovered in the other group. The dendrogram showed that the 39 proteins of Solanum Pennellii were generally closer to the proteins of Arabidopsis than rice, which indicated that they were evolutionarily more closely related.
Genetic structure analysis of SpCRK and SpCDPK genes
The genetic structure analysis of SpCRK and SpCDPK were carried out, and the results were shown in Figure 4. There were 4 groups in SpCRK and SpCDPK gene families, which were consistent with the respective corresponding phylogenetic relationships in Figure 3. As Figure 4 showed, genes within the same groups exhibited similar exon-intron organizations. The SpCDPK genes of the group I possessed 8-9 exons and 7-8 introns; group II genes included 7-8 exons and 6-7 introns; group III genes contained 7-9 exons and 6-8 introns. In particular, all SpCRK were part of the group IV. Except for SpCRK5 (involved 4 exons and 3 introns), all SpCRK introns and exons were 10 and 11, respectively. But SpCDPK of group IV contained 11 introns and 12 exons (Figure 4). These results showed that group IV exhibited more introns and exons than the other three groups, indicating that the gene structure of Group IV was more complex.
Conserved motif analysis of SpCRK and SpCDPK proteins
The results of motif analysis showed that the SpCRK and SpCDPK family proteins contained obvious structural characteristics. Motifs 9 and 10 were labeled as protein kinases ATP-binding region, motif 3 as Ser/Thr protein kinases region, and motifs 5, 6, 7 and 8 as EF-hand calcium-binding domain (Figure 6). As Figure 5 showed, interestingly, through analysis MEME online website, except for SpCRK5, the other 6 SpCRK included motifs 5, 6, 7 and 8; SpCDPK3, SpCDPK5 and SpCDPK27 included motif 5; and SpCDPK31 included motif 6. But these EF-hand calcium-binding region were not retrieved in Prosite (Table 1). It may be that these motifs were found only as recurring sequences, but they could not constitute EF-hand calcium-binding region structure. The SpCDPK of group I and II contained motifs 1-15, except for SpCDPK25 (Group II). The SpCDPK of Group III showed motifs 1-12 and 15, besides SpCDPK31. All proteins of group IV contained motifs 1-4, 10-11 and 15. These results indicated that all the identified proteins had typical ATP-binding region, Ser/Thr protein kinases region and EF-hand calcium-binding region, and each subgroup had similar motifs, which further supported the phylogenetic classification of SpCRK and SpCDPK families.
Promoter region analysis of SpCDPK and SpCRK genes
The cis elements in the upstream promoter region of SpCRK and SpCDPK genes were analyzed by PlantCARE. The results showed that the promoter region of SpCRK and SpCDPK genes contained 627 cis-acting elements, which can be divided into five types: hormone-related elements, growth-related elements, stress-related elements, secondary metabolite-related element and plant protein metabolism-related element (Figure 7). The hormone-related elements included auxin responsive element, gibberellin responsive element, salicylic acid responsive element, abscisic acid responsive element and MeJA responsive element. The growth-related elements included circadian control responsive element, palisade mesophyll cell differentiation responsive element, endosperm expression response element, seed-specific regulation responsive element, light responsive element and meristem expresses related elements. The stress-related elements included low temperature responsive element, drought responsive element, anoxic induction element, anaerobic induction element. The secondary metabolite-related element and plant protein metabolism-related element was flavonoid biosynthetic genes regulatory elements and zein metabolism regulatory elements, respectively. Among all the elements, the number of light responsive elements was the largest, with 356. Obviously, the number and distribution of cis-acting elements of SpCDPK10 and SpCDPK11 genes were basically the same, and SpCDPK16 and SpCDPK17 also showed similar phenomenon. All elements of SpCDPK20 were located between 900bp and 2000bp upstream sequence. The cis-acting elements of SpCDPK27 were mainly distributed at both ends of the promoter region (Figure 7). These results suggest that SpCRK and SpCDPK genes play an important role not only in the growth and development of Solanum Pennellii, but also in the response to biotic and abiotic stresses.