Inference of the transcriptional regulatory network in Staphylococcus aureus by integration of experimental and genomics-based evidence

Bioinformatics in Russia: history and present-day landscape

Bioinformatics has become an interdisciplinary subject due to its universal role in molecular biology research. The current status of Russia's bioinformatics research in Russia is not known. Here, we review the history of bioinformatics in Russia, present the current landscape, and highlight future directions and challenges. Bioinformatics research in Russia is driven by four major industries: information technology, pharmaceuticals, biotechnology, and agriculture. Over the past three decades, despite a delayed start, the field has gained momentum, especially in protein and nucleic acid research. Dedicated and shared centers for genomics, proteomics, and bioinformatics are active in different regions of Russia. Present-day bioinformatics in Russia is characterized by research issues related to genetics, metagenomics, OMICs, medical informatics, computational biology, environmental informatics, and structural bioinformatics. Notable developments are in the fields of software (tools, algorithms, and pipelines), use of high computation power (e.g. by the Siberian Supercomputer Center), and large-scale sequencing projects (the sequencing of 100 000 human genomes). Government funding is increasing, policies are being changed, and a National Genomic Information Database is being established. An increased focus on eukaryotic genome sequencing, the development of a common place for developers and researchers to share tools and data, and the use of biological modeling, machine learning, and biostatistics are key areas for future focus. Universities and research institutes have started to implement bioinformatics modules. A critical mass of bioinformaticians is essential to catch up with the global pace in the discipline.

LysR-type transcriptional regulators: state of the art

The LysR-type transcriptional regulators (LTTRs) are DNA-binding proteins present in bacteria, archaea, and in algae. Knowledge about their distribution, abundance, evolution, structural organization, transcriptional regulation, fundamental roles in free life, pathogenesis, and bacteria-plant interaction has been generated. This review focuses on these aspects and provides a current picture of LTTR biology.

Histidine transport is essential for the growth of Staphylococcus aureus at low pH

Staphylococcus aureus is an opportunistic pathogen capable of causing many different human diseases. During colonization and infection, S. aureus will encounter a range of hostile environments, including acidic conditions such as those found on the skin and within macrophages. However, little is known about the mechanisms that S. aureus uses to detect and respond to low pH. Here, we employed a transposon sequencing approach to determine on a genome-wide level the genes required or detrimental for growth at low pH. We identified 31 genes that were essential for the growth of S. aureus at pH 4.5 and confirmed the importance of many of them through follow up experiments using mutant strains inactivated for individual genes. Most of the genes identified code for proteins with functions in cell wall assembly and maintenance. These data suggest that the cell wall has a more important role than previously appreciated in promoting bacterial survival when under acid stress. We also identified several novel processes previously not linked to the acid stress response in S. aureus. These include aerobic respiration and histidine transport, the latter by showing that one of the most important genes, SAUSA300_0846, codes for a previously uncharacterized histidine transporter. We further show that under acid stress, the expression of the histidine transporter gene is increased in WT S. aureus. In a S. aureus SAUSA300_0846 mutant strain expression of the histidine biosynthesis genes is induced under acid stress conditions allowing the bacteria to maintain cytosolic histidine levels. This strain is, however, unable to maintain its cytosolic pH to the same extent as a WT strain, revealing an important function specifically for histidine transport in the acid stress response of S. aureus.

A genetic regulatory see-saw of biofilm and virulence in MRSA pathogenesis

Staphylococcus aureus is one of the most common opportunistic human pathogens causing several infectious diseases. Ever since the emergence of the first methicillin-resistant Staphylococcus aureus (MRSA) strain decades back, the organism has been a major cause of hospital-acquired infections (HA-MRSA). The spread of this pathogen across the community led to the emergence of a more virulent subtype of the strain, i.e., Community acquired Methicillin resistant Staphylococcus aureus (CA-MRSA). Hence, WHO has declared Staphylococcus aureus as a high-priority pathogen. MRSA pathogenesis is remarkable because of the ability of this "superbug" to form robust biofilm both in vivo and in vitro by the formation of polysaccharide intercellular adhesin (PIA), extracellular DNA (eDNA), wall teichoic acids (WTAs), and capsule (CP), which are major components that impart stability to a biofilm. On the other hand, secretion of a diverse array of virulence factors such as hemolysins, leukotoxins, enterotoxins, and Protein A regulated by agr and sae two-component systems (TCS) aids in combating host immune response. The up- and downregulation of adhesion genes involved in biofilm formation and genes responsible for synthesizing virulence factors during different stages of infection act as a genetic regulatory see-saw in the pathogenesis of MRSA. This review provides insight into the evolution and pathogenesis of MRSA infections with a focus on genetic regulation of biofilm formation and virulence factors secretion.

The gene regulatory network of Staphylococcus aureus ST239-SCCmecIII strain Bmb9393 and assessment of genes associated with the biofilm in diverse backgrounds

Introduction

Staphylococcus aureus is one of the most prevalent and relevant pathogens responsible for a wide spectrum of hospital-associated or community-acquired infections. In addition, methicillin-resistant Staphylococcus aureus may display multidrug resistance profiles that complicate treatment and increase the mortality rate. The ability to produce biofilm, particularly in device-associated infections, promotes chronic and potentially more severe infections originating from the primary site. Understanding the complex mechanisms involved in planktonic and biofilm growth is critical to identifying regulatory connections and ways to overcome the global health problem of multidrug-resistant bacteria.

Methods

In this work, we apply literature-based and comparative genomics approaches to reconstruct the gene regulatory network of the high biofilm-producing strain Bmb9393, belonging to one of the highly disseminating successful clones, the Brazilian epidemic clone. To the best of our knowledge, we describe for the first time the topological properties and network motifs for the Staphylococcus aureus pathogen. We performed this analysis using the ST239-SCCmecIII Bmb9393 strain. In addition, we analyzed transcriptomes available in the literature to construct a set of genes differentially expressed in the biofilm, covering different stages of the biofilms and genetic backgrounds of the strains.

Results and discussion

The Bmb9393 gene regulatory network comprises 1,803 regulatory interactions between 64 transcription factors and the non-redundant set of 1,151 target genes with the inclusion of 19 new regulons compared to the N315 transcriptional regulatory network published in 2011. In the Bmb9393 network, we found 54 feed-forward loop motifs, where the most prevalent were coherent type 2 and incoherent type 2. The non-redundant set of differentially expressed genes in the biofilm consisted of 1,794 genes with functional categories relevant for adaptation to the variable microenvironments established throughout the biofilm formation process. Finally, we mapped the set of genes with altered expression in the biofilm in the Bmb9393 gene regulatory network to depict how different growth modes can alter the regulatory systems. The data revealed 45 transcription factors and 876 shared target genes. Thus, the gene regulatory network model provided represents the most up-to-date model for Staphylococcus aureus, and the set of genes altered in the biofilm provides a global view of their influence on biofilm formation from distinct experimental perspectives and different strain backgrounds.

Evolution of transcriptional regulation of histidine metabolism in Gram-positive bacteria

BACKGROUND

The histidine metabolism and transport (his) genes are controlled by a variety of RNA-dependent regulatory systems among diverse taxonomic groups of bacteria including T-box riboswitches in Firmicutes and Actinobacteria and RNA attenuators in Proteobacteria. Using a comparative genomic approach, we previously identified a novel DNA-binding transcription factor (named HisR) that controls the histidine metabolism genes in diverse Gram-positive bacteria from the Firmicutes phylum.

RESULTS

Here we report the identification of HisR-binding sites within the regulatory regions of the histidine metabolism and transport genes in 395 genomes representing the Bacilli, Clostridia, Negativicutes, and Tissierellia classes of Firmicutes, as well as in 97 other HisR-encoding genomes from the Actinobacteria, Proteobacteria, and Synergistetes phyla. HisR belongs to the TrpR family of transcription factors, and their predicted DNA binding motifs have a similar 20-bp palindromic structure but distinct lineage-specific consensus sequences. The predicted HisR-binding motif was validated in vitro using DNA binding assays with purified protein from the human gut bacterium Ruminococcus gnavus. To fill a knowledge gap in the regulation of histidine metabolism genes in Firmicutes genomes that lack a hisR repressor gene, we systematically searched their upstream regions for potential RNA regulatory elements. As result, we identified 158 T-box riboswitches preceding the histidine biosynthesis and/or transport genes in 129 Firmicutes genomes. Finally, novel candidate RNA attenuators were identified upstream of the histidine biosynthesis operons in six species from the Bacillus cereus group, as well as in five Eubacteriales and six Erysipelotrichales species.

CONCLUSIONS

The obtained distribution of the HisR transcription factor and two RNA-mediated regulatory mechanisms for histidine metabolism genes across over 600 species of Firmicutes is discussed from functional and evolutionary points of view.

Homology-based reconstruction of regulatory networks for bacterial and archaeal genomes

Gene regulation is a key process for all microorganisms, as it allows them to adapt to different environmental stimuli. However, despite the relevance of gene expression control, for only a handful of organisms is there related information about genome regulation. In this work, we inferred the gene regulatory networks (GRNs) of bacterial and archaeal genomes by comparisons with six organisms with well-known regulatory interactions. The references we used are: Escherichia coli K-12 MG1655, Bacillus subtilis 168, Mycobacterium tuberculosis, Pseudomonas aeruginosa PAO1, Salmonella enterica subsp. enterica serovar typhimurium LT2, and Staphylococcus aureus N315. To this end, the inferences were achieved in two steps. First, the six model organisms were contrasted in an all-vs-all comparison of known interactions based on Transcription Factor (TF)-Target Gene (TG) orthology relationships and Transcription Unit (TU) assignments. In the second step, we used a guilt-by-association approach to infer the GRNs for 12,230 bacterial and 649 archaeal genomes based on TF-TG orthology relationships of the six bacterial models determined in the first step. Finally, we discuss examples to show the most relevant results obtained from these inferences. A web server with all the predicted GRNs is available at https://regulatorynetworks.unam.mx/ or http://132.247.46.6/.

From a large-scale genomic analysis of insertion sequences to insights into their regulatory roles in prokaryotes

Background

Insertion sequences (ISs) are mobile repeat sequences and most of them can copy themselves to new host genome locations, leading to genome plasticity and gene regulation in prokaryotes. In this study, we present functional and evolutionary relationships between IS and neighboring genes in a large-scale comparative genomic analysis.

Results

IS families were located in all prokaryotic phyla, with preferential occurrence of IS3, IS4, IS481, and IS5 families in Alpha-, Beta-, and Gammaproteobacteria, Actinobacteria and Firmicutes as well as in eukaryote host-associated organisms and autotrophic opportunistic pathogens. We defined the concept of the IS-Gene couple (IG), which allowed to highlight the functional and regulatory impacts of an IS on the closest gene. Genes involved in transcriptional regulation and transport activities were found overrepresented in IG. In particular, major facilitator superfamily (MFS) transporters, ATP-binding proteins and transposases raised as favorite neighboring gene functions of IS hotspots. Then, evolutionary conserved IS-Gene sets across taxonomic lineages enabled the classification of IS-gene couples into phylum, class-to-genus, and species syntenic IS-Gene couples. The IS5, IS21, IS4, IS607, IS91, ISL3 and IS200 families displayed two to four times more ISs in the phylum and/or class-to-genus syntenic IGs compared to other IS families. This indicates that those families were probably inserted earlier than others and then subjected to horizontal transfer, transposition and deletion events over time. In phylum syntenic IG category, Betaproteobacteria, Crenarchaeota, Calditrichae, Planctomycetes, Acidithiobacillia and Cyanobacteria phyla act as IS reservoirs for other phyla, and neighboring gene functions are mostly related to transcriptional regulators. Comparison of IS occurrences with predicted regulatory motifs led to ~ 26.5% of motif-containing ISs with 2 motifs per IS in average. These results, concomitantly with short IS-Gene distances, suggest that those ISs would interfere with the expression of neighboring genes and thus form strong candidates for an adaptive pairing.

Conclusions

All together, our large-scale study provide new insights into the IS genetic context and strongly suggest their regulatory roles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08678-3.

Limiting protease production plays a key role in the pathogenesis of the divergent clinical isolates of Staphylococcus aureus LAC and UAMS-1

Using the USA300, methicillin-resistant Staphylococcus aureus strain LAC, we previously examined the impact of regulatory mutations implicated in biofilm formation on protease production and virulence in a murine sepsis model. Here we examined the impact of these mutations in the USA200, methicillin-sensitive strain UAMS-1. Mutation of agr, mgrA, rot, sarA and sigB attenuated the virulence of UAMS-1. A common characteristic of codY, rot, sigB, and sarA mutants was increased protease production, with mutation of rot having the least impact followed by mutation of codY, sigB and sarA, respectively. Protein A was undetectable in conditioned medium from all four mutants, while extracellular nuclease was only present in the proteolytically cleaved NucA form. The abundance of high molecular weight proteins was reduced in all four mutants. Biofilm formation was reduced in codY, sarA and sigB mutants, but not in the rot mutant. Eliminating protease production partially reversed these phenotypes and enhanced biofilm formation. This was also true in LAC codY, rot, sarA and sigB mutants. Eliminating protease production enhanced the virulence of LAC and UAMS-1 sarA, sigB and rot mutants in a murine sepsis model but did not significantly impact the virulence of the codY mutant in either strain. Nevertheless, these results demonstrate that repressing protease production plays an important role in defining critical phenotypes in diverse clinical isolates of S. aureus and that Rot, SigB and SarA play critical roles in this regard.

Revealing 29 sets of independently modulated genes in Staphylococcus aureus, their regulators, and role in key physiological response

Staphylococcus aureus infections impose an immense burden on the healthcare system. To establish a successful infection in a hostile host environment, S. aureus must coordinate its gene expression to respond to a wide array of challenges. This balancing act is largely orchestrated by the transcriptional regulatory network. Here, we present a model of 29 independently modulated sets of genes that form the basis for a segment of the transcriptional regulatory network in clinical USA300 strains of S. aureus. Using this model, we demonstrate the concerted role of various cellular systems (e.g., metabolism, virulence, and stress response) underlying key physiological responses, including response during blood infection.

The ability of Staphylococcus aureus to infect many different tissue sites is enabled, in part, by its transcriptional regulatory network (TRN) that coordinates its gene expression to respond to different environments. We elucidated the organization and activity of this TRN by applying independent component analysis to a compendium of 108 RNA-sequencing expression profiles from two S. aureus clinical strains (TCH1516 and LAC). ICA decomposed the S. aureus transcriptome into 29 independently modulated sets of genes (i-modulons) that revealed: 1) High confidence associations between 21 i-modulons and known regulators; 2) an association between an i-modulon and σS, whose regulatory role was previously undefined; 3) the regulatory organization of 65 virulence factors in the form of three i-modulons associated with AgrR, SaeR, and Vim-3; 4) the roles of three key transcription factors (CodY, Fur, and CcpA) in coordinating the metabolic and regulatory networks; and 5) a low-dimensional representation, involving the function of few transcription factors of changes in gene expression between two laboratory media (RPMI, cation adjust Mueller Hinton broth) and two physiological media (blood and serum). This representation of the TRN covers 842 genes representing 76% of the variance in gene expression that provides a quantitative reconstruction of transcriptional modules in S. aureus, and a platform enabling its full elucidation.

Inhibition of Rho Activity Increases Expression of SaeRS-Dependent Virulence Factor Genes in Staphylococcus aureus, Showing a Link between Transcription Termination, Antibiotic Action, and Virulence

Staphylococcus aureus causes various diseases ranging from skin and soft tissue infections to life-threatening infections. Adaptation to the different host niches is controlled by a complex network of transcriptional regulators. Global profiling of condition-dependent transcription revealed adaptation of S. aureus HG001 at the levels of transcription initiation and termination. In particular, deletion of the gene encoding the Rho transcription termination factor triggered a remarkable overall increase in antisense transcription and gene expression changes attributable to indirect regulatory effects. The goal of the present study was a detailed comparative analysis of S. aureus HG001 and its isogenic rho deletion mutant. Proteome analysis revealed significant differences in cellular and extracellular protein profiles, most notably increased amounts of the proteins belonging to the SaeR regulon in the Rho-deficient strain. The SaeRS two-component system acts as a major regulator of virulence gene expression in staphylococci. Higher levels of SaeRS-dependent virulence factors such as adhesins, toxins, and immune evasion proteins in the rho mutant resulted in higher virulence in a murine bacteremia model, which was alleviated in a rho complemented strain. Inhibition of Rho activity by bicyclomycin, a specific inhibitor of Rho activity, also induced the expression of SaeRS-dependent genes, at both the mRNA and protein levels, to the same extent as observed in the rho mutant. Taken together, these findings indicate that activation of the Sae system in the absence of Rho is directly linked to Rho's transcription termination activity and establish a new link between antibiotic action and virulence gene expression in S. aureusIMPORTANCE The major human pathogen Staphylococcus aureus is a widespread commensal bacterium but also the most common cause of nosocomial infections. It adapts to the different host niches through a complex gene regulatory network. We show here that the Rho transcription termination factor, which represses pervasive antisense transcription in various bacteria, including S. aureus, plays a role in controlling SaeRS-dependent virulence gene expression. A Rho-deficient strain produces larger amounts of secreted virulence factors in vitro and shows increased virulence in mice. We also show that treatment of S. aureus with the antibiotic bicyclomycin, which inhibits Rho activity and is effective against Gram-negative bacteria, induces the same changes in the proteome as observed in the Rho-deficient strain. Our results reveal for the first time a link between transcription termination and virulence regulation in S. aureus, which implies a novel mechanism by which an antibiotic can modulate the expression of virulence factors.

Bacterial regulatory RNAs: complexity, function, and putative drug targeting

Over the past decade, RNA-deep sequencing has uncovered copious non-protein coding RNAs (npcRNAs) in bacteria. Many of them are key players in the regulation of gene expression, taking part in various regulatory circuits, such as metabolic responses to different environmental stresses, virulence, antibiotic resistance, and host-pathogen interactions. This has contributed to the high adaptability of bacteria to changing or even hostile environments. Their mechanisms include the regulation of transcriptional termination, modulation of translation, and alteration of messenger RNA (mRNA) stability, as well as protein sequestration. Here, the mechanisms of gene expression by regulatory bacterial npcRNAs are comprehensively reviewed and supplemented with well-characterized examples. This class of molecules and their mechanisms of action might be useful targets for the development of novel antibiotics.

Comparative Genomics Reveals the Regulatory Complexity of Bifidobacterial Arabinose and Arabino-Oligosaccharide Utilization

Members of the genus Bifidobacterium are common inhabitants of the human gastrointestinal tract. Previously it was shown that arabino-oligosaccharides (AOS) might act as prebiotics and stimulate the bifidobacterial growth in the gut. However, despite the rapid accumulation of genomic data, the precise mechanisms by which these sugars are utilized and associated transcription control still remain unclear. In the current study, we used a comparative genomic approach to reconstruct arabinose and AOS utilization pathways in over 40 bacterial species belonging to the Bifidobacteriaceae family. The results indicate that the gene repertoire involved in the catabolism of these sugars is highly diverse, and even phylogenetically close species may differ in their utilization capabilities. Using bioinformatics analysis we identified potential DNA-binding motifs and reconstructed putative regulons for the arabinose and AOS utilization genes in the Bifidobacteriaceae genomes. Six LacI-family transcriptional factors (named AbfR, AauR, AauU1, AauU2, BauR1 and BauR2) and a TetR-family regulator (XsaR) presumably act as local repressors for AOS utilization genes encoding various α- or β-L-arabinofuranosidases and predicted AOS transporters. The ROK-family regulator AraU and the LacI-family regulator AraQ control adjacent operons encoding putative arabinose transporters and catabolic enzymes, respectively. However, the AraQ regulator is universally present in all Bifidobacterium species including those lacking the arabinose catabolic genes araBDA, suggesting its control of other genes. Comparative genomic analyses of prospective AraQ-binding sites allowed the reconstruction of AraQ regulons and a proposed binary repression/activation mechanism. The conserved core of reconstructed AraQ regulons in bifidobacteria includes araBDA, as well as genes from the central glycolytic and fermentation pathways (pyk, eno, gap, tkt, tal, galM, ldh). The current study expands the range of genes involved in bifidobacterial arabinose/AOS utilization and demonstrates considerable variations in associated metabolic pathways and regulons. Detailed comparative and phylogenetic analyses allowed us to hypothesize how the identified reconstructed regulons evolved in bifidobacteria. Our findings may help to improve carbohydrate catabolic phenotype prediction and metabolic modeling, while it may also facilitate rational development of novel prebiotics.

Comparison of gene co-expression networks in Pseudomonas aeruginosa and Staphylococcus aureus reveals conservation in some aspects of virulence

Pseudomonas aeruginosa and Staphylococcus aureus are two evolutionary distant bacterial species that are frequently isolated from persistent infections such as chronic infectious wounds and severe lung infections in cystic fibrosis patients. To the best of our knowledge no comprehensive genome scale co-expression study has been already conducted on these two species and in most cases only the expression of very few genes has been the subject of investigation. In this study, in order to investigate the level of expressional conservation between these two species, using heterogeneous gene expression datasets the weighted gene co-expression network analysis (WGCNA) approach was applied to study both single and cross species genome scale co-expression patterns of these two species. Single species co-expression network analysis revealed that in P. aeruginosa, genes involved in quorum sensing (QS), iron uptake, nitrate respiration and type III secretion systems and in S. aureus, genes associated with the regulation of carbon metabolism, fatty acid-phospholipids metabolism and proteolysis represent considerable co-expression across a variety of experimental conditions. Moreover, the comparison of gene co-expression networks between P. aeruginosa and S. aureus was led to the identification of four co-expressed gene modules in both species totally consisting of 318 genes. Several genes related to two component signal transduction systems, small colony variants (SCVs) morphotype and protein complexes were found in the detected modules. We believe that targeting the key players among the identified co-expressed orthologous genes will be a potential intervention strategy to control refractory co-infections caused by these two bacterial species.

The SaeRS Two-Component System of Staphylococcus aureus

In the Gram-positive pathogenic bacterium Staphylococcus aureus, the SaeRS twocomponent system (TCS) plays a major role in controlling the production of over 20 virulence factors including hemolysins, leukocidins, superantigens, surface proteins, and proteases. The SaeRS TCS is composed of the sensor histidine kinase SaeS, response regulator SaeR, and two auxiliary proteins SaeP and SaeQ. Since its discovery in 1994, the sae locus has been studied extensively, and its contributions to staphylococcal virulence and pathogenesis have been well documented and understood; however, the molecular mechanism by which the SaeRS TCS receives and processes cognate signals is not. In this article, therefore, we review the literature focusing on the signaling mechanism and its interaction with other global regulators.

Heme Synthesis and Acquisition in Bacterial Pathogens

Bacterial pathogens require the iron-containing cofactor heme to cause disease. Heme is essential to the function of hemoproteins, which are involved in energy generation by the electron transport chain, detoxification of host immune effectors, and other processes. During infection, bacterial pathogens must synthesize heme or acquire heme from the host; however, host heme is sequestered in high-affinity hemoproteins. Pathogens have evolved elaborate strategies to acquire heme from host sources, particularly hemoglobin, and both heme acquisition and synthesis are important for pathogenesis. Paradoxically, excess heme is toxic to bacteria and pathogens must rely on heme detoxification strategies. Heme is a key nutrient in the struggle for survival between host and pathogen, and its study has offered significant insight into the molecular mechanisms of bacterial pathogenesis.

Regulatory Requirements for Staphylococcus aureus Nitric Oxide Resistance

The ability of Staphylococcus aureus to resist host innate immunity augments the severity and pervasiveness of its pathogenesis. Nitric oxide (NO˙) is an innate immune radical that is critical for the efficient clearance of a wide range of microbial pathogens. Exposure of microbes to NO˙ typically results in growth inhibition and induction of stress regulons. S. aureus, however, induces a metabolic state in response to NO˙ that allows for continued replication and precludes stress regulon induction. The regulatory factors mediating this distinctive response remain largely undefined. Here, we employ a targeted transposon screen and transcriptomics to identify and characterize five regulons essential for NO˙ resistance in S. aureus: three virulence regulons not formerly associated with NO˙ resistance, SarA, CodY, and Rot, as well as two regulons with established roles, Fur and SrrAB. We provide new insights into the contributions of Fur and SrrAB during NO˙ stress and show that the S. aureus ΔsarA mutant, the most sensitive of the newly identified mutants, exhibits metabolic dysfunction and widespread transcriptional dysregulation following NO˙ exposure. Altogether, our results broadly characterize the regulatory requirements for NO˙ resistance in S. aureus and suggest an intriguing overlap between the regulation of NO˙ resistance and virulence in this well-adapted human pathogen.

IMPORTANCE

The prolific human pathogen Staphylococcus aureus is uniquely capable of resisting the antimicrobial radical nitric oxide (NO˙), a crucial component of the innate immune response. However, a complete understanding of how S. aureus regulates an effective response to NO˙ is lacking. Here, we implicate three central virulence regulators, SarA, CodY, and Rot, as major players in the S. aureus NO˙ response. Additionally, we elaborate on the contribution of two regulators, SrrAB and Fur, already known to play a crucial role in S. aureus NO˙ resistance. Our study sheds light on a unique facet of S. aureus pathogenicity and demonstrates that the transcriptional response of S. aureus to NO˙ is highly pleiotropic and intrinsically tied to metabolism and virulence regulation.

Comparative genomics and evolution of transcriptional regulons in Proteobacteria

Comparative genomics approaches are broadly used for analysis of transcriptional regulation in bacterial genomes. In this work, we identified binding sites and reconstructed regulons for 33 orthologous groups of transcription factors (TFs) in 196 reference genomes from 21 taxonomic groups of Proteobacteria. Overall, we predict over 10 600 TF binding sites and identified more than 15 600 target genes for 1896 TFs constituting the studied orthologous groups of regulators. These include a set of orthologues for 21 metabolism-associated TFs from Escherichia coli and/or Shewanella that are conserved in five or more taxonomic groups and several additional TFs that represent non-orthologous substitutions of the metabolic regulators in some lineages of Proteobacteria. By comparing gene contents of the reconstructed regulons, we identified the core, taxonomy-specific and genome-specific TF regulon members and classified them by their metabolic functions. Detailed analysis of ArgR, TyrR, TrpR, HutC, HypR and other amino-acid-specific regulons demonstrated remarkable differences in regulatory strategies used by various lineages of Proteobacteria. The obtained genomic collection of in silico reconstructed TF regulons contains a large number of new regulatory interactions that await future experimental validation. The collection provides a framework for future evolutionary studies of transcriptional regulatory networks in Bacteria. It can be also used for functional annotation of putative metabolic transporters and enzymes that are abundant in the reconstructed regulons.

Staphylococcus aureus Transcriptome Architecture: From Laboratory to Infection-Mimicking Conditions

Staphylococcus aureus is a major pathogen that colonizes about 20% of the human population. Intriguingly, this Gram-positive bacterium can survive and thrive under a wide range of different conditions, both inside and outside the human body. Here, we investigated the transcriptional adaptation of S. aureus HG001, a derivative of strain NCTC 8325, across experimental conditions ranging from optimal growth in vitro to intracellular growth in host cells. These data establish an extensive repertoire of transcription units and non-coding RNAs, a classification of 1412 promoters according to their dependence on the RNA polymerase sigma factors SigA or SigB, and allow identification of new potential targets for several known transcription factors. In particular, this study revealed a relatively low abundance of antisense RNAs in S. aureus, where they overlap only 6% of the coding genes, and only 19 antisense RNAs not co-transcribed with other genes were found. Promoter analysis and comparison with Bacillus subtilis links the small number of antisense RNAs to a less profound impact of alternative sigma factors in S. aureus. Furthermore, we revealed that Rho-dependent transcription termination suppresses pervasive antisense transcription, presumably originating from abundant spurious transcription initiation in this A+T-rich genome, which would otherwise affect expression of the overlapped genes. In summary, our study provides genome-wide information on transcriptional regulation and non-coding RNAs in S. aureus as well as new insights into the biological function of Rho and the implications of spurious transcription in bacteria.

The major human pathogen Staphylococcus aureus can survive under a wide range of conditions, both inside and outside the human body. The goal of this study was to determine how S. aureus adapts to such different conditions and, additionally, we wanted to identify general factors governing the staphylococcal transcriptome architecture. Therefore, we performed a precise analysis of all RNA transcripts of S. aureus across experimental conditions ranging from in vitro growth in different media to internalization by eukaryotic host cells. We systematically mapped all transcription units, annotated non-coding RNAs, and assigned promoters controlled by particular RNA polymerase sigma factors and transcription factors. By a comparison with data available for the related Gram-positive bacterium Bacillus subtilis, we made key observations concerning the abundance and origin of antisense RNAs. Intriguingly, these findings support the view that many antisense RNAs in a bacterium like B. subtilis could be byproducts of spurious promoter recognition by condition-specific alternative sigma factors. We also report that the transcription termination factor Rho prevents widespread antisense transcription, presumably caused by pervasive transcription initiation in the A+T-rich genome of S. aureus. Altogether our study presents new perspectives on the biological significance of antisense and pervasive transcription in bacteria.

Transcriptional Regulation of Carbohydrate Utilization Pathways in the Bifidobacterium Genus

Bifidobacteria, which represent common commensals of mammalian gut, are believed to have positive effects on human health. The influence of certain non-digestible carbohydrates (and their use as so-called prebiotics) on growth and metabolic activity of bifidobacteria is of increasing interest; however, mechanisms of transcriptional control of carbohydrate metabolism are poorly understood in these species. We used a comparative genomics approach to reconstruct carbohydrate utilization pathways and transcriptional regulons in 10 Bifidobacterium genomes. Analysis of regulatory gene regions revealed candidate DNA motifs and reconstructed regulons for 268 transcription factors from the LacI, ROK, DeoR, AraC, GntR, and TetR families that form 64 orthologous groups of regulators. Most of the reconstructed regulons are local and control specific catabolic pathways for host- and diet-derived glycans and monosaccharides. Mosaic distributions of many of these local regulators across Bifidobacterium species correlate with distribution of corresponding catabolic pathways. In contrast, the maltose, galactose, sucrose, and fructose regulons, as well as a novel global LacI-family regulator that is predicted to control the central carbohydrate metabolism and arabinose catabolism genes, are universally present in all 10 studied bifidobacteria. A novel group of TetR-family regulators presumably controls the glucoside and galactoside utilization pathways. Paralogs of the ribose repressor RbsR control the pyrimidine nucleoside utilization genes. Multiple paralogs of the maltose regulator MalR co-regulate large sets of genes involved in maltodextrin utilization. The inferred metabolic regulons provide new insights on diverse carbohydrate utilization networks in bifidobacteria that can be employed in metabolic modeling, phenotype prediction and the rational development of novel prebiotics.

Comparative genomics of pyridoxal 5'-phosphate-dependent transcription factor regulons in Bacteria

The MocR-subfamily transcription factors (MocR-TFs) characterized by the GntR-family DNA-binding domain and aminotransferase-like sensory domain are broadly distributed among certain lineages of Bacteria. Characterized MocR-TFs bind pyridoxal 5'-phosphate (PLP) and control transcription of genes involved in PLP, gamma aminobutyric acid (GABA) and taurine metabolism via binding specific DNA operator sites. To identify putative target genes and DNA binding motifs of MocR-TFs, we performed comparative genomics analysis of over 250 bacterial genomes. The reconstructed regulons for 825 MocR-TFs comprise structural genes from over 200 protein families involved in diverse biological processes. Using the genome context and metabolic subsystem analysis we tentatively assigned functional roles for 38 out of 86 orthologous groups of studied regulators. Most of these MocR-TF regulons are involved in PLP metabolism, as well as utilization of GABA, taurine and ectoine. The remaining studied MocR-TF regulators presumably control genes encoding enzymes involved in reduction/oxidation processes, various transporters and PLP-dependent enzymes, for example aminotransferases. Predicted DNA binding motifs of MocR-TFs are generally similar in each orthologous group and are characterized by two to four repeated sequences. Identified motifs were classified according to their structures. Motifs with direct and/or inverted repeat symmetry constitute the majority of inferred DNA motifs, suggesting preferable TF dimerization in head-to-tail or head-to-head configuration. The obtained genomic collection of in silico reconstructed MocR-TF motifs and regulons in Bacteria provides a basis for future experimental characterization of molecular mechanisms for various regulators in this family.

Cautions about the reliability of pairwise gene correlations based on expression data

BACKGROUND

Rapid growth in the availability of genome-wide transcript abundance levels through gene expression microarrays and RNAseq promises to provide deep biological insights into the complex, genome-wide transcriptional behavior of single-celled organisms. However, this promise has not yet been fully realized.

RESULTS

We find that computation of pairwise gene associations (correlation; mutual information) across a set of 2782 total genome-wide expression samples from six diverse bacteria produces unexpectedly large variation in estimates of pairwise gene association-regardless of the metric used, the organism under study, or the number and source of the samples. We pinpoint the cause to sampling bias. In particular, in repositories of expression data (e.g., Gene Expression Omnibus, GEO), many individual genes show small differences in absolute gene expression levels across the set of samples. We demonstrate that these small differences are due mainly to "noise" instead of "signal" attributable to environmental or genetic perturbations. We show that downstream analysis using gene expression levels of genes with small differences yields biased estimates of pairwise association.

CONCLUSIONS

We propose flagging genes with small differences in absolute, RMA-normalized, expression levels (e.g., standard deviation less than 0.5), as potentially yielding biased pairwise association metrics. This strategy has the potential to substantially improve the confidence in genome-wide conclusions about transcriptional behavior in bacterial organisms. Further work is needed to further refine strategies to identify genes with small difference in expression levels prior to computing gene-gene association metrics.

The Role of Two-Component Signal Transduction Systems in Staphylococcus aureus Virulence Regulation

Staphylococcus aureus is a versatile, opportunistic human pathogen that can asymptomatically colonize a human host but can also cause a variety of cutaneous and systemic infections. The ability of S. aureus to adapt to such diverse environments is reflected in the presence of complex regulatory networks fine-tuning metabolic and virulence gene expression. One of the most widely distributed mechanisms is the two-component signal transduction system (TCS) which allows a pathogen to alter its gene expression profile in response to environmental stimuli. The simpler TCSs consist of only a transmembrane histidine kinase (HK) and a cytosolic response regulator. S. aureus encodes a total of 16 conserved pairs of TCSs that are involved in diverse signalling cascades ranging from global virulence gene regulation (e.g. quorum sensing by the Agr system), the bacterial response to antimicrobial agents, cell wall metabolism, respiration and nutrient sensing. These regulatory circuits are often interconnected and affect each other's expression, thus fine-tuning staphylococcal gene regulation. This manuscript gives an overview of the current knowledge of staphylococcal environmental sensing by TCS and its influence on virulence gene expression and virulence itself. Understanding bacterial gene regulation by TCS can give major insights into staphylococcal pathogenicity and has important implications for knowledge-based drug design and vaccine formulation.

Comparative genomics and evolution of regulons of the LacI-family transcription factors

DNA-binding transcription factors (TFs) are essential components of transcriptional regulatory networks in bacteria. LacI-family TFs (LacI-TFs) are broadly distributed among certain lineages of bacteria. The majority of characterized LacI-TFs sense sugar effectors and regulate carbohydrate utilization genes. The comparative genomics approaches enable in silico identification of TF-binding sites and regulon reconstruction. To study the function and evolution of LacI-TFs, we performed genomics-based reconstruction and comparative analysis of their regulons. For over 1300 LacI-TFs from over 270 bacterial genomes, we predicted their cognate DNA-binding motifs and identified target genes. Using the genome context and metabolic subsystem analyses of reconstructed regulons, we tentatively assigned functional roles and predicted candidate effectors for 78 and 67% of the analyzed LacI-TFs, respectively. Nearly 90% of the studied LacI-TFs are local regulators of sugar utilization pathways, whereas the remaining 125 global regulators control large and diverse sets of metabolic genes. The global LacI-TFs include the previously known regulators CcpA in Firmicutes, FruR in Enterobacteria, and PurR in Gammaproteobacteria, as well as the three novel regulators—GluR, GapR, and PckR—that are predicted to control the central carbohydrate metabolism in three lineages of Alphaproteobacteria. Phylogenetic analysis of regulators combined with the reconstructed regulons provides a model of evolutionary diversification of the LacI protein family. The obtained genomic collection of in silico reconstructed LacI-TF regulons in bacteria is available in the RegPrecise database (http://regprecise.lbl.gov). It provides a framework for future structural and functional classification of the LacI protein family and identification of molecular determinants of the DNA and ligand specificity. The inferred regulons can be also used for functional gene annotation and reconstruction of sugar catabolic networks in diverse bacterial lineages.

Comparative genomics and functional analysis of rhamnose catabolic pathways and regulons in bacteria

L-rhamnose (L-Rha) is a deoxy-hexose sugar commonly found in nature. L-Rha catabolic pathways were previously characterized in various bacteria including Escherichia coli. Nevertheless, homology searches failed to recognize all the genes for the complete L-Rha utilization pathways in diverse microbial species involved in biomass decomposition. Moreover, the regulatory mechanisms of L-Rha catabolism have remained unclear in most species. A comparative genomics approach was used to reconstruct the L-Rha catabolic pathways and transcriptional regulons in the phyla Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, and Thermotogae. The reconstructed pathways include multiple novel enzymes and transporters involved in the utilization of L-Rha and L-Rha-containing polymers. Large-scale regulon inference using bioinformatics revealed remarkable variations in transcriptional regulators for L-Rha utilization genes among bacteria. A novel bifunctional enzyme, L-rhamnulose-phosphate aldolase (RhaE) fused to L-lactaldehyde dehydrogenase (RhaW), which is not homologous to previously characterized L-Rha catabolic enzymes, was identified in diverse bacteria including Chloroflexi, Bacilli, and Alphaproteobacteria. By using in vitro biochemical assays we validated both enzymatic activities of the purified recombinant RhaEW proteins from Chloroflexus aurantiacus and Bacillus subtilis. Another novel enzyme of the L-Rha catabolism, L-lactaldehyde reductase (RhaZ), was identified in Gammaproteobacteria and experimentally validated by in vitro enzymatic assays using the recombinant protein from Salmonella typhimurium. C. aurantiacus induced transcription of the predicted L-Rha utilization genes when L-Rha was present in the growth medium and consumed L-Rha from the medium. This study provided comprehensive insights to L-Rha catabolism and its regulation in diverse Bacteria.

Polysaccharides utilization in human gut bacterium Bacteroides thetaiotaomicron: comparative genomics reconstruction of metabolic and regulatory networks

Background

Bacteroides thetaiotaomicron, a predominant member of the human gut microbiota, is characterized by its ability to utilize a wide variety of polysaccharides using the extensive saccharolytic machinery that is controlled by an expanded repertoire of transcription factors (TFs). The availability of genomic sequences for multiple Bacteroides species opens an opportunity for their comparative analysis to enable characterization of their metabolic and regulatory networks.

Results

A comparative genomics approach was applied for the reconstruction and functional annotation of the carbohydrate utilization regulatory networks in 11 Bacteroides genomes. Bioinformatics analysis of promoter regions revealed putative DNA-binding motifs and regulons for 31 orthologous TFs in the Bacteroides. Among the analyzed TFs there are 4 SusR-like regulators, 16 AraC-like hybrid two-component systems (HTCSs), and 11 regulators from other families. Novel DNA motifs of HTCSs and SusR-like regulators in the Bacteroides have the common structure of direct repeats with a long spacer between two conserved sites.

Conclusions

The inferred regulatory network in B. thetaiotaomicron contains 308 genes encoding polysaccharide and sugar catabolic enzymes, carbohydrate-binding and transport systems, and TFs. The analyzed TFs control pathways for utilization of host and dietary glycans to monosaccharides and their further interconversions to intermediates of the central metabolism. The reconstructed regulatory network allowed us to suggest and refine specific functional assignments for sugar catabolic enzymes and transporters, providing a substantial improvement to the existing metabolic models for B. thetaiotaomicron. The obtained collection of reconstructed TF regulons is available in the RegPrecise database (http://regprecise.lbl.gov).

RegPrecise 3.0--a resource for genome-scale exploration of transcriptional regulation in bacteria

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in prokaryotes is one of the critical tasks of modern genomics. Bacteria from different taxonomic groups, whose lifestyles and natural environments are substantially different, possess highly diverged transcriptional regulatory networks. The comparative genomics approaches are useful for in silico reconstruction of bacterial regulons and networks operated by both transcription factors (TFs) and RNA regulatory elements (riboswitches).

Description

RegPrecise (http://regprecise.lbl.gov) is a web resource for collection, visualization and analysis of transcriptional regulons reconstructed by comparative genomics. We significantly expanded a reference collection of manually curated regulons we introduced earlier. RegPrecise 3.0 provides access to inferred regulatory interactions organized by phylogenetic, structural and functional properties. Taxonomy-specific collections include 781 TF regulogs inferred in more than 160 genomes representing 14 taxonomic groups of Bacteria. TF-specific collections include regulogs for a selected subset of 40 TFs reconstructed across more than 30 taxonomic lineages. Novel collections of regulons operated by RNA regulatory elements (riboswitches) include near 400 regulogs inferred in 24 bacterial lineages. RegPrecise 3.0 provides four classifications of the reference regulons implemented as controlled vocabularies: 55 TF protein families; 43 RNA motif families; ~150 biological processes or metabolic pathways; and ~200 effectors or environmental signals. Genome-wide visualization of regulatory networks and metabolic pathways covered by the reference regulons are available for all studied genomes. A separate section of RegPrecise 3.0 contains draft regulatory networks in 640 genomes obtained by an conservative propagation of the reference regulons to closely related genomes.

Conclusions

RegPrecise 3.0 gives access to the transcriptional regulons reconstructed in bacterial genomes. Analytical capabilities include exploration of: regulon content, structure and function; TF binding site motifs; conservation and variations in genome-wide regulatory networks across all taxonomic groups of Bacteria. RegPrecise 3.0 was selected as a core resource on transcriptional regulation of the Department of Energy Systems Biology Knowledgebase, an emerging software and data environment designed to enable researchers to collaboratively generate, test and share new hypotheses about gene and protein functions, perform large-scale analyses, and model interactions in microbes, plants, and their communities.

Aureolib - a proteome signature library: towards an understanding of staphylococcus aureus pathophysiology

Gel-based proteomics is a powerful approach to study the physiology of Staphylococcus aureus under various growth restricting conditions. We analyzed 679 protein spots from a reference 2-dimensional gel of cytosolic proteins of S. aureus COL by mass spectrometry resulting in 521 different proteins. 4,692 time dependent protein synthesis profiles were generated by exposing S. aureus to nine infection-related stress and starvation stimuli (H2O2, diamide, paraquat, NO, fermentation, nitrate respiration, heat shock, puromycin, mupirocin). These expression profiles are stored in an online resource called Aureolib (http://www.aureolib.de). Moreover, information on target genes of 75 regulators and regulatory elements were included in the database. Cross-comparisons of this extensive data collection of protein synthesis profiles using the tools implemented in Aureolib lead to the identification of stress and starvation specific marker proteins. Altogether, 226 protein synthesis profiles showed induction ratios of 2.5-fold or higher under at least one of the tested conditions with 157 protein synthesis profiles specifically induced in response to a single stimulus. The respective proteins might serve as marker proteins for the corresponding stimulus. By contrast, proteins whose synthesis was increased or repressed in response to more than four stimuli are rather exceptional. The only protein that was induced by six stimuli is the universal stress protein SACOL1759. Most strikingly, cluster analyses of synthesis profiles of proteins differentially synthesized under at least one condition revealed only in rare cases a grouping that correlated with known regulon structures. The most prominent examples are the GapR, Rex, and CtsR regulon. In contrast, protein synthesis profiles of proteins belonging to the CodY and σ(B) regulon are widely distributed. In summary, Aureolib is by far the most comprehensive protein expression database for S. aureus and provides an essential tool to decipher more complex adaptation processes in S. aureus during host pathogen interaction.

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.

Genomic reconstruction of transcriptional regulatory networks in lactic acid bacteria

Background

Genome scale annotation of regulatory interactions and reconstruction of regulatory networks are the crucial problems in bacterial genomics. The Lactobacillales order of bacteria collates various microorganisms having a large economic impact, including both human and animal pathogens and strains used in the food industry. Nonetheless, no systematic genome-wide analysis of transcriptional regulation has been previously made for this taxonomic group.

Results

A comparative genomics approach was used for reconstruction of transcriptional regulatory networks in 30 selected genomes of lactic acid bacteria. The inferred networks comprise regulons for 102 orthologous transcription factors (TFs), including 47 novel regulons for previously uncharacterized TFs. Numerous differences between regulatory networks of the Streptococcaceae and Lactobacillaceae groups were described on several levels. The two groups are characterized by substantially different sets of TFs encoded in their genomes. Content of the inferred regulons and structure of their cognate TF binding motifs differ for many orthologous TFs between the two groups. Multiple cases of non-orthologous displacements of TFs that control specific metabolic pathways were reported.

Conclusions

The reconstructed regulatory networks substantially expand the existing knowledge of transcriptional regulation in lactic acid bacteria. In each of 30 studied genomes the obtained regulatory network contains on average 36 TFs and 250 target genes that are mostly involved in carbohydrate metabolism, stress response, metal homeostasis and amino acids biosynthesis. The inferred networks can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. All reconstructed regulons are captured within the Streptococcaceae and Lactobacillaceae collections in the RegPrecise database (http://regprecise.lbl.gov).

Functional diversification of ROK-family transcriptional regulators of sugar catabolism in the Thermotogae phylum

Large and functionally heterogeneous families of transcription factors have complex evolutionary histories. What shapes specificities toward effectors and DNA sites in paralogous regulators is a fundamental question in biology. Bacteria from the deep-branching lineage Thermotogae possess multiple paralogs of the repressor, open reading frame, kinase (ROK) family regulators that are characterized by carbohydrate-sensing domains shared with sugar kinases. We applied an integrated genomic approach to study functions and specificities of regulators from this family. A comparative analysis of 11 Thermotogae genomes revealed novel mechanisms of transcriptional regulation of the sugar utilization networks, DNA-binding motifs and specific functions. Reconstructed regulons for seven groups of ROK regulators were validated by DNA-binding assays using purified recombinant proteins from the model bacterium Thermotoga maritima. All tested regulators demonstrated specific binding to their predicted cognate DNA sites, and this binding was inhibited by specific effectors, mono- or disaccharides from their respective sugar catabolic pathways. By comparing ligand-binding domains of regulators with structurally characterized kinases from the ROK family, we elucidated signature amino acid residues determining sugar-ligand regulator specificity. Observed correlations between signature residues and the sugar-ligand specificities provide the framework for structure functional classification of the entire ROK family.

Transcription factor family-based reconstruction of singleton regulons and study of the Crp/Fnr, ArsR, and GntR families in Desulfovibrionales genomes

Accurate detection of transcriptional regulatory elements is essential for high-quality genome annotation, metabolic reconstruction, and modeling of regulatory networks. We developed a computational approach for reconstruction of regulons operated by transcription factors (TFs) from large protein families and applied this novel approach to three TF families in 10 Desulfovibrionales genomes. Phylogenetic analyses of 125 regulators from the ArsR, Crp/Fnr, and GntR families revealed that 65% of these regulators (termed reference TFs) are well conserved in Desulfovibrionales, while the remaining 35% of regulators (termed singleton TFs) are species specific and show a mosaic distribution. For regulon reconstruction in the group of singleton TFs, the standard orthology-based approach was inefficient, and thus, we developed a novel approach based on the simultaneous study of all homologous TFs from the same family in a group of genomes. As a result, we identified binding for 21 singleton TFs and for all reference TFs in all three analyzed families. Within each TF family we observed structural similarities between DNA-binding motifs of different reference and singleton TFs. The collection of reconstructed regulons is available at the RegPrecise database (http://regprecise.lbl.gov/RegPrecise/Desulfovibrionales.jsp).

Evolution of transcriptional regulation in closely related bacteria

BACKGROUND

The exponential growth of the number of fully sequenced genomes at varying taxonomic closeness allows one to characterize transcriptional regulation using comparative-genomics analysis instead of time-consuming experimental methods. A transcriptional regulatory unit consists of a transcription factor, its binding site and a regulated gene. These units constitute a graph which contains so-called "network motifs", subgraphs of a given structure. Here we consider genomes of closely related Enterobacteriales and estimate the fraction of conserved network motifs and sites as well as positions under selection in various types of non-coding regions.

RESULTS

Using a newly developed technique, we found that the highest fraction of positions under selection, approximately 50%, was observed in synvergon spacers (between consecutive genes from the same strand), followed by ~45% in divergon spacers (common 5'-regions), and ~10% in convergon spacers (common 3'-regions). The fraction of selected positions in functional regions was higher, 60% in transcription factor-binding sites and ~45% in terminators and promoters. Small, but significant differences were observed between Escherichia coli and Salmonella enterica. This fraction is similar to the one observed in eukaryotes.The conservation of binding sites demonstrated some differences between types of regulatory units. In E. coli, strains the interactions of the type "local transcriptional factor gene" turned out to be more conserved in feed-forward loops (FFLs) compared to non-motif interactions. The coherent FFLs tend to be less conserved than the incoherent FFLs. A natural explanation is that the former imply functional redundancy.

CONCLUSIONS

A naïve hypothesis that FFL would be highly conserved turned out to be not entirely true: its conservation depends on its status in the transcriptional network and also from its usage. The fraction of positions under selection in intergenic regions of bacterial genomes is roughly similar to that of eukaryotes. Known regulatory sites explain 20±5% of selected positions.

The auxiliary protein complex SaePQ activates the phosphatase activity of sensor kinase SaeS in the SaeRS two-component system of Staphylococcus aureus

In bacterial two-component regulatory systems (TCSs), dephosphorylation of phosphorylated response regulators is essential for resetting the activated systems to the pre-activation state. However, in the SaeRS TCS, a major virulence TCS of Staphylococcus aureus, the mechanism for dephosphorylation of the response regulator SaeR has not been identified. Here we report that two auxiliary proteins from the sae operon, SaeP and SaeQ, form a protein complex with the sensor kinase SaeS and activate the sensor kinase's phosphatase activity. Efficient activation of the phosphatase activity required the presence of both SaeP and SaeQ. When SaeP and SaeQ were ectopically expressed, the expression of coagulase, a sae target with low affinity for phosphorylated SaeR, was greatly reduced, while the expression of alpha-haemolysin, a sae target with high affinity for phosphorylated SaeR, was not, demonstrating a differential effect of SaePQ on sae target gene expression. When expression of SaePQ was abolished, most sae target genes were induced at an elevated level. Since the expression of SaeP and SaeQ is induced by the SaeRS TCS, these results suggest that the SaeRS TCS returns to the pre-activation state by a negative feedback mechanism.

Wisdom of crowds for robust gene network inference

Reconstructing gene regulatory networks from high-throughput data is a long-standing problem. Through the DREAM project (Dialogue on Reverse Engineering Assessment and Methods), we performed a comprehensive blind assessment of over thirty network inference methods on Escherichia coli, Staphylococcus aureus, Saccharomyces cerevisiae, and in silico microarray data. We characterize performance, data requirements, and inherent biases of different inference approaches offering guidelines for both algorithm application and development. We observe that no single inference method performs optimally across all datasets. In contrast, integration of predictions from multiple inference methods shows robust and high performance across diverse datasets. Thereby, we construct high-confidence networks for E. coli and S. aureus, each comprising ~1700 transcriptional interactions at an estimated precision of 50%. We experimentally test 53 novel interactions in E. coli, of which 23 were supported (43%). Our results establish community-based methods as a powerful and robust tool for the inference of transcriptional gene regulatory networks.

RegPrecise web services interface: programmatic access to the transcriptional regulatory interactions in bacteria reconstructed by comparative genomics

Web services application programming interface (API) was developed to provide a programmatic access to the regulatory interactions accumulated in the RegPrecise database (http://regprecise.lbl.gov), a core resource on transcriptional regulation for the microbial domain of the Department of Energy (DOE) Systems Biology Knowledgebase. RegPrecise captures and visualize regulogs, sets of genes controlled by orthologous regulators in several closely related bacterial genomes, that were reconstructed by comparative genomics. The current release of RegPrecise 2.0 includes >1400 regulogs controlled either by protein transcription factors or by conserved ribonucleic acid regulatory motifs in >250 genomes from 24 taxonomic groups of bacteria. The reference regulons accumulated in RegPrecise can serve as a basis for automatic annotation of regulatory interactions in newly sequenced genomes. The developed API provides an efficient access to the RegPrecise data by a comprehensive set of 14 web service resources. The RegPrecise web services API is freely accessible at http://regprecise.lbl.gov/RegPrecise/services.jsp with no login requirements.