Redefining fundamental concepts of transcription initiation in bacteria

Transcription activation in Escherichia coli and Salmonella

Promoter-specific activation of transcript initiation provides an important regulatory device in Escherichia coli and Salmonella. Here, we describe the different mechanisms that operate, focusing on how they have evolved to manage the "housekeeping" bacterial transcription machinery. Some mechanisms involve assisting the bacterial DNA-dependent RNA polymerase or replacing or remodeling one of its subunits. Others are directed to chromosomal DNA, improving promoter function, or relieving repression. We discuss how different activators work together at promoters and how the present complex network of transcription factors evolved.

Structure of the SigF1-dependent pilA1 gene promoter and characterization of the light-activated response in the cyanobacterium Synechococcus elongatus PCC 7942

In cyanobacteria that perform oxygenic photosynthesis, alternative sigma factors can play critical roles in environmental acclimation at the transcriptional initiation step. Here, we found in Synechococcus elongatus PCC 7942 that transcription of the pilA1 gene, encoding the type IV pilin, is dependent on one of the group 3 sigma factors, SigF1. We analyzed the promoter sequence determinants and proposed herein that the -10 and -35 boxes upstream of the transcriptional start site are critical for transcription. Interestingly, while the pilA1 promoter is activated by illumination, RNA polymerase containing SigF1 is already located on the promoter region under dark conditions, prior to illumination. This strongly suggests that promoter activation by light follows the recruitment of RNA polymerase during transcriptional initiation.

Oxidative stress response in avian pathogenic Escherichia coli

Avian pathogenic Escherichia coli (APEC) leads to significant economic losses in the poultry industry worldwide and restricts the development of the poultry industry. Oxidative stress, through the production of reactive oxygen species (ROS), damage iron‑sulfur (FeS) clusters, cysteine and methionine protein residues, and DNA, and then result in bacterial cells death. APEC has evolved a series of regulation systems to sense and quickly and appropriately respond to oxidative stress. Quorum sensing (QS), second messenger (SM), transcription factors (TFs), small regulatory RNAs (sRNAs), and two-component system (TCS) are important regulation systems ubiquitous in bacteria. It is of great significance to control APEC infection through investigating the molecular regulation mechanism on APEC adapting to oxidative stress. However, how the cross-talk among these regulation systems co-regulates transcription of oxidative stress-response genes in APEC has not been reported. This review suggests exploring connector proteins that co-regulate these regulation systems that co-activate transcription of oxidative stress-response genes to disrupt bacterial antioxidative defense mechanism in APEC, and then using these connector proteins as drug targets to control APEC infection. This review might contribute to illustrating the functional mechanism of APEC adapting to oxidative stress and exploring potential drug targets for the prevention and treatment of APEC infection.

Deciphering transcript architectural complexity in bacteria and archaea

RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6-1.7 kbp, while the median size of the 5'- and 3'-untranslated regions (UTRs) were 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community.IMPORTANCEOur understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

Assessing the transcriptional landscape of Pseudomonas phage 201ϕ2-1: Uncovering the small regulatory details of a giant phage

The transcriptional architecture of phages can deepen our understanding of the phage-host infection process and can be of key importance for phage engineering and biotechnological applications. Here, we applied ONT-cappable-sequencing, a long-read RNA-sequencing technique, to study the regulatory mechanisms of Pseudomonas infecting giant phage 201ϕ2-1. We identified 67 promoters and 132 terminators that together represent 92 transcriptional units. A full comparison of these data to the transcriptome of model Pseudomonas phage ϕKZ confirmed that the transcriptional programs of these prototypes of the Serwervirus and Phikzvirus genera are largely conserved, despite some subtle regulatory differences. Evidence supporting these shared mechanisms include the identification of highly similar sequence motifs for regulatory elements in both phages and the conservation of regulatory elements loci relative to homologous genes in each phage. Moreover, we discovered a sRNA in 201ϕ2-1 that is highly conserved among prototype members of different giant phage genera. Sequencing of the 201ϕ2-1 host genome resulted in its reclassification as Pseudomonas atacamensis, a close relative of the important agricultural biocontrol agent Pseudomonas chlororaphis. Finally, we conducted in vivo assays of eight 201ϕ2-1 terminators and found them to strongly terminate transcription in P. chlororaphis. Control elements from phage transcriptional programs have a rich history for applications in biotechnology. In these studies, we demonstrate new insight into the transcriptional program of 201ϕ2-1 and demonstrate the potential of its regulatory elements for novel and useful tools for synthetic biology circuitry.

Synonymous codon substitutions modulate transcription and translation of a divergent upstream gene by modulating antisense RNA production

Synonymous codons were originally viewed as interchangeable, with no phenotypic consequences. However, substantial evidence has now demonstrated that synonymous substitutions can perturb a variety of gene expression and protein homeostasis mechanisms, including translational efficiency, translational fidelity, and cotranslational folding of the encoded protein. To date, most studies of synonymous codon-derived perturbations have focused on effects within a single gene. Here, we show that synonymous codon substitutions made far within the coding sequence of Escherichia coli plasmid-encoded chloramphenicol acetyltransferase (cat) can significantly increase expression of the divergent upstream tetracycline resistance gene, tetR. In four out of nine synonymously recoded cat sequences tested, expression of the upstream tetR gene was significantly elevated due to transcription of a long antisense RNA (asRNA) originating from a transcription start site within cat. Surprisingly, transcription of this asRNA readily bypassed the native tet transcriptional repression mechanism. Even more surprisingly, accumulation of the TetR protein correlated with the level of asRNA, rather than total tetR RNA. These effects of synonymous codon substitutions on transcription and translation of a neighboring gene suggest that synonymous codon usage in bacteria may be under selection to both preserve the amino acid sequence of the encoded gene and avoid DNA sequence elements that can significantly perturb expression of neighboring genes. Avoiding such sequences may be especially important in plasmids and prokaryotic genomes, where genes and regulatory elements are often densely packed. Similar considerations may apply to the design of genetic circuits for synthetic biology applications.

Generating information-dense promoter sequences with optimal string packing

Dense arrangements of binding sites within nucleotide sequences can collectively influence downstream transcription rates or initiate biomolecular interactions. For example, natural promoter regions can harbor many overlapping transcription factor binding sites that influence the rate of transcription initiation. Despite the prevalence of overlapping binding sites in nature, rapid design of nucleotide sequences with many overlapping sites remains a challenge. Here, we show that this is an NP-hard problem, coined here as the nucleotide String Packing Problem (SPP). We then introduce a computational technique that efficiently assembles sets of DNA-protein binding sites into dense, contiguous stretches of double-stranded DNA. For the efficient design of nucleotide sequences spanning hundreds of base pairs, we reduce the SPP to an Orienteering Problem with integer distances, and then leverage modern integer linear programming solvers. Our method optimally packs sets of 20-100 binding sites into dense nucleotide arrays of 50-300 base pairs in 0.05-10 seconds. Unlike approximation algorithms or meta-heuristics, our approach finds provably optimal solutions. We demonstrate how our method can generate large sets of diverse sequences suitable for library generation, where the frequency of binding site usage across the returned sequences can be controlled by modulating the objective function. As an example, we then show how adding additional constraints, like the inclusion of sequence elements with fixed positions, allows for the design of bacterial promoters. The nucleotide string packing approach we present can accelerate the design of sequences with complex DNA-protein interactions. When used in combination with synthesis and high-throughput screening, this design strategy could help interrogate how complex binding site arrangements impact either gene expression or biomolecular mechanisms in varied cellular contexts.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Accelerating Genetic Sensor Development, Scale-up, and Deployment Using Synthetic Biology

Living cells are exquisitely tuned to sense and respond to changes in their environment. Repurposing these systems to create engineered biosensors has seen growing interest in the field of synthetic biology and provides a foundation for many innovative applications spanning environmental monitoring to improved biobased production. In this review, we present a detailed overview of currently available biosensors and the methods that have supported their development, scale-up, and deployment. We focus on genetic sensors in living cells whose outputs affect gene expression. We find that emerging high-throughput experimental assays and evolutionary approaches combined with advanced bioinformatics and machine learning are establishing pipelines to produce genetic sensors for virtually any small molecule, protein, or nucleic acid. However, more complex sensing tasks based on classifying compositions of many stimuli and the reliable deployment of these systems into real-world settings remain challenges. We suggest that recent advances in our ability to precisely modify nonmodel organisms and the integration of proven control engineering principles (e.g., feedback) into the broader design of genetic sensing systems will be necessary to overcome these hurdles and realize the immense potential of the field.

Molecular insights into the fine-tuning of pH-dependent ArsR-mediated regulation of the SabA adhesin in Helicobacter pylori

Adaptation to variations in pH is crucial for the ability of Helicobacter pylori to persist in the human stomach. The acid responsive two-component system ArsRS, constitutes the global regulon that responds to acidic conditions, but molecular details of how transcription is affected by the ArsR response regulator remains poorly understood. Using a combination of DNA-binding studies, in vitro transcription assays, and H. pylori mutants, we demonstrate that phosphorylated ArsR (ArsR-P) forms an active protein complex that binds DNA with high specificity in order to affect transcription. Our data showed that DNA topology is key for DNA binding. We found that AT-rich DNA sequences direct ArsR-P to specific sites and that DNA-bending proteins are important for the effect of ArsR-P on transcription regulation. The repression of sabA transcription is mediated by ArsR-P with the support of Hup and is affected by simple sequence repeats located upstream of the sabA promoter. Here stochastic events clearly contribute to the fine-tuning of pH-dependent gene regulation. Our results reveal important molecular aspects for how ArsR-P acts to repress transcription in response to acidic conditions. Such transcriptional control likely mediates shifts in bacterial positioning in the gastric mucus layer.

Escherichia coli non-coding regulatory regions are highly conserved

Microbial genome sequences are rapidly accumulating, enabling large-scale studies of sequence variation. Existing studies primarily focus on coding regions to study amino acid substitution patterns in proteins. However, non-coding regulatory regions also play a distinct role in determining physiologic responses. To investigate intergenic sequence variation on a large-scale, we identified non-coding regulatory region alleles across 2350 Escherichia coli strains. This 'alleleome' consists of 117 781 unique alleles for 1169 reference regulatory regions (transcribing 1975 genes) at single base-pair resolution. We find that 64% of nucleotide positions are invariant, and variant positions vary in a median of just 0.6% of strains. Additionally, non-coding alleles are sufficient to recover E. coli phylogroups. We find that core promoter elements and transcription factor binding sites are significantly conserved, especially those located upstream of essential or highly-expressed genes. However, variability in conservation of transcription factor binding sites is significant both within and across regulons. Finally, we contrast mutations acquired during adaptive laboratory evolution with wild-type variation, finding that the former preferentially alter positions that the latter conserves. Overall, this analysis elucidates the wealth of information found in E. coli non-coding sequence variation and expands pangenomic studies to non-coding regulatory regions at single-nucleotide resolution.

Predictive Biophysical Neural Network Modeling of a Compendium of in vivo Transcription Factor DNA Binding Profiles for Escherichia coli

The DNA binding of most Escherichia coli Transcription Factors (TFs) has not been comprehensively mapped, and few have models that can quantitatively predict binding affinity. We report the global mapping of in vivo DNA binding for 139 E. coli TFs using ChIP-Seq. We used these data to train BoltzNet, a novel neural network that predicts TF binding energy from DNA sequence. BoltzNet mirrors a quantitative biophysical model and provides directly interpretable predictions genome-wide at nucleotide resolution. We used BoltzNet to quantitatively design novel binding sites, which we validated with biophysical experiments on purified protein. We have generated models for 125 TFs that provide insight into global features of TF binding, including clustering of sites, the role of accessory bases, the relevance of weak sites, and the background affinity of the genome. Our paper provides new paradigms for studying TF-DNA binding and for the development of biophysically motivated neural networks.

Inhibition of bacterial RNA polymerase function and protein-protein interactions: a promising approach for next-generation antibacterial therapeutics

The increasing prevalence of multidrug-resistant pathogens necessitates the urgent development of new antimicrobial agents with innovative modes of action for the next generation of antimicrobial therapy. Bacterial transcription has been identified and widely studied as a viable target for antimicrobial development. The main focus of these studies has been the discovery of inhibitors that bind directly to the core enzyme of RNA polymerase (RNAP). Over the past two decades, substantial advancements have been made in understanding the properties of protein-protein interactions (PPIs) and gaining structural insights into bacterial RNAP and its associated factors. This has led to the crucial role of computational methods in aiding the identification of new PPI inhibitors to affect the RNAP function. In this context, bacterial transcriptional PPIs present promising, albeit challenging, targets for the creation of new antimicrobials. This review will succinctly outline the structural foundation of bacterial transcription networks and provide a summary of the known small molecules that target transcription PPIs.

5'UTR G-quadruplex structure enhances translation in size dependent manner

Translation initiation in bacteria is frequently regulated by various structures in the 5' untranslated region (5'UTR). Previously, we demonstrated that G-quadruplex (G4) formation in non-template DNA enhances transcription. In this study, we aim to explore how G4 formation in mRNA (RG4) at 5'UTR impacts translation using a T7-based in vitro translation system and in E. coli. We show that RG4 strongly promotes translation efficiency in a size-dependent manner. Additionally, inserting a hairpin upstream of the RG4 further enhances translation efficiency, reaching up to a 12-fold increase. We find that the RG4-dependent effect is not due to increased ribosome affinity, ribosome binding site accessibility, or mRNA stability. We propose a physical barrier model in which bulky structures in 5'UTR biases ribosome movement toward the downstream start codon, thereby increasing the translation output. This study provides biophysical insights into the regulatory role of 5'UTR structures in in vitro and bacterial translation, highlighting their potential applications in tuning gene expression.

Functional versatility of Zur in metal homeostasis, motility, biofilm formation, and stress resistance in Yersinia pseudotuberculosis

Zur (zinc uptake regulator) is a significant member of the Fur (ferric uptake regulator) superfamily, which is widely distributed in bacteria. Zur plays crucial roles in zinc homeostasis and influences cell development and environmental adaptation in various species. Yersinia pseudotuberculosis is a Gram-negative enteric that pathogen usually serves as a model organism in pathogenicity studies. The regulatory effects of Zur on the zinc transporter ZnuABC and the protein secretion system T6SS have been documented in Y. pseudotuberculosis. In this study, a comparative transcriptomics analysis between a ∆zur mutant and the wild-type (WT) strain of Y. pseudotuberculosis was conducted using RNA-seq. This analysis revealed global regulation by Zur across multiple functional categories, including membrane transport, cell motility, and molecular and energy metabolism. Additionally, Zur mediates the homeostasis not only of zinc but also ferric and magnesium in vivo. There was a notable decrease in 35 flagellar biosynthesis and assembly-related genes, leading to reduced swimming motility in the ∆zur mutant strain. Furthermore, Zur upregulated multiple simple sugar and oligopeptide transport system genes by directly binding to their promoters. The absence of Zur inhibited biofilm formation as well as reduced resistance to chloramphenicol and acidic stress. This study illustrates the comprehensive regulatory functions of Zur, emphasizing its importance in stress resistance and pathogenicity in Y. pseudotuberculosis.

IMPORTANCE

Bacteria encounter diverse stresses in the environment and possess essential regulators to modulate the expression of genes in responding to the stresses for better fitness and survival. Zur (zinc uptake regulator) plays a vital role in zinc homeostasis. Studies of Zur from multiple species reviewed that it influences cell development, stress resistance, and virulence of bacteria. Y. pseudotuberculosis is an enteric pathogen that serves a model organism in the study of pathogenicity, virulence factors, and mechanism of environmental adaptation. In this study, transcriptomics analysis of Zur's regulons was conducted in Y. pseudotuberculosis. The functions of Zur as a global regulator in metal homeostasis, motility, nutrient acquisition, glycan metabolism, and nucleotide metabolism, in turn, increasing the biofilm formation, stress resistance, and virulence were reviewed. The importance of Zur in environmental adaptation and pathogenicity of Y. pseudotuberculosis was emphasized.

Flexible gold standards for transcription factor regulatory interactions in Escherichia coli K-12: architecture of evidence types

Post-genomic implementations have expanded the experimental strategies to identify elements involved in the regulation of transcription initiation. Here, we present for the first time a detailed analysis of the sources of knowledge supporting the collection of transcriptional regulatory interactions (RIs) of Escherichia coli K-12. An RI groups the transcription factor, its effect (positive or negative) and the regulated target, a promoter, a gene or transcription unit. We improved the evidence codes so that specific methods are incorporated and classified into independent groups. On this basis we updated the computation of confidence levels, weak, strong, or confirmed, for the collection of RIs. These updates enabled us to map the RI set to the current collection of HT TF-binding datasets from ChIP-seq, ChIP-exo, gSELEX and DAP-seq in RegulonDB, enriching in this way the evidence of close to one-quarter (1329) of RIs from the current total 5446 RIs. Based on the new computational capabilities of our improved annotation of evidence sources, we can now analyze the internal architecture of evidence, their categories (experimental, classical, HT, computational), and confidence levels. This is how we know that the joint contribution of HT and computational methods increase the overall fraction of reliable RIs (the sum of confirmed and strong evidence) from 49% to 71%. Thus, the current collection has 3912 reliable RIs, with 2718 or 70% of them with classical evidence which can be used to benchmark novel HT methods. Users can selectively exclude the method they want to benchmark, or keep for instance only the confirmed interactions. The recovery of regulatory sites in RegulonDB by the different HT methods ranges between 33% by ChIP-exo to 76% by ChIP-seq although as discussed, many potential confounding factors limit their interpretation. The collection of improvements reported here provides a solid foundation to incorporate new methods and data, and to further integrate the diverse sources of knowledge of the different components of the transcriptional regulatory network. There is no other genomic database that offers this comprehensive high-quality architecture of knowledge supporting a corpus of transcriptional regulatory interactions.

Minimizing endogenous cryptic plasmids to construct antibiotic-free expression systems for Escherichia coli Nissle 1917

The probiotic bacterium Escherichia coli Nissle 1917 (EcN) holds significant promise for use in clinical and biological industries. However, the reliance on antibiotics to maintain plasmid-borne genes has overshadowed its benefits. In this study, we addressed this issue by engineering the endogenous cryptic plasmids pMUT1 and pMUT2. The non-essential elements were removed to create more stable derivatives pMUT1NR△ and pMUT2HBC△. Synthetic promoters by integrating binding motifs on sigma factors were further constructed and applied for expression of Bacteroides thetaiotaomicron heparinase III and the biosynthesis of ectoine. Compared to traditional antibiotic-dependent expression systems, our newly constructed antibiotic-free expression systems offer considerable advantages for clinical and synthetic biology applications.

Generating information-dense promoter sequences with optimal string packing

Dense arrangements of binding sites within nucleotide sequences can collectively influence downstream transcription rates or initiate biomolecular interactions. For example, natural promoter regions can harbor many overlapping transcription factor binding sites that influence the rate of transcription initiation. Despite the prevalence of overlapping binding sites in nature, rapid design of nucleotide sequences with many overlapping sites remains a challenge. Here, we show that this is an NP-hard problem, coined here as the nucleotide String Packing Problem (SPP). We then introduce a computational technique that efficiently assembles sets of DNA-protein binding sites into dense, contiguous stretches of double-stranded DNA. For the efficient design of nucleotide sequences spanning hundreds of base pairs, we reduce the SPP to an Orienteering Problem with integer distances, and then leverage modern integer linear programming solvers. Our method optimally packs libraries of 20-100 binding sites into dense nucleotide arrays of 50-300 base pairs in 0.05-10 seconds. Unlike approximation algorithms or meta-heuristics, our approach finds provably optimal solutions. We demonstrate how our method can generate large sets of diverse sequences suitable for library generation, where the frequency of binding site usage across the returned sequences can be controlled by modulating the objective function. As an example, we then show how adding additional constraints, like the inclusion of sequence elements with fixed positions, allows for the design of bacterial promoters. The nucleotide string packing approach we present can accelerate the design of sequences with complex DNA-protein interactions. When used in combination with synthesis and high-throughput screening, this design strategy could help interrogate how complex binding site arrangements impact either gene expression or biomolecular mechanisms in varied cellular contexts.

Author Summary

The way protein binding sites are arranged on DNA can control the regulation and transcription of downstream genes. Areas with a high concentration of binding sites can enable complex interplay between transcription factors, a feature that is exploited by natural promoters. However, designing synthetic promoters that contain dense arrangements of binding sites is a challenge. The task involves overlapping many binding sites, each typically about 10 nucleotides long, within a constrained sequence area, which becomes increasingly difficult as sequence length decreases, and binding site variety increases. We introduce an approach to design nucleotide sequences with optimally packed protein binding sites, which we call the nucleotide String Packing Problem (SPP). We show that the SPP can be solved efficiently using integer linear programming to identify the densest arrangements of binding sites for a specified sequence length. We show how adding additional constraints, like the inclusion of sequence elements with fixed positions, allows for the design of bacterial promoters. The presented approach enables the rapid design and study of nucleotide sequences with complex, dense binding site architectures.

Genomic basis of environmental adaptation in the widespread poly-extremophilic Exiguobacterium group

Delineating cohesive ecological units and determining the genetic basis for their environmental adaptation are among the most important objectives in microbiology. In the last decade, many studies have been devoted to characterizing the genetic diversity in microbial populations to address these issues. However, the impact of extreme environmental conditions, such as temperature and salinity, on microbial ecology and evolution remains unclear so far. In order to better understand the mechanisms of adaptation, we studied the (pan)genome of Exiguobacterium, a poly-extremophile bacterium able to grow in a wide range of environments, from permafrost to hot springs. To have the genome for all known Exiguobacterium type strains, we first sequenced those that were not yet available. Using a reverse-ecology approach, we showed how the integration of phylogenomic information, genomic features, gene and pathway enrichment data, regulatory element analyses, protein amino acid composition, and protein structure analyses of the entire Exiguobacterium pangenome allows to sharply delineate ecological units consisting of mesophilic, psychrophilic, halophilic-mesophilic, and halophilic-thermophilic ecotypes. This in-depth study clarified the genetic basis of the defined ecotypes and identified some key mechanisms driving the environmental adaptation to extreme environments. Our study points the way to organizing the vast microbial diversity into meaningful ecologically units, which, in turn, provides insight into how microbial communities adapt and respond to different environmental conditions in a changing world.

RegulonDB v12.0: a comprehensive resource of transcriptional regulation in E. coli K-12

RegulonDB is a database that contains the most comprehensive corpus of knowledge of the regulation of transcription initiation of Escherichia coli K-12, including data from both classical molecular biology and high-throughput methodologies. Here, we describe biological advances since our last NAR paper of 2019. We explain the changes to satisfy FAIR requirements. We also present a full reconstruction of the RegulonDB computational infrastructure, which has significantly improved data storage, retrieval and accessibility and thus supports a more intuitive and user-friendly experience. The integration of graphical tools provides clear visual representations of genetic regulation data, facilitating data interpretation and knowledge integration. RegulonDB version 12.0 can be accessed at https://regulondb.ccg.unam.mx.

Experimental promoter identification of a foodborne pathogen Salmonella enterica subsp. enterica serovar Typhimurium with near single base-pair resolution

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a common foodborne pathogen which is frequently used as the reference strain for Salmonella. Investigating the sigma factor network and protomers is crucial to understand the genomic and transcriptomic properties of the bacterium. Its promoters were identified using various methods such as dRNA-seq, ChIP-chip, or ChIP-Seq. However, validation using ChIP-exo, which exhibits higher-resolution performance compared to conventional ChIP, has not been conducted to date. In this study, using the representative strain S. Typhimurium LT2 (LT2), the ChIP-exo experiment was conducted to accurately determine the binding sites of catalytic RNA polymerase subunit RpoB and major sigma factors (RpoD, RpoN, RpoS, and RpoE) during exponential phase. Integrated with the results of RNA-Seq, promoters and sigmulons for the sigma factors and their association with RpoB have been discovered. Notably, the overlapping regions among binding sites of each alternative sigma factor were found. Furthermore, comparative analysis with Escherichia coli str. K-12 substr. MG1655 (MG1655) revealed conserved binding sites of RpoD and RpoN across different species. In the case of small RNAs (sRNAs), 50 sRNAs observed their expression during the exponential growth of LT2. Collectively, the integration of ChIP-exo and RNA-Seq enables genome-scale promoter mapping with high resolution and facilitates the characterization of binding events of alternative sigma factors, enabling a comprehensive understanding of the bacterial sigma factor network and condition-specific active promoters.

Targeting bacterial transcription factors for infection control: opportunities and challenges

The rising threat of antibiotic resistance in pathogenic bacteria emphasizes the need for new therapeutic strategies. This review focuses on bacterial transcription factors (TFs), which play crucial roles in bacterial pathogenesis. We discuss the regulatory roles of these factors through examples, and we outline potential therapeutic strategies targeting bacterial TFs. Specifically, we discuss the use of small molecules to interfere with TF function and the development of transcription factor decoys, oligonucleotides that compete with promoters for TF binding. We also cover peptides that target the interaction between the bacterial TF and other factors, such as RNA polymerase, and the targeting of sigma factors. These strategies, while promising, come with challenges, from identifying targets to designing interventions, managing side effects, and accounting for changing bacterial resistance patterns. We also delve into how Artificial Intelligence contributes to these efforts and how it may be exploited in the future, and we touch on the roles of multidisciplinary collaboration and policy to advance this research domain.Abbreviations: AI, artificial intelligence; CNN, convolutional neural networks; DTI: drug-target interaction; HTH, helix-turn-helix; IHF, integration host factor; LTTRs, LysR-type transcriptional regulators; MarR, multiple antibiotic resistance regulator; MRSA, methicillin resistant Staphylococcus aureus; MSA: multiple sequence alignment; NAP, nucleoid-associated protein; PROTACs, proteolysis targeting chimeras; RNAP, RNA polymerase; TF, transcription factor; TFD, transcription factor decoying; TFTRs, TetR-family transcriptional regulators; wHTH, winged helix-turn-helix.

Phosphorus starvation response and PhoB-independent utilization of organic phosphate sources by Salmonella enterica

IMPORTANCE

Phosphorus (P) is the fifth most abundant element in living cells. This element is acquired mainly as inorganic phosphate (Pi, PO4 3-). In enteric bacteria, P starvation activates a two-component signal transduction system which is composed of the membrane sensor protein PhoR and its cognate transcription regulator PhoB. PhoB, in turn, promotes the transcription of genes that help maintain Pi homeostasis. Here, we characterize the P starvation response of the bacterium Salmonella enterica. We determine the PhoB-dependent and independent transcriptional changes promoted by P starvation and identify proteins enabling the utilization of a range of organic substrates as sole P sources. We show that transcription and activity of a subset of these proteins are independent of PhoB and Pi availability. These results establish that Salmonella enterica can maintain Pi homeostasis and repress PhoB/PhoR activation even when cells are grown in medium lacking Pi.

The EcoCyc Database (2023)

EcoCyc is a bioinformatics database available online at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene product, metabolite, reaction, operon, and metabolic pathway. The database also includes information on the regulation of gene expression, E. coli gene essentiality, and nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for the analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc and can be executed online. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. Data generated from a whole-cell model that is parameterized from the latest data on EcoCyc are also available. This review outlines the data content of EcoCyc and of the procedures by which this content is generated.

A Gold Standard for Transcription Factor Regulatory Interactions in Escherichia coli K-12: Architecture of Evidence Types

Post-genomic implementations have expanded the experimental strategies to identify elements involved in the regulation of transcription initiation. As new methodologies emerge, a natural step is to compare their results with those from established methodologies, such as the classic methods of molecular biology used to characterize transcription factor binding sites, promoters, or transcription units. In the case of Escherichia coli K-12, the best-studied microorganism, for the last 30 years we have continuously gathered such knowledge from original scientific publications, and have organized it in two databases, RegulonDB and EcoCyc. Furthermore, since RegulonDB version 11.0 (1), we offer comprehensive datasets of binding sites from chromatin immunoprecipitation combined with sequencing (ChIP-seq), ChIP combined with exonuclease digestion and next-generation sequencing (ChIP-exo), genomic SELEX screening (gSELEX), and DNA affinity purification sequencing (DAP-seq) HT technologies, as well as additional datasets for transcription start sites, transcription units and RNA sequencing (RNA-seq) expression profiles. Here, we present for the first time an analysis of the sources of knowledge supporting the collection of transcriptional regulatory interactions (RIs) of E. coli K-12. An RI is formed by the transcription factor, its positive or negative effect on a promoter, a gene or transcription unit. We improved the evidence codes so that the specific methods are described, and we classified them into seven independent groups. This is the basis for our updated computation of confidence levels, weak, strong, or confirmed, for the collection of RIs. We compare the confidence levels of the RI collection before and after adding HT evidence illustrating how knowledge will change as more HT data and methods appear in the future. Users can generate subsets filtering out the method they want to benchmark and avoid circularity, or keep for instance only the confirmed interactions. The comparison of different HT methods with the available datasets indicate that ChIP-seq recovers the highest fraction (>70%) of binding sites present in RegulonDB followed by gSELEX, DAP-seq and ChIP-exo. There is no other genomic database that offers this comprehensive high-quality anatomy of evidence supporting a corpus of transcriptional regulatory interactions.

Selective RNA Processing and Stabilization are Multi-Layer and Stoichiometric Regulators of Gene Expression in Escherichia coli

Selective RNA processing and stabilization (SRPS) facilitates the differential expression of multiple genes in polycistronic operons. However, how the coordinated actions of SRPS-related enzymes affect stoichiometric regulation remains unclear. In the present study, the first genome-wide targetome analysis is reported of these enzymes in Escherichia coli, at a single-nucleotide resolution. A strictly linear relationship is observed between the RNA pyrophosphohydrolase processing ratio and scores assigned to the first three nucleotides of the primary transcript. Stem-loops associated with PNPase targetomes exhibit a folding free energy that is negatively correlated with the termination ratio of PNPase at the 3' end. More than one-tenth of the RNase E processing sites in the 5'-untranslated regions（UTR） form different stem-loops that affect ribosome-binding and translation efficiency. The effectiveness of the SRPS elements is validated using a dual-fluorescence reporter system. The findings highlight a multi-layer and quantitative regulatory method for optimizing the stoichiometric expression of genes in bacteria and promoting the application of SRPS in synthetic biology.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Genomic footprinting uncovers global transcription factor responses to amino acids in Escherichia coli

Our knowledge of transcriptional responses to changes in nutrient availability comes primarily from few well-studied transcription factors (TFs), often lacking an unbiased genome-wide perspective. Leveraging recent advances allowing bacterial genomic footprinting, we comprehensively mapped the genome-wide regulatory responses of Escherichia coli to exogenous leucine, methionine, alanine, and lysine. The global TF Lrp was found to individually sense three amino acids and mount three different target gene responses. Overall, 531 genes had altered RNA polymerase occupancy, and 32 TFs responded directly or indirectly to the presence of amino acids, including regulators of membrane and osmotic pressure homeostasis. About 70% of the detected TF-DNA interactions had not been reported before. We thus identified 682 previously unknown TF-binding locations, for a subset of which the involved TFs were identified by affinity purification. This comprehensive map of amino acid regulation illustrates the incompleteness of the known transcriptional regulation network, even in E. coli.

5'UTR G-quadruplex structure enhances translation in size dependent manner

Translation initiation in bacteria is frequently regulated by various structures in the 5' untranslated region (5'UTR). Previously, we demonstrated that G-quadruplex (G4) formation in non-template DNA enhances transcription. In this study, we aimed to explore how G4 formation in mRNA (RG4) at 5'UTR impacts translation using a T7-based in vitro translation system and in E. coli. We showed that RG4 strongly promotes translation efficiency in a size-dependent manner. Additionally, inserting a hairpin upstream of the RG4 further enhances translation efficiency, reaching up to a 12-fold increase. We found that the RG4-dependent effect is not due to increased ribosome affinity, ribosome binding site accessibility, or mRNA stability. We proposed a physical barrier model in which bulky structures in 5'UTR prevent ribosome dislodging and thereby increase the translation output. This study provides biophysical insights into the regulatory role of 5'UTR structures in bacterial translation, highlighting their potential applications in tuning gene expression.

Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Independent Component Analysis Reveals the Transcriptional Regulatory Modules in Bradyrhizobium diazoefficiens USDA110

The dynamic adaptation of bacteria to environmental changes is achieved through the coordinated expression of many genes, which constitutes a transcriptional regulatory network (TRN). Bradyrhizobium diazoefficiens USDA110 is an important model strain for the study of symbiotic nitrogen fixation (SNF), and its SNF ability largely depends on the TRN. In this study, independent component analysis was applied to 226 high-quality gene expression profiles of B. diazoefficiens USDA110 microarray datasets, from which 64 iModulons were identified. Using these iModulons and their condition-specific activity levels, we (1) provided new insights into the connection between the FixLJ-FixK2-FixK1 regulatory cascade and quorum sensing, (2) discovered the independence of the FixLJ-FixK2-FixK1 and NifA/RpoN regulatory cascades in response to oxygen, (3) identified the FixLJ-FixK2 cascade as a mediator connecting the FixK2-2 iModulon and the Phenylalanine iModulon, (4) described the differential activation of iModulons in B. diazoefficiens USDA110 under different environmental conditions, and (5) proposed a notion of active-TRN based on the changes in iModulon activity to better illustrate the relationship between gene regulation and environmental condition. In sum, this research offered an iModulon-based TRN for B. diazoefficiens USDA110, which formed a foundation for comprehensively understanding the intricate transcriptional regulation during SNF.

Localized modulation of DNA supercoiling, triggered by the Shigella anti-silencer VirB, is sufficient to relieve H-NS-mediated silencing

In Bacteria, nucleoid structuring proteins govern nucleoid dynamics and regulate transcription. In Shigella spp., at ≤30°C, the histone-like nucleoid structuring protein (H-NS) transcriptionally silences many genes on the large virulence plasmid. Upon a switch to 37°C, VirB, a DNA binding protein and key transcriptional regulator of Shigella virulence, is produced. VirB functions to counter H-NS-mediated silencing in a process called transcriptional anti-silencing. Here, we show that VirB mediates a loss of negative DNA supercoils from our plasmid-borne, VirB-regulated PicsP-lacZ reporter in vivo. The changes are not caused by a VirB-dependent increase in transcription, nor do they require the presence of H-NS. Instead, the VirB-dependent change in DNA supercoiling requires the interaction of VirB with its DNA binding site, a critical first step in VirB-dependent gene regulation. Using two complementary approaches, we show that VirB:DNA interactions in vitro introduce positive supercoils in plasmid DNA. Subsequently, by exploiting transcription-coupled DNA supercoiling, we reveal that a localized loss of negative supercoils is sufficient to alleviate H-NS-mediated transcriptional silencing independently of VirB. Together, our findings provide novel insight into VirB, a central regulator of Shigella virulence and, more broadly, a molecular mechanism that offsets H-NS-dependent silencing of transcription in bacteria.

An unexpected abundance of bidirectional promoters within Salmonella Typhimurium plasmids

Transcription of the DNA template, to generate an RNA message, is the first step in gene expression. The process initiates at DNA sequences called promoters. Conventionally, promoters have been considered to drive transcription in a specific direction. However, in recent work, we showed that many prokaryotic promoters can drive divergent transcription. This is a consequence of key DNA sequences for transcription initiation being inherently symmetrical. Here, we used global transcription start site mapping to determine the prevalence of such bidirectional promoters in Salmonella Typhimurium. Surprisingly, bidirectional promoters occur three times more frequently in plasmid components of the genome compared to chromosomal DNA. Implications for the evolution of promoter sequences are discussed.

Identification of Burkholderia cenocepacia non-coding RNAs expressed during Caenorhabditis elegans infection

Small non-coding RNAs (sRNAs) are key regulators of post-transcriptional gene expression in bacteria. Despite the identification of hundreds of bacterial sRNAs, their roles on bacterial physiology and virulence remain largely unknown, as is the case of bacteria of the Burkholderia cepacia complex (Bcc). Bcc is a group of opportunistic pathogens with relatively large genomes that can cause lethal lung infections amongst cystic fibrosis (CF) patients. To characterise sRNAs expressed by Bcc bacteria when infecting a host, the nematode Caenorhabditis elegans was used as an infection model by the epidemic CF strain B. cenocepacia J2315. A total of 108 new and 31 previously described sRNAs with a predicted Rho independent terminator were identified, most of them located on chromosome 1. RIT11b, a sRNA downregulated under C. elegans infection conditions, was shown to directly affect B. cenocepacia virulence, biofilm formation, and swimming motility. RIT11b overexpression reduced the expression of the direct targets dusA and pyrC, involved in biofilm formation, epithelial cell adherence, and chronic infections in other organisms. The in vitro direct interaction of RIT11b with the dusA and pyrC messengers was demonstrated by electrophoretic mobility shift assays. To the best of our knowledge this is the first report on the functional characterization of a sRNA directly involved in B. cenocepacia virulence. KEY POINTS: • 139 sRNAs expressed by B. cenocepacia during C. elegans infection were identified • The sRNA RIT11b affects B. cenocepacia virulence, biofilm formation, and motility • RIT11b directly binds to and regulates dusA and pyrC mRNAs.

An intracellular phosphorus-starvation signal activates the PhoB/PhoR two-component system in Salmonella enterica

Bacteria acquire P primarily as inorganic orthophosphate (Pi, PO43-). Once internalized, Pi is rapidly assimilated into biomass during the synthesis of ATP. Because Pi is essential, but excessive ATP is toxic, the acquisition of environmental Pi is tightly regulated. In the bacterium Salmonella enterica (Salmonella), growth in Pi-limiting environments activates the membrane sensor histidine kinase PhoR, leading to the phosphorylation of its cognate transcriptional regulator PhoB and subsequent transcription of genes involved in adaptations to low Pi. Pi limitation is thought to promote PhoR kinase activity by altering the conformation of a membrane signaling complex comprised by PhoR, the multicomponent Pi transporter system PstSACB and the regulatory protein PhoU. However, the identity of the low Pi signal and how it controls PhoR activity remain unknown. Here we characterize the PhoB-dependent and independent transcriptional changes elicited by Salmonella in response to P starvation, and identify PhoB-independent genes that are required for the utilization of several organic-P sources. We use this knowledge to identify the cellular compartment where the PhoR signaling complex senses the Pi-limiting signal. We demonstrate that the PhoB and PhoR signal transduction proteins can be maintained in an inactive state even when Salmonella is grown in media lacking Pi. Our results establish that PhoR activity is controlled by an intracellular signal resulting from P insufficiency.

No role for nuclear transcription regulators in mammalian mitochondria?

Although the mammalian mtDNA transcription machinery is simple and resembles bacteriophage systems, there are many reports that nuclear transcription regulators, as exemplified by MEF2D, MOF, PGC-1α, and hormone receptors, are imported into mammalian mitochondria and directly interact with the mtDNA transcription machinery. However, the supporting experimental evidence for this concept is open to alternate interpretations, and a main issue is the difficulty in distinguishing indirect regulation of mtDNA transcription, caused by altered nuclear gene expression, from direct intramitochondrial effects. We provide a critical discussion and experimental guidelines to stringently assess roles of intramitochondrial factors implicated in direct regulation of mammalian mtDNA transcription.

Natural variation in the hrpL promoter renders the phytopathogen Pseudomonas syringae pv. actinidiae nonpathogenic

The genetic basis underlying loss-of-virulence mutations that arise among natural phytopathogen populations is not well documented. In this study, we examined the virulence of 377 isolates of Pseudomonas syringae pv. actinidiae biovar 3 (Psa3) that were isolated from 76 kiwifruit orchards suffering from bacterial canker disease. Eighty-four nonpathogenic isolates were identified in 40 orchards. A nonpathogenic isolate G166 was found to be defective in hrpL transcription and the downstream type III secretion system (T3SS)-dependent phenotypes. Comparative genomics and complementary expression assay revealed that a single-base "G" insertion in the hrpL promoter blocks gene transcription by reducing promoter activity. The electrophoretic mobility shift assay showed that the genetic variation impairs σ⁵⁴ /promoter binding during gene transcription under hrp-inducing conditions, resulting in lower expression of hrpL. A PCR-restriction fragment length polymorphism assay was performed to trace the evolutionary history of this mutation, which revealed the independent onset of genetic variations in natural Psa3 populations. We also found that nonpathogenic variants outperformed virulent Psa3 bacteria for both epiphytic and apoplast colonization of kiwifruit leaves in mixed inoculations. Our study highlights a novel mechanism for loss of virulence in Psa3 and provides insight into bacterial adaptive evolution under natural settings.

Transcriptomics-Driven Characterization of LUZ100, a T7-like Pseudomonas Phage with Temperate Features

Autographiviridae is a diverse yet distinct family of bacterial viruses marked by a strictly lytic lifestyle and a generally conserved genome organization. Here, we characterized Pseudomonas aeruginosa phage LUZ100, a distant relative of type phage T7. LUZ100 is a podovirus with a limited host range which likely uses lipopolysaccharide (LPS) as a phage receptor. Interestingly, infection dynamics of LUZ100 indicated moderate adsorption rates and low virulence, hinting at temperate characteristics. This hypothesis was supported by genomic analysis, which showed that LUZ100 shares the conventional T7-like genome organization yet carries key genes associated with a temperate lifestyle. To unravel the peculiar characteristics of LUZ100, ONT-cappable-seq transcriptomics analysis was performed. These data provided a bird's-eye view of the LUZ100 transcriptome and enabled the discovery of key regulatory elements, antisense RNA, and transcriptional unit structures. The transcriptional map of LUZ100 also allowed us to identify new RNA polymerase (RNAP)-promoter pairs that can form the basis for biotechnological parts and tools for new synthetic transcription regulation circuitry. The ONT-cappable-seq data revealed that the LUZ100 integrase and a MarR-like regulator (proposed to be involved in the lytic/lysogeny decision) are actively cotranscribed in an operon. In addition, the presence of a phage-specific promoter transcribing the phage-encoded RNA polymerase raises questions on the regulation of this polymerase and suggests that it is interwoven with the MarR-based regulation. This transcriptomics-driven characterization of LUZ100 supports recent evidence that T7-like phages should not automatically be assumed to have a strictly lytic life cycle. IMPORTANCE Bacteriophage T7, considered the "model phage" of the Autographiviridae family, is marked by a strictly lytic life cycle and conserved genome organization. Recently, novel phages within this clade have emerged which display characteristics associated with a temperate life cycle. Screening for temperate behavior is of utmost importance in fields like phage therapy, where strictly lytic phages are generally required for therapeutic applications. In this study, we applied an omics-driven approach to characterize the T7-like Pseudomonas aeruginosa phage LUZ100. These results led to the identification of actively transcribed lysogeny-associated genes in the phage genome, pointing out that temperate T7-like phages are emerging more frequent than initially thought. In short, the combination of genomics and transcriptomics allowed us to obtain a better understanding of the biology of nonmodel Autographiviridae phages, which can be used to optimize the implementation of phages and their regulatory elements in phage therapy and biotechnological applications, respectively.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

Localized modulation of DNA supercoiling, triggered by the Shigella anti-silencer VirB, is sufficient to relieve H-NS-mediated silencing

In Bacteria, nucleoid structuring proteins govern nucleoid dynamics and regulate transcription. In Shigella spp ., at ≤ 30 °C, the histone-like nucleoid structuring protein (H-NS) transcriptionally silences many genes on the large virulence plasmid. Upon a switch to 37 °C, VirB, a DNA binding protein and key transcriptional regulator of Shigella virulence, is produced. VirB functions to counter H-NS-mediated silencing in a process called transcriptional anti-silencing. Here, we show that VirB mediates a loss of negative DNA supercoils from our plasmid-borne, VirB-regulated PicsP-lacZ reporter, in vivo . The changes are not caused by a VirB-dependent increase in transcription, nor do they require the presence of H-NS. Instead, the VirB-dependent change in DNA supercoiling requires the interaction of VirB with its DNA binding site, a critical first step in VirB-dependent gene regulation. Using two complementary approaches, we show that VirB:DNA interactions in vitro introduce positive supercoils in plasmid DNA. Subsequently, by exploiting transcription-coupled DNA supercoiling, we reveal that a localized loss of negative supercoils is sufficient to alleviate H-NS-mediated transcriptional silencing, independently of VirB. Together, our findings provide novel insight into VirB, a central regulator of Shigella virulence and more broadly, a molecular mechanism that offsets H-NS-dependent silencing of transcription in bacteria.

Variable DNA topology is an epigenetic generator of physiological heterogeneity in bacterial populations

Transcription is a noisy and stochastic process that produces sibling-to-sibling variations in physiology across a population of genetically identical cells. This pattern of diversity reflects, in part, the burst-like nature of transcription. Transcription bursting has many causes and a failure to remove the supercoils that accumulate in DNA during transcription elongation is an important contributor. Positive supercoiling of the DNA ahead of the transcription elongation complex can result in RNA polymerase stalling if this DNA topological roadblock is not removed. The relaxation of these positive supercoils is performed by the ATP-dependent type II topoisomerases DNA gyrase and topoisomerase IV. Interference with the action of these topoisomerases involving, inter alia, topoisomerase poisons, fluctuations in the [ATP]/[ADP] ratio, and/or the intervention of nucleoid-associated proteins with GapR-like or YejK-like activities, may have consequences for the smooth operation of the transcriptional machinery. Antibiotic-tolerant (but not resistant) persister cells are among the phenotypic outliers that may emerge. However, interference with type II topoisomerase activity can have much broader consequences, making it an important epigenetic driver of physiological diversity in the bacterial population.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, Emax, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with Emax but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and Emax values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Identifying the Gene Regulatory Network of the Starvation-Induced Transcriptional Activator Nla28

Myxococcus xanthus copes with starvation by producing fruiting bodies filled with dormant and stress-resistant spores. Here, we aimed to better define the gene regulatory network associated with Nla28, a transcriptional activator/enhancer binding protein (EBP) and a key regulator of the early starvation response. Previous work showed that Nla28 directly regulates EBP genes that are important for fruiting body development. However, the Nla28 regulatory network is likely to be much larger because hundreds of starvation-induced genes are downregulated in a nla28 mutant strain. To identify candidates for direct Nla28-mediated transcription, we analyzed the downregulated genes using a bioinformatics approach. Nine potential Nla28 target promoters (29 genes) were discovered. The results of in vitro promoter binding assays, coupled with in vitro and in vivo mutational analyses, suggested that the nine promoters along with three previously identified EBP gene promoters were indeed in vivo targets of Nla28. These results also suggested that Nla28 used tandem, imperfect repeats of an 8-bp sequence for promoter binding. Interestingly, eight of the new Nla28 target promoters were predicted to be intragenic. Based on mutational analyses, the newly identified Nla28 target loci contained at least one gene that was important for starvation-induced development. Most of these loci contained genes predicted to be involved in metabolic or defense-related functions. Using the consensus Nla28 binding sequence, bioinformatics, and expression profiling, 58 additional promoters and 102 genes were tagged as potential Nla28 targets. Among these putative Nla28 targets, functions, such as regulatory, metabolic, and cell envelope biogenesis, were assigned to many genes. IMPORTANCE In bacteria, starvation leads to profound changes in behavior and physiology. Some of these changes have economic and health implications because the starvation response has been linked to the formation of biofilms, virulence, and antibiotic resistance. To better understand how starvation contributes to changes in bacterial physiology and resistance, we identified the putative starvation-induced gene regulatory network associated with Nla28, a transcriptional activator from the bacterium Myxoccocus xanthus. We determined the mechanism by which starvation-responsive genes were activated by Nla28 and showed that several of the genes were important for the formation of a highly resistant cell type.

The Spatial Organization of Bacterial Transcriptional Regulatory Networks

The transcriptional regulatory network (TRN) is the central pivot of a prokaryotic organism to receive, process and respond to internal and external environmental information. However, little is known about its spatial organization so far. In recent years, chromatin interaction data of bacteria such as Escherichia coli and Bacillus subtilis have been published, making it possible to study the spatial organization of bacterial transcriptional regulatory networks. By combining TRNs and chromatin interaction data of E. coli and B. subtilis, we explored the spatial organization characteristics of bacterial TRNs in many aspects such as regulation directions (positive and negative), central nodes (hubs, bottlenecks), hierarchical levels (top, middle, bottom) and network motifs (feed-forward loops and single input modules) of the TRNs and found that the bacterial TRNs have a variety of stable spatial organization features under different physiological conditions that may be closely related with biological functions. Our findings provided new insights into the connection between transcriptional regulation and the spatial organization of chromosome in bacteria and might serve as a factual foundation for trying spatial-distance-based gene circuit design in synthetic biology.

MsmR1, a global transcription factor, regulates polymyxin synthesis and carbohydrate metabolism in Paenibacillus polymyxa SC2

The multiple-sugar metabolism regulator (MsmR), a transcription factor belonging to the AraC/XylS family, participates in polysaccharide metabolism and virulence. However, the transcriptional regulatory mechanisms of MsmR1 in Paenibacillus polymyxa remain unclear. In this study, knocking out msmR1 was found to reduce polymyxin synthesis by the SC2-M1 strain. Chromatin immunoprecipitation assay with sequencing (ChIP-seq) revealed that most enriched pathway was that of carbohydrate metabolism. Additionally, electromobility shift assays (EMSA) confirmed the direct interaction between MsmR1 and the promoter regions of oppC3, sucA, sdr3, pepF, yycN, PPSC2_23180, pppL, and ydfp. MsmR1 stimulates polymyxin biosynthesis by directly binding to the promoter regions of oppC3 and sdr3, while also directly regulating sucA and influencing the citrate cycle (TCA cycle). In addition, MsmR1 directly activates pepF and was beneficial for spore and biofilm formation. These results indicated that MsmR1 could regulate carbohydrate and amino acid metabolism, and indirectly affect biological processes such as polymyxin synthesis, biofilm formation, and motility. Moreover, MsmR1 could be autoregulated. Hence, this study expand the current knowledge of MsmR1 and will be beneficial for the application of P. polymyxa SC2 in the biological control against the certain pathogens in pepper.

Genetic context effects can override canonical cis regulatory elements in Escherichia coli

Recent experiments have shown that in addition to control by cis regulatory elements, the local chromosomal context of a gene also has a profound impact on its transcription. Although this chromosome-position dependent expression variation has been empirically mapped at high-resolution, the underlying causes of the variation have not been elucidated. Here, we demonstrate that 1 kb of flanking, non-coding synthetic sequences with a low frequency of guanosine and cytosine (GC) can dramatically reduce reporter expression compared to neutral and high GC-content flanks in Escherichia coli. Natural and artificial genetic context can have a similarly strong effect on reporter expression, regardless of cell growth phase or medium. Despite the strong reduction in the maximal expression level from the fully-induced reporter, low GC synthetic flanks do not affect the time required to reach the maximal expression level after induction. Overall, we demonstrate key determinants of transcriptional propensity that appear to act as tunable modulators of transcription, independent of regulatory sequences such as the promoter. These findings provide insight into the regulation of naturally occurring genes and an independent control for optimizing expression of synthetic biology constructs.

Automated model-predictive design of synthetic promoters to control transcriptional profiles in bacteria

Transcription rates are regulated by the interactions between RNA polymerase, sigma factor, and promoter DNA sequences in bacteria. However, it remains unclear how non-canonical sequence motifs collectively control transcription rates. Here, we combine massively parallel assays, biophysics, and machine learning to develop a 346-parameter model that predicts site-specific transcription initiation rates for any σ70 promoter sequence, validated across 22132 bacterial promoters with diverse sequences. We apply the model to predict genetic context effects, design σ70 promoters with desired transcription rates, and identify undesired promoters inside engineered genetic systems. The model provides a biophysical basis for understanding gene regulation in natural genetic systems and precise transcriptional control for engineering synthetic genetic systems.

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.