The MiST2 database: a comprehensive genomics resource on microbial signal transduction

Specificity residues determine binding affinity for two-component signal transduction systems

Two-component systems (TCS) comprise histidine kinases and their cognate response regulators and allow bacteria to sense and respond to a wide variety of signals. Histidine kinases (HKs) phosphorylate and dephosphorylate their cognate response regulators (RRs) in response to stimuli. In general, these reactions appear to be highly specific and require an appropriate association between the HK and RR proteins. The Myxococcus xanthus genome encodes one of the largest repertoires of signaling proteins in bacteria (685 open reading frames [ORFs]), including at least 127 HKs and at least 143 RRs. Of these, 27 are bona fide NtrC-family response regulators, 21 of which are encoded adjacent to their predicted cognate kinases. Using system-wide profiling methods, we determined that the HK-NtrC RR pairs display a kinetic preference during both phosphotransfer and phosphatase functions, thereby defining cognate signaling systems in M. xanthus. Isothermal titration calorimetry measurements indicated that cognate HK-RR pairs interact with dissociation constants (K_d) of approximately 1 µM, while noncognate pairs had no measurable binding. Lastly, a chimera generated between the histidine kinase, CrdS, and HK1190 revealed that residues conferring phosphotransfer and phosphatase specificity dictate binding affinity, thereby establishing discrete protein-protein interactions which prevent cross talk. The data indicate that binding affinity is a critical parameter governing system-wide signaling fidelity for bacterial signal transduction proteins.

Using in vitro phosphotransfer and phosphatase profiling assays and isothermal titration calorimetry, we have taken a system-wide approach to demonstrate specificity for a family of two-component signaling proteins in Myxococcus xanthus. Our results demonstrate that previously identified specificity residues dictate binding affinity and that phosphatase specificity follows phosphotransfer specificity for cognate HK-RR pairs. The data indicate that preferential binding affinity is the basis for signaling fidelity in bacterial two-component systems.

Canonical and non-canonical EcfG sigma factors control the general stress response in Rhizobium etli

A core component of the α-proteobacterial general stress response (GSR) is the extracytoplasmic function (ECF) sigma factor EcfG, exclusively present in this taxonomic class. Half of the completed α-proteobacterial genome sequences contain two or more copies of genes encoding σ^EcfG-like sigma factors, with the primary copy typically located adjacent to genes coding for a cognate anti-sigma factor (NepR) and two-component response regulator (PhyR). So far, the widespread occurrence of additional, non-canonical σ^EcfG copies has not satisfactorily been explained. This study explores the hierarchical relation between Rhizobium etli σ^EcfG1 and σ^EcfG2, canonical and non-canonical σ^EcfG proteins, respectively. Contrary to reports in other species, we find that σ^EcfG1 and σ^EcfG2 act in parallel, as nodes of a complex regulatory network, rather than in series, as elements of a linear regulatory cascade. We demonstrate that both sigma factors control unique yet also shared target genes, corroborating phenotypic evidence. σ^EcfG1 drives expression of rpoH2, explaining the increased heat sensitivity of an ecfG1 mutant, while katG is under control of σ^EcfG2, accounting for reduced oxidative stress resistance of an ecfG2 mutant. We also identify non-coding RNA genes as novel σ^EcfG targets. We propose a modified model for GSR regulation in R. etli, in which σ^EcfG1 and σ^EcfG2 function largely independently. Based on a phylogenetic analysis and considering the prevalence of α-proteobacterial genomes with multiple σ^EcfG copies, this model may also be applicable to numerous other species.

Automated design of bacterial genome sequences

Background

Organisms have evolved ways of regulating transcription to better adapt to varying environments. Could the current functional genomics data and models support the possibility of engineering a genome with completely rearranged gene organization while the cell maintains its behavior under environmental challenges? How would we proceed to design a full nucleotide sequence for such genomes?

Results

As a first step towards answering such questions, recent work showed that it is possible to design alternative transcriptomic models showing the same behavior under environmental variations than the wild-type model. A second step would require providing evidence that it is possible to provide a nucleotide sequence for a genome encoding such transcriptional model. We used computational design techniques to design a rewired global transcriptional regulation of Escherichia coli, yet showing a similar transcriptomic response than the wild-type. Afterwards, we “compiled” the transcriptional networks into nucleotide sequences to obtain the final genome sequence. Our computational evolution procedure ensures that we can maintain the genotype-phenotype mapping during the rewiring of the regulatory network. We found that it is theoretically possible to reorganize E. coli genome into 86% fewer regulated operons. Such refactored genomes are constituted by operons that contain sets of genes sharing around the 60% of their biological functions and, if evolved under highly variable environmental conditions, have regulatory networks, which turn out to respond more than 20% faster to multiple external perturbations.

Conclusions

This work provides the first algorithm for producing a genome sequence encoding a rewired transcriptional regulation with wild-type behavior under alternative environments.

Homology modeling of the CheW coupling protein of the chemotaxis signaling complex

Homology models of the E. coli and T. maritima chemotaxis protein CheW were constructed to assess the quality of structural predictions and their applicability in chemotaxis research: i) a model of E. coli CheW was constructed using the T. maritima CheW NMR structure as a template, and ii) a model of T. maritima CheW was constructed using the E. coli CheW NMR structure as a template. The conformational space accessible to the homology models and to the NMR structures was investigated using molecular dynamics and Monte Carlo simulations. The results show that even though static homology models of CheW may be partially structurally different from their corresponding experimentally determined structures, the conformational space they can access through their dynamic variations can be similar, for specific regions of the protein, to that of the experimental NMR structures. When CheW homology models are allowed to explore their local accessible conformational space, modeling can provide a rational path to predicting CheW interactions with the MCP and CheA proteins of the chemotaxis complex. Homology models of CheW (and potentially, of other chemotaxis proteins) should be seen as snapshots of an otherwise larger ensemble of accessible conformational space.

Freshwater bacterial lifestyles inferred from comparative genomics

While micro-organisms actively mediate and participate in freshwater ecosystem services, we know little about freshwater microbial genetic diversity. Genome sequences are available for many bacteria from the human microbiome and the ocean (over 800 and 200, respectively), but only two freshwater genomes are currently available: the streamlined genomes of Polynucleobacter necessarius ssp. asymbioticus and the Actinobacterium AcI-B1. Here, we sequenced and analysed draft genomes of eight phylogentically diverse freshwater bacteria exhibiting a range of lifestyle characteristics. Comparative genomics of these bacteria reveals putative freshwater bacterial lifestyles based on differences in predicted growth rate, capability to respond to environmental stimuli and diversity of useable carbon substrates. Our conceptual model based on these genomic characteristics provides a foundation on which further ecophysiological and genomic studies can be built. In addition, these genomes greatly expand the diversity of existing genomic context for future studies on the ecology and genetics of freshwater bacteria.

GoSynthetic database tool to analyse natural and engineered molecular processes

An essential topic for synthetic biologists is to understand the structure and function of biological processes and involved proteins and plan experiments accordingly. Remarkable progress has been made in recent years towards this goal. However, efforts to collect and present all information on processes and functions are still cumbersome. The database tool GoSynthetic provides a new, simple and fast way to analyse biological processes applying a hierarchical database. Four different search modes are implemented. Furthermore, protein interaction data, cross-links to organism-specific databases (17 organisms including six model organisms and their interactions), COG/KOG, GO and IntAct are warehoused. The built in connection to technical and engineering terms enables a simple switching between biological concepts and concepts from engineering, electronics and synthetic biology. The current version of GoSynthetic covers more than one million processes, proteins, COGs and GOs. It is illustrated by various application examples probing process differences and designing modifications.

Database URL:http://gosyn.bioapps.biozentrum.uni-wuerzburg.de.

Histidine phosphotransfer proteins in fungal two-component signal transduction pathways

The histidine phosphotransfer (HPt) protein Ypd1 is an important participant in the Saccharomyces cerevisiae multistep two-component signal transduction pathway and, unlike the expanded histidine kinase gene family, is encoded by a single gene in nearly all model and pathogenic fungi. Ypd1 is essential for viability in both S. cerevisiae and in Cryptococcus neoformans. These and other aspects of Ypd1 biology, combined with the availability of structural and mutational data in S. cerevisiae, suggest that the essential interactions between Ypd1 and response regulator domains would be a good target for antifungal drug development. The goal of this minireview is to summarize the wealth of data on S. cerevisiae Ypd1 and to consider the potential benefits of conducting related studies in pathogenic fungi.

A link between dimerization and autophosphorylation of the response regulator PhoB

Response regulator proteins within two-component signal transduction systems are activated by phosphorylation and can catalyze their own covalent phosphorylation using small molecule phosphodonors. To date, comprehensive kinetic characterization of response regulator autophosphorylation is limited to CheY, which follows a simple model of phosphodonor binding followed by phosphorylation. We characterized autophosphorylation of the response regulator PhoB, known to dimerize upon phosphorylation. In contrast to CheY, PhoB time traces exhibited an initial lag phase and gave apparent pseudo-first order rate constants that increased with protein concentration. Furthermore, plots of the apparent autophosphorylation rate constant versus phosphodonor concentration were sigmoidal, as were PhoB binding isotherms for the phosphoryl group analog BeF3(-). Successful mathematical modeling of the kinetic data necessitated inclusion of the formation of a PhoB heterodimer (one phosphorylated and one unphosphorylated monomer) with an enhanced rate of phosphorylation. Specifically, dimerization constants for the PhoB heterodimer and homodimer (two phosphorylated monomers) were similar, but the rate constant for heterodimer phosphorylation was ~10-fold higher than for the monomer. In a test of the model, disruption of the known PhoB(N) dimerization interface by mutation led to markedly slower and noncooperative autophosphorylation kinetics. Furthermore, phosphotransfer from the sensor kinase PhoR was enhanced by dimer formation. Phosphorylation-mediated dimerization allows many response regulators to bind to tandem DNA-binding sites and regulate transcription. Our data challenge the notion that response regulator dimers primarily form between two phosphorylated monomers and raise the possibility that response regulator heterodimers containing one phosphoryl group may participate in gene regulation.

The 3.2 Å resolution structure of a receptor: CheA:CheW signaling complex defines overlapping binding sites and key residue interactions within bacterial chemosensory arrays

Bacterial chemosensory arrays are composed of extended networks of chemoreceptors (also known as methyl-accepting chemotaxis proteins, MCPs), the histidine kinase CheA, and the adaptor protein CheW. Models of these arrays have been developed from cryoelectron microscopy, crystal structures of binary and ternary complexes, NMR spectroscopy, mutational, data and biochemical studies. A new 3.2 Å resolution crystal structure of a Thermotoga maritima MCP protein interaction region in complex with the CheA kinase-regulatory module (P4-P5) and adaptor protein CheW provides sufficient detail to define residue contacts at the interfaces formed among the three proteins. As in a previous 4.5 Å resolution structure, CheA-P5 and CheW interact through conserved hydrophobic surfaces at the ends of their β-barrels to form pseudo 6-fold symmetric rings in which the two proteins alternate around the circumference. The interface between P5 subdomain 1 and CheW subdomain 2 was anticipated from previous studies, whereas the related interface between CheW subdomain 1 and P5 subdomain 2 has only been observed in these ring assemblies. The receptor forms an unexpected structure in that the helical hairpin tip of each subunit has "unzipped" into a continuous α-helix; four such helices associate into a bundle, and the tetramers bridge adjacent P5-CheW rings in the lattice through interactions with both P5 and CheW. P5 and CheW each bind a receptor helix with a groove of conserved hydrophobic residues between subdomains 1 and 2. P5 binds the receptor helix N-terminal to the tip region (lower site), whereas CheW binds the same helix with inverted polarity near the bundle end (upper site). Sequence comparisons among different evolutionary classes of chemotaxis proteins show that the binding partners undergo correlated changes at key residue positions that involve the lower site. Such evolutionary analyses argue that both CheW and P5 bind to the receptor tip at overlapping positions. Computational genomics further reveal that two distinct CheW proteins in Thermotogae utilize the analogous recognition motifs to couple different receptor classes to the same CheA kinase. Important residues for function previously identified by mutagenesis, chemical modification and biophysical approaches also map to these same interfaces. Thus, although the native CheW-receptor interaction is not observed in the present crystal structure, the bioinformatics and previous data predict key features of this interface. The companion study of the P5-receptor interface in native arrays (accompanying paper Piasta et al. (2013) Biochemistry, DOI: 10.1021/bi400385c) shows that, despite the non-native receptor fold in the present crystal structure, the local helix-in-groove contacts of the crystallographic P5-receptor interaction are present in native arrays and are essential for receptor regulation of kinase activity.

Genomic reconstruction of the transcriptional regulatory network in Bacillus subtilis

The adaptation of microorganisms to their environment is controlled by complex transcriptional regulatory networks (TRNs), which are still only partially understood even for model species. Genome scale annotation of regulatory features of genes and TRN reconstruction are challenging tasks of microbial genomics. We used the knowledge-driven comparative-genomics approach implemented in the RegPredict Web server to infer TRN in the model Gram-positive bacterium Bacillus subtilis and 10 related Bacillales species. For transcription factor (TF) regulons, we combined the available information from the DBTBS database and the literature with bioinformatics tools, allowing inference of TF binding sites (TFBSs), comparative analysis of the genomic context of predicted TFBSs, functional assignment of target genes, and effector prediction. For RNA regulons, we used known RNA regulatory motifs collected in the Rfam database to scan genomes and analyze the genomic context of new RNA sites. The inferred TRN in B. subtilis comprises regulons for 129 TFs and 24 regulatory RNA families. First, we analyzed 66 TF regulons with previously known TFBSs in B. subtilis and projected them to other Bacillales genomes, resulting in refinement of TFBS motifs and identification of novel regulon members. Second, we inferred motifs and described regulons for 28 experimentally studied TFs with previously unknown TFBSs. Third, we discovered novel motifs and reconstructed regulons for 36 previously uncharacterized TFs. The inferred collection of regulons is available in the RegPrecise database (http://regprecise.lbl.gov/) and can be used in genetic experiments, metabolic modeling, and evolutionary analysis.

P2RP: a Web-based framework for the identification and analysis of regulatory proteins in prokaryotic genomes

Background

Regulatory proteins (RPs) such as transcription factors (TFs) and two-component system (TCS) proteins control how prokaryotic cells respond to changes in their external and/or internal state. Identification and annotation of TFs and TCSs is non-trivial, and between-genome comparisons are often confounded by different standards in annotation. There is a need for user-friendly, fast and convenient tools to allow researchers to overcome the inherent variability in annotation between genome sequences.

Results

We have developed the web-server P2RP (Predicted Prokaryotic Regulatory Proteins), which enables users to identify and annotate TFs and TCS proteins within their sequences of interest. Users can input amino acid or genomic DNA sequences, and predicted proteins therein are scanned for the possession of DNA-binding domains and/or TCS domains. RPs identified in this manner are categorised into families, unambiguously annotated, and a detailed description of their features generated, using an integrated software pipeline. P2RP results can then be outputted in user-specified formats.

Conclusion

Biologists have an increasing need for fast and intuitively usable tools, which is why P2RP has been developed as an interactive system. As well as assisting experimental biologists to interrogate novel sequence data, it is hoped that P2RP will be built into genome annotation pipelines and re-annotation processes, to increase the consistency of RP annotation in public genomic sequences. P2RP is the first publicly available tool for predicting and analysing RP proteins in users’ sequences. The server is freely available and can be accessed along with documentation at http://www.p2rp.org.

Nonconserved active site residues modulate CheY autophosphorylation kinetics and phosphodonor preference

In two-component signal transduction, response regulator proteins contain the catalytic machinery for their own covalent phosphorylation and can catalyze phosphotransfer from a partner sensor kinase or autophosphorylate using various small molecule phosphodonors. Although response regulator autophosphorylation is physiologically relevant and a powerful experimental tool, the kinetic determinants of the autophosphorylation reaction and how those determinants might vary for different response regulators and phosphodonors are largely unknown. We characterized the autophosphorylation kinetics of 21 variants of the model response regulator Escherichia coli CheY that contained substitutions primarily at nonconserved active site positions D + 2 (CheY residue 59) and T + 2 (CheY residue 89), two residues C-terminal to conserved D57 and T87, respectively. Overall, the CheY variants exhibited a >10(5)-fold range of rate constants (kphos/KS) for reaction with phosphoramidate, acetyl phosphate, or monophosphoimidazole, with the great majority of rates enhanced versus that of wild-type CheY. Although phosphodonor preference varied substantially, nearly all the CheY variants reacted faster with phosphoramidate than acetyl phosphate. Correlation between the increased positive charge of the D + 2 and T + 2 side chains and faster rates indicated electrostatic interactions are a kinetic determinant. Moreover, sensitivities of rate constants to ionic strength indicated that both long-range and localized electrostatic interactions influence autophosphorylation kinetics. The increased nonpolar surface area of the D + 2 and T + 2 side chains also correlated with an enhanced autophosphorylation rate, especially for reaction with phosphoramidate and monophosphoimidazole. Computer docking suggested that highly accelerated monophosphoimidazole autophosphorylation rates for CheY variants with a tyrosine at position T + 2 likely reflect structural mimicry of phosphotransfer from the sensor kinase histidyl phosphate.

Characterization of five ECF sigma factors in the genome of Pseudomonas syringae pv. syringae B728a

Pseudomonas syringae pv. syringae B728a, a bacterial pathogen of bean, utilizes large surface populations and extracellular signaling to initiate a fundamental change from an epiphytic to a pathogenic lifestyle. Extracytoplasmic function (ECF) sigma (σ) factors serve as important regulatory factors in responding to various environmental signals. Bioinformatic analysis of the B728a genome revealed 10 ECF sigma factors. This study analyzed deletion mutants of five previously uncharacterized ECF sigma factor genes in B728a, including three FecI-type ECF sigma factors (ECF5, ECF6, and ECF7) and two ECF sigma factors placed in groups ECF11 and ECF18. Transcriptional profiling by qRT-PCR analysis of ECF sigma factor mutants was used to measure expression of their associated anti-sigma and outer membrane receptor proteins, and expression of genes associated with production of extracellular polysaccharides, fimbriae, glycine betaine and syringomycin. Notably, the B728aΔecf7 mutant displayed reduced swarming and had decreased expression of CupC fimbrial genes. Growth and pathogenicity assays, using a susceptible bean host, revealed that none of the tested sigma factor genes are required for in planta growth and lesion formation.

Identification of disulphide stress-responsive extracytoplasmic function sigma factors in Rothia mucilaginosa

Rothia mucilaginosa is known as a member of commensal bacterial flora in the oral cavity and has received attention as a potential opportunistic pathogen. We previously determined the genomic sequence of R. mucilaginosa DY-18, a clinical strain with biofilm-like structures isolated from an infected root canal of a tooth with persistent apical periodontitis. We found that the DY-18 genome had only two sigma factor genes that encoded the primary and extracytoplasmic function (ECF) sigma factors. Genomic analysis on the available database of R. mucilaginosa ATCC 25296 (a type strain for R. mucilaginosa) revealed that ATCC 25296 has three sigma factors: one primary sigma factor and two ECF sigma factors, one of which was highly homologous to that of DY-18. ECF sigma factors play an important role in the response to environmental stress and to the production of virulence factors. Therefore, we first examined gene-encoding sigma factors on R. mucilaginosa genome in silico. The homologous ECF sigma factors found in strains DY-18 and ATCC 25296 formed a distinct SigH (SigR) clade in a phylogenetic tree and their cognate anti-sigma factor has a HXXXCXXC motif known to respond against disulphide stress. Quantitative reverse transcription polymerase chain reaction (PCR) and microarray analysis showed that the transcriptional levels of sigH were markedly up-regulated under disulphide stress in both strains. Microarray data also demonstrated that several oxidative-stress-related genes (thioredoxin, mycothione reductase, reductase and oxidoreductase) were significantly up-regulated under the diamide stress. On the basis of these results, we conclude that the alternative sigma factor SigH of R. mucilaginosa is a candidate regulator in the redox state.

Identification of proteins likely to be involved in morphogenesis, cell division, and signal transduction in Planctomycetes by comparative genomics

Members of the Planctomycetes clade share many unusual features for bacteria. Their cytoplasm contains membrane-bound compartments, they lack peptidoglycan and FtsZ, they divide by polar budding, and they are capable of endocytosis. Planctomycete genomes have remained enigmatic, generally being quite large (up to 9 Mb), and on average, 55% of their predicted proteins are of unknown function. Importantly, proteins related to the unusual traits of Planctomycetes remain largely unknown. Thus, we embarked on bioinformatic analyses of these genomes in an effort to predict proteins that are likely to be involved in compartmentalization, cell division, and signal transduction. We used three complementary strategies. First, we defined the Planctomycetes core genome and subtracted genes of well-studied model organisms. Second, we analyzed the gene content and synteny of morphogenesis and cell division genes and combined both methods using a "guilt-by-association" approach. Third, we identified signal transduction systems as well as sigma factors. These analyses provide a manageable list of candidate genes for future genetic studies and provide evidence for complex signaling in the Planctomycetes akin to that observed for bacteria with complex life-styles, such as Myxococcus xanthus.

Genetic and biochemical dissection of a HisKA domain identifies residues required exclusively for kinase and phosphatase activities

Two-component signal transduction systems, composed of histidine kinases (HK) and response regulators (RR), allow bacteria to respond to diverse environmental stimuli. The HK can control both phosphorylation and subsequent dephosphorylation of its cognate RR. The majority of HKs utilize the HisKA subfamily of dimerization and histidine phosphotransfer (DHp) domains, which contain the phospho-accepting histidine and directly contact the RR. Extensive genetics, biochemistry, and structural biology on several prototypical TCS systems including NtrB-NtrC and EnvZ-OmpR have provided a solid basis for understanding the function of HK–RR signaling. Recently, work on NarX, a HisKA_3 subfamily protein, indicated that two residues in the highly conserved region of the DHp domain are responsible for phosphatase activity. In this study we have carried out both genetic and biochemical analyses on Myxococcus xanthus CrdS, a member of the HisKA subfamily of bacterial HKs. CrdS is required for the regulation of spore formation in response to environmental stress. Following alanine-scanning mutagenesis of the α1 helix of the DHp domain of CrdS, we determined the role for each mutant protein for both kinase and phosphatase activity. Our results indicate that the conserved acidic residue (E372) immediately adjacent to the site of autophosphorylation (H371) is specifically required for kinase activity but not for phosphatase activity. Conversely, we found that the conserved Thr/Asn residue (N375) was required for phosphatase activity but not for kinase activity. We extended our biochemical analyses to two CrdS homologs from M. xanthus, HK1190 and HK4262, as well as Thermotoga maritima HK853. The results were similar for each HisKA family protein where the conserved acidic residue is required for kinase activity while the conserved Thr/Asn residue is required for phosphatase activity. These data are consistent with conserved mechanisms for kinase and phosphatase activities in the broadly occurring HisKA family of sensor kinases in bacteria.

Bacterial histidine kinases (HK) serve as bifunctional enzymes capable of both phosphorylation and dephosphorylation of their cognate response regulators (RR). The majority of HKs (77%) belong to the HisKA subfamily. While both kinase and phosphatase functions have been assayed for HisKA proteins, relatively few examples have been studied to determine which residues are required for kinase and phosphatase activity. Recent studies on NarX, a HisKA_3 family protein, and the dedicated phosphatases CheZ and CheX illustrate requirements for two amino acids for phosphatase function. In this study, we undertook saturating mutagenesis of the proposed interaction surface between the HK and its cognate RR and conclude that only one residue (T/N) is required exclusively for phosphatase activity for HisKA family proteins in evolutionarily distant organisms Myxococcus xanthus and Thermotoga maritima. In addition, we identified only one residue (E/D), adjacent to the conserved site of phosphorylation, required exclusively for kinase activity within the highly conserved motif H-E/D-x-x-T/N. Because similar sequences are found in nearly all HisKA kinases, these residues provide excellent targets for dissection of kinase and phosphatase activities within this broadly occurring family of bacterial kinases.

A connecter-like factor, CacA, links RssB/RpoS and the CpxR/CpxA two-component system in Salmonella

Background

Bacteria integrate numerous environmental stimuli when generating cellular responses. Increasing numbers of examples describe how one two-component system (TCS) responds to signals detected by the sensor of another TCS. However, the molecular mechanisms underlying this phenomenon remain poorly defined.

Results

Here, we report a connector-like factor that affects the activity of the CpxR/CpxA two-component system in Salmonella enterica serovar Typhimurium. We isolated a clone that induced the expression of a cpxP-lac gene fusion from a high-copy-number plasmid pool of random Salmonella genomic fragments. A 63-amino acid protein, CacA, was responsible for the CpxA/CpxR-dependent activation of the cpxP gene. The CpxR-activated genes cpxP and spy exhibited approximately 30% and 50% reductions in transcription, respectively, in a clean cacA deletion mutant strain in comparison to wild-type. From 33 response regulator (RR) deletion mutants, we identified that the RssB regulator represses cacA transcription. Substitution mutations in a conserved -10 region harboring the RNA polymerase recognition sequence, which is well conserved with a known RpoS -10 region consensus sequence, rendered the cacA promoter RpoS-independent. The CacA-mediated induction of cpxP transcription was affected in a trxA deletion mutant, which encodes thioredoxin 1, suggesting a role for cysteine thiol-disulfide exchange(s) in CacA-dependent Cpx activation.

Conclusions

We identified CacA as an activator of the CpxR/CpxA system in the plasmid clone. We propose that CacA may integrate the regulatory status of RssB/RpoS into the CpxR/CpxA system. Future investigations are necessary to thoroughly elucidate how CacA activates the CpxR/CpxA system.

Extracytoplasmic function σ factors of the widely distributed group ECF41 contain a fused regulatory domain

Bacteria need signal transducing systems to respond to environmental changes. Next to one- and two-component systems, alternative σ factors of the extra-cytoplasmic function (ECF) protein family represent the third fundamental mechanism of bacterial signal transduction. A comprehensive classification of these proteins identified more than 40 phylogenetically distinct groups, most of which are not experimentally investigated. Here, we present the characterization of such a group with unique features, termed ECF41. Among analyzed bacterial genomes, ECF41 σ factors are widely distributed with about 400 proteins from 10 different phyla. They lack obvious anti-σ factors that typically control activity of other ECF σ factors, but their structural genes are often predicted to be cotranscribed with carboxymuconolactone decarboxylases, oxidoreductases, or epimerases based on genomic context conservation. We demonstrate for Bacillus licheniformis and Rhodobacter sphaeroides that the corresponding genes are preceded by a highly conserved promoter motif and are the only detectable targets of ECF41-dependent gene regulation. In contrast to other ECF σ factors, proteins of group ECF41 contain a large C-terminal extension, which is crucial for σ factor activity. Our data demonstrate that ECF41 σ factors are regulated by a novel mechanism based on the presence of a fused regulatory domain.

An orphan chemotaxis sensor regulates virulence and antibiotic tolerance in the human pathogen Pseudomonas aeruginosa

The synthesis of virulence factors by pathogenic bacteria is highly regulated and occurs in response to diverse environmental cues. An array of two component systems (TCSs) serves to link perception of different cues to specific changes in gene expression and/or bacterial behaviour. Those TCSs that regulate functions associated with virulence represent attractive targets for interference in anti-infective strategies for disease control. We have previously identified PA2572 as a putative response regulator required for full virulence of Pseudomonas aeruginosa, the opportunistic human pathogen, to Galleria mellonella (Wax moth) larvae. Here we have investigated the involvement of candidate sensors for signal transduction involving PA2572. Mutation of PA2573, encoding a probable methyl-accepting chemotaxis protein, gave rise to alterations in motility, virulence, and antibiotic resistance, functions which are also controlled by PA2572. Comparative transcriptome profiling of mutants revealed that PA2572 and PA2573 regulate expression of a common set of 49 genes that are involved in a range of biological functions including virulence and antibiotic resistance. Bacterial two-hybrid analysis indicated a REC-dependent interaction between PA2572 and PA2573 proteins. Finally expression of PA2572 in the PA2573 mutant background restored virulence to G. mellonella towards wild-type levels. The findings indicate a role for the orphan chemotaxis sensor PA2573 in the regulation of virulence and antibiotic tolerance in P. aeruginosa and indicate that these effects are exerted in part through signal transduction involving PA2572.

Evolution of two-component signal transduction systems

To exist in a wide range of environmental niches, bacteria must sense and respond to a variety of external signals. A primary means by which this occurs is through two-component signal transduction pathways, typically composed of a sensor histidine kinase that receives the input stimuli and then phosphorylates a response regulator that effects an appropriate change in cellular physiology. Histidine kinases and response regulators have an intrinsic modularity that separates signal input, phosphotransfer, and output response; this modularity has allowed bacteria to dramatically expand and diversify their signaling capabilities. Recent work has begun to reveal the molecular basis by which two-component proteins evolve. How and why do orthologous signaling proteins diverge? How do cells gain new pathways and recognize new signals? What changes are needed to insulate a new pathway from existing pathways? What constraints are there on gene duplication and lateral gene transfer? Here, we review progress made in answering these questions, highlighting how the integration of genome sequence data with experimental studies is providing major new insights.

Predicting inter-species cross-talk in two-component signalling systems

Phosphosignalling pathways are an attractive option for the synthetic biologist looking for a wide repertoire of modular components from which to build. We demonstrate that two-component systems can be used in synthetic biology. However, their potential is limited by the fact that host cells contain many of their own phosphosignalling pathways and these may interact with, and cross-talk to, the introduced synthetic components. In this paper we also demonstrate a simple bioinformatic tool that can help predict whether interspecies cross-talk between introduced and native two-component signalling pathways will occur and show both in vitro and in vivo that the predicted interactions do take place. The ability to predict potential cross-talk prior to designing and constructing novel pathways or choosing a host organism is essential for the promise that phosphosignalling components hold for synthetic biology to be realised.

Different evolutionary modifications as a guide to rewire two-component systems

Two-component systems (TCS) are short signalling pathways generally occurring in prokaryotes. They frequently regulate prokaryotic stimulus responses and thus are also of interest for engineering in biotechnology and synthetic biology. The aim of this study is to better understand and describe rewiring of TCS while investigating different evolutionary scenarios. Based on large-scale screens of TCS in different organisms, this study gives detailed data, concrete alignments, and structure analysis on three general modification scenarios, where TCS were rewired for new responses and functions: (i) exchanges in the sequence within single TCS domains, (ii) exchange of whole TCS domains; (iii) addition of new components modulating TCS function. As a result, the replacement of stimulus and promotor cassettes to rewire TCS is well defined exploiting the alignments given here. The diverged TCS examples are non-trivial and the design is challenging. Designed connector proteins may also be useful to modify TCS in selected cases.

DNA binding site analysis of Burkholderia thailandensis response regulators

Bacterial response regulators (RR) that function as transcription factors in two component signaling pathways are crucial for ensuring tight regulation and coordinated expression of the genome. Currently, consensus DNA binding sites in the promoter for very few bacterial RRs have been identified. A systematic method to characterize these DNA binding sites for RRs would enable prediction of specific gene expression patterns in response to extracellular stimuli. To identify RR DNA binding sites, we functionally activated RRs using beryllofluoride and applied them to a protein-binding microarray (PBM) to discover DNA binding motifs for RRs expressed in Burkholderia, a Gram-negative bacterial genus. We identified DNA binding motifs for conserved RRs in Burkholderia thailandensis, including KdpE, RisA, and NarL, as well as for a previously uncharacterized RR at locus BTH_II2335 and its ortholog in the human pathogen Burkholderia pseudomallei at locus BPSS2315. We further demonstrate RR binding of predicted genomic targets for the two orthologs using gel shift assays and reveal a pattern of RR regulation of expression of self and other two component systems. Our studies illustrate the use of PBMs to identify DNA binding specificities for bacterial RRs and enable prediction of gene regulatory networks in response to two component signaling.

The passive yet successful way of planktonic life: genomic and experimental analysis of the ecology of a free-living polynucleobacter population

Background

The bacterial taxon Polynucleobacter necessarius subspecies asymbioticus represents a group of planktonic freshwater bacteria with cosmopolitan and ubiquitous distribution in standing freshwater habitats. These bacteria comprise <1% to 70% (on average about 20%) of total bacterioplankton cells in various freshwater habitats. The ubiquity of this taxon was recently explained by intra-taxon ecological diversification, i.e. specialization of lineages to specific environmental conditions; however, details on specific adaptations are not known. Here we investigated by means of genomic and experimental analyses the ecological adaptation of a persistent population dwelling in a small acidic pond.

Findings

The investigated population (F10 lineage) contributed on average 11% to total bacterioplankton in the pond during the vegetation periods (ice-free period, usually May to November). Only a low degree of genetic diversification of the population could be revealed. These bacteria are characterized by a small genome size (2.1 Mb), a relatively small number of genes involved in transduction of environmental signals, and the lack of motility and quorum sensing. Experiments indicated that these bacteria live as chemoorganotrophs by mainly utilizing low-molecular-weight substrates derived from photooxidation of humic substances.

Conclusions

Evolutionary genome streamlining resulted in a highly passive lifestyle so far only known among free-living bacteria from pelagic marine taxa dwelling in environmentally stable nutrient-poor off-shore systems. Surprisingly, such a lifestyle is also successful in a highly dynamic and nutrient-richer environment such as the water column of the investigated pond, which was undergoing complete mixis and pronounced stratification in diurnal cycles. Obviously, metabolic and ecological versatility is not a prerequisite for long-lasting establishment of abundant bacterial populations under highly dynamic environmental conditions. Caution should be exercised when generalizing the obtained insights into the ecology and adaptation of the investigated lineage to other Polynucleobacter lineages.

Extra cytoplasmic function σ factor activation

The bacterial cell envelope is essential for cell viability and is a target for numerous antibiotics and host immune defenses. Thus bacteria must sense and respond to damage to the cell envelope. Many bacteria utilize alternative σ factors such as extracytoplasmic function (ECF) σ factors to respond to cell envelope stress. Although ECF σ factors are utilized by both Gram negative and Gram positive bacteria to respond to cell envelope stress, the mechanisms of sensing differ. In this review, we examine the events and proteins that are required for activation of two model extracytoplasmic function σ factors, σ(E) in E. coli and σ(W) in B. subtilis.

Structural basis of response regulator inhibition by a bacterial anti-activator protein

The complex interplay between the response regulator ComA, the anti-activator RapF, and the signaling peptide PhrF controls competence development in Bacillus subtilis. More specifically, ComA drives the expression of genetic competence genes, while RapF inhibits the interaction of ComA with its target promoters. The signaling peptide PhrF accumulates at high cell density and upregulates genetic competence by antagonizing the interaction of RapF and ComA. How RapF functions mechanistically to inhibit ComA activity and how PhrF in turn antagonizes the RapF-ComA interaction were unknown. Here we present the X-ray crystal structure of RapF in complex with the ComA DNA binding domain. Along with biochemical and genetic studies, the X-ray crystal structure reveals how RapF mechanistically regulates ComA function. Interestingly, we found that a RapF surface mimics DNA to block ComA binding to its target promoters. Furthermore, RapF is a monomer either alone or in complex with PhrF, and it undergoes a conformational change upon binding to PhrF, which likely causes the dissociation of ComA from the RapF-ComA complex. Finally, we compare the structure of RapF complexed with the ComA DNA binding domain and the structure of RapH complexed with Spo0F. This comparison reveals that RapF and RapH have strikingly similar overall structures, and that they have evolved different, non-overlapping surfaces to interact with diverse cellular targets. To our knowledge, the data presented here reveal the first atomic level insight into the inhibition of response regulator DNA binding by an anti-activator. Compounds that affect the interaction of Rap and Rap-like proteins with their target domains could serve to regulate medically and commercially important phenotypes in numerous Bacillus species, such as sporulation in B. anthracis and sporulation and the production of Cry protein endotoxin in B. thuringiensis.

Sinorhizobium meliloti CheA complexed with CheS exhibits enhanced binding to CheY1, resulting in accelerated CheY1 dephosphorylation

Retrophosphorylation of the histidine kinase CheA in the chemosensory transduction chain is a widespread mechanism for efficient dephosphorylation of the activated response regulator. First discovered in Sinorhizobium meliloti, the main response regulator CheY2-P shuttles its phosphoryl group back to CheA, while a second response regulator, CheY1, serves as a sink for surplus phosphoryl groups from CheA-P. We have identified a new component in this phospho-relay system, a small 97-amino-acid protein named CheS. CheS has no counterpart in enteric bacteria but revealed distinct similarities to proteins of unknown function in other members of the α subgroup of proteobacteria. Deletion of cheS causes a phenotype similar to that of a cheY1 deletion strain. Fluorescence microscopy revealed that CheS is part of the polar chemosensory cluster and that its cellular localization is dependent on the presence of CheA. In vitro binding, as well as coexpression and copurification studies, gave evidence of CheA/CheS complex formation. Using limited proteolysis coupled with mass spectrometric analyses, we defined CheA(163-256) to be the CheS binding domain, which overlaps with the N-terminal part of the CheY2 binding domain (CheA(174-316)). Phosphotransfer experiments using isolated CheA-P showed that dephosphorylation of CheY1-P but not CheY2-P is increased in the presence of CheS. As determined by surface plasmon resonance spectroscopy, CheY1 binds ∼100-fold more strongly to CheA/CheS than to CheA. We propose that CheS facilitates signal termination by enhancing the interaction of CheY1 and CheA, thereby promoting CheY1-P dephosphorylation, which results in a more efficient drainage of the phosphate sink.

Azospirillum genomes reveal transition of bacteria from aquatic to terrestrial environments

Fossil records indicate that life appeared in marine environments ∼3.5 billion years ago (Gyr) and transitioned to terrestrial ecosystems nearly 2.5 Gyr. Sequence analysis suggests that “hydrobacteria” and “terrabacteria” might have diverged as early as 3 Gyr. Bacteria of the genus Azospirillum are associated with roots of terrestrial plants; however, virtually all their close relatives are aquatic. We obtained genome sequences of two Azospirillum species and analyzed their gene origins. While most Azospirillum house-keeping genes have orthologs in its close aquatic relatives, this lineage has obtained nearly half of its genome from terrestrial organisms. The majority of genes encoding functions critical for association with plants are among horizontally transferred genes. Our results show that transition of some aquatic bacteria to terrestrial habitats occurred much later than the suggested initial divergence of hydro- and terrabacterial clades. The birth of the genus Azospirillum approximately coincided with the emergence of vascular plants on land.

Genome sequencing and analysis of plant-associated beneficial soil bacteria Azospirillum spp. reveals that these organisms transitioned from aquatic to terrestrial environments significantly later than the suggested major Precambrian divergence of aquatic and terrestrial bacteria. Separation of Azospirillum from their close aquatic relatives coincided with the emergence of vascular plants on land. Nearly half of the Azospirillum genome has been acquired horizontally, from distantly related terrestrial bacteria. The majority of horizontally acquired genes encode functions that are critical for adaptation to the rhizosphere and interaction with host plants.

Direct-coupling analysis of residue coevolution captures native contacts across many protein families

The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced direct-coupling analysis (DCA). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intradomain residue contacts, arising, e.g., from alternative protein conformations, ligand-mediated residue couplings, and interdomain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, contingent on the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.

Activated chemoreceptor arrays remain intact and hexagonally packed

Bacterial chemoreceptors cluster into exquisitively sensitive, tunable, highly ordered, polar arrays. While these arrays serve as paradigms of cell signalling in general, it remains unclear what conformational changes transduce signals from the periplasmic tips, where attractants and repellents bind, to the cytoplasmic signalling domains. Conflicting reports support and contest the hypothesis that activation causes large changes in the packing arrangement of the arrays, up to and including their complete disassembly. Using electron cryotomography, here we show that in Caulobacter crescentus, chemoreceptor arrays in cells grown in different media and immediately after exposure to the attractant galactose all exhibit the same 12 nm hexagonal packing arrangement, array size and other structural parameters. ΔcheB and ΔcheR mutants mimicking attractant- or repellent-bound states prior to adaptation also show the same lattice structure. We conclude that signal transduction and amplification must be accomplished through only small, nanoscale conformational changes.

Two-component systems and their co-option for eukaryotic signal transduction

Two-component signaling pathways involve histidine kinases, response regulators, and sometimes histidine-containing phosphotransfer proteins. Prevalent in prokaryotes, these signaling elements have also been co-opted to meet the needs of signal transduction in eukaryotes such as fungi and plants. Here we consider the evolution of such regulatory systems, with a particular emphasis on the roles they play in signaling by the plant hormones cytokinin and ethylene, in phytochrome-mediated perception of light, and as integral components of the circadian clock.

CrdS and CrdA comprise a two-component system that is cooperatively regulated by the Che3 chemosensory system in Myxococcus xanthus

Myxococcus xanthus serves as a model organism for development and complex signal transduction. Regulation of developmental aggregation and sporulation is controlled, in part, by the Che3 chemosensory system. The Che3 pathway consists of homologs to two methyl-accepting chemotaxis proteins (MCPs), CheA, CheW, CheB, and CheR but not CheY. Instead, the output for Che3 is the NtrC homolog CrdA, which functions to regulate developmental gene expression. In this paper we have identified an additional kinase, CrdS, which directly regulates the phosphorylation state of CrdA. Both epistasis and in vitro phosphotransfer assays indicate that CrdS functions as part of the Che3 pathway and, in addition to CheA3, serves to regulate CrdA phosphorylation in M. xanthus. We provide kinetic data for CrdS autophosphorylation and demonstrate specificity for phosphotransfer from CrdS to CrdA. We further demonstrate that CheA3 destabilizes phosphorylated CrdA (CrdA~P), indicating that CheA3 likely acts as a phosphatase. Both CrdS and CheA3 control developmental progression by regulating the phosphorylation state of CrdA~P in the cell. These results support a model in which a classical two-component system and a chemosensory system act synergistically to control the activity of the response regulator CrdA.

While phosphorylation-mediated signal transduction is well understood in prototypical chemotaxis and two-component systems (TCS), chemosensory regulation of alternative cellular functions (ACF) has not been clearly defined. The Che3 system in Myxococcus xanthus is a member of the ACF class of chemosensory systems and regulates development via the transcription factor CrdA (chemosensory regulator of development) (K. Wuichet and I. B. Zhulin, Sci. Signal. 3:ra50, 2010; J. R. Kirby and D. R. Zusman, Proc. Natl. Acad. Sci. U. S. A. 100:2008–2013, 2003). We have identified and characterized a homolog of NtrB, designated CrdS, capable of specifically phosphorylating the NtrC homolog CrdA in M. xanthus. Additionally, we demonstrate that the CrdSA two-component system is negatively regulated by CheA3, the central processor within the Che3 system of M. xanthus. To our knowledge, this study provides the first example of an ACF chemosensory system regulating a prototypical two-component system and extends our understanding of complex regulation of developmental signaling pathways.

The atypical two-component sensor kinase Lpl0330 from Legionella pneumophila controls the bifunctional diguanylate cyclase-phosphodiesterase Lpl0329 to modulate bis-(3'-5')-cyclic dimeric GMP synthesis

A significant part of bacterial two-component system response regulators contains effector domains predicted to be involved in metabolism of bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP), a second messenger that plays a key role in many physiological processes. The intracellular level of c-di-GMP is controlled by diguanylate cyclase and phosphodiesterases activities associated with GGDEF and EAL domains, respectively. The Legionella pneumophila Lens genome displays 22 GGDEF/EAL domain-encoding genes. One of them, lpl0329, encodes a protein containing a two-component system receiver domain and both GGDEF and EAL domains. Here, we demonstrated that the GGDEF and EAL domains of Lpl0329 are both functional and lead to simultaneous synthesis and hydrolysis of c-di-GMP. Moreover, these two opposite activities are finely regulated by Lpl0329 phosphorylation due to the atypical histidine kinase Lpl0330. Indeed, Lpl0330 was found to autophosphorylate on a histidine residue in an atypical H box, which is conserved in various bacteria species and thus defines a new histidine kinase subfamily. Lpl0330 also catalyzes the phosphotransferase to Lpl0329, which results in a diguanylate cyclase activity decrease whereas phosphodiesterase activity remains efficient. Altogether, these data present (i) a new histidine kinase subfamily based on the conservation of an original H box that we named HGN H box, and (ii) the first example of a bifunctional enzyme that modulates synthesis and turnover of c-di-GMP in response to phosphorylation of its receiver domain.

Comparative genomics of the dormancy regulons in mycobacteria

In response to stresses, Mycobacterium cells become dormant. This process is regulated by the DosR transcription factor. In Mycobacterium tuberculosis, the dormancy regulon is well characterized and contains the dosR gene itself and dosS and dosT genes encoding DosR kinases, nitroreductases (acg; Rv3131), diacylglycerol acyltransferase (DGAT) (Rv3130c), and many universal stress proteins (USPs). In this study, we apply comparative genomic analysis to characterize the DosR regulons in nine Mycobacterium genomes, Rhodococcus sp. RHA1, Nocardia farcinica, and Saccharopolyspora erythraea. The regulons are highly labile, containing eight core gene groups (regulators, kinases, USPs, DGATs, nitroreductases, ferredoxins, heat shock proteins, and the orthologs of the predicted kinase [Rv2004c] from M. tuberculosis) and 10 additional genes with more restricted taxonomic distribution that are mostly involved in anaerobic respiration. The largest regulon is observed in M. marinum and the smallest in M. abscessus. Analysis of large gene families encoding USPs, nitroreductases, and DGATs demonstrates a mosaic distribution of regulated and nonregulated members, suggesting frequent acquisition and loss of DosR-binding sites.

Functional characterisation of Lp_2714, an EAL-domain protein from Lactobacillus plantarum

Bioinformatic analysis of lp_2714 from Lactobacillus plantarum WCFS1 demonstrates that it encodes an EAL-domain protein associated with a membrane targeting signal-sequence. Comparison of the predicted primary amino-acid sequence of Lp_2714 shows that it lacks critical catalytic residues and heterologous expression has determined that it does not encode a functional phosphodiesterase. We designate Lp_2714 as a class-3 EAL domain protein probably involved in regulating polysaccharide synthesis on the cell surface the cell.

Coevolution of ABC transporters and two-component regulatory systems as resistance modules against antimicrobial peptides in Firmicutes Bacteria

In Firmicutes bacteria, ATP-binding cassette (ABC) transporters have been recognized as important resistance determinants against antimicrobial peptides. Together with neighboring two-component systems (TCSs), which regulate their expression, they form specific detoxification modules. Both the transport permease and sensor kinase components show unusual domain architecture: the permeases contain a large extracellular domain, while the sensor kinases lack an obvious input domain. One of the best-characterized examples is the bacitracin resistance module BceRS-BceAB of Bacillus subtilis. Strikingly, in this system, the ABC transporter and TCS have an absolute mutual requirement for each other in both sensing of and resistance to bacitracin, suggesting a novel mode of signal transduction in which the transporter constitutes the actual sensor. We identified over 250 such BceAB-like ABC transporters in the current databases. They occurred almost exclusively in Firmicutes bacteria, and 80% of the transporters were associated with a BceRS-like TCS. Phylogenetic analyses of the permease and sensor kinase components revealed a tight evolutionary correlation. Our findings suggest a direct regulatory interaction between the ABC transporters and TCSs, mediating communication between both components. Based on their observed coclustering and conservation of response regulator binding sites, we could identify putative corresponding two-component systems for transporters lacking a regulatory system in their immediate neighborhood. Taken together, our results show that these types of ABC transporters and TCSs have coevolved to form self-sufficient detoxification modules against antimicrobial peptides, widely distributed among Firmicutes bacteria.

Dissecting the specificity of protein-protein interaction in bacterial two-component signaling: orphans and crosstalks

Predictive understanding of the myriads of signal transduction pathways in a cell is an outstanding challenge of systems biology. Such pathways are primarily mediated by specific but transient protein-protein interactions, which are difficult to study experimentally. In this study, we dissect the specificity of protein-protein interactions governing two-component signaling (TCS) systems ubiquitously used in bacteria. Exploiting the large number of sequenced bacterial genomes and an operon structure which packages many pairs of interacting TCS proteins together, we developed a computational approach to extract a molecular interaction code capturing the preferences of a small but critical number of directly interacting residue pairs. This code is found to reflect physical interaction mechanisms, with the strongest signal coming from charged amino acids. It is used to predict the specificity of TCS interaction: Our results compare favorably to most available experimental results, including the prediction of 7 (out of 8 known) interaction partners of orphan signaling proteins in Caulobacter crescentus. Surveying among the available bacterial genomes, our results suggest 15∼25% of the TCS proteins could participate in out-of-operon “crosstalks”. Additionally, we predict clusters of crosstalking candidates, expanding from the anecdotally known examples in model organisms. The tools and results presented here can be used to guide experimental studies towards a system-level understanding of two-component signaling.

c-di-GMP turn-over in Clostridium difficile is controlled by a plethora of diguanylate cyclases and phosphodiesterases

Clostridium difficile infections have become a major healthcare concern in the last decade during which the emergence of new strains has underscored this bacterium's capacity to cause persistent epidemics. c-di-GMP is a bacterial second messenger regulating diverse bacterial phenotypes, notably motility and biofilm formation, in proteobacteria such as Vibrio cholerae, Pseudomonas aeruginosa, and Salmonella. c-di-GMP is synthesized by diguanylate cyclases (DGCs) that contain a conserved GGDEF domain. It is degraded by phosphodiesterases (PDEs) that contain either an EAL or an HD-GYP conserved domain. Very little is known about the role of c-di-GMP in the regulation of phenotypes of Gram-positive or fastidious bacteria. Herein, we exposed the main components of c-di-GMP signalling in 20 genomes of C. difficile, revealed their prevalence, and predicted their enzymatic activity. Ectopic expression of 31 of these conserved genes was carried out in V. cholerae to evaluate their effect on motility and biofilm formation, two well-characterized phenotype alterations associated with intracellular c-di-GMP variation in this bacterium. Most of the predicted DGCs and PDEs were found to be active in the V. cholerae model. Expression of truncated versions of CD0522, a protein with two GGDEF domains and one EAL domain, suggests that it can act alternatively as a DGC or a PDE. The activity of one purified DGC (CD1420) and one purified PDE (CD0757) was confirmed by in vitro enzymatic assays. GTP was shown to be important for the PDE activity of CD0757. Our results indicate that, in contrast to most Gram-positive bacteria including its closest relatives, C. difficile encodes a large assortment of functional DGCs and PDEs, revealing that c-di-GMP signalling is an important and well-conserved signal transduction system in this human pathogen.

The genetic organisation of prokaryotic two-component system signalling pathways

Background

Two-component systems (TCSs) are modular and diverse signalling pathways, involving a stimulus-responsive transfer of phosphoryl groups from transmitter to partner receiver domains. TCS gene and domain organisation are both potentially informative regarding biological function, interaction partnerships and molecular mechanisms. However, there is currently little understanding of the relationships between domain architecture, gene organisation and TCS pathway structure.

Results

Here we classify the gene and domain organisation of TCS gene loci from 1405 prokaryotic replicons (>40,000 TCS proteins). We find that 200 bp is the most appropriate distance cut-off for defining whether two TCS genes are functionally linked. More than 90% of all TCS gene loci encode just one or two transmitter and/or receiver domains, however numerous other geometries exist, often with large numbers of encoded TCS domains. Such information provides insights into the distribution of TCS domains between genes, and within genes. As expected, the organisation of TCS genes and domains is affected by phylogeny, and plasmid-encoded TCS exhibit differences in organisation from their chromosomally-encoded counterparts.

Conclusions

We provide here an overview of the genomic and genetic organisation of TCS domains, as a resource for further research. We also propose novel metrics that build upon TCS gene/domain organisation data and allow comparisons between genomic complements of TCSs. In particular, 'percentage orphaned TCS genes' (or 'Dissemination') and 'percentage of complex loci' (or 'Sophistication') appear to be useful discriminators, and to reflect mechanistic aspects of TCS organisation not captured by existing metrics.

P2CS: a database of prokaryotic two-component systems

P2CS (http://www.p2cs.org) is a specialized database for prokaryotic two-component systems (TCSs), virtually ubiquitous signalling proteins which regulate a wide range of physiological processes. The primary aim of the database is to annotate and classify TCS proteins from completely sequenced prokaryotic genomes and metagenomes. Information within P2CS can be accessed through a variety of routes-TCS complements can be browsed by metagenome, replicon or sequence cluster (and these genesets are available for download by users). Alternatively a variety of database-wide or taxon-specific searches are supported. Each TCS protein is fully annotated with sequence-feature information including replicon context, while properties of the predicted proteins can be queried against several external prediction servers to suggest homologues, interaction networks, sub-cellular localization and domain complements. Another unique feature of P2CS is the analysis of ORFeomes to identify TCS genes missed during genome annotation. Recent innovations for P2CS include a CGView representation of the distribution of TCS genes around a replicon, categorization of TCS genes based on gene organization, an expanded domain-based classification scheme, a P2CS 'gene cart' and categorization on the basis of sequence clusters.

Transcriptional responses of Leptospira interrogans to host innate immunity: significant changes in metabolism, oxygen tolerance, and outer membrane

Background

Leptospira interrogans is the major causative agent of leptospirosis. Phagocytosis plays important roles in the innate immune responses to L. interrogans infection, and L. interrogans can evade the killing of phagocytes. However, little is known about the adaptation of L. interrogans during this process.

Methodology/Principal Findings

To better understand the interaction of pathogenic Leptospira and innate immunity, we employed microarray and comparative genomics analyzing the responses of L. interrogans to macrophage-derived cells. During this process, L. interrogans altered expressions of many genes involved in carbohydrate and lipid metabolism, energy production, signal transduction, transcription and translation, oxygen tolerance, and outer membrane proteins. Among them, the catalase gene expression was significantly up-regulated, suggesting it may contribute to resisting the oxidative pressure of the macrophages. The expressions of several major outer membrane protein (OMP) genes (e.g., ompL1, lipL32, lipL41, lipL48 and ompL47) were dramatically down-regulated (10–50 folds), consistent with previous observations that the major OMPs are differentially regulated in vivo. The persistent down-regulations of these major OMPs were validated by immunoblotting. Furthermore, to gain initial insight into the gene regulation mechanisms in L. interrogans, we re-defined the transcription factors (TFs) in the genome and identified the major OmpR TF gene (LB333) that is concurrently regulated with the major OMP genes, suggesting a potential role of LB333 in OMPs regulation.

Conclusions/Significance

This is the first report on global responses of pathogenic Leptospira to innate immunity, which revealed that the down-regulation of the major OMPs may be an immune evasion strategy of L. interrogans, and a putative TF may be involved in governing these down-regulations. Alterations of the leptospiral OMPs up interaction with host antigen-presenting cells (APCs) provide critical information for selection of vaccine candidates. In addition, genome-wide annotation and comparative analysis of TFs set a foundation for further studying regulatory networks in Leptospira spp.

Leptospirosis is an important tropical disease around the world, particularly in humid tropical and subtropical countries. As a major pathogen of this disease, Leptospira interrogans can be shed from the urine of reservoir hosts, survive in soil and water, and infect humans through broken skin or mucous membranes. Recently, host adaptability and immune evasion of L. interrogans to host innate immunity was partially elucidated in infection or animal models. A better understanding of the molecular mechanisms of L. interrogans in response to host innate immunity is required to learn the nature of early leptospirosis. This study focused on the transcriptome of L. interrogans during host immune cells interaction. Significant changes in energy metabolism, oxygen tolerance and outer membrane protein profile were identified as potential immune evasion strategies by pathogenic Leptospira during the early stage of infection. The major outer membrane proteins (OMPs) of L. interrogans may be regulated by the major OmpR specific transcription factor (LB333). These results provide a foundation for further studying the pathogenesis of leptospirosis, as well as identifying gene regulatory networks in Leptospira spp.

Identification of YsrT and evidence that YsrRST constitute a unique phosphorelay system in Yersinia enterocolitica

Two-component systems (TCS) and phosphorelay systems are mechanisms used by bacteria and fungi to quickly adapt to environmental changes to produce proteins necessary for survival in new environments. Bacterial pathogens use TCS and phosphorelay systems to regulate genes necessary to establish infection within their hosts, including type III secretion systems (T3SS). The Yersinia enterocolitica ysa T3SS is activated in response to NaCl by YsrS and YsrR, a putative hybrid sensor kinase and a response regulator, respectively. Hybrid TCS consist of a sensor kinase that typically has three well-conserved sites of phosphorylation: autophosphorylation site H1, D1 within a receiver domain, and H2 in the histidine phosphotransferase (HPt) domain. From H2, the phosphoryl group is transferred to D2 on the response regulator. A curious feature of YsrS is that it lacks the terminal HPt domain. We report here the identification of the HPt-containing protein (YsrT) that provides this activity for the Ysr system. YsrT is an 82-residue protein predicted to be cytosolic and α-helical in nature and is encoded by a gene adjacent to ysrS. To demonstrate predicted functions of YsrRST as a phosphorelay system, we introduced alanine substitutions at H1, D1, H2, and D2 and tested the mutant proteins for the ability to activate a ysaE-lacZ reporter. As expected, substitutions at H1, H2, and D2 resulted in a loss of activation of ysaE expression. This indicates an interruption of normal protein function, most likely from loss of phosphorylation. A similar result was expected for D1; however, an intriguing "constitutive on" phenotype was observed. In addition, the unusual feature of a separate HPt domain led us to compare the sequences surrounding the ysrS-ysrT junction in several Yersinia strains. In every strain examined, ysrT is a separate gene, leading to speculation that there is a functional advantage to YsrT being an independent protein.

Origins and diversification of a complex signal transduction system in prokaryotes

The molecular machinery that controls chemotaxis in bacteria is substantially more complex than any other signal transduction system in prokaryotes, and its origins and variability among living species are unknown. We found that this multiprotein "chemotaxis system" is present in most prokaryotic species and evolved from simpler two-component regulatory systems that control prokaryotic transcription. We discovered, through genomic analysis, signaling systems intermediate between two-component systems and chemotaxis systems. Evolutionary genomics established central and auxiliary components of the chemotaxis system. While tracing its evolutionary history, we also developed a classification scheme that revealed more than a dozen distinct classes of chemotaxis systems, enabling future predictive modeling of chemotactic behavior in unstudied species.

Coupling metabolism and chemotaxis-dependent behaviours by energy taxis receptors

Bacteria have evolved the ability to monitor changes in various physico-chemical parameters and to adapt their physiology and metabolism by implementing appropriate cellular responses to these changes. Energy taxis is a metabolism-dependent form of taxis and is the directed movement of motile bacteria in gradients of physico-chemical parameters that affect metabolism. Energy taxis has been described in diverse bacterial species and several dedicated energy sensors have been identified. The molecular mechanism of energy taxis has not been studied in as much detail as chemotaxis, but experimental evidence indicates that this behaviour differs from metabolism-independent taxis only by the presence of dedicated energy taxis receptors. Energy taxis receptors perceive changes in energy-related parameters, including signals related to the redox and/or intracellular energy status of the cell. The best-characterized energy taxis receptors are those that sense the redox state of the electron transport chain via non-covalently bound FAD cofactors. Other receptors shown to mediate energy taxis lack any recognizable redox cofactor or conserved energy-sensing motif, and some have been suggested to monitor changes in the proton motive force. The exact energy-sensing mechanism(s) involved are yet to be elucidated for most of these energy sensors. By monitoring changes in energy-related parameters, energy taxis receptors allow cells to couple motility behaviour with metabolism under diverse environmental conditions. Energy taxis receptors thus provide fruitful models to decipher how cells integrate sensory behaviours with metabolic activities.