Direct and indirect effects of H-NS and Fis on global gene expression control in Escherichia coli

Bacterial epigenetics and its implication for agriculture, probiotics development, and biotechnology design

In their natural habitats, organisms encounter numerous external stimuli and must be able to sense and adapt to those stimuli to survive. Unlike mutations, epigenetic changes do not alter the underlying DNA sequence. Instead, they create modifications that promote or silence gene expression. Bacillus subtilis has long been a model organism in studying genetics and development. It is beneficial for numerous biotechnological applications where it is included as a probiotic, in fermentation, or in bio-concrete design. This bacterium has also emerged recently as a model organism for studying bacterial epigenetic adaptation. In this review, we examine the evolving knowledge of epigenetic regulation (restriction-modification systems (RM), orphan methyltransferases, and chromosome condensation) in B. subtilis and related bacteria, and utilize it as a case study to test their potential roles and future applications in genetic engineering and microbial biotechnology. Finally, we suggest how the implementation of these fundamental findings promotes the design of synthetic epigenetic memory circuits and their future applications in agriculture, medicine, and biotechnology.

Structure of the E. coli nucleoid-associated protein YejK reveals a novel DNA binding clamp

Nucleoid-associated proteins (NAPs) play central roles in bacterial chromosome organization and DNA processes. The Escherichia coli YejK protein is a highly abundant, yet poorly understood NAP. YejK proteins are conserved among Gram-negative bacteria but show no homology to any previously characterized DNA-binding protein. Hence, how YejK binds DNA is unknown. To gain insight into YejK structure and its DNA binding mechanism we performed biochemical and structural analyses on the E. coli YejK protein. Biochemical assays demonstrate that, unlike many NAPs, YejK does not show a preference for AT-rich DNA and binds non-sequence specifically. A crystal structure revealed YejK adopts a novel fold comprised of two domains. Strikingly, each of the domains harbors an extended arm that mediates dimerization, creating an asymmetric clamp with a 30 Å diameter pore. The lining of the pore is electropositive and mutagenesis combined with fluorescence polarization assays support DNA binding within the pore. Finally, our biochemical analyses on truncated YejK proteins suggest a mechanism for YejK clamp loading. Thus, these data reveal YejK contains a newly described DNA-binding motif that functions as a novel clamp.

Bacterial chromatin proteins, transcription, and DNA topology: Inseparable partners in the control of gene expression

DNA in bacterial chromosomes is organized into higher-order structures by DNA-binding proteins called nucleoid-associated proteins (NAPs) or bacterial chromatin proteins (BCPs). BCPs often bind to or near DNA loci transcribed by RNA polymerase (RNAP) and can either increase or decrease gene expression. To understand the mechanisms by which BCPs alter transcription, one must consider both steric effects and the topological forces that arise when DNA deviates from its fully relaxed double-helical structure. Transcribing RNAP creates DNA negative (-) supercoils upstream and positive (+) supercoils downstream whenever RNAP and DNA are unable to rotate freely. This (-) and (+) supercoiling generates topological forces that resist forward translocation of DNA through RNAP unless the supercoiling is constrained by BCPs or relieved by topoisomerases. BCPs also may enhance topological stress and overall can either inhibit or aid transcription. Here, we review current understanding of how RNAP, BCPs, and DNA topology interplay to control gene expression.

Spatio-temporal organization of the E. coli chromosome from base to cellular length scales

Escherichia coli has been a vital model organism for studying chromosomal structure, thanks, in part, to its small and circular genome (4.6 million base pairs) and well-characterized biochemical pathways. Over the last several decades, we have made considerable progress in understanding the intricacies of the structure and subsequent function of the E. coli nucleoid. At the smallest scale, DNA, with no physical constraints, takes on a shape reminiscent of a randomly twisted cable, forming mostly random coils but partly affected by its stiffness. This ball-of-spaghetti-like shape forms a structure several times too large to fit into the cell. Once the physiological constraints of the cell are added, the DNA takes on overtwisted (negatively supercoiled) structures, which are shaped by an intricate interplay of many proteins carrying out essential biological processes. At shorter length scales (up to about 1 kb), nucleoid-associated proteins organize and condense the chromosome by inducing loops, bends, and forming bridges. Zooming out further and including cellular processes, topological domains are formed, which are flanked by supercoiling barriers. At the megabase-scale both large, highly self-interacting regions (macrodomains) and strong contacts between distant but co-regulated genes have been observed. At the largest scale, the nucleoid forms a helical ellipsoid. In this review, we will explore the history and recent advances that pave the way for a better understanding of E. coli chromosome organization and structure, discussing the cellular processes that drive changes in DNA shape, and what contributes to compaction and formation of dynamic structures, and in turn how bacterial chromatin affects key processes such as transcription and replication.

Elucidation of Spatial Positioning of Ribosomes around Chromosome in Escherichia coli Cytoplasm via a Data-Informed Polymer-Based Model

The spatial arrangement of ribosomes and chromosome in Escherichia coli's cytoplasm challenges conventional wisdom. Contrary to the notion of ribosomes acting as inert crowders to the chromosome in the cytoplasm, here we propose a nuanced view by integrating a wide array of experimental data sets into a polymer-based computer model. A set of data-informed computer simulations determines that a delicate balance of attractive and repulsive interactions between ribosomes and the chromosome is required in order to reproduce experimentally obtained linear densities and brings forth the view that ribosomes are not mere inert crowders in the cytoplasm. The model finds that the ribosomes represent themselves as a poor solvent for the chromosome with a 50 nm mesh size, consistent with previous experimental analysis. Our multidimensional analysis of ribosome distribution, both free (30S and 50S) and bound (70S polysome), uncovers a relatively less pronounced segregation pattern than previously thought. Notably, we identify a ribosome-rich central region within the innermost core of the nucleoid. Moreover, our exploration of the chromosome mesh size and the conformation of bound ribosomes suggests that these ribosomes maintain elongated shapes, enabling them to navigate through the chromosome mesh and access the central core. This dynamic localization challenges the static segregation model and underscores the pivotal role of ribosome-chromosome interactions in cellular media.

Comprehensive Characterization of fucAO Operon Activation in Escherichia coli

Wildtype Escherichia coli cells cannot grow on L-1,2-propanediol, as the fucAO operon within the fucose (fuc) regulon is thought to be silent in the absence of L-fucose. Little information is available concerning the transcriptional regulation of this operon. Here, we first confirm that fucAO operon expression is highly inducible by fucose and is primarily attributable to the upstream operon promoter, while the fucO promoter within the 3'-end of fucA is weak and uninducible. Using 5'RACE, we identify the actual transcriptional start site (TSS) of the main fucAO operon promoter, refuting the originally proposed TSS. Several lines of evidence are provided showing that the fucAO locus is within a transcriptionally repressed region on the chromosome. Operon activation is dependent on FucR and Crp but not SrsR. Two Crp-cAMP binding sites previously found in the regulatory region are validated, where the upstream site plays a more critical role than the downstream site in operon activation. Furthermore, two FucR binding sites are identified, where the downstream site near the first Crp site is more important than the upstream site. Operon transcription relies on Crp-cAMP to a greater degree than on FucR. Our data strongly suggest that FucR mainly functions to facilitate the binding of Crp to its upstream site, which in turn activates the fucAO promoter by efficiently recruiting RNA polymerase.

Bacterial nucleoid is a riddle wrapped in a mystery inside an enigma

Bacterial chromosome, the nucleoid, is traditionally modeled as a rosette of DNA mega-loops, organized around proteinaceous central scaffold by nucleoid-associated proteins (NAPs), and mixed with the cytoplasm by transcription and translation. Electron microscopy of fixed cells confirms dispersal of the cloud-like nucleoid within the ribosome-filled cytoplasm. Here, I discuss evidence that the nucleoid in live cells forms DNA phase separate from riboprotein phase, the "riboid." I argue that the nucleoid-riboid interphase, where DNA interacts with NAPs, transcribing RNA polymerases, nascent transcripts, and ssRNA chaperones, forms the transcription zone. An active part of phase separation, transcription zone enforces segregation of the centrally positioned information phase (the nucleoid) from the surrounding action phase (the riboid), where translation happens, protein accumulates, and metabolism occurs. I speculate that HU NAP mostly tiles up the nucleoid periphery-facilitating DNA mobility but also supporting transcription in the interphase. Besides extruding plectonemically supercoiled DNA mega-loops, condensins could compact them into solenoids of uniform rings, while HU could support rigidity and rotation of these DNA rings. The two-phase cytoplasm arrangement allows the bacterial cell to organize the central dogma activities, where (from the cell center to its periphery) DNA replicates and segregates, DNA is transcribed, nascent mRNA is handed over to ribosomes, mRNA is translated into proteins, and finally, the used mRNA is recycled into nucleotides at the inner membrane. The resulting information-action conveyor, with one activity naturally leading to the next one, explains the efficiency of prokaryotic cell design-even though its main intracellular transportation mode is free diffusion.

Lsr2, a pleiotropic regulator at the core of the infectious strategy of Mycobacterium abscessus

Mycobacterium abscessus is a non-tuberculous mycobacterium, causing lung infections in cystic fibrosis patients. During pulmonary infection, M. abscessus switches from smooth (Mabs-S) to rough (Mabs-R) morphotypes, the latter being hyper-virulent. Previously, we isolated the lsr2 gene as differentially expressed during S-to-R transition. lsr2 encodes a pleiotropic transcription factor that falls under the superfamily of nucleoid-associated proteins. Here, we used two functional genomic methods, RNA-seq and chromatin immunoprecipitation-sequencing (ChIP-seq), to elucidate the molecular role of Lsr2 in the pathobiology of M. abscessus. Transcriptomic analysis shows that Lsr2 differentially regulates gene expression across both morphotypes, most of which are involved in several key cellular processes of M. abscessus, including host adaptation and antibiotic resistance. These results were confirmed through quantitative real-time PCR, as well as by minimum inhibitory concentration tests and infection tests on macrophages in the presence of antibiotics. ChIP-seq analysis revealed that Lsr2 extensively binds the M. abscessus genome at AT-rich sequences and appears to form long domains that participate in the repression of its target genes. Unexpectedly, the genomic distribution of Lsr2 revealed no distinctions between Mabs-S and Mabs-R, implying more intricate mechanisms at play for achieving target selectivity.IMPORTANCELsr2 is a crucial transcription factor and chromosome organizer involved in intracellular growth and virulence in the smooth and rough morphotypes of Mycobacterium abscessus. Using RNA-seq and chromatin immunoprecipitation-sequencing (ChIP-seq), we investigated the molecular role of Lsr2 in gene expression regulation along with its distribution on M. abscessus genome. Our study demonstrates the pleiotropic regulatory role of Lsr2, regulating the expression of many genes coordinating essential cellular and molecular processes in both morphotypes. In addition, we have elucidated the role of Lsr2 in antibiotic resistance both in vitro and in vivo, where lsr2 mutant strains display heightened sensitivity to antibiotics. Through ChIP-seq, we reported the widespread distribution of Lsr2 on M. abscessus genome, revealing a direct repressive effect due to its extensive binding on promoters or coding sequences of its targets. This study unveils the significant regulatory role of Lsr2, intricately intertwined with its function in shaping the organization of the M. abscessus genome.

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.

Deciphering Sequence-Specific DNA Binding by H-NS Using Molecular Simulation

H-NS is a DNA organizing protein that occurs in Gram-negative bacteria. It can form long filaments between two DNA duplexes by first binding to a high-affinity AT-rich nucleotide sequence and extending from there. Using molecular dynamics simulations and steered molecular dynamics, we are able to determine the free energy of formation and dissociation of a protein-DNA complex comprising an H-NS DNA-binding domain and a specific nucleotide sequence. The molecular dynamics simulations allow detailed characterization of the interactions between the protein and a specific nucleotide sequence. To quantify the strength of the interaction, we employ an additional potential based on protein-DNA contacts to speed up dissociation of the protein-DNA complex. The work required for the dissociation results in an estimate of the free energy of dissociation/complex formation. Our protocol can provide quantitative prediction of protein-DNA complex stability, while also providing high-resolution insights into recognition mechanisms. In this chapter, we have used this approach to quantify the sequence specificity of H-NS DNA-binding domains to various nucleotide sequences, thus elucidating the mechanism with which H-NS can specifically bind to AT-rich DNA.

GeF-seq: A Simple Procedure for Base-Pair Resolution ChIP-seq

Nucleotide sequences recognized and bound by DNA-binding proteins (DBPs) are critical to controlling and maintaining gene expression, replication, chromosome segregation, cell division, and nucleoid structure in bacterial cells. Therefore, determination of the binding sequences of DBPs is important not only to study DBP recognition mechanisms but also to understand the fundamentals of cell homeostasis. While ChIP-seq analysis appears to be an effective way to determine DBP binding sites on the genome, the resolution is sometimes not sufficient to identify the sites precisely. Here we introduce a simple and effective method named Genome Footprinting with high-throughput sequencing (GeF-seq) to determine binding sites of DBPs with single base-pair resolution. GeF-seq detects binding sites of DBPs as sharp peaks and thus makes it possible to identify the recognition sequence in each "binding peak" more easily and accurately compared to the common ChIP-seq.

Shedding Light on Bacterial Chromosome Structure: Exploring the Significance of 3C-Based Approaches

The complex architecture of DNA within living organisms is essential for maintaining the genetic information that dictates their functions and characteristics. Among the many complexities of genetic material organization, the folding and arrangement of DNA into chromosomes play a critical role in regulating gene expression, replication, and other essential cellular processes. Bacteria, despite their apparently simple cellular structure, exhibit a remarkable level of chromosomal organization that influences their adaptability and survival in diverse environments. Understanding the three-dimensional arrangement of bacterial chromosomes has long been a challenge due to technical limitations, but the development of Chromosome Conformation Capture (3C) methods revolutionized our ability to explore the hierarchical structure and the dynamics of bacterial genomes. Here, we review the major advances in the field of bacterial chromosome structure using 3C technology over the past decade.

Macromolecular Crowding and DNA: Bridging the Gap between In Vitro and In Vivo

The cellular environment is highly crowded, with up to 40% of the volume fraction of the cell occupied by various macromolecules. Most laboratory experiments take place in dilute buffer solutions; by adding various synthetic or organic macromolecules, researchers have begun to bridge the gap between in vitro and in vivo measurements. This is a review of the reported effects of macromolecular crowding on the compaction and extension of DNA, the effect of macromolecular crowding on DNA kinetics, and protein-DNA interactions. Theoretical models related to macromolecular crowding and DNA are briefly reviewed. Gaps in the literature, including the use of biologically relevant crowders, simultaneous use of multi-sized crowders, empirical connections between macromolecular crowding and liquid-liquid phase separation of nucleic materials are discussed.

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.

The heparin-binding hemagglutinin protein of Mycobacterium tuberculosis is a nucleoid-associated protein

Nucleoid-associated proteins (NAPs) regulate multiple cellular processes such as gene expression, virulence, and dormancy throughout bacterial species. NAPs help in the survival and adaptation of Mycobacterium tuberculosis (Mtb) within the host. Fourteen NAPs have been identified in Escherichia coli; however, only seven NAPs are documented in Mtb. Given its complex lifestyle, it is reasonable to assume that Mtb would encode for more NAPs. Using bioinformatics tools and biochemical experiments, we have identified the heparin-binding hemagglutinin (HbhA) protein of Mtb as a novel sequence-independent DNA-binding protein which has previously been characterized as an adhesion molecule required for extrapulmonary dissemination. Deleting the carboxy-terminal domain of HbhA resulted in a complete loss of its DNA-binding activity. Atomic force microscopy showed HbhA-mediated architectural modulations in the DNA, which may play a regulatory role in transcription and genome organization. Our results showed that HbhA colocalizes with the nucleoid region of Mtb. Transcriptomics analyses of a hbhA KO strain revealed that it regulates the expression of ∼36% of total and ∼29% of essential genes. Deletion of hbhA resulted in the upregulation of ∼73% of all differentially expressed genes, belonging to multiple pathways suggesting it to be a global repressor. The results show that HbhA is a nonessential NAP regulating gene expression globally and acting as a plausible transcriptional repressor.

The environmentally-regulated interplay between local three-dimensional chromatin organisation and transcription of proVWX in E. coli

Nucleoid associated proteins (NAPs) maintain the architecture of bacterial chromosomes and regulate gene expression. Thus, their role as transcription factors may involve three-dimensional chromosome re-organisation. While this model is supported by in vitro studies, direct in vivo evidence is lacking. Here, we use RT-qPCR and 3C-qPCR to study the transcriptional and architectural profiles of the H-NS (histone-like nucleoid structuring protein)-regulated, osmoresponsive proVWX operon of Escherichia coli at different osmolarities and provide in vivo evidence for transcription regulation by NAP-mediated chromosome re-modelling in bacteria. By consolidating our in vivo investigations with earlier in vitro and in silico studies that provide mechanistic details of how H-NS re-models DNA in response to osmolarity, we report that activation of proVWX in response to a hyperosmotic shock involves the destabilization of H-NS-mediated bridges anchored between the proVWX downstream and upstream regulatory elements (DRE and URE), and between the DRE and ygaY that lies immediately downstream of proVWX. The re-establishment of these bridges upon adaptation to hyperosmolarity represses the operon. Our results also reveal additional structural features associated with changes in proVWX transcript levels such as the decompaction of local chromatin upstream of the operon, highlighting that further complexity underlies the regulation of this model operon. H-NS and H-NS-like proteins are wide-spread amongst bacteria, suggesting that chromosome re-modelling may be a typical feature of transcriptional control in bacteria.

The roles of nucleoid-associated proteins and topoisomerases in chromosome structure, strand segregation, and the generation of phenotypic heterogeneity in bacteria

How to adapt to a changing environment is a fundamental, recurrent problem confronting cells. One solution is for cells to organize their constituents into a limited number of spatially extended, functionally relevant, macromolecular assemblies or hyperstructures, and then to segregate these hyperstructures asymmetrically into daughter cells. This asymmetric segregation becomes a particularly powerful way of generating a coherent phenotypic diversity when the segregation of certain hyperstructures is with only one of the parental DNA strands and when this pattern of segregation continues over successive generations. Candidate hyperstructures for such asymmetric segregation in prokaryotes include those containing the nucleoid-associated proteins (NAPs) and the topoisomerases. Another solution to the problem of creating a coherent phenotypic diversity is by creating a growth-environment-dependent gradient of supercoiling generated along the replication origin-to-terminus axis of the bacterial chromosome. This gradient is modulated by transcription, NAPs, and topoisomerases. Here, we focus primarily on two topoisomerases, TopoIV and DNA gyrase in Escherichia coli, on three of its NAPs (H-NS, HU, and IHF), and on the single-stranded binding protein, SSB. We propose that the combination of supercoiling-gradient-dependent and strand-segregation-dependent topoisomerase activities result in significant differences in the supercoiling of daughter chromosomes, and hence in the phenotypes of daughter cells.

Search for Dispersed Repeats in Bacterial Genomes Using an Iterative Procedure

We have developed a de novo method for the identification of dispersed repeats based on the use of random position-weight matrices (PWMs) and an iterative procedure (IP). The created algorithm (IP method) allows detection of dispersed repeats for which the average number of substitutions between any two repeats per nucleotide (x) is less than or equal to 1.5. We have shown that all previously developed methods and algorithms (RED, RECON, and some others) can only find dispersed repeats for x ≤ 1.0. We applied the IP method to find dispersed repeats in the genomes of E. coli and nine other bacterial species. We identify three families of approximately 1.09 × 10⁶, 0.64 × 10⁶, and 0.58 × 10⁶ DNA bases, respectively, constituting almost 50% of the complete E. coli genome. The length of the repeats is in the range of 400 to 600 bp. Other analyzed bacterial genomes contain one to three families of dispersed repeats with a total number of 10³ to 6 × 10³ copies. The existence of such highly divergent repeats could be associated with the presence of a single-type triplet periodicity in various genes or with the packing of bacterial DNA into a nucleoid.

Genome-Wide Mapping of the Escherichia coli PhoB Regulon Reveals Many Transcriptionally Inert, Intragenic Binding Sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of, genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert and hence are tolerated as genomic "noise." IMPORTANCE Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of nonregulatory, intragenic binding sites for transcription factors and that these binding sites are not under selective pressure.

The morphology and metabolic changes of Actinobacillus pleuropneumoniae during its growth as a biofilm

Actinobacillus pleuropneumoniae is an important swine respiratory pathogen. Previous studies have suggested that growth as a biofilm is a natural state of A. pleuropneumoniae infection. To understand the survival features involved in the biofilm state, the growth features, morphology and gene expression profiles of planktonic and biofilm A. pleuropneumoniae were compared. A. pleuropneumoniae in biofilms showed reduced viability but maintained the presence of extracellular polymeric substances (EPS) after late log-phase. Under the microscope, bacteria in biofilms formed dense aggregated structures that were connected by abundant EPS, with reduced condensed chromatin. By construction of Δpga and ΔdspB mutants, polymeric β-1,6-linked N-acetylglucosamine and dispersin B were confirmed to be critical for normal biofilm formation. RNA-seq analysis indicated that, compared to their planktonic counterparts, A. pleuropneumoniae in biofilms had an extensively altered transcriptome. Carbohydrate metabolism, energy metabolism and translation were significantly repressed, while fermentation and genes contributing to EPS synthesis and translocation were up-regulated. The regulators Fnr (HlyX) and Fis were found to be up-regulated and their binding motifs were identified in the majority of the differentially expressed genes, suggesting their coordinated global role in regulating biofilm metabolism. By comparing the transcriptome of wild-type biofilm and Δpga, the utilization of oligosaccharides, iron and sulfur and fermentation were found to be important in adhesion and aggregation during biofilm formation. Additionally, when used as inocula, biofilm bacteria showed reduced virulence in mouse, compared with planktonic grown cells. Thus, these results have identified new facets of A. pleuropneumoniae biofilm maintenance and regulation.

Morphological Diversity of Dps Complex with Genomic DNA

In response to adverse environmental factors, Escherichia coli cells actively produce Dps proteins which form ordered complexes (biocrystals) with bacterial DNA to protect the genome. The effect of biocrystallization has been described extensively in the scientific literature; furthermore, to date, the structure of the Dps-DNA complex has been established in detail in vitro using plasmid DNA. In the present work, for the first time, Dps complexes with E. coli genomic DNA were studied in vitro using cryo-electron tomography. We demonstrate that genomic DNA forms one-dimensional crystals or filament-like assemblies which transform into weakly ordered complexes with triclinic unit cells, similar to what is observed for plasmid DNA. Changing such environmental factors as pH and KCl and MgCl₂ concentrations leads to the formation of cylindrical structures.

Revealing Causes for False-Positive and False-Negative Calling of Gene Essentiality in Escherichia coli Using Transposon Insertion Sequencing

The massive sequencing of transposon insertion mutant libraries (Tn-Seq) represents a commonly used method to determine essential genes in bacteria. Using a hypersaturated transposon mutant library consisting of 400,096 unique Tn insertions, 523 genes were classified as essential in Escherichia coli K-12 MG1655. This provided a useful genome-wide gene essentiality landscape for rapidly identifying 233 of 301 essential genes previously validated by a knockout study. However, there was a discrepancy in essential gene sets determined by conventional gene deletion methods and Tn-Seq, although different Tn-Seq studies reported different extents of discrepancy. We have elucidated two causes of this discrepancy. First, 68 essential genes not detected by Tn-Seq contain nonessential subgenic domains that are tolerant to transposon insertion, which leads to the false assignment of an essential gene as a nonessential or dispensable gene. These genes exhibited a high level of transposon insertion in their subgenic nonessential domains. In contrast, 290 genes were additionally categorized as essential by Tn-Seq, although their knockout mutants were available. The comparative analysis of Tn-Seq and high-resolution footprinting of nucleoid-associated proteins (NAPs) revealed that a protein-DNA interaction hinders transposon insertion. We identified 213 false-positive genes caused by NAP-genome interactions. These two limitations have to be considered when addressing essential bacterial genes using Tn-Seq. Furthermore, a comparative analysis of high-resolution Tn-Seq with other data sets is required for a more accurate determination of essential genes in bacteria. IMPORTANCE Transposon mutagenesis is an efficient way to explore gene essentiality of a bacterial genome. However, there was a discrepancy between the essential gene set determined by transposon mutagenesis and that determined using single-gene knockout strains. In this study, we generated a hypersaturated Escherichia coli transposon mutant library comprising approximately 400,000 different mutants. Determination of transposon insertion sites using next-generation sequencing provided a high-resolution essentiality landscape of the E. coli genome. We identified false negatives of essential gene discovery due to the permissive insertion of transposons in the C-terminal region. Comparisons between the transposon insertion landscape with binding profiles of DNA-binding proteins revealed interference of nucleoid-associated proteins to transposon insertion, generating false positives of essential gene discovery. Consideration of these findings is required to avoid the misinterpretation of transposon mutagenesis results.

An Attempt of Seeking Favorable Binding Free Energy Prediction Schemes Considering the Entropic Effect on Fis-DNA Binding

Protein-DNA binding mechanisms in a complex manner are essential for understanding many biological processes. Over the past decades, numerous experiments and calculations have analyzed the specificity of protein-DNA binding. However, the accuracy of binding free energy prediction for multi-base DNA systems still needs to be improved. Fis is a DNA-binding protein that regulates various transcription and recombination reactions. In the present work, we tested several methods of predict binding free energy based on this system to find a favorable prediction scheme and explore the binding mechanism of Fis protein and DNA. Two solvent models (explicit and implicit solvent models) were chosen for the dynamics process, and the predicted binding free energy was more accurate under the explicit solvent model. When different Poisson-Boltzmann/Generalized Born (PB/GB) models were tested for DNA force fields (BSC1 and OL15), it was found that the binding free energy predicted by the selected OL15 force field performed better and the correlation between predicted and experimental values was improved with the increasing interior dielectric constant (Dk). Finally, using Dk = 8, the GB^OBC1 model combined with interaction entropy (IE), which was calculated for entropic contribution (GB^OBC1_IE_8), was screened out for the binding free energy prediction and analysis of the Fis-DNA system, and the validity of the method was further verified by testing the Cren7-DNA system. By performing conformational analysis of the minor groove, it was found that mutation of the DNA central sequence A/T to C/G and deletion of the guanine 2-amino group would change the minor groove width and thus affect the formation of the major groove, altering the interaction and atomic contact between the protein and the major groove, thus changing the binding affinity of Fis and DNA. Hopefully, the series of tests in this work can shed some light on the related studies of protein and DNA systems.

The Hypersaline Archaeal Histones HpyA and HstA Are DNA Binding Proteins That Defy Categorization According to Commonly Used Functional Criteria

Histone proteins are found across diverse lineages of Archaea, many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum, HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii. hstA, like hpyA, is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species-from yeast to humans-possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea, the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others.

Genome-wide mapping of the Escherichia coli PhoB regulon reveals many transcriptionally inert, intragenic binding sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert, and hence are tolerated as genomic "noise".

IMPORTANCE

Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of non-regulatory, intragenic binding sites for transcription factors, and that these binding sites are not under selective pressure.

A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli

Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.

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, E_max, 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 E_max 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 E_max 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.

RfaH Counter-Silences Inhibition of Transcript Elongation by H-NS-StpA Nucleoprotein Filaments in Pathogenic Escherichia coli

Expression of virulence genes in pathogenic Escherichia coli is controlled in part by the transcription silencer H-NS and its paralogs (e.g., StpA), which sequester DNA in multi-kb nucleoprotein filaments to inhibit transcription initiation, elongation, or both. Some activators counter-silence initiation by displacing H-NS from promoters, but how H-NS inhibition of elongation is overcome is not understood. In uropathogenic E. coli (UPEC), elongation regulator RfaH aids expression of some H-NS-silenced pathogenicity operons (e.g., hlyCABD encoding hemolysin). RfaH associates with elongation complexes (ECs) via direct contacts to a transiently exposed, nontemplate DNA strand sequence called operon polarity suppressor (ops). RfaH-ops interactions establish long-lived RfaH-EC contacts that allow RfaH to recruit ribosomes to the nascent mRNA and to suppress transcriptional pausing and termination. Using ChIP-seq, we mapped the genome-scale distributions of RfaH, H-NS, StpA, RNA polymerase (RNAP), and σ⁷⁰ in the UPEC strain CFT073. We identify eight RfaH-activated operons, all of which were bound by H-NS and StpA. Four are new additions to the RfaH regulon. Deletion of RfaH caused premature termination, whereas deletion of H-NS and StpA allowed elongation without RfaH. Thus, RfaH is an elongation counter-silencer of H-NS. Consistent with elongation counter-silencing, deletion of StpA alone decreased the effect of RfaH. StpA increases DNA bridging, which inhibits transcript elongation via topological constraints on RNAP. Residual RfaH effect when both H-NS and StpA were deleted was attributable to targeting of RfaH-regulated operons by a minor H-NS paralog, Hfp. These operons have evolved higher levels of H-NS-binding features, explaining minor-paralog targeting. IMPORTANCE Bacterial pathogens adapt to hosts and host defenses by reprogramming gene expression, including by H-NS counter-silencing. Counter-silencing turns on transcription initiation when regulators bind to promoters and rearrange repressive H-NS nucleoprotein filaments that ordinarily block transcription. The specialized NusG paralog RfaH also reprograms virulence genes but regulates transcription elongation. To understand how elongation regulators might affect genes silenced by H-NS, we mapped H-NS, StpA (an H-NS paralog), RfaH, σ⁷⁰, and RNA polymerase (RNAP) locations on DNA in the uropathogenic E. coli strain CFT073. Although H-NS-StpA filaments bind only 18% of the CFT073 genome, all loci at which RfaH binds RNAP are also bound by H-NS-StpA and are silenced when RfaH is absent. Thus, RfaH represents a distinct class of counter-silencer that acts on elongating RNAP to enable transcription through repressive nucleoprotein filaments. Our findings define a new mechanism of elongation counter-silencing and explain how RfaH functions as a virulence regulator.

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.

Bacterial Transcription Factors Bind to Coding Regions and Regulate Internal Cryptic Promoters

Transcription factors (TFs) regulate transcription by binding to the specific sequences at the promoter region. However, the mechanisms and functions of TFs binding within the coding sequences (CDS) remain largely elusive in prokaryotes. To this end, we collected 409 data sets for bacterial TFs, including 104 chromatin immunoprecipitation sequencing (ChIP-seq) assays and 305 data sets from the systematic evolution of ligands by exponential enrichment (SELEX) in seven model bacteria. Interestingly, these TFs displayed the same binding capabilities for both coding and intergenic regions. Subsequent biochemical and genetic experiments demonstrated that several TFs bound to the coding regions and regulated the transcription of the binding or adjacent genes. Strand-specific RNA sequencing revealed that these CDS-binding TFs regulated the activity of the cryptic promoters, resulting in the altered transcription of the corresponding antisense RNA. TF RhpR hindered the transcriptional elongation of a subgenic transcript within a CDS. A ChIP-seq and Ribo-seq coanalysis revealed that RhpR influenced the translational efficiency of binding genes. Taken together, the present study reveals three regulatory mechanisms of CDS-bound TFs within individual genes, operons, and antisense RNAs, which demonstrate the variability of the regulatory mechanisms of TFs and expand upon the complexity of bacterial transcriptomes.

Interaction between transcribing RNA polymerase and topoisomerase I prevents R-loop formation in E. coli

Bacterial topoisomerase I (TopoI) removes excessive negative supercoiling and is thought to relax DNA molecules during transcription, replication and other processes. Using ChIP-Seq, we show that TopoI of Escherichia coli (EcTopoI) is colocalized, genome-wide, with transcribing RNA polymerase (RNAP). Treatment with transcription elongation inhibitor rifampicin leads to EcTopoI relocation to promoter regions, where RNAP also accumulates. When a 14 kDa RNAP-binding EcTopoI C-terminal domain (CTD) is overexpressed, colocalization of EcTopoI and RNAP along the transcription units is reduced. Pull-down experiments directly show that the two enzymes interact in vivo. Using ChIP-Seq and Topo-Seq, we demonstrate that EcTopoI is enriched upstream (within up to 12-15 kb) of highly-active transcription units, indicating that EcTopoI relaxes negative supercoiling generated by transcription. Uncoupling of the RNAP:EcTopoI interaction by either overexpression of EcTopoI competitor (CTD or inactive EcTopoI Y319F mutant) or deletion of EcTopoI domains involved in the interaction is toxic for cells and leads to excessive negative plasmid supercoiling. Moreover, uncoupling of the RNAP:EcTopoI interaction leads to R-loops accumulation genome-wide, indicating that this interaction is required for prevention of R-loops formation.

In E. coli, disruption of TopoI and RNAP interaction decreases cells viability and leads to hypernegative DNA supercoiling and R loops accumulation. TopoI and DNA gyrase bind around transcription units and TopoI recognizes cleavage sites by a specific motif and negative supercoiling.

Bacterial H-NS contacts DNA at the same irregularly spaced sites in both bridged and hemi-sequestered linear filaments

Gene silencing in bacteria is mediated by chromatin proteins, of which Escherichia coli H-NS is a paradigmatic example. H-NS forms nucleoprotein filaments with either one or two DNA duplexes. However, the structures, arrangements of DNA-binding domains (DBDs), and positions of DBD-DNA contacts in linear and bridged filaments are uncertain. To characterize the H-NS DBD contacts that silence transcription by RNA polymerase, we combined ·OH footprinting, molecular dynamics, statistical modeling, and DBD mapping using a chemical nuclease (Fe²⁺-EDTA) tethered to the DBDs (TEN-map). We find that H-NS DBDs contact DNA at indistinguishable locations in bridged or linear filaments and that the DBDs vary in orientation and position with ∼10-bp average spacing. Our results support a hemi-sequestration model of linear-to-bridged H-NS switching. Linear filaments able to inhibit only transcription initiation switch to bridged filaments able to inhibit both initiation and elongation using the same irregularly spaced DNA contacts.

Spatiotemporal Coupling of DNA Supercoiling and Genomic Sequence Organization-A Timing Chain for the Bacterial Growth Cycle?

In this article we describe the bacterial growth cycle as a closed, self-reproducing, or autopoietic circuit, reestablishing the physiological state of stationary cells initially inoculated in the growth medium. In batch culture, this process of self-reproduction is associated with the gradual decline in available metabolic energy and corresponding change in the physiological state of the population as a function of "travelled distance" along the autopoietic path. We argue that this directional alteration of cell physiology is both reflected in and supported by sequential gene expression along the chromosomal OriC-Ter axis. We propose that during the E. coli growth cycle, the spatiotemporal order of gene expression is established by coupling the temporal gradient of supercoiling energy to the spatial gradient of DNA thermodynamic stability along the chromosomal OriC-Ter axis.

Psoralen mapping reveals a bacterial genome supercoiling landscape dominated by transcription

DNA supercoiling is a key regulator of all DNA metabolic processes including replication, transcription, and recombination, yet a reliable genomic assay for supercoiling is lacking. Here, we present a robust and flexible method (Psora-seq) to measure whole-genome supercoiling at high resolution. Using this tool in Escherichia coli, we observe a supercoiling landscape that is well correlated to transcription. Supercoiling twin-domains generated by RNA polymerase complexes span 25 kb in each direction - an order of magnitude farther than previous measurements in any organism. Thus, ribosomal and many other highly expressed genes strongly affect the topology of about 40 neighboring genes each, creating highly integrated gene circuits. Genomic patterns of supercoiling revealed by Psora-seq could be aptly predicted from modeling based on gene expression levels alone, indicating that transcription is the major determinant of chromosome supercoiling. Large-scale supercoiling patterns were highly symmetrical between left and right chromosome arms (replichores), indicating that DNA replication also strongly influences supercoiling. Skew in the axis of symmetry from the natural ori-ter axis supports previous indications that the rightward replication fork is delayed several minutes after initiation. Implications of supercoiling on DNA replication and chromosome domain structure are discussed.

Regulatory role of CsuR (YiaU) in determination of cell surface properties of Escherichia coli K-12

Genomic SELEX screening was performed to identify the binding sites of YiaU, an uncharacterized LysR family transcription factor, on the Escherichia coli K-12 genome. Five high-affinity binding targets of YiaU were identified, all of which were involved in the structures of the bacterial cell surface such as outer and inner membrane proteins, and lipopolysaccharides. Detailed in vitro and in vivo analyses suggest that YiaU activates these target genes. To gain insight into the effects of YiaU in vivo on physiological properties, we used phenotype microarrays, biofilm screening assays and the sensitivity against serum complement analysed using a yiaU deletion mutant or YiaU expression strain. Together, these results suggest that the YiaU regulon confers resistance to some antibiotics, and increases biofilm formation and complement sensitivity. We propose renaming YiaU as CsuR (regulator of cell surface).

RfaH May Oppose Silencing by H-NS and YmoA Proteins during Transcription Elongation

Nucleoid-associated proteins (NAPs) silence xenogenes by blocking RNA polymerase binding to promoters and hindering transcript elongation. In Escherichia coli, H-NS and its homolog SptA interact with YmoA proteins Hha and YdgT to assemble nucleoprotein filaments that facilitate transcription termination by Rho, which acts in synergy with NusG. Countersilencing during initiation is facilitated by proteins that exclude NAPs from promoter regions, but auxiliary factors that alleviate silencing during elongation are not known. A specialized NusG paralog, RfaH, activates lipopolysaccharide core biosynthesis operons, enabling survival in the presence of detergents and antibiotics. RfaH strongly inhibits Rho-dependent termination by reducing RNA polymerase pausing, promoting translation, and competing with NusG. We hypothesize that RfaH also acts as a countersilencer of NAP/YmoA filaments. We show that deletions of hns and hha+ydgT suppress the growth defects of ΔrfaH by alleviating Rho-mediated polarity within the waa operon. The absence of YmoA proteins exacerbates cellular defects caused by reduced Rho levels or Rho inhibition by bicyclomycin but has negligible effects at a strong model Rho-dependent terminator. Our findings that the distribution of Hha and RfaH homologs is strongly correlated supports a model in which they comprise a silencing/countersilencing pair that controls expression of chromosomal and plasmid-encoded xenogenes. IMPORTANCE Horizontally acquired DNA drives bacterial evolution, but its unregulated expression may harm the recipient. Xenogeneic silencers recognize foreign genes and inhibit their transcription. However, some xenogenes, such as those encoding lipo- and exopolysaccharides, confer resistance to antibiotics, bile salts, and detergents, necessitating the existence of countersilencing fitness mechanisms. Here, we present evidence that Escherichia coli antiterminator RfaH alleviates silencing of the chromosomal waa operon and propose that plasmid-encoded RfaH homologs promote dissemination of antibiotic resistance genes through conjugation.

Xenogeneic silencing strategies in bacteria are dictated by RNA polymerase promiscuity

Horizontal gene transfer facilitates dissemination of favourable traits among bacteria. However, foreign DNA can also reduce host fitness: incoming sequences with a higher AT content than the host genome can misdirect transcription. Xenogeneic silencing proteins counteract this by modulating RNA polymerase binding. In this work, we compare xenogeneic silencing strategies of two distantly related model organisms: Escherichia coli and Bacillus subtilis. In E. coli, silencing is mediated by the H-NS protein that binds extensively across horizontally acquired genes. This prevents spurious non-coding transcription, mostly intragenic in origin. By contrast, binding of the B. subtilis Rok protein is more targeted and mostly silences expression of functional mRNAs. The difference reflects contrasting transcriptional promiscuity in E. coli and B. subtilis, largely attributable to housekeeping RNA polymerase σ factors. Thus, whilst RNA polymerase specificity is key to the xenogeneic silencing strategy of B. subtilis, transcriptional promiscuity must be overcome to silence horizontally acquired DNA in E. coli.

Bacteria use specific silencing proteins to prevent spurious transcription of horizontally acquired DNA. Here, Forrest et al. describe differences in silencing strategies between E. coli and Bacillus subtilis, driven by the respective specificities of the silencing protein and the RNA polymerase in each organism.

Facilitated Dissociation of Nucleoid Associated Proteins from DNA in the Bacterial Confinement

Transcription machinery depends on the temporal formation of protein-DNA complexes. Recent experiments demonstrated that not only the formation but also the lifetime of such complexes can affect the transcriptional machinery. In parallel, in vitro single-molecule studies showed that nucleoid-associated proteins (NAPs) leave the DNA rapidly as the bulk concentration of the protein increases via facilitated dissociation (FD). Nevertheless, whether such a concentration-dependent mechanism is functional in a bacterial cell, in which NAP levels and the 3d chromosomal structure are often coupled, is not clear a priori. Here, by using extensive coarse-grained molecular simulations, we model the unbinding of specific and nonspecific dimeric NAPs from a high-molecular-weight circular DNA molecule in a cylindrical structure mimicking the cellular confinement of a bacterial chromosome. Our simulations confirm that physiologically relevant peak protein levels (tens of micromolar) lead to highly compact chromosomal structures. This compaction results in rapid off rates (shorter DNA residence times) for specifically DNA-binding NAPs, such as the factor for inversion stimulation, which mostly dissociate via a segmental jump mechanism. Contrarily, for nonspecific NAPs, which are more prone to leave their binding sites via 1d sliding, the off rates decrease as the protein levels increase. The simulations with restrained chromosome models reveal that chromosome compaction is in favor of faster dissociation but only for specific proteins, and nonspecific proteins are not affected by the chromosome compaction. Overall, our results suggest that the cellular concentration level of a structural DNA-binding protein can be highly intermingled with its DNA residence time.

Global pleiotropic effects in adaptively evolved Escherichia coli lacking CRP reveal molecular mechanisms that define the growth physiology

Evolution facilitates emergence of fitter phenotypes by efficient allocation of cellular resources in conjunction with beneficial mutations. However, system-wide pleiotropic effects that redress the perturbations to the apex node of the transcriptional regulatory networks remain unclear. Here, we elucidate that absence of global transcriptional regulator CRP in Escherichia coli results in alterations in key metabolic pathways under glucose respiratory conditions, favouring stress- or hedging-related functions over growth-enhancing functions. Further, we disentangle the growth-mediated effects from the CRP regulation-specific effects on these metabolic pathways. We quantitatively illustrate that the loss of CRP perturbs proteome efficiency, as evident from metabolic as well as ribosomal proteome fractions, that corroborated with intracellular metabolite profiles. To address how E. coli copes with such systemic defect, we evolved Δcrp mutant in the presence of glucose. Besides acquiring mutations in the promoter of glucose transporter ptsG, the evolved populations recovered the metabolic pathways to their pre-perturbed state coupled with metabolite re-adjustments, which altogether enabled increased growth. By contrast to Δcrp mutant, the evolved strains remodelled their proteome efficiency towards biomass synthesis, albeit at the expense of carbon efficiency. Overall, we comprehensively illustrate the genetic and metabolic basis of pleiotropic effects, fundamental for understanding the growth physiology.

Distinct heterochromatin-like domains promote transcriptional memory and silence parasitic genetic elements in bacteria

There is increasing evidence that prokaryotes maintain chromosome structure, which in turn impacts gene expression. We recently characterized densely occupied, multi-kilobase regions in the E. coli genome that are transcriptionally silent, similar to eukaryotic heterochromatin. These extended protein occupancy domains (EPODs) span genomic regions containing genes encoding metabolic pathways as well as parasitic elements such as prophages. Here, we investigate the contributions of nucleoid-associated proteins (NAPs) to the structuring of these domains, by examining the impacts of deleting NAPs on EPODs genome-wide in E. coli and B. subtilis. We identify key NAPs contributing to the silencing of specific EPODs, whose deletion opens a chromosomal region for RNA polymerase binding at genes contained within that region. We show that changes in E. coli EPODs facilitate an extra layer of transcriptional regulation, which prepares cells for exposure to exotic carbon sources. Furthermore, we distinguish novel xenogeneic silencing roles for the NAPs Fis and Hfq, with the presence of at least one being essential for cell viability in the presence of domesticated prophages. Our findings reveal previously unrecognized mechanisms through which genomic architecture primes bacteria for changing metabolic environments and silences harmful genomic elements.

Genomic inversions in Escherichia coli alter gene expression and are associated with nucleoid protein binding sites

Genomic reorganization, like rearrangements and inversions, influence how genetic information is organized within bacterial genomes. Inversions in particular, facilitate genome evolution through gene gain and loss, and can alter gene expression. Previous studies investigating the impact inversions have on gene expression induced inversions targeting specific genes or examine inversions between distantly related species. This fails to encompass a genome wide perspective on naturally occurring inversions and their post adaptation impact on gene expression. Here we use bioinformatic techniques and multiple RNA-seq datasets to investigate the short- and long-range impact inversions have on genomic gene expression within <i>Escherichia coli</i>. We observed differences in gene expression between homologous inverted and non-inverted genes, even after long term exposure to adaptive selection. In 4% of inversions representing 33 genes, differential gene expression between inverted and non-inverted homologs was detected, with nearly two thirds (71%) of differentially expressed inverted genes having 9.4-85.6 fold higher gene expression. The identified inversions had more overlap than expected with nucleoid associated protein binding sites, which assist in genomic gene expression regulation. Some inversions can drastically impact gene expression even between different strains of <i>E.coli</i>, and could provide a mechanism for the diversification of genetic content through controlled expression changes.

Fis Connects Two Sensory Pathways, Quorum Sensing and Surface Sensing, to Control Motility in Vibrio parahaemolyticus

Factor for inversion stimulation (Fis) is a global regulator that is highly expressed during exponential phase growth and undetectable in stationary phase growth. Quorum sensing (QS) is a global regulatory mechanism that controls gene expression in response to changes in cell density and growth phase. In Vibrio parahaemolyticus, a marine species and a significant human pathogen, the QS regulatory sRNAs, Qrr1 to Qrr5, are expressed during exponential growth and negatively regulate the high cell density QS master regulator OpaR. OpaR is a positive regulator of capsule polysaccharide (CPS) formation, which is required for biofilm formation, and is a repressor of lateral flagella required for swarming motility. In V. parahaemolyticus, we show that Fis is a positive regulator of the qrr sRNAs expression. In an in-frame fis deletion mutant, qrr expression was repressed and opaR expression was induced. The Δfis mutant produced CPS and biofilm, but swarming motility was abolished. Also, the fis deletion mutant was more sensitive to polymyxin B. Swarming motility requires expression of both the surface sensing scrABC operon and lateral flagella laf operon. Our data showed that in the Δfis mutant both laf and scrABC genes were repressed. Fis controlled swarming motility indirectly through the QS pathway and directly through the surface sensing pathway. To determine the effects of Fis on cellular metabolism, we performed in vitro growth competition assays, and found that Δfis was outcompeted by wild type in minimal media supplemented with intestinal mucus as a sole nutrient source. The data showed that Fis positively modulated mucus components L-arabinose, D-gluconate and N-acetyl-D-glucosamine catabolism gene expression. In an in vivo colonization competition assay, Δfis was outcompeted by wild type, indicating Fis is required for fitness. Overall, these data demonstrate a global regulatory role for Fis in V. parahaemolyticus that includes QS, motility, and metabolism.

Testing mechanisms of DNA sliding by architectural DNA-binding proteins: dynamics of single wild-type and mutant protein molecules in vitro and in vivo

Architectural DNA-binding proteins (ADBPs) are abundant constituents of eukaryotic or bacterial chromosomes that bind DNA promiscuously and function in diverse DNA reactions. They generate large conformational changes in DNA upon binding yet can slide along DNA when searching for functional binding sites. Here we investigate the mechanism by which ADBPs diffuse on DNA by single-molecule analyses of mutant proteins rationally chosen to distinguish between rotation-coupled diffusion and DNA surface sliding after transient unbinding from the groove(s). The properties of yeast Nhp6A mutant proteins, combined with molecular dynamics simulations, suggest Nhp6A switches between two binding modes: a static state, in which the HMGB domain is bound within the minor groove with the DNA highly bent, and a mobile state, where the protein is traveling along the DNA surface by means of its flexible N-terminal basic arm. The behaviors of Fis mutants, a bacterial nucleoid-associated helix-turn-helix dimer, are best explained by mobile proteins unbinding from the major groove and diffusing along the DNA surface. Nhp6A, Fis, and bacterial HU are all near exclusively associated with the chromosome, as packaged within the bacterial nucleoid, and can be modeled by three diffusion modes where HU exhibits the fastest and Fis the slowest diffusion.

Multiscale Dynamic Structuring of Bacterial Chromosomes

Since the nucleoid was isolated from bacteria in the 1970s, two fundamental questions emerged and are still in the spotlight: how bacteria organize their chromosomes to fit inside the cell and how nucleoid organization enables essential biological processes. During the last decades, knowledge of bacterial chromosome organization has advanced considerably, and today, such chromosomes are considered to be highly organized and dynamic structures that are shaped by multiple factors in a multiscale manner. Here we review not only the classical well-known factors involved in chromosome organization but also novel components that have recently been shown to dynamically shape the 3D structuring of the bacterial genome. We focus on the different functional elements that control short-range organization and describe how they collaborate in the establishment of the higher-order folding and disposition of the chromosome. Recent advances have opened new avenues for a deeper understanding of the principles and mechanisms of chromosome organization in bacteria. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Hierarchy of transcription factor network in Escherichia coli K-12: H-NS-mediated silencing and Anti-silencing by global regulators

Transcriptional regulation for genome expression determines growth and adaptation of single-cell bacteria that are directly exposed to environment. The transcriptional apparatus in Escherichia coli K-12 is composed of RNA polymerase core enzyme and two groups of its regulatory proteins, seven species of promoter-recognition subunit sigma and about 300 species of transcription factors. The identification of regulatory targets for all these regulatory proteins is critical toward understanding the genome regulation as a whole. For this purpose, we performed a systematic search in vitro of the whole set of binding sites for each factor by gSELEX system. This review summarizes the accumulated knowledge of regulatory targets for more than 150 TFs from E. coli K-12. Overall TFs could be classified into four families: nucleoid-associated bifunctional TFs; global regulators; local regulators; and single-target regulators, in which the regulatory functions remain uncharacterized for the nucleoid-associated TFs. Here we overview the regulatory targets of two nucleoid-associated TFs, H-NS and its paralog StpA, both together playing the silencing role of a set of non-essential genes. Participation of LeuO and other global regulators have been indicated for the anti-silencing. Finally, we propose the hierarchy of TF network as a key framework of the bacterial genome regulation.

Single-Target Regulators Constitute the Minority Group of Transcription Factors in Escherichia coli K-12

The identification of regulatory targets of all transcription factors (TFs) is critical for understanding the entire network of genome regulation. A total of approximately 300 TFs exist in the model prokaryote Escherichia coli K-12, but the identification of whole sets of their direct targets is impossible with use of in vivo approaches. For this end, the most direct and quick approach is to identify the TF-binding sites in vitro on the genome. We then developed and utilized the gSELEX screening system in vitro for identification of more than 150 E. coli TF-binding sites along the E. coli genome. Based on the number of predicted regulatory targets, we classified E. coli K-12 TFs into four groups, altogether forming a hierarchy ranging from a single-target TF (ST-TF) to local TFs, global TFs, and nucleoid-associated TFs controlling as many as 1,000 targets. Using the collection of purified TFs and a library of genome DNA segments from a single and the same E. coli K-12, we identified here a total of 11 novel ST-TFs, CsqR, CusR, HprR, NorR, PepA, PutA, QseA, RspR, UvrY, ZraR, and YqhC. The regulation of single-target promoters was analyzed in details for the hitherto uncharacterized QseA and RspR. In most cases, the ST-TF gene and its regulatory target genes are adjacently located on the E. coli K-12 genome, implying their simultaneous transfer in the course of genome evolution. The newly identified 11 ST-TFs and the total of 13 hitherto identified altogether constitute the minority group of TFs in E. coli K-12.

Dynamic landscape of protein occupancy across the Escherichia coli chromosome

Free-living bacteria adapt to environmental change by reprogramming gene expression through precise interactions of hundreds of DNA-binding proteins. A predictive understanding of bacterial physiology requires us to globally monitor all such protein-DNA interactions across a range of environmental and genetic perturbations. Here, we show that such global observations are possible using an optimized version of in vivo protein occupancy display technology (in vivo protein occupancy display-high resolution, IPOD-HR) and present a pilot application to Escherichia coli. We observe that the E. coli protein-DNA interactome organizes into 2 distinct prototypic features: (1) highly dynamic condition-dependent transcription factor (TF) occupancy; and (2) robust kilobase scale occupancy by nucleoid factors, forming silencing domains analogous to eukaryotic heterochromatin. We show that occupancy dynamics across a range of conditions can rapidly reveal the global transcriptional regulatory organization of a bacterium. Beyond discovery of previously hidden regulatory logic, we show that these observations can be utilized to computationally determine sequence specificity models for the majority of active TFs. Our study demonstrates that global observations of protein occupancy combined with statistical inference can rapidly and systematically reveal the transcriptional regulatory and structural features of a bacterial genome. This capacity is particularly crucial for non-model bacteria that are not amenable to routine genetic manipulation.

Widespread divergent transcription from bacterial and archaeal promoters is a consequence of DNA-sequence symmetry

Transcription initiates at promoters, DNA regions recognized by a DNA-dependent RNA polymerase. We previously identified horizontally acquired Escherichia coli promoters from which the direction of transcription was unclear. In the present study, we show that more than half of these promoters are bidirectional and drive divergent transcription. Using genome-scale approaches, we demonstrate that 19% of all transcription start sites detected in E. coli are associated with a bidirectional promoter. Bidirectional promoters are similarly common in diverse bacteria and archaea, and have inherent symmetry: specific bases required for transcription initiation are reciprocally co-located on opposite DNA strands. Bidirectional promoters enable co-regulation of divergent genes and are enriched in both intergenic and horizontally acquired regions. Divergent transcription is conserved among bacteria, archaea and eukaryotes, but the underlying mechanisms for bidirectionality are different.

Prespacers formed during primed adaptation associate with the Cas1-Cas2 adaptation complex and the Cas3 interference nuclease-helicase

Primed adaptation allows rapid acquisition of protective spacers derived from foreign mobile genetic elements into CRISPR arrays of the host. Primed adaptation requires ongoing CRISPR interference that destroys foreign genetic elements, but the nature of this requirement is unknown. Using the Escherichia coli I-E CRISPR-Cas as a model, we show that prespacers, short fragments of foreign DNA on their way to become incorporated into CRISPR arrays as spacers, are associated with both the adaptation integrase Cas1 and the interference nuclease Cas3, implying physical association of the interference and adaptation machineries during priming.

For Type I CRISPR-Cas systems, a mode of CRISPR adaptation named priming has been described. Priming allows specific and highly efficient acquisition of new spacers from DNA recognized (primed) by the Cascade-crRNA (CRISPR RNA) effector complex. Recognition of the priming protospacer by Cascade-crRNA serves as a signal for engaging the Cas3 nuclease–helicase required for both interference and primed adaptation, suggesting the existence of a primed adaptation complex (PAC) containing the Cas1–Cas2 adaptation integrase and Cas3. To detect this complex in vivo, we here performed chromatin immunoprecipitation with Cas3-specific and Cas1-specific antibodies using cells undergoing primed adaptation. We found that prespacers are bound by both Cas1 (presumably, as part of the Cas1–Cas2 integrase) and Cas3, implying direct physical association of the interference and adaptation machineries as part of PAC.