- Research article
- Open Access
Genome-wide survey of putative Serine/Threonine protein kinases in cyanobacteria
BMC Genomics volume 8, Article number: 395 (2007)
Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution.
286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains.
The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.
Cyanobacteria, dating back 2.5–3.5 billion years and constituting a single but large taxonomic and phylogenetic group within the domain Eubacteria , are characterized by their ability to carry out oxygenic photosynthesis. Moreover, fossilized cyanobacteria appear similar in form to extant species . Cyanobacteria have a pronounced variation in genome size from 1.6 Mb to 9.2 Mb and exhibit extraordinary diversity in terms of morphology and cell activity. They also exhibit the widest range of diversity in ecological habitats of all photosynthetic organisms, including environments that are extremely hot, extremely cold, alkaline and acidic, marine, freshwater, saline, terrestrial, and symbiotic . Prochlorococcus marinus, which has the smallest genome size and can be divided into two distinct ecotypes (high-light adapted and low-light adapted), is the dominant photosynthetic prokaryote in the open ocean . The diazotrophic filamentous cyanobacteria have the largest genome sizes and include strains isolated from fresh water (Anabaena PCC7120), from a plant-cyanobacterial symbionsis (Nostoc punctiforme PCC73102), or from tropical and subtropical oceans (Trichodesmium erythraeum IMS101). Crocosphaera, a novel genus of marine unicellular diazotrophic cyanobacterium, and Gloeobacter, a rod-shaped unicellular cyanobacterium isolated from calcareous rocks, have larger genome sizes (6.3 Mb and 4.6 Mb) than other unicellular cyanobacteria.
The diversity of cyanobacteria is also reflected in the complexity of their signal transduction systems. To cope with changing environmental conditions, cyanobacteria have developed a variety of adaptive mechanisms to respond to external or internal changes. Two-component signal transduction systems, characterized by the transfer of phosphate by a sensor kinase from a His residue on the enzyme to an Asp residue on the response regulator, are widely distributed among bacteria [5, 6]. One-component systems, defined as proteins that contain known or predicted input and output domains in a single protein molecule but lack histidine kinase and receiver domains, are considered to be the pre-eminent mechanism for signal transduction in bacteria and archaea, except for cyanobacteria . In contrast, the Ser/Thr-specific protein kinases (STKs) serve as the backbone of the eukaryotes transduction network. However, with the first identification of an STK in Myxococcus xanthus in 1991 , regulatory STKs have been repeatedly identified in prokaryotes. Protein phosphorylation on serine/threonine residues in cyanobacteria was first revealed by radioactive labeling of proteins in 1994 . Numerous bacterial STK genes have since been predicted within genome sequences [10–12], and they have been associated primarily with three different processes, namely regulation of development, stress responses, and pathogenicity.
According to Hanks and Hunter, canonical Ser/Thr protein kinases contain 12 conserved subdomains  that fold into a common catalytic core structure, as revealed by the 3-dimensional structures of several protein-serine kinases. These 12 conserved sequence motifs, about 280 amino acids in length, are named as subdomains I-V, VIa, VIb and VII-XI [13, 14]. There are three structural subdomains with separate roles: N-terminal nucleotide-binding domains containing subdomains I-IV, C-terminal phosphotransfer and protein-substrate-binding domains containing subdomains VIa-XI, and the intervening linker containing subdomain V . In addition to the conserved catalytic domains, some STKs contain at least one additional domain, such as FHA, WD40, PAS and GAF , which may endow STKs with more complicated functions. For example, the FHA domain can participate in a wide range of processes in bacteria, such as intracellular signal transduction, transcription, protein transport, DNA repair, and protein degradation. Moreover, there is a near-perfect correlation between the presence of FHA-containing proteins and Ser/Thr kinases and phosphatases in bacterial genomes .
In cyanobacteria, mutation analyses have revealed the functions of some STKs, such as SpkA and SpkB involved in the cellular motility in Synechocystis PCC6803 [18, 19], and PknD involved in the regulation of nitrogen metabolism in Anabaena PCC7120 . A complete signaling pathway involving a bacteria STK remains to be described, but one may anticipate that some STKs are important or even essential in regulating bacterial activities.
As of January 2007, 19 cyanobacterial genomes have been fully sequenced and 6 are in the draft stages. In addition, more than 20 genomes are in the process of being sequenced. The availability of multiple sequenced genomes has been very helpful in phylogenetic and functional studies of cyanobacterial genes, such as those in restriction-modification systems  and two-component systems . Comparative genome analysis has been employed in the study on STKs in archaea , mycobacteria  streptomyces  and cyanobacteria [25, 26]. Here we present a detailed analysis of the repertoire of STKs in sequenced cyanobacterial genomes. An attempt has been made to identify all STK sequences encoded in the 19 fully sequenced genomes and 2 draft sequences (Crocosphaera watsonii WH8501 and Nostoc punctiforme PCC73102). An analysis of STK sequences according to their classifications, conservation, domain organization, phylogeny and evolution are presented in this paper, aiming towards a deep understanding of the biological role of STKs in cyanobacteria.
Identification of STK proteins
The 21 cyanobacterial genomes fully available from IMG database (Version 2.0, as of Jan 2007) were considered in our analysis (Fig. 1, Table 1). Synechococcus elongatus PCC6301 and Synechococcus elongatus PCC7942 are virtually identical except for a chromosomal inversion, but both were included in our analysis. Using BlastP and TBlastN programs to look for proteins similar to proven cyanobacterial STKs, we obtained 303 protein sequences from the 21 cyanobacterial genomes, 284 of which were originally annotated as protein kinase or Serine/Threonine Protein Kinase. The remaining 19 proteins were accepted as STKs in this study after sequence alignment and SMART analysis . Seventeen out of these were originally annotated not by their STK domains but by other additional domains, such as ATPase and WD40 repeats, and the last two were annotated as hypothetical proteins .
17 proteins, all originally annotated as STKs, from the 303 proteins identified by sequence similarity lack all or part of at least one important catalytic domains , and these were excluded from further consideration. Some slightly truncated proteins, such as CwatDRAFT_6607 and CwatDRAFT_1979, lacking structural motifs I, II or XI but possessing all catalytic motifs were retained. Nine proteins (including four STKs in marine unicellular Prochlorococcus and Synechococcus, 6803_slr1443, CwatDRAFT_0901, Tery_0088, Tery_4781 and Ava_3867) exhibit some deviations from the canonical motifs (Fig. 2) that were deemed tolerable. In all 286 putative STK sequences were considered in this study (Table 1).
The number of STK genes within the different cyanobacterial genomes varies considerably, from 0 to 56. Only one STK gene was found within the 5 fully sequenced Prochlorococcus marinus strains, in strain MIT9313. Eight out of nine sequenced Synechococcus strains contain STKs, with the marine Synechococcus strain WH8102 being the only exception. Three marine Synechococcus strains CC9902, CC9605, and CC9311 contain only one STK gene for each. Two Synechococcus elongatus strains PCC6301 and PCC7942, which are virtually identical in sequence, have five pairs of orthologous STK genes. More STK genes were found in the thermophilic Synechococcus strains JA-3-3Ab (8), JA-2-3Ba (8), and BP-1 (10). The other two unicellular cyanobacterial strains, Gloeobacter violaceus PCC7421 and Crocosphaera watsonii WH8501, have 14 and 28 STK genes respectively.
Filamentous diazotrophic cyanobacteria have the largest number of STK genes (40 for Trichodesmium erythraeum IMS101, 48 for Anabaena PCC7120, 53 for Anabaena variabilis ATCC 29413, and 56 for Nostoc punctiforme PCC73102). Among the unicellular cyanobacteria, the nitrogen-fixing unicellular strain Crocosphaera WH8501 has the biggest genome, the largest number of STKs and the highest percentage (0.52%) of STKs within protein-encoding genes (Fig. 1). STKs in Prochlorococcus MIT9313 represent only 0.04% of the total proteins, while in Anabaena ATCC29413, the percentage is twenty-fold higher at 0.85%. It is evident from these findings that filamentous diazotrophic cyanobacteria contain more STK genes than unicellular species (T-test, p < 0.001), and the number of STK genes is overrepresented in their genomes even after allowing for their larger genome sizes (Fig. 1).
Conserved domain features
The core catalytic domain contains 12 conserved subdomains, within which a total of 12 amino acids are invariant or nearly invariant in eukaryotes: two G residues in subdomain I, K in subdomain II, E in subdomain III, D and N in subdomain VIb, D and G in subdomain VII, E in subdomain VIII, D and G in subdomain IX, and R in subdomain XI . In cyanobacteria, 12 similar conserved subdomains were found and 10 residues are conserved well, with the exception of the two G residues in subdomain I, which are common but not invariant. It should be noted that additional conserved amino acids appear in cyanobacterial subdomains I, III, V, VIa, VIb, VII, VIII and IX (Fig. 2).
In eukaryotes, subdomain I serves as the ATP binding site (P-loop) and is characterized by the sequence GXGXXG. Cyanobacterial STKs are divergent in this motif with only 152 (53.3%) sequences exactly conserved in this domain, but 94.4% of them contain at least one Gly. It is possible that STKs lacking a conventional subdomain I may employ high energy phosphate compounds other than ATP, as does a PK from Sulfolobus acidocaldarius that was reported to use polyphosphate as phosphoryl donor [28, 29]. Besides the conserved Gly, cyanobacterial subdomain I differs from its eukaryotic counterpart in that it starts at Tyr (93% of all cyanobacterial sequences) and ends with a non-polar residues (99%).
The Lys (96.1%) in subdomain II is important for maximal catalytic activity and helping anchor and orient ATP, and it is well conserved in cyanobacteria. But the upstream Ala described by Hanks and Hunter  is generally replaced by a different non-polar residue (97.2%). In subdomain III, Glu (86.0%) is not conserved well and is substituted by Gln in 27 sequences belonging to cyanobacterial Family III (cbSTKIII, details in section Structure and function). Upstream of the Glu, a conserved non-polar residue (99.3%) was found and most of them were Phe (87.7%).
Cyanobacterial subdomain V is much more conserved than the eukaryotic counterpart, which has no invariant residues. In contrast, the cyanobacterial subdomain possesses a highly conserved Gly (94.1%). This conserved residue may help anchor ATP by forming hydrogen bonds. Gly is also conserved in other bacteria (e.g. Streptomyces coelicolor ). In eukaryotes, subdomain VIa contributes to structural stabilization and possesses only a relatively conserved His. In cyanobacteria, besides this His (88.1%), a more conserved Leu (97.2%) has been found, which is also conserved well in Mycobacterium tuberculosis .
Subdomain VIb, a candidate for direct catalysis, is the most conserved motif in eukaryotes and is characterized as "hrDxkxxN". Cyanobacterial subdomain VIb is more conserved and is characterized as "HrDikPxN". His (100%) is the most conserved residue, followed by Asn (98.9%), Pro (98.6%), Asp (97.2%), Lys (86.7%), Ile (82.5%) and Arg (73.8%). Pro is also conserved in other bacteria, such as Streptomyces coelicolor  and Mycobacterium tuberculosis . STKs are also referred as "RD" kinases because their activation requires the Arg and catalytic Asp in this subdomain. In cyanobacteria, however, Arg is not conserved well and it can be substituted by Lys, Gly, Cys or Gln. All members of cyanobacterial Family III (cbSTKIII) belong to the family of "KD" kinases. There are five cyanobacterial STKs lacking the highly conserved D, a candidate for the catalytic base. Among these five putative STKs, SpkE (Slr1443), which also lack of several other conserved residues, did not show any kinase activity . Whether the other four could phosphorylate or not may need more experimental proofs.
Subdomain VII helps to orient the γ phosphate of ATP for transfer and is characterized by a conserved DFG triplet. The invariant Asp in subdomain VII, as well as the invariant Lys in subdomain II, was found to anchor and orient ATP . In cyanobacteria, Asp (97.2%) and Phe (97.9%) are highly conserved but Gly (87.4%) less so. In 30 STK proteins from cyanobacterial Family III (cbSTKIII), Gly is substituted by Ser. The Lys residue often found in eukaryotes before the DFG triplet is not conserved in cyanobacteria, and the residues with near neutral polarity were replaced by non-polar ones (95.1%).
The typical subdomain VIII includes a highly conserved PE motif and plays a major role in the recognition of peptide or protein substrates. In cyanobacteria, the Pro residue (62.2%) of this motif is not well conserved, and the Glu residue (97.5%) is followed by a conserved Gln residue (90.5%) of unknown function. The PE motif in archaea, eubacteria and eukaryotes may represent an ancient and universal substrate recognition mechanism.
In subdomain IX, besides the Asp (99.3%), Gly (84.6%) and Pro (90.2%) residues conserved in eukaryotes, another three conserved amino acids were found: Thr (90.9%), Leu (87.0%) and Gly (80.1%). Gly is also conserved in Mycobacterium tuberculosis  and Streptomyces coelicolor .
In cyanobacteria there are six STKs, SpkA (Sll1574) , SpkB (Slr1697) , SpkC (Slr0599), SpkD (Sll0776), SpkF (Slr1225) , PknC (All4813) , which have been shown to phosphorylate themselves or other substrate proteins. From the sequence alignment it can seen that all six of these STKs have sequences that are in accordance with the conserved features found in cyanobacteria (Fig. 2).
Structure and function
Cyanobacterial STKs were classified according to structural characteristics into three major families: cbSTKI, cbSTKII and cbSTKIII (Fig. 3).
Cyanobacterial STK Family I (cbSTKI) groups together the 155 (54.2%) STK proteins that possess no identifiable domain besides that associated with core catalysis. This family was further divided into three subfamilies. Subfamily I (cbSTKI-I) contains 35 proteins with fewer than 400 amino acids long and may be considered to possess the basic structure supplemented in other families by additional domains. Proteins in this family are found only in Gloeobacter PCC7421, Crocosphaera WH8501 and the four filamentous diazotrophic strains. Subfamily II (cbSTKI-TM) contains 13 proteins and is named for their 3 to 6 transmemberane (TM) domains behind the STK domains. These 13 proteins are distributed amongst unicellular Synechocystis PCC6803, Thermosynechococcus elongatus BP1, Crocosphaera watsonii WH 8501 and the four filamentous diazotrophic strains. Subfamily III (cbSTKI-other) consists of 107 STK proteins distributed amongst all 16 strains, each protein with more than 400 amino acids and possessing a long unidentified C-terminal domains. The functions of two cbSTKI proteins in Synechocystis PCC6803 and three in Anabaena PCC7120 have been confirmed by experiment. In Synechocystis PCC6803, Sll1574 (SpkA, I-other) was found to be involved in regulating cellular motility via phosphorylation of membrane proteins . No phenotype defect could be discerned in a knockout mutant of Slr0599 (SpkC, I-other); however the protein was demonstrated to phosphorylate itself and some general substrate proteins . In Anabaena PCC7120, Alr4366 (PknA, I-other) showed a complex pattern of regulation during heterocyst development , and Alr4368 (PknD, I-other) appeared essential for normal growth under diazotrophic conditions . The function of Alr3732 (PknE, I-other) is known to be required for the formation of heterocyst envelope structures and nitrogen fixation .
Cyanobacterial STK Family II (cbSTKII) includes 92 (32.2%) proteins from 12 cyanobacterial species and are defined as those STK proteins that possess at least one additional domain in addition to their STK domains but do not possess Hik kinase domains. The additional domains with distinct functions are prevalent in protein kinases and phosphatases of prokaryotic and eukaryotic signal transduction systems. In total, 14 types of additional functional domains were identified in cyanobacterial STKs: PAS, PAC, GAF, ANF, WD40, FHA, GUN4, TPR, DUF323, PbH1, CHASE2, Pentapeptide, RDD, and SH3b. Most are involved in known signaling proteins and play various functions in the signal transduction process. For example, TPR (tetratricopeptide repeat), which are repeated in chaperone, cell-cycle, transcription, and protein transport complexes, form signaling domains in higher eukaryotes and may participate in the aggregation of proteins into multi-protein complexes . GAF domains (cyclic GMP, adenylyl cyclases, FhlA) are linked to the binding of small-molecules, in particular the second messengers cyclic AMP (cAMP) and cGMP . WD40 repeats, mainly in eukaryotic proteins, cover a wide variety of functions, for example, in adaptor/regulatory modules in signal transduction [36, 37]. Finally, GUN4 is thought to participate in plastid-to-nucleus signaling by regulating magnesium-protoporphyrin IX synthesis or trafficking . There are also some domains with unknown functions. These are found mainly in bacteria and include DUF323, Pentapeptide, and RDD. The putative functions of most of these additional domains have been described in detail by Krupa . Most of these additional domains are located at C-terminal, with the exceptions of FHA and CHASE2.
Unicellular non-nitrogen-fixing strains have none or a limited number of additional domains. TPR is identified exclusively in Gloeobacter violaceus PCC7421, and DUF323 is the most prevalent in Trichodesmium erythraeum IMS101. Crocosphaera watsonii WH8501 has the lowest percent of additional domains except for marine unicellular non-nitrogen-fixing strains (Fig. 1). Filamentous diazotrophic cyanobacteria not only contain vast number of STKs but their STKs also recruit a diversity of additional domains. The STKs of Anabaena PCC7120, for example, contain 12 types of different additional domains. The marine filamentous diazotrophic strain Trichodesmium shows an exceptionally high percent (77.5%) and the greatest number of additional domains. Pentapeptide repeats, which have been found to participate in the accumulation of glycolipids into the heterocysts , are prevalent in both unicellular and filamentous cyanobacteria. WD40, containing 6 or 7 repeats, are found exclusively in nitrogen-fixing strains. Functions of three proteins in this family have been identified. Alr2502 (Pkn22, PbH1) was found to be associated with both iron-depletion and oxidative stress in Anabaena PCC7120 . In Synechocystis PCC6803, Slr1697 (SpkB, Pentapeptide) was shown to be involved in the control of cell motility exclusive of positive phototaxis , and Sll0776 (SpkD, SH3b) could not be knocked out completely, indicating that it is essential for survival .
Kinases of STK Family III (cbSTKIII) are also named as dual protein kinases, as they contain both N-terminal STK domains and C-terminal His kinase domains, with at least one GAF domain in between. They are encountered in quite large numbers (38) in filamentous nitrogen-fixing strains, with the exception of Trichodesmium erythraeum IMS101, which possesses only one. This observation is consistent with an earlier comparative study reporting that GAF and PAS domains, possibly involved in signal recognition, are extremely abundant in Anabaena: 87 GAF domains in 62 ORFs and 140 PAS domains in 59 ORFs . The expression of one cbSTKIII protein, Alr2258 (HstK) from Anabaena PCC7120, has been shown to depend on the type of available nitrogen source , and the expression of two others, Alr0709 (Pkn40) and Alr0710 (Pkn41), whose genes are adjacent along the chromosome and co-transcribed, are induced by iron deprivation and are under the control of the global nitrogen-regulator NtcA .
Phylogenetic analysis was performed using the conserved catalytic domains of STKs rather than their whole sequences, as the additional domains with their possibly separate evolutionary histories might confuse the phylogeny. The STK catalytic domains, about 280 amino acids in length, were identified using SMART and CDD databases [27, 44].
The STK phylogenetic tree was rooted in the archaeal STKs, deviating from the canonical view of the universal phylogenetic tree rooted in eubacteria . This STK phylogenetic tree indicates that these STK proteins display complicated relationships, at odds with the phylogeny of the species. STK proteins from the unicellular marine Synechococcus and Prochlorococcus cluster together and have a close relationship with Slr1443 from Synechocystis PCC6803. STKs from the unicellular cyanobacteria, Synechocystis PCC6803, Synechococcus elongatus PCC6301, Synechococcus elongatus PCC 7942, and Thermosynechococcus elongatus BP1 are dispersed throughout the phylogenetic tree. Seven pairs of STK orthologs are found in two Synechococcus strains, JA-2-3Ba and JA-3-3Ab, that share a close evolutionary relationship and are located at the bottom of 16S rRNA tree, along with Gloeobacter violaceus PCC7421 (Fig. 1). STKs from Trichodesmium erythraeum IMS101 and Gloeobacter violaceus PCC7421 each form several separate clusters indicating obvious lineage-specific duplication events in these strains. Most kinase sequences from Anabaena and Nostoc are orthologs for their close evolutionary relationships, while the nonorthologous STK genes in Anabaena and Nostoc may have been produced by recent gene duplication or lateral transfer. This phylogenetic tree shows also that cyanobacterial STKs do not cluster strictly according to their structural characteristics, except for members of cbSTKIII which are clustered together in a big clade. Most STKs of cbSTKII do not cluster according to their additional domains, while some with specific additional domains, such as Pentapeptide and WD40, do cluster together (Fig. 4).
Thirty-eight STKs in cbSTKIII belong to a large cluster that also includes four STK proteins that do not contain HisKA and ATPase domains and therefore have been placed in other families. Tery_2857 and Tery_3863 have GAF domains and belong to cbSTKII, while Tery_2107 and Npun02008276 have only the STK domains and belong to cbSTKI-other (Fig. 4). To test the evolutionary relationships amongst the catalytic domains of dual protein kinases, their HiskA domains, and their GAF domains, a detailed analysis was performed by constructing three separate phylogenetic trees (Fig. 5). Only two groups of orthologs as defined by the catalytic domain (highlighted in Fig. 5) clustered together in the other two phylogenetic trees, indicating that two original genes were formed and appeared in the lineage prior to the divergence of Nostoc and Anabaena. Genes ancestral to orthologs that are conserved only in Anabaena PCC7120 and Anabaena ATCC29413 are assumed to have formed before the divergence of these two strains. Apparent incongruities may be explained by new recruitment or domain shuffling. For example, the HiskA domain of N. punctiforme gene Npun02001148 appears to have a distinct history from the HiskA domains of Anabaena PCC7120 All0323 and Anabaena ATCC29413 Ava_4716, but the three proteins share a common history with respect to their catalytic and GAF domains (see proteins highlighted in blue in Fig. 5) The second GAF domains in All3691, Ava_3596, Npun02004457, and Npun02004756 and the third GAF domain in Npun02004457 bear only a distant relationships with other STK GAF domains in cyanobacteria, including those in the same genes, indicating a recently recruitment. Alr0709 and Alr0710 lie adjacent to each other on the chromosome of Anabaena PCC7120 and show 75% similarities, indicating a recently gene duplication.
Ser/Thr protein kinases, members of signal transduction systems, play important roles in responses to the environmental changes and intra-cellular signals . This may be vital for their host cyanobacteria, which originated 2.5–3.5 billion years ago and exhibit the widest range of diversity in growth habitats of all photosynthetic organisms. Some of these kinases have demonstrated, in at least one case, possibly essential roles, illustrated by Sll0776 (SpkD) in Synechocystis PCC6803 which could not be completely knocked out .
The STK proteins used in this study were identified by Local Blast rather than from the COG database in IMG, and were manually checked to avoid false-negative and false-positive hits that commonly arise during large-scale automated analyses. Sequence alignment of conserved catalytic domains show that cyanobacterial STKs contain nearly the same subdomains as eukaryotic STKs but seem to be more conserved except in subdomain I. Eukaryotic STK proteins shared 12 nearly identical residues. Whereas a total of 15 highly conserved (shared by >90% sequences) sites were found in cyanobacteria, seven of them were the same in eukaryotes. Such conserved motifs and amino acids may provide useful targets to evaluate their functions, and the high proportion of conserved residues shared by both eukaryotic organisms and cyanobacteria indicate a quite similar mechanism for domain organization and catalyzation.
Signal transduction systems provide basic mechanisms for cellular responses to environmental changes and connect environmental signals to cellular activities . Thus different distributions of STKs through gene gain or loss may reflect various environmental selective pressures. Such as Synechococcus strains, which have similar genome size, hold different numbers of STK genes from zero to ten mainly for the different environmental conditions. Moreover, the number of STK genes in cyanobacteria is a function of the genome size, ecophysiology, and physiological properties of organism, as has been noted in general for bacteria [46, 47]. STK genes are not ubiquitous in cyanobacteria. They are absent in four Prochlorococcus strains and one marine Synechococcus WH8102. Previous study showed that signal transduction and environmental stress response systems are dramatically reduced in marine unicellular Prochlorococcus and Synechococcus , both of whom live in the oligotrophic open ocean, and the major driving force behind might be a selective process favoring the adaptation of these cyanobacteria to adapt to the oligotrophic environment . In contrast, filamentous heterocystous cyanobacteria, which differentiate heterocysts in response to the absence of combined nitrogen, and which display physiological and ecological properties including broad symbiotic competence with plants and fungi , showed a disproportionate number of STK genes.
Although Ser/Thr protein kinases were first discovered in eukaryotes and are evidently more widely dispersed amongst them, this is an insufficient basis to conclude that STKs originated in eukaryotes and were obtained by prokaryotes through lateral gene transfer. Zhang  and Leonard  have demonstrated that an ancestral protein kinase existed prior to the divergence of eucarya, bacteria, and archaea. The phylogenetic tree from our analysis shows STKs from archaea clustering at the root and cyanobacterial STKs clustering together with those from eukaryotes. The phylogenetic tree of cyanobacterial STKs, which was constructed by their conserved catalytic domains, is very complicated, and the distribution of STKs along the tree does not follow their domain structures or species phylogenies, indicating that frequent gene gain-and-loss events in cyanobacteria may obscure the actual relationships of STKs.
Besides gene gain-and-loss, novel proteins could also be generated by inserting or shuffling their domains. In these 21 cyanobacteria genomes, a total of 14 types of additional functional domains were found, much more than in other bacteria such as Mycobacterium and Streptomyces [23, 24]. Additional domains are usually employed as sensor response modules [52–55] and can help organisms to assemble a sophisticated signal transduction apparatus. Their modularity may facilitate their spread into new genes and the rapid creation of novel functions . Some additional domains (e.g. TPR) may have existed in the putative common ancestor of all life forms, and some (e.g. FHA) are thought to have been transferred from eukaryotes to bacteria . Some prevalent additional domains, such as pentapeptides that are widely distributed in bacteria, may have been obtained from the ancestor at an early evolutional stage. The species-specific additional domains might have been recruited after speciation and duplicated recently, as in the example of the DUF323 domain in Trichodesmium erythraeum IMS101. Four STKs with DUF323 domains of T. erythraeum _0640–0643, which have high sequence similarities and neighbor on the chromosome, are closest paralogues and might have arisen from a recent duplication. STKs that cluster together according to their conserved catalytic domains do not always do so according to their additional domains, indicating the lateral recruitment of additional domains by STKs. A similar phenomenon was previously observed regarding the GAF and PAS domains found in PPM-family protein phosphatases in Streptomyces . The two clusters with FHA and Pentapeptide domains are indicative of an early recruitment. GAF domains are prevalent in cyanobacterial two-component system and exist in all sequenced cyanobacterial genomes. But different STK orthologues did not recruit the same GAF and HiskA orthologs. It is assumed that HiskA and GAF domains were added to STK domains together and followed by domain shuffling. Orthologs that are conserved in three strains indicate their existence before the divergence of Anabaena and Nostoc, and the orthologs conserved in only two Anabaena strains suggest that they appeared before the divergence of the two Anabaena strains and after the split of Anabaena and Nostoc. Some sequences appear to have lost the additional domains, and some appear to have gained PAC and PAS domains, indicating that gene shuffling, deletion and insertion, as well as gene duplication, occurred during evolution.
Results from PPM-family protein phosphatases in Streptomyces, showed that PAS and GAF domains were clustered by physiological functions rather than taxonomic relationship . So functions of substrate recognition in cyanobacterial STKs may be conferred largely by their additional domains too. For example, the CHASE (cyclase/histidine kinase-associated sensing extracellular) domain, always followed by three transmembrane regions, is present in more than one type of sensory proteins and can recognize several signals, such as cytokines and short peptides that are important for the development of an organism . STKs, lacking transmembrane regions and acting as a receptor-like protein kinase in the cytoplasm, also appeared to have recruitted some additional domains such as FHA and GAF. The FHA domain is thought to play roles in a wide range of processes, including intracellular signal transduction, transcription, protein transport, DNA repair, and protein degradation . The GAF domain family, linked to the binding of second messengers, is involved in cyclic nucleotide signaling, transcription, phototransduction, and probably many more unidentified processes . Therefore, besides the conserved catalytic domain, the additional domains also play important roles in functions of kinases.
The availability of genome sequences provides a good opportunity for comparative analysis of gene families. Ser/Thr protein kinases have important effect on cyanobacteria living, as well as two component system. A total of 286 putative STK genes have been identified from 21 species of cyanobacteria using BlastP, TblastN and ClustalW, and 19 were not annotated originally as STKs. The number of STKs varied as a function of the genome size, ecophysiology and physiological properties. Fourteen types and 131 additional domains embedded in STK genes may assist their signal recognition functions. A similar catalytic mechanism was inherited with conventional conserved motifs, amino acids, and some uniquely cyanobacterial conserved residues. Gene duplication, loss, shuffling, insertion, and/or horizontal transfer appear to have played important roles during the evolution of cyanobacterial STKs.
Twenty-one species of cyanobacteria, including Prochlorococcus, Synechococcus, Synechocystis, Crocosphaera, Gloeobacter, Trichodesmium, Anabaena and Nostoc were used in this analysis. Since sequences of three Synechococcus strains RS9917, WH5701, WH7805 and one Prochlorococcus marinus strain MIT9312 are not yet fully released, they were not considered in our comparisons. All 21 genome sequences (as of Jan 2007) were accessed from IMG  in FASTA format. Four archaeal STK proteins in Methanopyrus kandleri AV19 were downloaded from the IMG database, and five eukaryotic STK proteins from human and mouse were downloaded from NCBI .
In order to identify genes that may encode STKs, a set of proven cyanobacterial STKs was used to search individual cyanobacterial genomes. The proven cyanobacterial STKs used were alr2502 (Pkn22) , alr4366 (PknA) , all4813 (PknC) , alr4368 (PknD) , alr3732 (PknE)  from Anabaena PCC 7120, and sll1574 (SpkA) , slr1697 (SpkB) , slr0599 (SpkC) , sll0776 (SpkD)  from Synechocystis PCC6803. Sll1574 (SpkA) is interrupted by an insertion element in the sequenced strain, and in this paper "Sll1574" refers to the reconstituted Sll1574-1575. BlastP and TblastN were conducted locally to search all proteins from each of the 21 cyanobacteria, using a threshold e-value of 1e-10. Proteins found by this method that fit the criteria for a genuine STK (see below) were added to the query set for another round of Blast searches. The procedure was continued until no new proteins were found.
Proteins identified by Blast were aligned by using ClustalW program  with a gap opening penalty of 10, a gap extension penalty of 0.2, and Gonnet as the weight matrix. Alignment of query sets were examined by inspection of apparent Motifs I through XI per Hanks and Hunter  (Figure 2). A protein was accepted as an STK if it was possible to recognize the most conserved subdomains (i.e. II, V, VIb, VII, VIII, IX, and XI) and if those conserved amino acid residues known to participate in the function of STKs were present. Some minor deviations were tolerated, considered on a case-by-case basis. The ultimate decision was in the end subjective but in most cases quite clear cut.
Structure analyses of the obtained STKs were performed using the SMART [27, 61] and CDD (Conserved Domains Database) [44, 62] databases, relying on hidden Markov models and Reverse Position-Specific BLAST separately. Sequences of the catalytic domains (about 280 aa in length) and additional domains used for phylogenetic tree construction were obtained from the SMART database and sequences with poor instances of Motif I and XI were checked manually. Trees based on 16s rRNA, STK conserved catalytic domains, GAF domains (mean of 150–200 aa), and HiskA domains (mean of 65–100 aa) were constructed using NJ methods of the MEGA (Version 3.0) package , and the reliability of each branch was tested by 1000 bootstrap replications.
Castenholz RW, Waterbury JB: Preface. Bergey's manual of systematic bacteriology. Williams and Wilkins, Baltimore. Edited by: Staley JT, Bryant MP, Pfennig N, Holt JG. 1989, 3: 1710-1727.
Castenholz R: Species usage, concept, and evolution in the cyanobacteria (blue-green algae). J Phycol. 1992, 28: 737-745. 10.1111/j.0022-3646.1992.00737.x.
Stanier RY, Cohen-Bazire G: Phototrophic prokaryotes: the cyanobacteria. Annu Rev Microbiol. 1977, 31: 225-274. 10.1146/annurev.mi.31.100177.001301.
Partensky F, Hess WR, Vaulot D: Prochlorococcus, a marine photosynthetic prokaryote of global significance. Microbiol Mol Biol Rev. 1999, 63: 106-127.
Parkinson JS: Signal transduction schemes of bacteria. Cell. 1993, 73: 857-871. 10.1016/0092-8674(93)90267-T.
Stock JB, Stock AM, Motten JM: Signal transduction in bacteria. Nature. 1990, 344: 395-400. 10.1038/344395a0.
Ulrich LE, Koonin EV, Zhulin IB: One-component systems dominate signal transduction in prokaryotes. Trends Microbiol. 2005, 13: 52-56. 10.1016/j.tim.2004.12.006.
Munoz-Dorado J, Inouye S, Inouye M: A gene encoding a protein serine/threonine kinase is required for normal development of M. xanthus, a gram-negative bacterium. Cell. 1991, 67: 995-1006. 10.1016/0092-8674(91)90372-6.
Mann NH: Protein phosphorylation in cyanobacteria. Microbiology. 1994, 140 (Pt 12): 3207-3215.
Kennelly PJ, Potts M: Fancy meeting you here: a fresh look at 'prokaryotic' protein phosphorylation. J Bacteriol. 1996, 178: 4759-4764.
Kennelly PJ, Potts M: Life among the primitives: protein O-phosphatases in prokaryotes. Front Biosci. 1999, 4: D372-D385. 10.2741/Kennelly.
Zhang CC: Bacterial signalling involving eukaryotic-type protein kinases. Mol Microbiol. 1996, 20: 9-15. 10.1111/j.1365-2958.1996.tb02483.x.
Hanks SK, Hunter T: The eukaryotic protein kinase superfamily: kinase (catalytic) domain structure and classification. FASEB J. 1995, 9: 576-596.
Taylor SS, Knighton DR, Zheng J, Sowadski JM, Gibbs CS, Zoller MJ: A template for the protein kinase family. Trends Biochem Sci. 1993, 18: 84-89. 10.1016/0968-0004(93)80001-R.
Kennelly PJ: Archaeal protein kinases and protein phosphatases : insights from genomics and biochemistry. Biochem J. 2003, 370: 373-389. 10.1042/BJ20021547.
Wang L, Sun YP, Chen WL, Li JH, Zhang CC: Genomic analysis of protein kinases, protein phosphatases and two-component regulatory systems of the cyanobacterium Anabaena sp. strain PCC 7120. FEMS Microbiol Lett. 2002, 217: 155-165. 10.1111/j.1574-6968.2002.tb11469.x.
Durochera D, Jacksonc SP: The FHA domain. FEBS Lett. 2002, 513: 58-66. 10.1016/S0014-5793(01)03294-X.
Kamei A, Yuasa T, Orikawa K, Geng XX, Ikeuchi M: A eukaryotic-type protein kinase, SpkA, is required for normal motility of the unicellular cyanobacterium Synechocystis sp. strain PCC 6803. J Bacteriol. 2001, 183: 1505-1510. 10.1128/JB.183.5.1505-1510.2001.
Kamei A, Yoshihara S, Yuasa T, Geng X, Ikeuchi M: Biochemical and functional characterization of a eukaryotic-type protein kinase, SpkB, in the cyanobacterium Synechocystis sp. PCC 6803. Curr Microbiol. 2003, 46: 296-301. 10.1007/s00284-002-3887-2.
Zhang CC, Libs L: Cloning and characterisation of the pknD gene encoding a eukaryotic-type protein kinase in the cyanobacterium Anabaena sp. PCC7120. Mol Gen Genet. 1998, 258: 26-33. 10.1007/s004380050703.
Zhao F, Zhang X, Liang C, Wu J, Bao Q, Qin S: Genome-wide Analysis of Restriction-Modification System in Unicellular and Filamentous Cyanobacteria. Physiol Genomics. 2006, 24: 181-190.
Ashby MK, Houmard J: Cyanobacterial Two-Component Proteins: Structure, Diversity, Distribution, and Evolution. Microbiol Mol Biol R. 2006, 70 (2): 472-509. 10.1128/MMBR.00046-05.
Av-Gay Y, Everett M: The eukaryotic-like Ser/Thr protein kinases of Mycobacterium tuberculosis. Trends Microbiol. 2002, 8 (5): 238-244. 10.1016/S0966-842X(00)01734-0.
Petrickova K, Petricek M: Review Eukaryotic-type protein kinases in Streptomyces coelicolor: variations on a common theme. Microbiol. 2003, 149: 1609-1621. 10.1099/mic.0.26275-0.
Krupa A, Srinivasan N: Diversity in domain architectures of Ser/Thr kinases and their homologues in prokaryotes. BMC Genomics. 2005, 6: 129-149. 10.1186/1471-2164-6-129.
Zhang CC, Jiang JC, Sakr S, Wang L: Protein Phosphorylation on Ser, Thr and Tyr Residues In Cyanobacteria. J Mol Microbiol Biotechnol. 2005, 9: 154-166. 10.1159/000089644.
Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: Identification of signaling domains. Proc Natl Acad Sci USA. 1998, 95: 5857-5864. 10.1073/pnas.95.11.5857.
Skorko R: Polyphosphate as a source of phosphoryl group in protein modification in the archaebacterium Sulfolobus acidocaldarius. Biochimie. 1989, 71: 1089-1093. 10.1016/0300-9084(89)90115-6.
Shi L, Potts M, Kennelly PJ: The serine, threonine, and/or tyrosine-specific protein kinases and protein phosphatases of prokaryotic organisms: a family portrait. FEMS Microbiol Rev. 1998, 22: 229-253. 10.1111/j.1574-6976.1998.tb00369.x.
Kamei A, Yuasa T, Geng XX, Ikeuchi M: Biochemical Examination of the Potential Eukaryotic-type Protein Kinase Genes in the Complete Genome of the Unicellular Cyanobacterium Synechocystis sp. PCC 6803. DNA Res. 2002, 9: 71-78. 10.1093/dnares/9.3.71.
Gonzalez L, Phalip V, Zhang CC: Characterization of PknC, a Ser/Thr kinase with broad substrate specificity from the cyanobacterium Anabaena sp. strain PCC 7120. Eur J Biochem. 2001, 268: 1869-1875. 10.1046/j.1432-1327.2001.02057.x.
Zhang CC: A gene encoding a protein related to eukaryotic protein kinases from the filamentous heterocystous cyanobacterium Anabaena PCC 7120. Proc Natl Acad Sci USA. 1993, 90: 11840-11844. 10.1073/pnas.90.24.11840.
Zhang CC, Friry A, Peng L: Molecular and genetic analysis of two closely linked genes that encode, respectively, a protein phosphatase 1/2A/2B homolog and a protein kinase homolog in the cyanobacterium Anabaena sp. Strain PCC 7120. J Bacteriol. 1998, 180: 2616-2622.
Pratt WB: The hsp90-based chaperone system: involvement in signal transduction from a variety of hormone and growth factor receptors. Proc Soc Exp Biol Med. 1998, 217: 420-434.
Ho YJ, Burden LM, Hurley JH: Structure of the GAF domain, a ubiquitous signaling motif and a new class of cyclic GMP receptor. EMBO J. 2000, 19: 5288-5299. 10.1093/emboj/19.20.5288.
Joshi B, Janda L, Stoytcheva Z, Tichy P: PkwA, a WD-repeat protein, is expressed inspore-derived mycelium of Thermomonospora curvata and phosphorylation of its WD domain could act as a molecular switch. Microbiology. 2000, 146: 3259-3267.
Ponting CP, Aravind L, Schultz J, Bork P, Koonin EV: Eukaryotic signalling domain homologues in archaea and bacteria. Ancient ancestry and horizontal gene transfer. J Mol Biol. 1999, 289: 729-745. 10.1006/jmbi.1999.2827.
Larkin RM, Alonso JM, Ecker JR, Chory J: GUN4, a regulator of chlorophyll synthesis and intracellular signaling. Science. 2003, 299: 902-906. 10.1126/science.1079978.
Black K, Buikema WJ, Haselkorn R: The hglK gene is required for localization of heterocyst-specific glycolipids in the cyanobacterium Anabaena sp. strain PCC 7120. J Bacteriol. 1995, 177: 6440-6448.
Xu WL, Jeanjean R, Liu YD, Zhang CC: Pkn22 (alr2502) encoding a putative Ser/Thr kinase in the cyanobacterium Anabaena sp. PCC 7120 is induced by both iron starvation and oxidative stress and regulates the expression of isiA. FEBS Lett. 2003, 553: 179-182. 10.1016/S0014-5793(03)01019-6.
Ohmori M, Ikeuchi M, Sato N, other 15 authors: Characterization of Genes Encoding Multi-domain Proteins in the Genome of the Filamentous Nitrogen-fixing Cyanobacterium Anabaena sp. Strain PCC 7120. DNA Res. 2001, 8: 271-284. 10.1093/dnares/8.6.271.
Phalip V, Li JH, Zhang CC: HstK, a cyanobacterial protein with both a serine/threonine kinase domain and a histidine kinase domain: implication for the mechanism of signal transduction. Biochem J. 2001, 360: 639-644. 10.1042/0264-6021:3600639.
Cheng Y, Li JH, Shi L, Wang L, Latifi A, Zhang CC: A Pairof Iron-Responsive Genes Encoding Protein Kinases with a Ser/Thr Kinase Domain and a His Kinase Domain Are Regulated by NtcA in the Cyanobacterium Anabaena sp. Strain PCC 7120. J Bacteriol. 2006, 6: 4822-4829. 10.1128/JB.00258-06.
Marchler BA, Bryant SH: CD-Search: protein domain annotations on the fly. Nucleic Acids Res. 2004, 32 (W): 327-331. 10.1093/nar/gkh454.
Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci USA. 1990, 87: 4576-4579. 10.1073/pnas.87.12.4576.
Ashby MK: Survey of the number of two-component response regulator genes in the complete and annotated genome sequences of prokaryotes. FEMS Microbiol Lett. 2004, 231: 277-281. 10.1016/S0378-1097(04)00004-7.
Galperin MY: A census of membrane-bound and intracellular signal transduction proteins in bacteria: bacterial IQ, extroverts and introverts. BMC Microbiol. 2005, 5: 35-55. 10.1186/1471-2180-5-35.
Dufresne A, Salanoubat M, Partensky F, other 18 authors: Genome sequence of the cyanobacterium Prochlorococcus marinus SS120, a nearly minimal oxyphototrophic genome. PNAS. 2003, 100 (17): 10020-10025. 10.1073/pnas.1733211100.
Dufresne A, Garczarek L, Partensky F: Accelerated evolution associated with genome reduction in a free-living prokaryote. Genome Biol. 2005, 6: R14-10.1186/gb-2005-6-2-r14.
Han GN, Zhang CC: On the origin of Ser/Thr kinases in a prokaryote. FEMS Microbiol Let. 2001, 200: 79-84. 10.1111/j.1574-6968.2001.tb10696.x.
Leonard CJ, Aravind L, Koonin EV: Novel families of putative protein kinases in bacteria and archaea: evolution of the "Eukaryotic" protein kinase superfamily. Genome Res. 1998, 8: 1038-1047.
Aravind L, Ponting CP: The GAF domain: an evolutionary link between diverse phototransducing proteins. Trends Biochem Sci. 1997, 22: 458-459. 10.1016/S0968-0004(97)01148-1.
Galperin MY, Nikolskaya AN, Koonin EV: Novel domains of the prokaryotic two-component signal transduction systems. FEMS Microbiol Lett. 2001, 203: 11-21. 10.1111/j.1574-6968.2001.tb10814.x.
Hofmann K, Bucher P: The FHA domain: a putative nuclear signalling domain found in protein kinases and transcription factors. Trends Biochem Sci. 1995, 20: 347-349. 10.1016/S0968-0004(00)89072-6.
Ponting CP, Aravind L: PAS: a multifunctional domain family comes to light. Curr Biol. 1997, 7: R674-677. 10.1016/S0960-9822(06)00352-6.
Zhang WW, Shi L: Evolution of the PPM-family protein phosphatases in Streptomyces: duplication of catalytic domain and lateral recruitment of additional sensory domains. Microbiol. 2004, 150: 4189-4197. 10.1099/mic.0.27480-0.
Zhulin IB, Nikolskaya AN, Galperin MY: Common Extracellular Sensory Domains in Transmembrane Receptors for Diverse Signal Transduction Pathways in Bacteria and Archaea. J Bacteriol. 2003, 185: 285-294. 10.1128/JB.185.1.285-294.2003.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
IMG Home. [http://img.jgi.doe.gov/cgi-bin/pub/main.cgi]
NCBI Home. [http://www.ncbi.nlm.nih.gov]
SMART Home. [http://smart.embl-heidelberg.de/]
NCBI Conserved Domain Database (CDD). [http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml]
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.
This work was supported by the Key Innovative Project of Chinese Academy of Science (KZCX2-YW-209, KZCX3-SW-223), Hi-Tech Research and Development Program (2006AA090303) of China, and the CAS/SAFEA International Partnership Program for Creative Research Teams (Research and Applications of Marine Functional Genomics).
XWZ, FQZ and SQ devised the overall strategy for these studies. XWZ performed all database searches, acquired the sequence data, and prepared all figures and tables. XWZ and FQZ performed the conserved domain, structure, alignments, and phylogenetic analyses. XWZ, FQZ, YY, CWL and XYG jointly wrote the paper, and all authors have read and accepted the final version of the manuscript.
Xiaowen Zhang, Fangqing Zhao contributed equally to this work.
About this article
Cite this article
Zhang, X., Zhao, F., Guan, X. et al. Genome-wide survey of putative Serine/Threonine protein kinases in cyanobacteria. BMC Genomics 8, 395 (2007). https://doi.org/10.1186/1471-2164-8-395
- Signal Transduction System
- Additional Domain
- Unicellular Cyanobacterium
- Cyanobacterial Genome
- Domain Shuffling