Skip to main content
  • Research article
  • Open access
  • Published:

Comparative genomic analysis of two-component regulatory proteins in Pseudomonas syringae

Abstract

Background

Pseudomonas syringae is a widespread bacterial plant pathogen, and strains of P. syringae may be assigned to different pathovars based on host specificity among different plant species. The genomes of P. syringae pv. syringae (Psy) B728a, pv. tomato (Pto) DC3000 and pv. phaseolicola (Pph) 1448A have been recently sequenced providing a major resource for comparative genomic analysis. A mechanism commonly found in bacteria for signal transduction is the two-component system (TCS), which typically consists of a sensor histidine kinase (HK) and a response regulator (RR). P. syringae requires a complex array of TCS proteins to cope with diverse plant hosts, host responses, and environmental conditions.

Results

Based on the genomic data, pattern searches with Hidden Markov Model (HMM) profiles have been used to identify putative HKs and RRs. The genomes of Psy B728a, Pto DC3000 and Pph 1448A were found to contain a large number of genes encoding TCS proteins, and a core of complete TCS proteins were shared between these genomes: 30 putative TCS clusters, 11 orphan HKs, 33 orphan RRs, and 16 hybrid HKs. A close analysis of the distribution of genes encoding TCS proteins revealed important differences in TCS proteins among the three P. syringae pathovars.

Conclusion

In this article we present a thorough analysis of the identification and distribution of TCS proteins among the sequenced genomes of P. syringae. We have identified differences in TCS proteins among the three P. syringae pathovars that may contribute to their diverse host ranges and association with plant hosts. The identification and analysis of the repertoire of TCS proteins in the genomes of P. syringae pathovars constitute a basis for future functional genomic studies of the signal transduction pathways in this important bacterial phytopathogen.

Background

Bacterial signal transduction pathways sense the cellular external environment and regulate cellular functions in response to environmental signals. A mechanism commonly found in bacteria for signal transduction is the two-component system (TCS). Bacterial TCSs are common components of complex regulatory networks and cascades, often associated with global regulation as well as with regulation of virulence. TCS genes are typically located within the same operon encoding two signalling proteins: a transmembrane sensor histidine kinase (HK) and a cytoplasmic response regulator (RR), which may sometimes be carried by a single polypeptide to form the hybrid HKs [1]. The mechanism of signal transduction by TCS proteins is based on phosphotransfer reactions between histidine (H) and aspartate (D) residues in highly conserved signalling domains of the HKs and their cognate RRs. TCS proteins have a modular organization, which may give rise to highly complex structures, but the core structures and activities are maintained [2]. HKs are typically organized as homodimers with two functionally and structurally distinct domains: a highly variable N-terminal extracytoplasmic sensory domain, and a more conserved C-terminal cytoplasmic transmitter domain, also known as the dimerization/phosphoacceptor domain [2, 3]. The sensor domain varies in length and amino acid sequence from one HK to another, conferring specificity for different environmental stimuli. In most HKs, the transmitter domain shows high sequence conservation, especially within a set of six recognizable motifs or boxes designated H, N, F, G1, G2, and G3. In particular, the H box contains an invariant H residue that is autophosphorylated in an ATP-dependent manner [4]. In contrast, CheA-like HKs that function in chemotaxis lack the sensor domain and differ from other HKs in their domain constitution and organization, where the H box of the transmitter domain resides at the N-terminal end of the protein [58]. LytS-like HKs also differ significantly in their domain architecture from other HKs [911]. RRs generally contain at least two functional domains: a conserved N-terminal receiver domain (REC domain) that is phosphorylated by the HK at a strictly conserved D residue, and one or more variable C-terminal output domains [12]. Modulation of the phosphorylated state of the RR controls either expression of the target genes or cellular behaviour. The principal type of bacterial RRs are transcription factors that regulate gene-expression with DNA-binding helix-turn-helix (HTH) output domains [1, 3, 12, 13]. Hybrid HKs contain both a HK transmitter domain and a REC domain in a single large polypeptide, and are characterized by multi-step phosphotransfer reactions [1, 7, 14].

The availability of complete genome sequences for a continually growing number of bacteria has allowed the definitive assessment that TCS proteins are present in almost all bacterial species [1, 8, 12]. Genomic analyses demonstrate the enormous impact of TCSs on environmental adaptation of bacteria, and reveal a wide variation of HK and RR numbers between different bacterial species [7, 8, 12, 1520].

The bacterial plant pathogen Pseudomonas syringae causes disease on a variety of plant species, and strains of P. syringae have been classified into different pathovars depending on their host range among different plant species [21]. Infection of host plants by P. syringae involves growth on leaf surfaces as an epiphyte, that enters plant leaves through stomata, multiplies to large populations in the apoplast and produces disease symptoms [21, 22]. P. syringae injects effector proteins into the cytoplasm of plant cells by means of the Hrp type III secretion system [21]. Genome comparisons indicate that P. syringae is significantly different from other Pseudomonas species [23, 24], suggesting that in the adaptation to the phytopathogenic lifestyle its genome must have undergone fundamental changes without a reduction in size. The complete genomic sequences of three economically important pathovars of this plant pathogenic bacteria have been determined: P. syringae pv. tomato (Pto) DC3000, pv. syringae (Psy) B728a and pv. phaseolicola (Pph) 1448A [2527]. In these genomes, over 10 to 12 % of the genes are dedicated to regulation, which may reflect the need for rapid adaptation to the diverse environments encountered during epiphytic growth, plant colonization and pathogenesis [2527]. Genome analyses of these P. syringae pathovars revealed fewer extracytoplasmic function (ECF) sigma factors (10 ECF sigma factors) than in related Pseudomonas with different lifestyles [24]. Recently, analysis of the Pto DC3000 genome sequence allowed the identification of 69 HKs [28, 29] and 71 RRs, 21 of which were hybrid HKs [12]. In a different study not including CheA-like HKs, 64 HKs were identified in Pto DC3000, 20 of which were hybrid HKs [30]. Hence, P. syringae requires a complex array of TCS proteins to cope with diverse plant hosts, host responses, and environmental conditions. The availability of complete genomic sequences of three different P. syringae pathovars makes it possible to conduct this comparative genomic study to identify and analyse the TCS proteins of P. syringae.

Results and Discussion

Distribution of TCS proteins in P. syringae

The putative HKs and RRs in Psy B728a, Pto DC3000 and Pph 1448A were identified by searching the complete genome sequences for proteins containing HK and RR domains using Pfam HMM profiles. Four CheA-like HKs in each P. syringae genome were identified in BLASTP searches using as template the CheA HK of E. coli [31] (Table 1). In addition, BLASTP searches of the HKs and RRs found in each P. syringae pathovar against the genomes of the other two pathovars allowed the identification of additional HKs and RRs. The genomes of P. syringae pathovars were found to contain large numbers of genes encoding TCS proteins: 68 HKs and 93 RRs in Psy B728a, 69 HKs and 95 RRs in Pto DC3000, and 70 HKs and 92 RRs in Pph 1448A (Table 1; see Additional File 1 and 2). The number of genes encoding hybrid HKs (REC-HKs) was 20 in Psy B728a, 22 in Pto DC3000 and 24 in Pph 1448A (Tables 1 and 4). The HMM search method used in this work retrieved hybrid HKs as well as RRs (Table 1). No TCS proteins were identified on any of the plasmids of Pto DC3000 and Pph 1448A. In recent studies, similar numbers of TCS proteins for Pto DC3000 have been reported: 69 HKs [28, 29] and 71 RRs, 21 of which were hybrid HKs [12]; or 64 HKs in a study not including CheA-like HKs, 20 of which where hybrid HKs [30]. Although the number of ECF sigma factors in all three P. syringae genomes (10 ECF sigma factors) is only about half that found in other Pseudomonas species [24, 32], the number of TCS proteins is close to that found in other Pseudomonas genomes [33].

Table 1 Distribution of HKs and RRs found in the genomes of P. syringae pv. syringae B728a, pv. tomato DC3000 and pv. phaseolicola 1448A.
Table 2 Putative TCS gene clusters in the genomes of P. syringae pv. syringae B728a, pv. tomato DC3000 and pv. phaseolicola 1448A.
Table 3 Orphan HK and RR genes in the genomes of P. syringae pv. syringae B728a, pv. tomato DC3000 and pv. phaseolicola 1448A.
Table 4 Hybrid HK genes in the genomes of P. syringae pv. syringae B728a, pv. tomato DC3000 and pv. phaseolicola 1448A.

HK and RR genes were scattered over the entire chromosomes of the three P. syringae pathovars. Conservation of the genetic organization between HK and RR genes was analysed in the genomes of Psy B728a, Pto DC3000 and Pph 1448A allowing the identification of gene clusters containing HKs and RRs that constitute putative TCSs (Table 2). Like in other bacterial species, many P. syringae HKs and RRs were encoded by clusters of adjacent genes: 37 putative clusters of complete TCS genes in Psy B728a, 34 in Pto DC3000, and 33 in Pph 1448A (Table 2). For the remaining HK or RR genes, their partner genes could not be predicted from genetic organization and, therefore, they were considered as orphan HKs or RRs. The orphan HKs were 11 in each P. syringae genome, and the number of genes encoding orphan RRs was very high: 36 in Psy B728a, 38 in Pto DC3000 and 35 in Pph 1448A (Table 3). Finally, the comparative genomic analysis allowed the identification of a core of complete TCS protein orthologues among the three P. syringae pathovars, that is composed by 30 putative TCS clusters (HK and RR) (Table 2), 11 orphan HKs, 33 orphan RRs (Table 3), and 16 hybrid HKs (Table 4).

Classification of HKs

HKs have been classified on the basis of phylogenetic analyses and the sequence relationships of the residues surrounding the H-box [7, 8, 17, 34]. Furthermore, several new domains with putative biological functions have been described in HKs, and domain architecture has proven particularly informative for analysing multi-domain proteins involved in signal transduction [2, 11, 12, 35]. The phylogenetic analysis and examination of the region around the H box of P. syringae HKs showed that three of the five major HK types found in E. coli [8] were present in P. syringae: Type I (IA, IB, IC), III, and CheA-like HKs (Table 1; see Additional File 1). In contrast, Type II and IV HKs were totally absent from P. syringae. However, the LytS-like HK FimS/AlgZ and HKs containing GAF domains did not cluster within any of the defined HK types of E. coli [8], and formed two separate HK groups: LytS-like HKs and GAF-HKs. GAF sensor domains are commonly found cytoplasmic signalling domains in the N-terminal region of HKs [2, 34], and appear to act as binding sites for small ligands, such as cyclic nucleotides (cAMP and cGMP) and small molecules, which modulate the catalytic activity of the target protein [36, 37]. In addition, analysis of domain architecture of P. syringae HKs showed a conserved core structure for each HK type in P. syringae (Figure 1). The conserved core of Type III HKs and LytS-like HKs only had a HK-like ATPase (HATPase_c) catalytic domain and a His_kinase domain, respectively. The conserved core of CheA-like HKs contained a C-terminal CheA regulatory domain but lacked the HisKA domain. The conserved core of Type I HKs and GAF-HKs had a central region with HisKA and HATPase_c domains fused to additional domains on the N-terminal end: a HAMP domain in Type IA, a PAS domain in Type IC, and GAF plus phytochrome (PHY) binding domains in GAF-HKs (Figure 1).

Figure 1
figure 1

Schematic representation of the conserved core structures found in P. syringae HK types. The domains are not drawn to scale. HAMP, domain found in H Ks, A denylyl cyclases, M ethyl binding proteins and P hosphatases (PF00672); HisKA, HK dimerization/phosphoaceptor domain (PF00512); HATPase_c, HK-type ATPase catalytic domain (PF02518); REC, receiver domain (PF00072); PAS, signal sensor domain (PF00989); HPt, H istidine-containing P hosphot ransfer domain (PF01627); H-kinase_dim, HK homodimeric domain (PF02895); GAF, signal sensor domain (PF01590); PHY, phytochrome domain (PF00360); His_kinase, region within bacterial HKs (PF06850).

Orphan HKs fell into two HK types: Type I (IA, IB and IC), and GAF-HKs (Table 3); and hybrid HKs of P. syringae belong either to the Type I (IB and IC) or CheA-like HKs (Table 4). PSYR3504 (BphP1) and PSYR2385 (BphP2) HKs have been previously described as bacteriophytochromes (BphPs) that belong to the HWE_HK family [4, 38]. Similar to other BphPs, the bphP1 (PSYR3504) gene of P. syringae pathovars is located in an operon downstream from a bphO gene, encoding a putative heme oxigenase.

Classification of RRs

RRs show a great variety of output domains and domain combinations. Recently, bacterial and archaeal RRs have been classified into families based in their domain architectures [12]. RRs typically consist of an N-terminal REC domain fused to a C-terminal HTH DNA-binding output domain (OmpR, NarL, NtrC, LytR, AraC, Spo0A, Fis, YcbB, RpoE, and MerR) that activates or represses transcription of specific target genes [2, 12]. In addition, prokaryotic genomes encode a variety of RRs with unusual domain organization: RRs with enzymatic output domains (GGDEF, EAL, HD-GYP, CheB, CheC, PP2C, and HisKA), RRs with RNA-binding output domains (ANTAR and CsrA), RRs with protein- or ligand-binding output domains (CheW, PAS, GAF, TPR, CAP_ED, and Hpt), RRs with the REC domain as a stand-alone module, and RRs with domains of unknown function [12]. The RRs identified from the genomes of P. syringae pathovars were assigned to these different RR families [12] according to the domain architecture and phylogenetic analysis (Table 1; see Additional File 2).

Bacterial RRs without a REC domain are extremely rare, but a number of enhancer-binding proteins (EBPs) lack the REC domain and normally function as RRs [39]. EBPs are involved in the activation of the bacterial transcription by interaction with the sigma-54 RNA polymerase holoenzyme [40]. In P. syringae, the HrpR and HrpS proteins show a high sequence similarity to the NtrC family of transcriptional RRs and have been previously identified as unusual EBPs lacking the N-terminal REC domain; however, similar to other EBPs, they retain the domain that interacts with the sigma-54 RNA polymerase holoenzyme plus the C-terminal DNA-binding domain [3942]. In addition, the NarL-like RR CorP of Pto DC3000 that is involved in the regulation of coronatine biosynthesis [43, 44] also lacks the REC domain. Thus, HrpR, HrpS and CorP proteins were not identified during the search of RRs in P. syringae genomes with the HMM profile that targets the RR REC domain, nevertheless these proteins were considered orphan RRs (Table 3).

Differences in TCS genes among pathovars that may contribute to plant host specificity

A close analysis of the distribution of genes encoding TCS proteins revealed that there are important differences in TCS proteins among the three pathovars of P. syringae that may contribute to their diverse host ranges and association with particular host plants. A number of the identified TCS genes were unique to each P. syringae pathovar without counterparts in the other two pathovars. The corRSP regulatory region (PSPTO4704-4706) of coronatine biosynthesis and the copRS TCS (PSYR1497/1498) regulating copper resistance were only present in Pto DC3000 and Psy B728a, respectively. Other TCS genes unique to each P. syringae pathovar were: PSYR2114, PSYR2939, PSYR2940 and PSYR3084 in Psy B728a; PSPTO0785/0786, PSPTO2329, PSPTO4079, PSPTO4080 and PSPTO5573/5574 in Pto DC3000; PSPPH0770, PSPPH0778, PSPPH0944 and PSPPH1261 in Pph 1448A. The unique hybrid HKs PSPPH0770 and PSPPH0944 were flanked by transposases. However, the unique RRs PSPTO2329 and PSPTO5574 were disrupted by transposases [25, 27], and it is unlikely that these genes encode functional products. Finally, 11 TCS proteins were only shared between two of these P. syringae pathovars.

Variations among P. syringae pathovars were also produced by the insertion of mobile genetic elements or point mutations in TCS genes resulting in disrupted reading frames. PSPTO2326 and PSPPH2083 encoded truncated hybrid HKs by comparison with the length of their orthologue PSYR2113 (Table 4) that is located next to the unique RR PSYR2114. PSPTO2326 and PSPPH2083 were located adjacent to a transposase and to a site-specific recombinase, respectively. Probably these elements caused the disrupted hybrid HKs and the lack of PSYR2114 orthologues in Pto DC3000 and Pph 1448A. Similarly, PSPTO2983 (baeS2) and PSPPH2510 encoded truncated HKs compared to the length of their P. syringae orthologues, and PSPPH2980 was interrupted by an ISPsy18 transposase. PSPTO2983, PSPPH2510 and PSPPH2980 HKs were unpaired without a RR gene in its vicinity, whereas their P. syringae orthologues are located on TCS gene clusters with adjacent RRs (Table 2).

Although the PSPPH1362 gene was disrupted by an authentic frameshift, Psy B728a (PSYR1292) and Pto DC3000 (PSPTO1482) orthologues encoded intact hybrid HKs with similarity to BvgS of Bordetella species that controls the regulation of many virulence factors [45]. In each pathovar, these hybrid HK genes were adjacent to orphan RR genes transcribed in the same direction (PSYR1293, PSPTO1482 and PSPPH1363), and their encoded proteins exhibited significant homology to the PvrR RR of P. aeruginosa PA14 which controls antibiotic susceptibility and biofilm formation [46], and to the virulence related protein VieA of Vibrio cholerae [47].

Conclusion

In this article we present a thorough analysis of the identification and distribution of TCS proteins among the sequenced genomes of P. syringae. A large set of TCS proteins is required for the capacity of P. syringae to detect and adapt to changing environments during plant association and pathogenesis. Moreover, P. syringae has been isolated from non-plant environments such as river epilithon (rock-attached biofilms) [48] in which TCS proteins may have also important regulatory roles. P. syringae pathovars posses between 68–70 HKs and 92–95 RRs (Table 1), however there is little information describing their regulatory functions and the major part of these TCS proteins is uncharacterized. Many of the TCS proteins investigated so far in P. syringae have been shown to be involved in plant pathogenicity and association with host plants. The orphan RRs HrpR and HrpS are involved in a complex regulatory cascade that activates the transcription of the Hrp type III secretion genes and all known effector genes [42, 49]. Expression of the type III secretion genes and effector genes is also regulated by the particular TCS GacA/GacS [50] and the RhpRS system [51]. Furthermore, the GacA/GacS system controls the expression of a variety of virulence factors, including protease and syringomycin biosynthesis [52]. The TCS CopRS and the modified CorRSP system regulate resistance to copper [53] and coronatine synthesis [43, 44], respectively. Finally, the hybrid HK PSPTO2896 contains an N-terminal LOV (light, oxygen, or voltage) domain and is blue-light-activated [54].

Bacteria with large genomes are disproportionately enriched in regulatory proteins involved in transcription control and signal transduction compared to medium and small-size genomes, and typically have complex regulatory networks relative to bacteria with smaller genomes [5557]. The existence of large numbers of HKs and RRs in P. syringae strongly suggests that TCS proteins play important regulatory roles in the adaptation of this bacterium to different plant and non-plant environments. Comparative genomics of closely related species of pathogenic bacteria represents a powerful tool for the identification of genes potentially involved in host specificity and pathogenesis. The availability of the genome sequences of Pto DC3000, Psy B728a and Pph 1448A provides us with the unique capability of comparing the complement of TCS proteins in these P. syringae pathovars that differ in host range and other interactions with plants. This comparative genomic analysis reveals a core of orthologues and important differences in TCS genes between P. syringae pathovars. It is especially worth noting the high number of genes encoding orphan HKs and RRs in these genomes. Moreover, differences in the repertoires of TCS proteins are likely to facilitate the adaptation of P. syringae pathovars to different plant hosts and/or could be responsible for the different disease characteristics induced. Consequently, the TCS proteins unique to each P. syringae pathovar are interesting targets for future investigations to identify TCS proteins involved in the different host ranges and/or plant pathogenesis. However, the challenge remains to associate these differences in TCS proteins to specific traits of P. syringae pathovars. Additionally, pathovar-specific differences in gene content might be used to design targeted approaches for disease control and could allow the precise PCR-based diagnosis of bacterial diseases [58].

Analysis of the regulatory functions, molecular mechanisms and signal transduction pathways of TCS proteins should contribute to the understanding of the complex events that occur in P. syringae during pathogenesis and adaptation to different plant hosts and different non-plant environments. Rapid progress in the study of TCS proteins is being made by the combination of molecular genetic approaches with genome-scale analysis [59]. Genetic and biochemical studies are necessary to further explore the signal transduction pathways mediated by some of these TCS proteins at the molecular level: construction and analysis of deletion mutants in TCS genes in order to determine the signals sensed by the HK and the targets for the RR of each system. In addition, the application of more extensive analysis with global methods, such as DNA microarray studies reported for B. subtilis [60] and S. pneumoniae [61], might allow defining the regulons and the potential regulatory functions of TCS proteins in response to environmental signals. Furthermore, unravelling these signal transduction pathways could potentially lead to the design of innovative strategies to control P. syringae. In conclusion, this comparative genomic analysis constitutes a basis for future functional genomic analysis of P. syringae to establish which TCS proteins participate in the pathogenesis and the adaptation to different plant and non-plant environments.

Methods

Identification of TCS proteins in P. syringae genomes

The identification of HKs and RRs is based on the computational domain analysis of protein sequences. The approach used to identify putative HKs and RRs from the complete genome sequences of Psy B728a, Pto DC3000 and Pph 1448A was similar to that described previously [33] with slight modification. Briefly, five different HMM profiles (accession numbers PF00512, PF07568, PF07730, PF07536 and PF06580) were found in Pfam database that target different families of HKs (HisKA, HisKA_2, HisKA_3, HWE_HK and His_kinase). The HWE_HK domain is defined by the absence of a recognizable F box, and the presence of a highly conserved H residue and a WxE motif within the N and G1 boxes of the C-terminal transmitter domain, respectively [4]. These five different HMM profiles were used to recognize the different HKs in the P. syringae genomes, and hits with a E-value below a selected cut-off (10-6) were extracted. A profile HMM downloaded from Pfam protein families database [62], which targets the RR REC domain (accession number PF00072), was used to recognize the RRs in each P. syringae genome. Hits with an E-value below a selected cut-off (10-12) were extracted. Additionally, the CheA HK of Escherichia coli [31] was used as template in BLASTP searches to identify CheA-like HKs in the P. syringae genomes and hits with an E-value below a selected cut-off (10-10) were extracted. Hybrid HKs (REC-HKs) were determined by the presence of complete HK transmitter and REC domains in a single protein. Detection of orthologues of the identified HKs and RRs between the genomes of Psy B728a, Pto DC3000 and Pph 1448A was determined by BLASTP [63] based on the reciprocal best hits of each P. syringae genome against each other genome, completed by the phylogenetic analyses. Finally, functional domains of the HKs and RRs were identified by search the Conserved Domain Databases (CDD) with Reverse Specific Position BLAST [64].

Sequence alignment and phylogenetic analysis

Multiple sequence alignments and phylogenetic trees of HKs and RRs were constructed using the ClustalW program [65], and aligned sequences were imported into the MEGA 3.1 program [66] where phylogenetic trees were inferred. Default parameters were used. Phylogenetic trees were subdivided into groups of orthologues, and co-clustering with members of specific TCS proteins allowed a definitive assignation to a given HK type or RR family.

Abbreviations

TCS:

two-component system

HK:

histidine kinase

RR:

response regulator

HMM:

Hidden Markov Model

HTH:

helix-turn-helix

ECF:

extracytoplasmic function

EBP:

enhancer-binding protein

Psy :

pv. syringae Pseudomonas syringae

Pto :

pv. tomato P. syringae

Pph :

pv. phaseolicola P. syringae

REC:

receiver

PHY:

phytochrome

LOV:

light, oxygen, and voltage

References

  1. Stock AM, Robinson VL, Goudreau PN: Two-component signal transduction. Annu Rev Biochem. 2000, 69: 183-215. 10.1146/annurev.biochem.69.1.183.

    Article  CAS  PubMed  Google Scholar 

  2. Galperin MY: Bacterial signal transduction network in a genomic perspective. Environ Microbiol. 2004, 6: 552-567. 10.1111/j.1462-2920.2004.00633.x.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Hoch JA: Two-component and phosphorelay signal transduction. Curr Opin Microbiol. 2000, 3: 165-170. 10.1016/S1369-5274(00)00070-9.

    Article  CAS  PubMed  Google Scholar 

  4. Karniol B, Vierstra RD: The HWE histidine kinases, a new family of bacterial two-component sensor kinases with potentially diverse roles in environmental signalling. J Bacteriol. 2004, 186: 445-453. 10.1128/JB.186.2.445-453.2004.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Garzon A, Parkinson JS: Chemotactic signalling by the P1 phosphorylation domain liberated from the CheA histidine kinase of Escherichia coli. J Bacteriol. 1996, 178: 6752-6758.

    CAS  PubMed Central  PubMed  Google Scholar 

  6. Bilwes AM, Alex LA, Crane BR, Simon MI: Structure of CheA, a signal transducing histidine kinase. Cell. 1999, 96: 131-141. 10.1016/S0092-8674(00)80966-6.

    Article  CAS  PubMed  Google Scholar 

  7. Grebe TW, Stock JB: The histidine protein kinase superfamily. Adv Microb Physiol. 1999, 41: 139-227.

    Article  CAS  PubMed  Google Scholar 

  8. Kim D, Forst S: Genomic analysis of the histidine kinase family in bacteria and archaea. Microbiology. 2001, 147 (Pt 5): 1197-1212.

    Article  CAS  PubMed  Google Scholar 

  9. Witchurch CB, Alm RA, Mattick JS: The alginate regulator AlgR and an associated FimS are required for twitching motility in Pseudomonas aeruginosa. Proc Natl Acad Sci USA. 1996, 93: 9839-9843. 10.1073/pnas.93.18.9839.

    Article  Google Scholar 

  10. Anantharaman V, Aravind L: Application of comparative genomics in the identification and analysis of novel families of membrane-associated receptors in bacteria. BMC Genomics. 2003, 4: 34-10.1186/1471-2164-4-34.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Mascher T, Helmann JD, Unden G: Stimulus perception in bacterial signal-transducing histidine kinases. Microbiol Mol Biol Rev. 2006, 70: 910-938. 10.1128/MMBR.00020-06.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. Galperin MY: Structural classification of bacterial response regulators: diversity of output domains and domain combinations. J Bacteriol. 2006, 188: 4169-4182. 10.1128/JB.01887-05.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. 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.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Appleby JL, Parkinson JS, Bourret RB: Signal transduction via the multistep phosphorelay: not necessarily a road less traveled. Cell. 1996, 86: 845-848. 10.1016/S0092-8674(00)80158-0.

    Article  CAS  PubMed  Google Scholar 

  15. Mizuno T, Kaneko T, Tabata S: Compilation of all genes encoding bacterial two-component signal transducers in the genome of the cyanobacterium, Synechocystis sp. Strain PCC 6803. DNA Res. 1996, 3: 407-414. 10.1093/dnares/3.6.407.

    Article  CAS  PubMed  Google Scholar 

  16. Mizuno T: Compilation of all genes encoding two-component phosphotransfer signal transducers in the genome of Escherichia coli. DNA Res. 1997, 4: 161-168. 10.1093/dnares/4.2.161.

    Article  CAS  PubMed  Google Scholar 

  17. Fabret C, Feher VA, Hoch JA: Two-component signal transduction in Bacillus subtilis: how one organims sees its world. J Bacteriol. 1999, 181: 1975-1983.

    CAS  PubMed Central  PubMed  Google Scholar 

  18. 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 (2): 277-281. 10.1016/S0378-1097(04)00004-7.

    Article  CAS  PubMed  Google Scholar 

  19. Ashby MK, Houmard J: Cyanobacterial two-component proteins: structure, diversity, distribution, and evolution. Microbiol Mol Biol Rev. 2006, 70: 472-509. 10.1128/MMBR.00046-05.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. de Been M, Francke C, Moezelaar R, Abee T, Siezen RJ: Comparative analysis of two-component signal transduction systems of Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis. Microbiology UK. 2006, 152: 3035-3048. 10.1099/mic.0.29137-0.

    Article  CAS  Google Scholar 

  21. Hirano SS, Upper CD: Bacteria in the leaf ecosystem with emphasis on Pseudomonas syringae-a pathogen, ice nucleus and epiphyte. Microbiol Mol Biol Rev. 2000, 64: 624-653. 10.1128/MMBR.64.3.624-653.2000.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Boch J, Joardar V, Gao L, Robertson TL, Lim M, Kunkel BN: Identification of Pseudomonas syringae pv. tomato genes induced during infection of Arabidopsis thaliana. Mol Microbiol. 2002, 44: 73-88. 10.1046/j.1365-2958.2002.02877.x.

    Article  CAS  PubMed  Google Scholar 

  23. Jensen LJ, Skovgaard M, Sicheritz-Pontén T, Hansen NT, Johansson H, Jørgensen MK, Kiil K, Hallin PF, Ussery DW: Comparative genomics of four Pseudomonas species. The Pseudomonas. Genomics, Life Style and Molecular Architecture. Edited by: Ramos JL. 2004, New York: Kluwer Academic/Plenum Publishers, 1: 139-164.

    Google Scholar 

  24. Oguiza JA, Kiil K, Ussery DW: Extracytoplasmic function sigma factors in Pseudomonas syringae. Trends Microbiol. 2005, 13: 565-568.

    Article  CAS  PubMed  Google Scholar 

  25. Buell CR, Joardar V, Lindeberg M, Selengut J, Paulsen IT, Gwinn ML, Dodson RJ, Deboy RT, Durkin AS, Kolonay JF, Madupu R, Daugherty S, Brinkac L, Beanan MJ, Haft DH, Nelson WC, Davidsen T, Zafar N, Zhou L, Liu J, Yuan Q, Khouri H, Fedorova N, Tran B, Russell D, Berry K, Utterback T, Van Aken SE, Feldblyum TV, D'Ascenzo M, Deng WL, Ramos AR, Alfano JR, Cartinhour S, Chatterjee AK, Delaney TP, Lazarowitz SG, Martin GB, Schneider DJ, Tang X, Bender CL, White O, Fraser CM, Collmer A: The complete genome sequence of the Arabidopsis and tomato pathogen Pseudomonas syringae pv. tomato DC3000. Proc Natl Acad Sci USA. 2003, 100: 10181-10186. 10.1073/pnas.1731982100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Feil H, Feil WS, Chain P, Larimer F, DiBartolo G, Copeland A, Lykidis A, Trong S, Nolan M, Goltsman E, Thiel J, Malfatti S, Loper JE, Lapidus A, Detter JC, Land M, Richardson PM, Kyrpides NC, Ivanova N, Lindow SE: Comparison of the complete genome sequences of Pseudomonas syringae pv. syringae B728a and pv. tomato DC3000. Proc Natl Acad Sci USA. 2005, 102: 11064-11069. 10.1073/pnas.0504930102.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Joardar V, Lindeberg M, Jackson RW, Selengut J, Dodson R, Brinkac LM, Daugherty SC, Deboy R, Durkin AS, Giglio MG, Madupu R, Nelson WC, Rosovitz MJ, Sullivan S, Crabtree J, Creasy T, Davidsen T, Haft DH, Zafar N, Zhou L, Halpin R, Holley T, Khouri H, Feldblyum T, White O, Fraser CM, Chatterjee AK, Cartinhour S, Schneider DJ, Mansfield J, Collmer A, Buell CR: Whole-genome sequence analysis of Pseudomonas syringae pv. phaseolicola 1448A reveals divergence among pathovars in genes involved in virulence and transposition. J Bacteriol. 2005, 187: 6488-6498. 10.1128/JB.187.18.6488-6498.2005.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Galperin MY: A census of membrane-bound and intracellular signal transduction proteins in bacteria: bacterial IQ, extroverts and introverts. BMC Microbiol. 2005, 5: 35-10.1186/1471-2180-5-35.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Alm E, Huang K, Arkin A: The evolution of two-component systems in bacteria reveals different strategies for niche adaptation. PLoS Comput Biol. 2006, 2: e143-10.1371/journal.pcbi.0020143.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Zhang W, Shi L: Distribution and evolution of multiple-step phosphorelay in prokaryotes: lateral domain recruitment involved in the formation of hybrid-type histidine kinases. Microbiology UK. 2005, 151: 2159-2173. 10.1099/mic.0.27987-0.

    Article  CAS  Google Scholar 

  31. Kofoid EC, Parkinson JS: Tandem translation in the cheA locus of Escherichia coli. J Bacteriol. 1991, 173: 2116-2119.

    CAS  PubMed Central  PubMed  Google Scholar 

  32. Kill K, Binnewies TT, Sicheritz-Pontén T, Willenbrock H, Hallin PF, Wassenaar TM, Ussery DW: Genome update: sigma factors in 240 bacterial genomes. Microbiology UK. 2005, 151: 3147-3150. 10.1099/mic.0.28339-0.

    Article  CAS  Google Scholar 

  33. Kiil K, Ferchaud JB, David C, Binnewies TT, Wu H, Sicheritz-Pontén T, Willenbrock H, Ussery DW: Genome update: distribution of two-component transduction systems in 250 bacterial genomes. Microbiology UK. 2005, 151: 3447-3452. 10.1099/mic.0.28423-0.

    Article  Google Scholar 

  34. Hutchings MI, Hoskisson PA, Chandra G, Buttner MJ: Sensing and responding to diverse extracellular signals? Analysis of the sensor kinases and response regulators of Streptomyces coelicolor A3(2). Microbiology UK. 2004, 150: 2795-2806. 10.1099/mic.0.27181-0.

    Article  CAS  Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. Aravind L, Ponting CP: The GAF domain: and evolutionary link between diverse phototransducing proteins. Trends Biochem Sci. 1997, 22 (12): 458-459. 10.1016/S0968-0004(97)01148-1.

    Article  CAS  PubMed  Google Scholar 

  37. Ho YS, Burden LM, Hurley JH: Structure of the GAF domain, a ubiquitous signalling motif and a new class of cyclic GMP receptor. EMBO J. 2000, 19: 5288-5299. 10.1093/emboj/19.20.5288.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Karniol B, Wagner JR, Walker JM, Vierstra RD: Phylogenetic analysis of the phytochrome superfamily reveals distinct microbial subfamilies of photoreceptors. Biochem J. 2005, 392: 103-116. 10.1042/BJ20050826.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Studholme DJ, Dixon R: Domain Architectures of σ54-Dependent Transcriptional Activators. J Bacteriol. 2003, 185: 1757-1767. 10.1128/JB.185.6.1757-1767.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Morett E, Segovia L: The σ54 bacterial enhancer-binding protein family: mechanism of action and phylogenetic relationship of their functional domains. J Bacteriol. 1993, 175: 6067-6074.

    CAS  PubMed Central  PubMed  Google Scholar 

  41. Grimm C, Panopoulos NJ: The predicted protein product of a pathogenicity locus from Pseudomonas syringae pv. phaseolicola is homologous to a highly conserved domain of several prokaryotic regulatory proteins. J Bacteriol. 1989, 171: 5031-5038.

    CAS  PubMed Central  PubMed  Google Scholar 

  42. Xiao Y, Heu S, Yi J, Lu Y, Hutcheson SW: Identification of a putative alternative sigma factor and characterization of a multicomponent regulatory cascade controlling the expression of Pseudomonas syringae pv. syringae Pss61 hrp and hrmA genes. J Bacteriol. 1994, 176: 1025-1036.

    CAS  PubMed Central  PubMed  Google Scholar 

  43. Liyanage H, Palmer DA, Ullrich M, Bender CL: Characterization and transcriptional analysis of the gene cluster for coronafacic acid, the polyketide component of the phytotoxin coronatine. Appl Environ Microbiol. 1995, 61: 3843-3848.

    CAS  PubMed Central  PubMed  Google Scholar 

  44. Ullrich M, Peñaloza-Vázquez A, Bailey A, Bender CL: A modified two-component regulatory system is involved in temperature-dependent biosynthesis of the Pseudomonas syringae phytotoxin coronatine. J Bacteriol. 1995, 177: 6160-6169.

    CAS  PubMed Central  PubMed  Google Scholar 

  45. Cotter PA, Jones AM: Phosphorelay control of virulence gene expression in Bordetella. Trends Microbiol. 2003, 11: 367-373. 10.1016/S0966-842X(03)00156-2.

    Article  CAS  PubMed  Google Scholar 

  46. Drenkard E, Ausubel FM: Pseudomonas biofilm formation and antibiotic resistance are linked to phenotypic variation. Nature. 2002, 416: 740-743. 10.1038/416740a.

    Article  CAS  PubMed  Google Scholar 

  47. Lee SH, Angelichio MJ, Mekalanos JJ, Camilli A: Nucleotide sequence and spatiotemporal expression of the Vibrio cholerae vieSAB genes during infection. J Bacteriol. 1998, 180: 2298-2305.

    CAS  PubMed Central  PubMed  Google Scholar 

  48. Morris CE, Kinkel LL, Xiao T, Prior P, Sands DC: Surprising niche for the plant pathogen Pseudomonas syringae. Infect Genet Evol. 2007, 7: 84-92. 10.1016/j.meegid.2006.05.002.

    Article  CAS  PubMed  Google Scholar 

  49. Hutcheson SW, Bretz J, Sussan T, Jin S, Pak K: Enhancer-binding HrpR and HrpS interact to regulate hrp-encoded type III protein secretion in Pseudomonas syringae strains. J Bacteriol. 2001, 183: 5589-5598. 10.1128/JB.183.19.5589-5598.2001.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. Chatterjee A, Cul Y, Yang H, Collmer A, Alfano JR, Chatterjee AK: GacA, the response regulator of a two-component system, acts as a master regulator in Pseudomonas syringae pv. tomato DC3000 by controlling regulatory RNA, transcriptional activators, and alternate sigma factors. Mol Plant-Microbe Interact. 2003, 16: 1106-1117. 10.1094/MPMI.2003.16.12.1106.

    Article  CAS  PubMed  Google Scholar 

  51. Xiao Y, Lan L, Yin C, Deng X, Baker D, Zhou JM, Tang X: Two-component sensor RhpS promotes induction of Pseudomonas syringae type III secretion system by repressing negative regulator RhpR. Mol Plant Microbe Interact. 2007, 20: 223-234. 10.1094/MPMI-20-3-0223.

    Article  CAS  PubMed  Google Scholar 

  52. Rich JJ, Kinscherf TG, Kitten T, Willis DK: Genetic evidence that the gacA gene encodes the cognate response regulator for the lemA sensor in Pseudomonas syringae. J Bacteriol. 1994, 176: 7468-7475.

    CAS  PubMed Central  PubMed  Google Scholar 

  53. Mills SD, Jasalavich CA, Cooksey DA: A two-component regulatory system required for copper-inducible expression of the copper resistance operon of Pseudomonas syringae. J Bacteriol. 1993, 175: 1656-1664.

    CAS  PubMed Central  PubMed  Google Scholar 

  54. Swartz TE, Tseng TS, Frederickson MA, Paris G, Comerci DJ, Rajashekara G, Kim JG, Mudgett MB, Splitter GA, Ugalde RA, Goldbaum FA, Briggs WR, Bogomolni RA: Blue-light-activated histidine kinases: two-component sensors in bacteria. Science. 2007, 317: 1090-1093. 10.1126/science.1144306.

    Article  CAS  PubMed  Google Scholar 

  55. Stover CK, Pham XQ, Erwin AL, Mizoguchi SD, Warrener P, Hickey MJ, Brinkma FS, Hufnagle WO, Kowalik DJ, Lagrou M, Garber RL, Goltry L, Tolentino E, Westbrock-Wadman S, Yuan Y, Brody LL, Coulter SN, Folger KR, Kas A, Larbig K, Lim R, Smith K, Spencer D, Wong GK, Wu Z, Paulsen IT, Reizer J, Saier MH, Hancock RE, Lory S, Olson MV: Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen. Nature. 2000, 406: 959-964. 10.1038/35023079.

    Article  CAS  PubMed  Google Scholar 

  56. van Nimwegen E: Scaling laws in the functional content of genomes. Trends Genet. 2003, 19: 479-484. 10.1016/S0168-9525(03)00203-8.

    Article  CAS  PubMed  Google Scholar 

  57. Kostantinidis KT, Tiedje JM: Trends between gene content and genome size in prokaryotic species with larger genomes. Proc Natl Acad Sci USA. 2004, 101: 3160-3165. 10.1073/pnas.0308653100.

    Article  Google Scholar 

  58. Büttner D, Noël L, Thieme F, Bonas U: Genomic approaches in Xanthomonas campestris pv. vesicatoria allow fishing for virulence genes. J Biotechnol. 2003, 106: 203-214. 10.1016/j.jbiotec.2003.07.012.

    Article  PubMed  Google Scholar 

  59. Seshasayee AS, Bertone P, Fraser GM, Luscombe NM: Transcriptional regulatory networks in bacteria: from input signals to output responses. Curr Opin Microbiol. 2006, 9: 511-519. 10.1016/j.mib.2006.08.007.

    Article  CAS  PubMed  Google Scholar 

  60. Kobayashi K, Ogura M, Yamaguchi H, Yoshida KI, Ogasawara N, Tanaka T, Fujita Y: Comprehensive DNA microarray analysis of Bacillus subtilis two-component regulatory 2001 systems. J Bacteriol. 2001, 183: 7365-7370. 10.1128/JB.183.24.7365-7370.2001.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  61. Throup JP, Koretke KK, Bryant AP, Ingraham KA, Chalker AF, Ge Y, Marra A, Wallis NG, Brown JR, Holmes DJ, Rosenberg M, Burnham MKR: A genomic analysis of two-component signal transduction in Streptococcus pneumoniae. Mol Microbiol. 2000, 35: 566-576. 10.1046/j.1365-2958.2000.01725.x.

    Article  CAS  PubMed  Google Scholar 

  62. Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, 32: D138-D141. 10.1093/nar/gkh121.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  63. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  64. Marchler-Bauer A, Anderson JB, DeWeese-Scott C, Fedorova ND, Geer LY, He S, Hurwitz DI, Jackson JD, Jacobs AR, Lanczycki CJ, Liebert CA, Liu C, Madej T, Marchler GH, Mazumder R, Nikolskaya AN, Panchenko AR, Rao BS, Shoemaker BA, Simonyan V, Song JS, Thiessen PA, Vasudevan S, Wang Y, Yamashita RA, Yin JJ, Bryant SH: CDD: a curated Entrez database of conserved domain alignments. Nucleic Acids Res. 2003, 31: 383-387. 10.1093/nar/gkg087.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. 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.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  66. 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.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

JLL was a recipient of a predoctoral fellowship from the Public University of Navarra. JAO was supported by the Ramón y Cajal Programme and Complementary Action Grant BIO2006-28484-E of the Spanish Ministerio de Educación y Ciencia (MEC). JLL and JAO thank Antonio G Pisabarro and Lucía Ramírez for continued support. DWU and KK would like to acknowledge funding by the Danish Center for Scientific Computing.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to José A Oguiza.

Additional information

Authors' contributions

DWU and JAO designed and coordinated the project. JLL, KK and OR performed the bioinformatics studies and interpreted the results. JAO wrote the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

12864_2007_1110_MOESM1_ESM.xls

Additional file 1: HKs in the genomes of P. syringae pv. syringae B728a, pv. tomato DC3000 and pv. phaseolicola 1448A. (XLS 118 KB)

12864_2007_1110_MOESM2_ESM.xls

Additional file 2: RRs in the genomes of P. syringae pv. syringae B728a, pv. tomato DC3000 and pv. phaseolicola 1448A. (XLS 149 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lavín, J.L., Kiil, K., Resano, O. et al. Comparative genomic analysis of two-component regulatory proteins in Pseudomonas syringae. BMC Genomics 8, 397 (2007). https://doi.org/10.1186/1471-2164-8-397

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-8-397

Keywords