Evolution of transcriptional regulatory circuits in bacteria

Inter-species Transcriptomic Analysis Reveals a Constitutive Adaptation Against Oxidative Stress for the Highly Virulent Leptospira Species

Transcriptomic analyses across large scales of evolutionary distance have great potential to shed light on regulatory evolution but are complicated by difficulties in establishing orthology and limited availability of accessible software. We introduce here a method and a graphical user interface wrapper, called Annotator-RNAtor, for performing interspecies transcriptomic analysis and studying intragenus evolution. The pipeline uses third-party software to infer homologous genes in various species and highlight differences in the expression of the core-genes. To illustrate the methodology and demonstrate its usefulness, we focus on the emergence of the highly virulent Leptospira subclade known as P1+, which includes the causative agents of leptospirosis. Here, we expand on the genomic study through the comparison of transcriptomes between species from P1+ and their related P1- counterparts (low-virulent pathogens). In doing so, we shed light on differentially expressed pathways and focused on describing a specific example of adaptation based on a differential expression of PerRA-controlled genes. We showed that P1+ species exhibit higher expression of the katE gene, a well-known virulence determinant in pathogenic Leptospira species correlated with greater tolerance to peroxide. Switching PerRA alleles between P1+ and P1- species demonstrated that the lower repression of katE and greater tolerance to peroxide in P1+ species was solely controlled by PerRA and partly caused by a PerRA amino-acid permutation. Overall, these results demonstrate the strategic fit of the methodology and its ability to decipher adaptive transcriptomic changes, not observable by comparative genome analysis, that may have been implicated in the emergence of these pathogens.

The multifarious MerR family of transcriptional regulators

The MerR family of transcriptional regulators includes a variety of bacterial cytoplasmic proteins that respond to a wide range of signals, including toxins, metal ions, and endogenous metabolites. Its best-characterized members share similar structural and functional features with the family founder, the mercury sensor MerR, although most of them do not respond to metal ions. The group of "canonical" MerR homologs displays common molecular mechanisms for controlling the transcriptional activation of their target genes in response to inducer signals. This includes the recognition of distinctive operator sequences located at suboptimal σ70 -dependent promoters. Interestingly, an increasing number of proteins assigned to the MerR family based on their DNA-binding domain do not match in structure, sequence, or mode of action with any of the canonical MerR-like regulators. Here, we analyzed several members of the family, including this last group. Based on a phylogenetic analysis, and similarities in structural/functional features and position of their target operators relative to the promoter elements, we propose to assign these "atypical/divergent" MerR regulators to a phylogenetically separated group. These atypical/divergent homologs represent a new class of transcriptional regulators with novel regulatory mechanisms.

Nitrous oxide reduction by two partial denitrifying bacteria requires denitrification intermediates that cannot be respired

Denitrification is a form of anaerobic respiration wherein nitrate (NO3-) is sequentially reduced via nitrite (NO2-), nitric oxide, and nitrous oxide (N2O) to dinitrogen gas (N2) by four reductase enzymes. Partial denitrifying bacteria possess only one or some of these four reductases and use them as independent respiratory modules. However, it is unclear if partial denitrifiers sense and respond to denitrification intermediates outside of their reductase repertoire. Here, we tested the denitrifying capabilities of two purple nonsulfur bacteria, Rhodopseudomonas palustris CGA0092 and Rhodobacter capsulatus SB1003. Each had denitrifying capabilities that matched their genome annotation; CGA0092 reduced NO2- to N2, and SB1003 reduced N2O to N2. For each bacterium, N2O reduction could be used both for electron balance during growth on electron-rich organic compounds in light and for energy transformation via respiration in darkness. However, N2O reduction required supplementation with a denitrification intermediate, including those for which there was no associated denitrification enzyme. For CGA0092, NO3- served as a stable, non-catalyzable molecule that was sufficient to activate N2O reduction. Using a β-galactosidase reporter, we found that NO3- acted, at least in part, by stimulating N2O reductase gene expression. In SB1003, NO2- but not NO3- activated N2O reduction, but NO2- was slowly removed, likely by a promiscuous enzyme activity. Our findings reveal that partial denitrifiers can still be subject to regulation by denitrification intermediates that they cannot use.IMPORTANCEDenitrification is a form of microbial respiration wherein nitrate is converted via several nitrogen oxide intermediates into harmless dinitrogen gas. Partial denitrifying bacteria, which individually have some but not all denitrifying enzymes, can achieve complete denitrification as a community by cross-feeding nitrogen oxide intermediates. However, the last intermediate, nitrous oxide (N2O), is a potent greenhouse gas that often escapes, motivating efforts to understand and improve the efficiency of denitrification. Here, we found that at least some partial denitrifying N2O reducers can sense and respond to nitrogen oxide intermediates that they cannot otherwise use. The regulatory effects of nitrogen oxides on partial denitrifiers are thus an important consideration in understanding and applying denitrifying bacterial communities to combat greenhouse gas emissions.

Distinct ecological fitness factors coordinated by a conserved Escherichia coli regulator during systemic bloodstream infection

The ability of bacterial pathogens to adapt to host niches is driven by the carriage and regulation of genes that benefit pathogenic lifestyles. Genes that encode virulence or fitness-enhancing factors must be regulated in response to changing host environments to allow rapid response to challenges presented by the host. Furthermore, this process can be controlled by preexisting transcription factors (TFs) that acquire new roles in tailoring regulatory networks, specifically in pathogens. However, the mechanisms underlying this process are poorly understood. The highly conserved Escherichia coli TF YhaJ exhibits distinct genome-binding dynamics and transcriptome control in pathotypes that occupy different host niches, such as uropathogenic E. coli (UPEC). Here, we report that this important regulator is required for UPEC systemic survival during murine bloodstream infection (BSI). This advantage is gained through the coordinated regulation of a small regulon comprised of both virulence and metabolic genes. YhaJ coordinates activation of both Type 1 and F1C fimbriae, as well as biosynthesis of the amino acid tryptophan, by both direct and indirect mechanisms. Deletion of yhaJ or the individual genes under its control leads to attenuated survival during BSI. Furthermore, all three systems are up-regulated in response to signals derived from serum or systemic host tissue, but not urine, suggesting a niche-specific regulatory trigger that enhances UPEC fitness via pleiotropic mechanisms. Collectively, our results identify YhaJ as a pathotype-specific regulatory aide, enhancing the expression of key genes that are collectively required for UPEC bloodstream pathogenesis.

Optimal transcriptional regulation of dynamic bacterial responses to sudden drug exposures

Cellular responses to the presence of toxic compounds in their environment require prompt expression of the correct levels of the appropriate enzymes, which are typically regulated by transcription factors that control gene expression for the duration of the response. The characteristics of each response dictate the choice of regulatory parameters such as the affinity of the transcription factor to its binding sites and the strength of the promoters it regulates. Although much is known about the dynamics of cellular responses, we still lack a framework to understand how different regulatory strategies evolved in natural systems relate to the selective pressures acting in each particular case. Here, we analyze a dynamical model of a typical antibiotic response in bacteria, where a transcriptionally repressed enzyme is induced by a sudden exposure to the drug that it processes. We identify strategies of gene regulation that optimize this response for different types of selective pressures, which we define as a set of costs associated with the drug, enzyme, and repressor concentrations during the response. We find that regulation happens in a limited region of the regulatory parameter space. While responses to more costly (toxic) drugs favor the usage of strongly self-regulated repressors, responses where expression of enzyme is more costly favor the usage of constitutively expressed repressors. Only a very narrow range of selective pressures favor weakly self-regulated repressors. We use this framework to determine which costs and benefits are most critical for the evolution of a variety of natural cellular responses that satisfy the approximations in our model and to analyze how regulation is optimized in new environments with different demands.

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

Auxotrophic and prototrophic conditional genetic networks reveal the rewiring of transcription factors in Escherichia coli

Bacterial transcription factors (TFs) are widely studied in Escherichia coli. Yet it remains unclear how individual genes in the underlying pathways of TF machinery operate together during environmental challenge. Here, we address this by applying an unbiased, quantitative synthetic genetic interaction (GI) approach to measure pairwise GIs among all TF genes in E. coli under auxotrophic (rich medium) and prototrophic (minimal medium) static growth conditions. The resulting static and differential GI networks reveal condition-dependent GIs, widespread changes among TF genes in metabolism, and new roles for uncharacterized TFs (yjdC, yneJ, ydiP) as regulators of cell division, putrescine utilization pathway, and cold shock adaptation. Pan-bacterial conservation suggests TF genes with GIs are co-conserved in evolution. Together, our results illuminate the global organization of E. coli TFs, and remodeling of genetic backup systems for TFs under environmental change, which is essential for controlling the bacterial transcriptional regulatory circuits.

The bacterium E. coli has around 300 transcriptional factors, but the functions of many of them, and the interactions between their respective regulatory networks, are unclear. Here, the authors study genetic interactions among all transcription factor genes in E. coli, revealing condition-dependent interactions and roles for uncharacterized transcription factors.

Evolution of a cis-Acting SNP That Controls Type VI Secretion in Vibrio cholerae

Mutations in regulatory mechanisms that control gene expression contribute to phenotypic diversity and thus facilitate the adaptation of microbes and other organisms to new niches. Comparative genomics can be used to infer rewiring of regulatory architecture based on large effect mutations like loss or acquisition of transcription factors but may be insufficient to identify small changes in noncoding, intergenic DNA sequence of regulatory elements that drive phenotypic divergence. In human-derived Vibrio cholerae, the response to distinct chemical cues triggers production of multiple transcription factors that can regulate the type VI secretion system (T6), a broadly distributed weapon for interbacterial competition. However, to date, the signaling network remains poorly understood because no regulatory element has been identified for the major T6 locus. Here we identify a conserved cis-acting single nucleotide polymorphism (SNP) controlling T6 transcription and activity. Sequence alignment of the T6 regulatory region from diverse V. cholerae strains revealed conservation of the SNP that we rewired to interconvert V. cholerae T6 activity between chitin-inducible and constitutive states. This study supports a model of pathogen evolution through a noncoding cis-regulatory mutation and preexisting, active transcription factors that confers a different fitness advantage to tightly regulated strains inside a human host and unfettered strains adapted to environmental niches. IMPORTANCE Organisms sense external cues with regulatory circuits that trigger the production of transcription factors, which bind specific DNA sequences at promoters ("cis" regulatory elements) to activate target genes. Mutations of transcription factors or their regulatory elements create phenotypic diversity, allowing exploitation of new niches. Waterborne pathogen Vibrio cholerae encodes the type VI secretion system "nanoweapon" to kill competitor cells when activated. Despite identification of several transcription factors, no regulatory element has been identified in the promoter of the major type VI locus, to date. Combining phenotypic, genetic, and genomic analysis of diverse V. cholerae strains, we discovered a single nucleotide polymorphism in the type VI promoter that switches its killing activity between a constitutive state beneficial outside hosts and an inducible state for constraint in a host. Our results support a role for noncoding DNA in adaptation of this pathogen.

Disrupting autorepression circuitry generates "open-loop lethality" to yield escape-resistant antiviral agents

Across biological scales, gene-regulatory networks employ autorepression (negative feedback) to maintain homeostasis and minimize failure from aberrant expression. Here, we present a proof of concept that disrupting transcriptional negative feedback dysregulates viral gene expression to therapeutically inhibit replication and confers a high evolutionary barrier to resistance. We find that nucleic-acid decoys mimicking cis-regulatory sites act as "feedback disruptors," break homeostasis, and increase viral transcription factors to cytotoxic levels (termed "open-loop lethality"). Feedback disruptors against herpesviruses reduced viral replication >2-logs without activating innate immunity, showed sub-nM IC50, synergized with standard-of-care antivirals, and inhibited virus replication in mice. In contrast to approved antivirals where resistance rapidly emerged, no feedback-disruptor escape mutants evolved in long-term cultures. For SARS-CoV-2, disruption of a putative feedback circuit also generated open-loop lethality, reducing viral titers by >1-log. These results demonstrate that generating open-loop lethality, via negative-feedback disruption, may yield a class of antimicrobials with a high genetic barrier to resistance.

Regulatory Evolution of the phoH Ancestral Gene in Salmonella enterica Serovar Typhimurium

One important event for the divergence of Salmonella from Escherichia coli was the acquisition by horizontal transfer of the Salmonella pathogenicity island 1 (SPI-1), containing genes required for the invasion of host cells by Salmonella. HilD is an AraC-like transcriptional regulator in SPI-1 that induces the expression of the SPI-1 and many other acquired virulence genes located in other genomic regions of Salmonella. Additionally, HilD has been shown to positively control the expression of some ancestral genes (also present in E. coli and other bacteria), including phoH. In this study, we determined that both the gain of HilD and cis-regulatory evolution led to the integration of the phoH gene into the HilD regulon. Our results indicate that a HilD-binding sequence was generated in the regulatory region of the S. enterica serovar Typhimurium phoH gene, which mediates the activation of promoter 1 of this gene under SPI-1-inducing conditions. Furthermore, we found that repression by H-NS, a histone-like protein, was also adapted on the S. Typhimurium phoH gene and that HilD activates the expression of this gene in part by antagonizing H-NS. Additionally, our results revealed that the expression of the S. Typhmurium phoH gene is also activated in response to low phosphate but independently of the PhoB/R two-component system, known to regulate the E. coli phoH gene in response to low phosphate. Thus, our results indicate that cis-regulatory evolution has played a role in the expansion of the HilD regulon and illustrate the phenomenon of differential regulation of ortholog genes. IMPORTANCE Two mechanisms mediating differentiation of bacteria are well known: acquisition of genes by horizontal transfer events and mutations in coding DNA sequences. In this study, we found that the phoH ancestral gene is differentially regulated between Salmonella Typhimurium and Escherichia coli, two closely related bacterial species. Our results indicate that this differential regulation was generated by mutations in the regulatory sequence of the S. Typhimurium phoH gene and by the acquisition by S. Typhimurium of foreign DNA encoding the transcriptional regulator HilD. Thus, our results, together with those from an increasing number of studies, indicate that cis-regulatory evolution can lead to the rewiring and reprogramming of transcriptional regulation, which also plays an important role in the divergence of bacteria through time.

The ancestral stringent response potentiator, DksA has been adapted throughout Salmonella evolution to orchestrate the expression of metabolic, motility, and virulence pathways

DksA is a conserved RNA polymerase-binding protein known to play a key role in the stringent response of proteobacteria species, including many gastrointestinal pathogens. Here, we used RNA-sequencing of Escherichia coli, Salmonella bongori and Salmonella enterica serovar Typhimurium, together with phenotypic comparison to study changes in the DksA regulon, during Salmonella evolution. Comparative RNA-sequencing showed that under non-starved conditions, DksA controls the expression of 25%, 15%, and 20% of the E. coli, S. bongori, and S. enterica genes, respectively, indicating that DksA is a pleiotropic regulator, expanding its role beyond the canonical stringent response. We demonstrate that DksA is required for the growth of these three enteric bacteria species in minimal medium and controls the expression of the TCA cycle, glycolysis, pyrimidine biosynthesis, and quorum sensing. Interestingly, at multiple steps during Salmonella evolution, the type I fimbriae and various virulence genes encoded within SPIs 1, 2, 4, 5, and 11 have been transcriptionally integrated under the ancestral DksA regulon. Consequently, we show that DksA is necessary for host cells invasion by S. Typhimurium and S. bongori and for intracellular survival of S. Typhimurium in bone marrow-derived macrophages (BMDM). Moreover, we demonstrate regulatory inversion of the conserved motility-chemotaxis regulon by DksA, which acts as a negative regulator in E. coli, but activates this pathway in S. bongori and S. enterica. Overall, this study demonstrates the regulatory assimilation of multiple horizontally acquired virulence genes under the DksA regulon and provides new insights into the evolution of virulence genes regulation in Salmonella spp.

Development of a bacterial regulatory motif database

O b j e c t i v e s . The amount of data generated by modern methods of high-throughput sequencing is such that their analysis is performed mainly in automatic mode. In particular, the use of newly decoded genomic sequences is possible only after the annotation of functional elements of the genome, which, as a rule, is performed by automatic pipelines. Such annotation pipelines do a good job to identify the genes, but none of them annotate regulatory elements. Without these elements it is not possible to understand when and how genes can be expressed. Information on the regulatory elements of bacteria is collected in several specialized databases (RegulonDB, CollecTF, Prodoric2, etc.), however, only a part of this information can be used for annotation of regulatory elements, and only for a very limited range of bacteria. Previously, we proposed a clear formal criterion for applying regulatory information to any bacterial genome. Such a criterion is the CR tag, a sequence of amino acid residues of a transcriptional regulator that specifically contacts the nitrogenous bases of regulatory element in genomic DNA. The mathematical model of a regulatory element (motif) associated with a CR tag can be correctly applied to annotate similar elements in any genomes encoding a transcriptional regulator with an identical CR tag. The accumulation of motifs associated with CR tags raised the question of their ordered storage for the convenience of subsequent use in the annotation of genomic sequences. Since no one of well-known databases uses the concept of CR tags, a new database ought to be developed. Thus, the goal of this work is to create a database with information about bacterial transcription factors and DNA sequences recognized by them, suitable for annotation of regulatory sequences in bacterial genomes.

M e t h o d s .  Infological  modeling  of  the  subject  area  was  carried  out  using  the  IDEF1X  methodology. The database was developed using the Microsoft SQL Server DBMS. A cross-platform application for importing data into a database is written in C++ using Qt technology.

Re s u l t s . As a result of the study of the subject area, a relational data model was developed and implemented in the Microsoft SQL Server DBMS, which allows holistic storage of information about accumulated transcription regulation motifs in bacteria, including information about the publications confirming their correctness. To automate the process of entering accumulated data, a cross-platform application was developed for importing structured data on transcription factors.

Co n c l u s i o n .  The  main difference of  the  developed database is  the  use  of  CR-tag  concept. Records of mathematical models of regulatory elements (motifs) in the database are associated with a CR tag and, therefore, can be correctly used to annotate similar elements in any genomes encoding a transcriptional regulator with an identical CR tag. The developed database will provide structured and holistic data storage, as well as their quick search when used in the pipeline for automatic annotation of regulatory elements in bacterial genomic sequences.

Coenzyme B12 -dependent and independent photoregulation of carotenogenesis across Myxococcales

Light-induced carotenogenesis in Myxococcus xanthus is controlled by the B₁₂ -based CarH repressor and photoreceptor, and by a separate intricate pathway involving singlet oxygen, the B₁₂ -independent CarH paralog CarA and various other proteins, some eukaryotic-like. Whether other myxobacteria conserve these pathways and undergo photoregulated carotenogenesis is unknown. Here, comparative analyses across 27 Myxococcales genomes identified carotenogenic genes, albeit arranged differently, with carH often in their genomic vicinity, in all three Myxococcales suborders. However, CarA and its associated factors were found exclusively in suborder Cystobacterineae, with carA-carH invariably in tandem in a syntenic carotenogenic operon, except for Cystobacter/Melittangium, which lack CarA but retain all other factors. We experimentally show B₁₂ -mediated photoregulated carotenogenesis in representative myxobacteria, and a remarkably plastic CarH operator design and DNA binding across Myxococcales. Unlike the two characterized CarH from other phyla, which are tetrameric, Cystobacter CarH (the first myxobacterial homolog amenable to analysis in vitro) is a dimer that combines direct CarH-like B₁₂ -based photoregulation with CarA-like DNA-binding and inhibition by an antirepressor. This study provides new molecular insights into B₁₂ -dependent photoreceptors. It further establishes the B₁₂ -dependent pathway for photoregulated carotenogenesis as broadly prevalent across myxobacteria and its evolution, exclusively in one suborder, into a parallel complex B₁₂ -independent circuit. This article is protected by copyright. All rights reserved.

Inhibitors of bacterial RNA polymerase transcription complex

A series of hybrid compounds that incorporated anthranilic acid with activated 1H-indoles through a glyoxylamide linker were designed to target bacterial RNA polymerase holoenzyme formation using computational docking. Synthesis, in vitro transcription inhibition assays, and biological testing of the hybrids identified a range of potent anti-transcription inhibitors with activity against a range of pathogenic bacteria with MICs as low as 3.1 μM. A structure activity relationship study identified the key structural components necessary for inhibition of both bacterial growth and transcription. Correlation of in vitro transcription inhibition activity with in vivo mechanism of action was established using fluorescence microscopy and resistance passaging using Gram-positive bacteria showed no resistance development over 30 days. Furthermore, no toxicity was observed from the compounds in a wax moth larvae model, establishing a platform for the development of a series of new antibacterial drugs with an established mode of action.

Determination of the two-component systems regulatory network reveals core and accessory regulations across Pseudomonas aeruginosa lineages

Pseudomonas aeruginosa possesses one of the most complex bacterial regulatory networks, which largely contributes to its success as a pathogen. However, most of its transcription factors (TFs) are still uncharacterized and the potential intra-species variability in regulatory networks has been mostly ignored so far. Here, we used DAP-seq to map the genome-wide binding sites of all 55 DNA-binding two-component systems (TCSs) response regulators (RRs) across the three major P. aeruginosa lineages. The resulting networks encompass about 40% of all genes in each strain and contain numerous new regulatory interactions across most major physiological processes. Strikingly, about half of the detected targets are specific to only one or two strains, revealing a previously unknown large functional diversity of TFs within a single species. Three main mechanisms were found to drive this diversity, including differences in accessory genome content, as exemplified by the strain-specific plasmid in IHMA87 outlier strain which harbors numerous binding sites of conserved chromosomally-encoded RRs. Additionally, most RRs display potential auto-regulation or RR-RR cross-regulation, bringing to light the vast complexity of this network. Overall, we provide the first complete delineation of the TCSs regulatory network in P. aeruginosa that will represent an important resource for future studies on this pathogen.

The CovR regulatory network drives the evolution of Group B Streptococcus virulence

Virulence of the neonatal pathogen Group B Streptococcus is under the control of the master regulator CovR. Inactivation of CovR is associated with large-scale transcriptome remodeling and impairs almost every step of the interaction between the pathogen and the host. However, transcriptome analyses suggested a plasticity of the CovR signaling pathway in clinical isolates leading to phenotypic heterogeneity in the bacterial population. In this study, we characterized the CovR regulatory network in a strain representative of the CC-17 hypervirulent lineage responsible of the majority of neonatal meningitis. Transcriptome and genome-wide binding analysis reveal the architecture of the CovR network characterized by the direct repression of a large array of virulence-associated genes and the extent of co-regulation at specific loci. Comparative functional analysis of the signaling network links strain-specificities to the regulation of the pan-genome, including the two specific hypervirulent adhesins and horizontally acquired genes, to mutations in CovR-regulated promoters, and to variability in CovR activation by phosphorylation. This regulatory adaptation occurs at the level of genes, promoters, and of CovR itself, and allows to globally reshape the expression of virulence genes. Overall, our results reveal the direct, coordinated, and strain-specific regulation of virulence genes by the master regulator CovR and suggest that the intra-species evolution of the signaling network is as important as the expression of specific virulence factors in the emergence of clone associated with specific diseases.

Streptococcus agalactiae, commonly known as the Group B Streptococcus (GBS), is a commensal bacterium of the intestinal and vaginal tracts found in approximately 30% of healthy adults. However, GBS is also an opportunistic pathogen and the leading cause of neonatal invasive infections. Epidemiologic data have identified a particular GBS clone, designated the CC-17 hypervirulent clonal complex, as responsible for the overwhelming majority of neonatal meningitis. The hypervirulence of CC-17 has been linked to the expression of two specific surface proteins increasing their abilities to cross epithelial and endothelial barriers. In this study, we characterized the role of the major regulator of virulence gene expression, the CovR response regulator, in a representative hypervirulent strain. Transcriptome and genome-wide binding analysis reveal the architecture of the CovR signaling network characterized by the direct repression of a large array of virulence-associated genes, including the specific hypervirulent adhesins. Comparative analysis in a non-CC-17 wild type strain demonstrates a high level of plasticity of the regulatory network, allowing to globally reshape pathogen-host interaction. Overall, our results suggest that the intra-species evolution of the regulatory network is an important factor in the emergence of GBS clones associated with specific pathologies.

How the PhoP/PhoQ System Controls Virulence and Mg2+ Homeostasis: Lessons in Signal Transduction, Pathogenesis, Physiology, and Evolution

The PhoP/PhoQ two-component system governs virulence, Mg2+ homeostasis, and resistance to a variety of antimicrobial agents, including acidic pH and cationic antimicrobial peptides, in several Gram-negative bacterial species. Best understood in Salmonella enterica serovar Typhimurium, the PhoP/PhoQ system consists o-regulated gene products alter PhoP-P amounts, even under constant inducing conditions. PhoP-P controls the abundance of hundreds of proteins both directly, by having transcriptional effects on the corresponding genes, and indirectly, by modifying the abundance, activity, or stability of other transcription factors, regulatory RNAs, protease regulators, and metabolites. The investigation of PhoP/PhoQ has uncovered novel forms of signal transduction and the physiological consequences of regulon evolution.

Plastic rewiring of Sef1 transcriptional networks and the potential of non-functional transcription factor binding in facilitating adaptive evolution

Prior and extensive plastic rewiring of a transcriptional network, followed by a functional switch of the conserved transcriptional regulator, can shape the evolution of a new network with diverged functions. The presence of three distinct iron regulatory systems in fungi that use orthologous transcriptional regulators suggests that these systems evolved in that manner. Orthologs of the transcriptional activator Sef1 are believed to be central to how iron regulatory systems developed in fungi, involving gene gain, plastic network rewiring, and switches in regulatory function. We show that, in the protoploid yeast Lachancea kluyveri, plastic rewiring of the L. kluyveri Sef1 (Lk-Sef1) network, together with a functional switch, enabled Lk-Sef1 to regulate TCA cycle genes, unlike Candida albicans Sef1 that mainly regulates iron-uptake genes. Moreover, we observed pervasive nonfunctional binding of Sef1 to its target genes. Enhancing Lk-Sef1 activity resuscitated the corresponding transcriptional network, providing immediate adaptive benefits in changing environments. Our study not only sheds light on the evolution of Sef1-centered transcriptional networks but also shows the adaptive potential of nonfunctional transcription factor binding for evolving phenotypic novelty and diversity.

Genome Features of Asaia sp. W12 Isolated from the Mosquito Anopheles stephensi Reveal Symbiotic Traits

Asaia bacteria commonly comprise part of the microbiome of many mosquito species in the genera Anopheles and Aedes, including important vectors of infectious agents. Their close association with multiple organs and tissues of their mosquito hosts enhances the potential for paratransgenesis for the delivery of antimalaria or antivirus effectors. The molecular mechanisms involved in the interactions between Asaia and mosquito hosts, as well as Asaia and other bacterial members of the mosquito microbiome, remain underexplored. Here, we determined the genome sequence of Asaia strain W12 isolated from Anopheles stephensi mosquitoes, compared it to other Asaia species associated with plants or insects, and investigated the properties of the bacteria relevant to their symbiosis with mosquitoes. The assembled genome of strain W12 had a size of 3.94 MB, the largest among Asaia spp. studied so far. At least 3585 coding sequences were predicted. Insect-associated Asaia carried more glycoside hydrolase (GH)-encoding genes than those isolated from plants, showing their high plant biomass-degrading capacity in the insect gut. W12 had the most predicted regulatory protein components comparatively among the selected Asaia, indicating its capacity to adapt to frequent environmental changes in the mosquito gut. Two complete operons encoding cytochrome bo₃-type ubiquinol terminal oxidases (cyoABCD-1 and cyoABCD-2) were found in most Asaia genomes, possibly offering alternative terminal oxidases and allowing the flexible transition of respiratory pathways. Genes involved in the production of 2,3-butandiol and inositol have been found in Asaia sp. W12, possibly contributing to biofilm formation and stress tolerance.

A stress-induced block in dicarboxylate uptake and utilization in Salmonella

Bacteria have evolved to sense and respond to their environment by altering gene expression and metabolism to promote growth and survival. In this work we demonstrate that Salmonella displays an extensive (>30 hour) lag in growth when subcultured into media where dicarboxylates such as succinate are the sole carbon source. This growth lag is regulated in part by RpoS, the RssB anti-adaptor IraP, translation elongation factor P, and to a lesser degree the stringent response. We also show that small amounts of proline or citrate can trigger early growth in succinate media and that, at least for proline, this effect requires the multifunctional enzyme/regulator PutA. We demonstrate that activation of RpoS results in the repression of dctA, encoding the primary dicarboxylate importer, and that constitutive expression of dctA induced growth. This dicarboxylate growth lag phenotype is far more severe across multiple Salmonella isolates than in its close relative E. coli Replacing 200 nt of the Salmonella dctA promoter region with that of E. coli was sufficient to eliminate the observed lag in growth. We hypothesized that this cis-regulatory divergence might be an adaptation to Salmonella's virulent lifestyle where levels of phagocyte-produced succinate increase in response to bacterial LPS, however we found that impairing dctA repression had no effect on Salmonella's survival in acidified succinate or in macrophages.Importance Bacteria have evolved to sense and respond to their environment to maximize their chance of survival. By studying differences in the responses of pathogenic bacteria and closely related non-pathogens, we can gain insight into what environments they encounter inside of an infected host. Here we demonstrate that Salmonella diverges from its close relative E. coli in its response to dicarboxylates such as the metabolite succinate. We show that this is regulated by stress response proteins and ultimately can be attributed to Salmonella repressing its import of dicarboxylates. Understanding this phenomenon may reveal a novel aspect of the Salmonella virulence cycle, and our characterization of its regulation yields a number of mutant strains that can be used to further study it.

The pursuit of mechanism of action: uncovering drug complexity in TB drug discovery

Whole cell-based phenotypic screens have become the primary mode of hit generation in tuberculosis (TB) drug discovery during the last two decades. Different drug screening models have been developed to mirror the complexity of TB disease in the laboratory. As these culture conditions are becoming more and more sophisticated, unraveling the drug target and the identification of the mechanism of action (MOA) of compounds of interest have additionally become more challenging. A good understanding of MOA is essential for the successful delivery of drug candidates for TB treatment due to the high level of complexity in the interactions between Mycobacterium tuberculosis (Mtb) and the TB drug used to treat the disease. There is no single "standard" protocol to follow and no single approach that is sufficient to fully investigate how a drug restrains Mtb. However, with the recent advancements in -omics technologies, there are multiple strategies that have been developed generally in the field of drug discovery that have been adapted to comprehensively characterize the MOAs of TB drugs in the laboratory. These approaches have led to the successful development of preclinical TB drug candidates, and to a better understanding of the pathogenesis of Mtb infection. In this review, we describe a plethora of efforts based upon genetic, metabolomic, biochemical, and computational approaches to investigate TB drug MOAs. We assess these different platforms for their strengths and limitations in TB drug MOA elucidation in the context of Mtb pathogenesis. With an emphasis on the essentiality of MOA identification, we outline the unmet needs in delivering TB drug candidates and provide direction for further TB drug discovery.

Transcriptional rewiring of the GcrA/CcrM bacterial epigenetic regulatory system in closely related bacteria

Transcriptional rewiring is the regulation of different target genes by orthologous regulators in different organisms. While this phenomenon has been observed, it has not been extensively studied, particularly in core regulatory systems. Several global cell cycle regulators are conserved in the Alphaproteobacteria, providing an excellent model to study this phenomenon. First characterized in Caulobacter crescentus, GcrA and CcrM compose a DNA methylation-based regulatory system that helps coordinate the complex life cycle of this organism. These regulators are well-conserved across Alphaproteobacteria, but the extent to which their regulatory targets are conserved is not known. In this study, the regulatory targets of GcrA and CcrM were analyzed by SMRT-seq, RNA-seq, and ChIP-seq technologies in the Alphaproteobacterium Brevundimonas subvibrioides, and then compared to those of its close relative C. crescentus that inhabits the same environment. Although the regulators themselves are highly conserved, the genes they regulate are vastly different. GcrA directly regulates 204 genes in C. crescentus, and though B. subvibrioides has orthologs to 147 of those genes, only 48 genes retained GcrA binding in their promoter regions. Additionally, only 12 of those 48 genes demonstrated significant transcriptional change in a gcrA mutant, suggesting extensive transcriptional rewiring between these organisms. Similarly, out of hundreds of genes CcrM regulates in each of these organisms, only 2 genes were found in common. When multiple Alphaproteobacterial genomes were analyzed bioinformatically for potential GcrA regulatory targets, the regulation of genes involved in DNA replication and cell division was well conserved across the Caulobacterales but not outside this order. This work suggests that significant transcriptional rewiring can occur in cell cycle regulatory systems even over short evolutionary distances.

The degree to which genetic or physiological systems evolve over evolutionary distance is often untested. One can assume that the same system in different organisms will change very little if 1) the evolutionary distance between the organisms is small, 2) the systems perform critical functions, and 3) the organisms have been under similar selective pressures (i.e. the organisms inhabited the same ecological niche). The Alphaproteobacteria offer an excellent opportunity to test this assertion as several critical global transcriptional regulators are conserved throughout this clade. In this study, the regulons of two such global regulators, GcrA and CcrM, in two closely related Alphaproteobacteria that inhabit the same ecological niche were compared and it was found that they regulate vastly different genes. In many cases, genes were present in both organisms, but targeted by a regulator in one organism and not in the other. These results suggest that significant transcriptional rewiring can occur even in a core regulatory system over small evolutionary distances and indicate that conservation of genes and genetic regulators may not be a complete indicator of their physiological function in an organism.

Transcriptional regulation of organohalide pollutant utilisation in bacteria

Organohalides are organic molecules formed biotically and abiotically, both naturally and through industrial production. They are usually toxic and represent a health risk for living organisms, including humans. Bacteria capable of degrading organohalides for growth express dehalogenase genes encoding enzymes that cleave carbon-halogen bonds. Such bacteria are of potential high interest for bioremediation of contaminated sites. Dehalogenase genes are often part of gene clusters that may include regulators, accessory genes and genes for transporters and other enzymes of organohalide degradation pathways. Organohalides and their degradation products affect the activity of regulatory factors, and extensive genome-wide modulation of gene expression helps dehalogenating bacteria to cope with stresses associated with dehalogenation, such as intracellular increase of halides, dehalogenase-dependent acid production, organohalide toxicity and misrouting and bottlenecks in metabolic fluxes. This review focuses on transcriptional regulation of gene clusters for dehalogenation in bacteria, as studied in laboratory experiments and in situ. The diversity in gene content, organization and regulation of such gene clusters is highlighted for representative organohalide-degrading bacteria. Selected examples illustrate a key, overlooked role of regulatory processes, often strain-specific, for efficient dehalogenation and productive growth in presence of organohalides.

Widespread Strain-Specific Distinctions in Chromosomal Binding Dynamics of a Highly Conserved Escherichia coli Transcription Factor

Bacterial gene regulation is governed by often hundreds of transcription factors (TFs) that bind directly to targets on the chromosome. Global studies of TFs usually make assumptions that regulatory targets within model strains will be conserved between members of the same species harboring common genetic targets. We recently discovered that YhaJ of Escherichia coli is crucial for virulence in two different pathotypes but binds to distinct regions of their genomes and regulates no common genes. This surprising result leads to strain-specific mechanisms of virulence regulation, but the implications for other E. coli pathotypes or commensals were unclear. Here, we report that heterogenous binding of YhaJ is widespread within the E. coli species. We analyzed the global YhaJ binding dynamics of four evolutionarily distinct E. coli isolates under two conditions, revealing 78 significant sites on the core genome as well as horizontally acquired loci. Condition-dependent dosage of YhaJ correlated with the number of occupied sites in vivo but did not significantly alter its enrichment at regions bound in both conditions, explaining the availability of this TF to occupy accessory sites in response to the environment. Strikingly, only ∼15% of YhaJ binding sites were common to all strains. Furthermore, differences in enrichment of uncommon sites were observed largely in chromosomal regions found in all strains and not explained exclusively by binding to strain-specific horizontally acquired elements or mutations in the DNA binding sequence. This observation suggests that intraspecies distinctions in TF binding dynamics are a widespread phenomenon and represent strain-specific gene regulatory potential.IMPORTANCE In bacterial cells, hundreds of transcription factors coordinate gene regulation and thus are a major driver of cellular processes. However, the immense diversity in bacterial genome structure and content makes deciphering regulatory networks challenging. This is particularly apparent for the model organism Escherichia coli as evolution has driven the emergence of species members with highly distinct genomes, which occupy extremely different niches in nature. While it is well-known that transcription factors must integrate horizontally acquired DNA into the regulatory network of the cell, the extent of regulatory diversity beyond single model strains is unclear. We have explored this concept in four evolutionarily distinct E. coli strains and show that a highly conserved transcription factor displays unprecedented diversity in chromosomal binding sites. Importantly, this diversity is not restricted to strain-specific DNA or mutation in binding sites. This observation suggests that strain-specific regulatory networks are potentially widespread within individual bacterial species.

Sensing and responding to diverse extracellular signals: an updated analysis of the sensor kinases and response regulators of Streptomyces species

Streptomyces venezuelae is a Gram-positive, filamentous actinomycete with a complex developmental life cycle. Genomic analysis revealed that S. venezuelae encodes a large number of two-component systems (TCSs): these consist of a membrane-bound sensor kinase (SK) and a cognate response regulator (RR). These proteins act together to detect and respond to diverse extracellular signals. Some of these systems have been shown to regulate antimicrobial biosynthesis in Streptomyces species, making them very attractive to researchers. The ability of S. venezuelae to sporulate in both liquid and solid cultures has made it an increasingly popular model organism in which to study these industrially and medically important bacteria. Bioinformatic analysis identified 58 TCS operons in S. venezuelae with an additional 27 orphan SK and 18 orphan RR genes. A broader approach identified 15 of the 58 encoded TCSs to be highly conserved in 93 Streptomyces species for which high-quality and complete genome sequences are available. This review attempts to unify the current work on TCS in the streptomycetes, with an emphasis on S. venezuelae.

Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

When is a transcription factor a NAP?

Proteins that regulate transcription often also play an architectural role in the genome. Thus, it has been difficult to define with precision the distinctions between transcription factors and nucleoid-associated proteins (NAPs). Anachronistic descriptions of NAPs as 'histone-like' implied an organizational function in a bacterial chromatin-like complex. Definitions based on protein abundance, regulatory mechanisms, target gene number, or the features of their DNA-binding sites are insufficient as marks of distinction, and trying to distinguish transcription factors and NAPs based on their ranking within regulatory hierarchies or positions in gene-control networks is also unsatisfactory. The terms 'transcription factor' and 'NAP' are ad hoc operational definitions with each protein lying along a spectrum of structural and functional features extending from highly specific actors with few gene targets to those with a pervasive influence on the transcriptome. The Streptomyces BldC protein is used to illustrate these issues.

Plastic Circuits: Regulatory Flexibility in Fine Tuning Pathogen Success

Bacterial pathogens employ diverse fitness and virulence mechanisms to gain an advantage in competitive niches. These lifestyle-specific traits require integration into the regulatory network of the cell and are often controlled by pre-existing transcription factors. In this review, we highlight recent advances that have been made in characterizing this regulatory flexibility in prominent members of the Enterobacteriaceae. We focus on the direct global interactions between transcription factors and their target genes in pathogenic Escherichia coli and Salmonella revealed using chromatin immunoprecipitation coupled with next-generation sequencing. Furthermore, the implications and advantages of such regulatory adaptations in benefiting distinct pathogenic lifestyles are discussed.

Horizontally Acquired Quorum-Sensing Regulators Recruited by the PhoP Regulatory Network Expand the Host Adaptation Repertoire in the Phytopathogen Pectobacterium brasiliense

In this study, we examine the impact of transcriptional network rearrangements driven by horizontal gene acquisition in PhoP and SlyA regulons using as a case study a phytopathosystem comprised of potato tubers and the soft-rot pathogen Pectobacterium brasiliense 1692 (Pb1692). Genome simulations and statistical analyses uncovered the tendency of PhoP and SlyA networks to mobilize lineage-specific traits predicted as horizontal gene transfer at late infection, highlighting the prominence of regulatory network rearrangements in this stage of infection. The evidence further supports the circumscription of two horizontally acquired quorum-sensing regulators (carR and expR1) by the PhoP network. By recruiting carR and expR1, the PhoP network also impacts certain host adaptation- and bacterial competition-related systems, seemingly in a quorum sensing-dependent manner, such as the type VI secretion system, carbapenem biosynthesis, and plant cell wall-degrading enzymes (PCWDE) like cellulases and pectate lyases. Conversely, polygalacturonases and the type III secretion system (T3SS) exhibit a transcriptional pattern that suggests quorum-sensing-independent regulation by the PhoP network. This includes an uncharacterized novel phage-related gene family within the T3SS gene cluster that has been recently acquired by two Pectobacterium species. The evidence further suggests a PhoP-dependent regulation of carbapenem- and PCWDE-encoding genes based on the synthesized products' optimum pH. The PhoP network also controls slyA expression in planta, which seems to impact carbohydrate metabolism regulation, especially at early infection, when 76.2% of the SlyA-regulated genes from that category also require PhoP to achieve normal expression levels.IMPORTANCE Exchanging genetic material through horizontal transfer is a critical mechanism that drives bacteria to efficiently adapt to host defenses. In this report, we demonstrate that a specific plant-pathogenic species (from the Pectobacterium genus) successfully integrated a population density-based behavior system (quorum sensing) acquired through horizontal transfer into a resident stress-response gene regulatory network controlled by the PhoP protein. Evidence found here underscores that subsets of bacterial weaponry critical for colonization, typically known to respond to quorum sensing, are also controlled by PhoP. Some of these traits include different types of enzymes that can efficiently break down plant cell walls depending on the environmental acidity level. Thus, we hypothesize that PhoP's ability to elicit regulatory responses based on acidity and nutrient availability fluctuations has strongly impacted the fixation of its regulatory connection with quorum sensing. In addition, another global gene regulator, known as SlyA, was found under the PhoP regulatory network. The SlyA regulator controls a series of carbohydrate metabolism-related traits, which also seem to be regulated by PhoP. By centralizing quorum sensing and slyA under PhoP scrutiny, Pectobacterium cells added an advantageous layer of control over those two networks that potentially enhances colonization efficiency.

Distinct intraspecies virulence mechanisms regulated by a conserved transcription factor

Bacterial pathogens emerge by adapting mechanisms of virulence, differentiating them from their nonpathogenic progenitor. Virulence factors are often encoded on accessory genomic elements not part of the core genome and therefore must be integrated into the regulatory architecture of the cell. Here, we show that a highly conserved transcription factor in Escherichia coli has been relieved of a common purpose and adapted to regulate virulence pleiotropically in 2 distinct genetic backgrounds. This leads to enhanced virulence of both intestinal enterohemorrhagic E. coli and extraintestinal uropathogenic E. coli by exclusive mechanisms. These findings challenge the assumption that conserved transcription factors regulate common pathways maintained within a species and suggest that transcriptional repurposing creates new primary roles on an individual basis.

Tailoring transcriptional regulation to coordinate the expression of virulence factors in tandem with the core genome is a hallmark of bacterial pathogen evolution. Bacteria encode hundreds of transcription factors forming the base-level control of gene regulation. Moreover, highly homologous regulators are assumed to control conserved genes between members within a species that harbor the same genetic targets. We have explored this concept in 2 Escherichia coli pathotypes that employ distinct virulence mechanisms that facilitate specification of a different niche within the host. Strikingly, we found that the transcription factor YhaJ actively regulated unique gene sets between intestinal enterohemorrhagic E. coli (EHEC) and extraintestinal uropathogenic E. coli (UPEC), despite being very highly conserved. In EHEC, YhaJ directly activates expression of type 3 secretion system components and effectors. Alternatively, YhaJ enhances UPEC virulence regulation by binding directly to the phase-variable type 1 fimbria promoter, driving its expression. Additionally, YhaJ was found to override the universal GAD acid tolerance system but exclusively in EHEC, thereby indirectly enhancing type 3 secretion pleiotropically. These results have revealed that within a species, conserved regulators are actively repurposed in a “personalized” manner to benefit particular lifestyles and drive virulence via multiple distinct mechanisms.

The antitoxin MqsA homologue in Pseudomonas fluorescens 2P24 has a rewired regulatory circuit through evolution

The mqsRA operon encodes a toxin-antitoxin pair that was characterized to participate in biofilm and persister cell formation in Escherichia coli. Notably, the antitoxin MqsA possesses a C-terminal DNA-binding domain that recognizes the [5'-AACCT(N)_2-4 AGGTT-3'] motif and acts as a transcriptional regulator controlling multiple genes including the general stress response regulator RpoS. However, it is unknown how the transcriptional circuits of MqsA homologues have changed in bacteria over evolutionary time. Here, we found mqsA in Pseudomonas fluorescens (PfmqsA) is acquired through horizontal gene transfer and binds to a slightly different motif [5'-TACCCT(N)₃ AGGGTA-3'], which exists upstream of the PfmqsRA operon. Interestingly, an adjacent GntR-type transcriptional regulator, which was termed AgtR, is under negative control of PfMqsA. It was further demonstrated that PfMqsA reduces production of biofilm components through AgtR, which directly regulates the pga and fap operons involved in the synthesis of extracellular polymeric substances. Moreover, through quantitative proteomics analysis, we showed AgtR is a highly pleiotropic regulator that influences up to 252 genes related to diverse processes including chemotaxis, oxidative phosphorylation and carbon and nitrogen metabolism. Taken together, our findings suggest the rewired regulatory circuit of PfMqsA influences diverse physiological aspects of P. fluorescens 2P24 via the newly characterized AgtR.

Genome-Wide Transcriptional Responses of Mycobacterium to Antibiotics

Antibiotics can stimulate or depress gene expression in bacteria. The analysis of transcriptional responses of Mycobacterium to antimycobacterial compounds has improved our understanding of the mode of action of various drug classes and the efficacy and effect of such compounds on the global metabolism of Mycobacterium. This approach can provide new insights for known antibiotics, for example those currently used for tuberculosis treatment, as well as help to identify the mode of action and predict the targets of new compounds identified by whole-cell screening assays. In addition, changes in gene expression profiles after antimycobacterial treatment can provide information about the adaptive ability of bacteria to escape the effects of antibiotics and allow monitoring of the physiology of the bacteria during treatment. Genome-wide expression profiling also makes it possible to pinpoint genes differentially expressed between drug sensitive Mycobacterium and multidrug-resistant clinical isolates. Finally, genes involved in adaptive responses and drug tolerance could become new targets for improving the efficacy of existing antibiotics.

MarR Family Transcription Factors from Burkholderia Species: Hidden Clues to Control of Virulence-Associated Genes

Species within the genus Burkholderia exhibit remarkable phenotypic diversity. Genomic plasticity, including genome reduction and horizontal gene transfer, has been correlated with virulence traits in several species. However, the conservation of virulence genes in species otherwise considered to have limited potential for infection suggests that phenotypic diversity may not be explained solely on the basis of genetic diversity. Instead, differential organization and control of gene regulatory networks may underlie many phenotypic differences. In this review, we evaluate how regulation of gene expression by members of the multiple antibiotic resistance regulator (MarR) family of transcription factors may contribute to shaping the physiological diversity of Burkholderia species, with a focus on the clinically relevant human pathogens. All Burkholderia species encode a relatively large number of MarR proteins, a feature common to bacteria that must respond to environmental changes such as those associated with host invasion. However, evolution of gene regulatory networks has likely resulted in orthologous transcription factors controlling disparate sets of genes. Adaptation to, and survival in, diverse habitats, including a human or plant host, is key to the success of Burkholderia species as (opportunistic) pathogens, and recent reports suggest that control of virulence-associated genes by MarR proteins features prominently among the survival strategies employed by these species. We suggest that identification of MarR regulons will contribute significantly to clarification of virulence determinants and phenotypic diversity.

Driving the expression of the Salmonella enterica sv Typhimurium flagellum using flhDC from Escherichia coli results in key regulatory and cellular differences

The flagellar systems of Escherichia coli and Salmonella enterica exhibit a significant level of genetic and functional synteny. Both systems are controlled by the flagellar specific master regulator FlhD₄C₂. Since the early days of genetic analyses of flagellar systems it has been known that E. coli flhDC can complement a ∆flhDC mutant in S. enterica. The genomic revolution has identified how genetic changes to transcription factors and/or DNA binding sites can impact the phenotypic outcome across related species. We were therefore interested in asking: using modern tools to interrogate flagellar gene expression and assembly, what would the impact be of replacing the flhDC coding sequences in S. enterica for the E. coli genes at the flhDC S. entercia chromosomal locus? We show that even though all strains created are motile, flagellar gene expression is measurably lower when flhDC_EC are present. These changes can be attributed to the impact of FlhD₄C₂ DNA recognition and the protein-protein interactions required to generate a stable FlhD₄C₂ complex. Furthermore, our data suggests that in E. coli the internal flagellar FliT regulatory feedback loop has a marked difference with respect to output of the flagellar systems. We argue due diligence is required in making assumptions based on heterologous expression of regulators and that even systems showing significant synteny may not behave in exactly the same manner.

Dissecting the Repertoire of DNA-Binding Transcription Factors of the Archaeon Pyrococcus furiosus DSM 3638

In recent years, there has been a large increase in the amount of experimental evidence for diverse archaeal organisms, and these findings allow for a comprehensive analysis of archaeal genetic organization. However, studies about regulatory mechanisms in this cellular domain are still limited. In this context, we identified a repertoire of 86 DNA-binding transcription factors (TFs) in the archaeon Pyrococcus furiosus DSM 3638, that are clustered into 32 evolutionary families. In structural terms, 45% of these proteins are composed of one structural domain, 41% have two domains, and 14% have three structural domains. The most abundant DNA-binding domain corresponds to the winged helix-turn-helix domain; with few alternative DNA-binding domains. We also identified seven regulons, which represent 13.5% (279 genes) of the total genes in this archaeon. These analyses increase our knowledge about gene regulation in P. furiosus DSM 3638 and provide additional clues for comprehensive modeling of transcriptional regulatory networks in the Archaea cellular domain.

Transcriptional Regulation in Archaea: From Individual Genes to Global Regulatory Networks

Archaea are major contributors to biogeochemical cycles, possess unique metabolic capabilities, and resist extreme stress. To regulate the expression of genes encoding these unique programs, archaeal cells use gene regulatory networks (GRNs) composed of transcription factor proteins and their target genes. Recent developments in genetics, genomics, and computational methods used with archaeal model organisms have enabled the mapping and prediction of global GRN structures. Experimental tests of these predictions have revealed the dynamical function of GRNs in response to environmental variation. Here, we review recent progress made in this area, from investigating the mechanisms of transcriptional regulation of individual genes to small-scale subnetworks and genome-wide global networks. At each level, archaeal GRNs consist of a hybrid of bacterial, eukaryotic, and uniquely archaeal mechanisms. We discuss this theme from the perspective of the role of individual transcription factors in genome-wide regulation, how these proteins interact to compile GRN topological structures, and how these topologies lead to emergent, high-level GRN functions. We conclude by discussing how systems biology approaches are a fruitful avenue for addressing remaining challenges, such as discovering gene function and the evolution of GRNs.

The evolutionary impact of intragenic FliA promoters in proteobacteria

In Escherichia coli, one sigma factor recognizes the majority of promoters, and six 'alternative' sigma factors recognize specific subsets of promoters. The alternative sigma factor FliA (σ28 ) recognizes promoters upstream of many flagellar genes. We previously showed that most E. coli FliA binding sites are located inside genes. However, it was unclear whether these intragenic binding sites represent active promoters. Here, we construct and assay transcriptional promoter-lacZ fusions for all 52 putative FliA promoters previously identified by ChIP-seq. These experiments, coupled with integrative analysis of published genome-scale transcriptional datasets, strongly suggest that most intragenic FliA binding sites are active promoters that transcribe highly unstable RNAs. Additionally, we show that widespread intragenic FliA-dependent transcription may be a conserved phenomenon, but that specific promoters are not themselves conserved. We conclude that intragenic FliA-dependent promoters and the resulting RNAs are unlikely to have important regulatory functions. Nonetheless, one intragenic FliA promoter is broadly conserved and constrains evolution of the overlapping protein-coding gene. Thus, our data indicate that intragenic regulatory elements can influence bacterial protein evolution and suggest that the impact of intragenic regulatory sequences on genome evolution should be considered more broadly.

Evolutionary Plasticity of AmrZ Regulation in Pseudomonas

amrZ encodes a master regulator protein conserved across pseudomonads, which can be either a positive or negative regulator of swimming motility depending on the species examined. To better understand plasticity in the regulatory function of AmrZ, we characterized the mode of regulation for this protein for two different motility-related phenotypes in Pseudomonas stutzeri As in Pseudomonas syringae, AmrZ functions as a positive regulator of swimming motility within P. stutzeri, which suggests that the functions of this protein with regard to swimming motility have switched at least twice across pseudomonads. Shifts in mode of regulation cannot be explained by changes in AmrZ sequence alone. We further show that AmrZ acts as a positive regulator of colony spreading within this strain and that this regulation is at least partially independent of swimming motility. Closer investigation of mechanistic shifts in dual-function regulators like AmrZ could provide unique insights into how transcriptional pathways are rewired between closely related species.IMPORTANCE Microbes often display finely tuned patterns of gene regulation across different environments, with major regulatory changes controlled by a small group of "master" regulators within each cell. AmrZ is a master regulator of gene expression across pseudomonads and can be either a positive or negative regulator for a variety of pathways depending on the strain and genomic context. Here, we demonstrate that the phenotypic outcomes of regulation of swimming motility by AmrZ have switched at least twice independently in pseudomonads, so that AmrZ promotes increased swimming motility in P. stutzeri and P. syringae but represses this phenotype in Pseudomonas fluorescens and Pseudomonas aeruginosa Since examples of switches in regulatory mode are relatively rare, further investigation into the mechanisms underlying shifts in regulator function for AmrZ could provide unique insights into the evolution of bacterial regulatory proteins.

Systematic and synthetic approaches to rewire regulatory networks

Microbial gene regulatory networks are composed of cis- and trans-components that in concert act to control essential and adaptive cellular functions. Regulatory components and interactions evolve to adopt new configurations through mutations and network rewiring events, resulting in novel phenotypes that may benefit the cell. Advances in high-throughput DNA synthesis and sequencing have enabled the development of new tools and approaches to better characterize and perturb various elements of regulatory networks. Here, we highlight key recent approaches to systematically dissect the sequence space of cis-regulatory elements and trans-regulators as well as their inter-connections. These efforts yield fundamental insights into the architecture, robustness, and dynamics of gene regulation and provide models and design principles for building synthetic regulatory networks for a variety of practical applications.

The Transcriptional Regulators of the CRP Family Regulate Different Essential Bacterial Functions and Can Be Inherited Vertically and Horizontally

One of the best-studied transcriptional regulatory proteins in bacteria is the Escherichia coli catabolite repressor protein (CRP) that when complexed with 3′-5′-cyclic AMP (cAMP) changes its conformation and interacts with specific DNA-sequences. CRP DNA-binding can result in positive or negative regulation of gene expression depending on the position of its interaction with respect to RNA polymerase binding site. The aim of this work is to review the biological role and phylogenetic relations that some members of the CRP family of transcriptional regulators (also known as cAMP receptor protein family) have in different bacterial species. This work is not intended to give an exhaustive revision of bacterial CRP-orthologs, but to provide examples of the role that these proteins play in the expression of genes that are fundamental for the life style of some bacterial species. We highlight the conservation of their structural characteristics and of their binding to conserved-DNA sequences, in contrast to their very diverse repertoire of gene activation. CRP activates a wide variety of fundamental genes for the biological characteristic of each bacterial species, which in several instances form part of their core-genome (defined as the gene sequences present in all members of a bacterial species). We present evidence that support the fact that some of the transcriptional regulators that belong to the CRP family in different bacterial species, and some of the genes that are regulated by them, can be inherited by horizontal gene transfer. These data are discussed in the framework of bacterial evolution models.

An Essential Regulatory System Originating from Polygenic Transcriptional Rewiring of PhoP-PhoQ of Xanthomonas campestris

How essential, regulatory genes originate and evolve is intriguing because mutations of these genes not only lead to lethality in organisms, but also have pleiotropic effects since they control the expression of multiple downstream genes. Therefore, the evolution of essential, regulatory genes is not only determined by genetic variations of their own sequences, but also by the biological function of downstream genes and molecular mechanisms of regulation. To understand the origin of essential, regulatory genes, experimental dissection of the complete regulatory cascade is needed. Here, we provide genetic evidences to reveal that PhoP-PhoQ is an essential two-component signal transduction system in the gram-negative bacterium Xanthomonas campestris, but that its orthologs in other bacteria belonging to Proteobacteria are nonessential. Mutational, biochemical, and chromatin immunoprecipitation together with high-throughput sequencing analyses revealed that phoP and phoQ of X. campestris and its close relative Pseudomonas aeruginosa are replaceable, and that the consensus binding motifs of the transcription factor PhoP are also highly conserved. PhoP _Xcc in X. campestris regulates the transcription of a number of essential, structural genes by directly binding to cis-regulatory elements (CREs); however, these CREs are lacking in the orthologous essential, structural genes in P. aeruginosa, and thus the regulatory relationships between PhoP _Pae and these downstream essential genes are disassociated. Our findings suggested that the recruitment of regulatory proteins by critical structural genes via transcription factor-CRE rewiring is a driving force in the origin and functional divergence of essential, regulatory genes.

Transcriptomic buffering of cryptic genetic variation contributes to meningococcal virulence

Background

Commensal bacteria like Neisseria meningitidis sometimes cause serious disease. However, genomic comparison of hyperinvasive and apathogenic lineages did not reveal unambiguous hints towards indispensable virulence factors. Here, in a systems biological approach we compared gene expression of the invasive strain MC58 and the carriage strain α522 under different ex vivo conditions mimicking commensal and virulence compartments to assess the strain-specific impact of gene regulation on meningococcal virulence.

Results

Despite indistinguishable ex vivo phenotypes, both strains differed in the expression of over 500 genes under infection mimicking conditions. These differences comprised in particular metabolic and information processing genes as well as genes known to be involved in host-damage such as the nitrite reductase and numerous LOS biosynthesis genes. A model based analysis of the transcriptomic differences in human blood suggested ensuing metabolic flux differences in energy, glutamine and cysteine metabolic pathways along with differences in the activation of the stringent response in both strains. In support of the computational findings, experimental analyses revealed differences in cysteine and glutamine auxotrophy in both strains as well as a strain and condition dependent essentiality of the (p)ppGpp synthetase gene relA and of a short non-coding AT-rich repeat element in its promoter region.

Conclusions

Our data suggest that meningococcal virulence is linked to transcriptional buffering of cryptic genetic variation in metabolic genes including global stress responses. They further highlight the role of regulatory elements for bacterial virulence and the limitations of model strain approaches when studying such genetically diverse species as N. meningitidis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3616-7) contains supplementary material, which is available to authorized users.

The dynamics of mobile promoters: Enhanced stability in promoter regions

Mobile promoters are emerging as a new class of mobile genetic elements, first identified by examining prokaryote genome sequences, and more recently confirmed by experimental observations in bacteria. Recent datasets have identified over 40,000 putative mobile promoters in sequenced prokaryote genomes, however only one-third of these are in regions of the genome directly upstream from coding sequences, that is, in promoter regions. The presence of many promoter sequences in non-promoter regions is unexplained. Here we develop a general mathematical model for the dynamics of mobile promoters, extending previous work to capture the dynamics both within and outside promoter regions. From this general model, we apply rigorous model selection techniques to identify which parameters are statistically justified in describing the available mobile promoter data, and find best-fit values of these parameters. Our results suggest that high rates of horizontal gene transfer maintain the population of mobile promoters in promoter regions, and that once established at these sites, mobile promoters are rarely lost, but are commonly copied to other genomic regions. In contrast, mobile promoter copies in non-promoter regions are more numerous and more volatile, experiencing substantially higher rates of duplication, loss and diversification.

The Tip of the Iceberg: On the Roles of Regulatory Small RNAs in the Virulence of Enterohemorrhagic and Enteropathogenic Escherichia coli

Enterohemorrhagic and enteropathogenic Escherichia coli are gastrointestinal pathogens that disrupt the intestinal microvilli to form attaching and effacing (A/E) lesions on infected cells and cause diarrhea. This pathomorphological trait is encoded within the pathogenicity island locus of enterocyte effacement (LEE). The LEE houses a type 3 secretion system (T3SS), which upon assembly bridges the bacterial cytosol to that of the host and enables the bacterium to traffic dozens of effectors into the host where they hijack regulatory and signal transduction pathways and contribute to bacterial colonization and disease. Owing to the importance of the LEE to EHEC and EPEC pathogenesis, much of the research on these pathogens has centered on its regulation. To date, over 40 proteinaceous factors have been identified that control the LEE at various hierarchical levels of gene expression. In contrast, RNA-based regulatory mechanisms that converge on the LEE have only just begun to be unraveled. In this minireview, we highlight major breakthroughs in small RNAs (sRNAs)-dependent regulation of the LEE, with an emphasis on their mechanisms of action and/or LEE-encoded targets.

The Impact of Gene Silencing on Horizontal Gene Transfer and Bacterial Evolution

The H-NS family of DNA-binding proteins is the subject of intense study due to its important roles in the regulation of horizontally acquired genes critical for virulence, antibiotic resistance, and metabolism. Xenogeneic silencing proteins, typified by the H-NS protein of Escherichia coli, specifically target and downregulate expression from AT-rich genes by selectively recognizing specific structural features unique to the AT-rich minor groove. In doing so, these proteins facilitate bacterial evolution; enabling these cells to engage in horizontal gene transfer while buffering potential any detrimental fitness consequences that may result from it. Xenogeneic silencing and counter-silencing explain how bacterial cells can evolve effective gene regulatory strategies in the face of rampant gene gain and loss and it has extended our understanding of bacterial gene regulation beyond the classic operon model. Here we review the structures and mechanisms of xenogeneic silencers as well as their impact on bacterial evolution. Several H-NS-like proteins appear to play a role in facilitating gene transfer by other mechanisms including by regulating transposition, conjugation, and participating in the activation of virulence loci like the locus of enterocyte effacement pathogenicity island of pathogenic strains of E. coli. Evidence suggests that the critical determinants that dictate whether an H-NS-like protein will be a silencer or will perform a different function do not lie in the DNA-binding domain but, rather, in the domains that control oligomerization. This suggests that H-NS-like proteins are transcription factors that both recognize and alter the shape of DNA to exert specific effects that include but are not limited to gene silencing.

Multi-host environments select for host-generalist conjugative plasmids

Background

Conjugative plasmids play an important role in bacterial evolution by transferring ecologically important genes within and between species. A key limit on interspecific horizontal gene transfer is plasmid host range. Here, we experimentally test the effect of single and multi-host environments on the host-range evolution of a large conjugative mercury resistance plasmid, pQBR57. Specifically, pQBR57 was conjugated between strains of a single host species, either P. fluorescens or P. putida, or alternating between P. fluorescens and P. putida. Crucially, the bacterial hosts were not permitted to evolve allowing us to observe plasmid evolutionary responses in isolation.

Results

In all treatments plasmids evolved higher conjugation rates over time. Plasmids evolved in single-host environments adapted to their host bacterial species becoming less costly, but in the case of P. fluorescens-adapted plasmids, became costlier in P. putida, suggesting an evolutionary trade-off. When evolved in the multi-host environment plasmids adapted to P. fluorescens without a higher cost in P. putida.

Conclusion

Whereas evolution in a single-host environment selected for host-specialist plasmids due to a fitness trade-off, this trade-off could be circumvented in the multi-host environment, leading to the evolution of host-generalist plasmids.