Cells are complex systems composed of hundreds of genes whose products interact to produce elaborated behaviors. To control such behaviors, cells rely on transcription factors to regulate gene expression, and gene regulatory networks (GRNs) are employed to describe and understand such behavior. However, GRNs are static models, and dynamic models are difficult to obtain due to their size, complexity, stochastic dynamics and interactions with other cell processes.

Antibiotic-resistant bacteria are an increasingly serious public health concern, as strains emerge that demonstrate resistance to almost all available treatments. One factor that contributes to the crisis is the adaptive ability of bacteria, which exhibit remarkable phenotypic and gene expression heterogeneity in order to gain a survival advantage in damaging environments. This high degree of variability in gene expression across biological populations makes it a challenging task to identify key regulators of bacterial adaptation. Here, we research the regulation of adaptive resistance by investigating transcriptome profiles of Escherichia coli upon adaptation to disparate toxins, including antibiotics and biofuels. We locate potential target genes via conventional gene expression analysis as well as using a new analysis technique examining differential gene expression variability. By investigating trends across the diverse adaptation conditions, we identify a focused set of genes with conserved behavior, including those involved in cell motility, metabolism, membrane structure, and transport, and several genes of unknown function. To validate the biological relevance of the observed changes, we synthetically perturb gene expression using clustered regularly interspaced short palindromic repeat (CRISPR)-dCas9. Manipulation of select genes in combination with antibiotic treatment promotes adaptive resistance as demonstrated by an increased degree of antibiotic tolerance and heterogeneity in MICs. We study the mechanisms by which identified genes influence adaptation and find that select differentially variable genes have the potential to impact metabolic rates, mutation rates, and motility. Overall, this work provides evidence for a complex nongenetic response, encompassing shifts in gene expression and gene expression variability, which underlies adaptive resistance. IMPORTANCE Even initially sensitive bacteria can rapidly thwart antibiotic treatment through stress response processes known as adaptive resistance. Adaptive resistance fosters transient tolerance increases and the emergence of mutations conferring heritable drug resistance. In order to extend the applicable lifetime of new antibiotics, we must seek to hinder the occurrence of bacterial adaptive resistance; however, the regulation of adaptation is difficult to identify due to immense heterogeneity emerging during evolution. This study specifically seeks to generate heterogeneity by adapting bacteria to different stresses and then examines gene expression trends across the disparate populations in order to pinpoint key genes and pathways associated with adaptive resistance. The targets identified here may eventually inform strategies for impeding adaptive resistance and prolonging the effectiveness of antibiotic treatment.

Classic cadherins are cell adhesion molecules whose expression patterns are dynamically modulated in association with their diverse functions during morphogenesis. The large size and complexity of cadherin loci have made it a challenge to investigate the organization of cis-regulatory modules that control their spatiotemporal patterns of expression. Towards this end, we utilized bacterial artificial chromosomes (BACs) containing the Cdh6 gene, a mouse type II classic cadherin, to systematically identify cis-regulatory modules that govern its expression. By inserting a lacZ reporter gene into the Cdh6 BAC and generating a series of modified variants via homologous recombination or transposon insertions that have been examined in transgenic mice, we identified an array of genomic regions that contribute to specific regulation of the gene. These regions span approximately 350 kb of the locus between 161-kb upstream and 186-kb downstream of the Cdh6 transcription start site. Distinct modules independently regulate compartmental expression (i.e. forebrain, hindbrain rhombomeres, and spinal cord) and/or cell lineage-specific expression patterns (i.e. neural crest subpopulations such as Schwann cells) of Cdh6 at the early developmental stages. With respect to regulation of expression in neural crest cells, we have found that distinct regions contribute to different aspects of expression and have identified a short 79-bp region that is implicated in regulating expression in cells once they have emigrated from the neural tube. These results build a picture of the complex organization of Cdh6 cis-regulatory modules and highlight the diverse inputs that contribute to its dynamic expression during early mouse embryonic development.

Although it is often tacitly assumed that gene regulatory interactions are finely tuned, how accurate gene regulation could evolve from a state without regulation is unclear. Moreover, gene expression noise would seem to impede the evolution of accurate gene regulation, and previous investigations have provided circumstantial evidence that natural selection has acted to lower noise levels. By evolving synthetic Escherichia coli promoters de novo, we here show that, contrary to expectations, promoters exhibit low noise by default. Instead, selection must have acted to increase the noise levels of highly regulated E. coli promoters. We present a general theory of the interplay between gene expression noise and gene regulation that explains these observations. The theory shows that propagation of expression noise from regulators to their targets is not an unwanted side-effect of regulation, but rather acts as a rudimentary form of regulation that facilitates the evolution of more accurate regulation.

Pulmonary exposure to Francisella tularensis is associated with severe lung pathology and a high mortality rate. The lack of induction of classical inflammatory mediators, including IL1-β and TNF-α, during early infection has led to the suggestion that F. tularensis evades detection by host innate immune surveillance and/or actively suppresses inflammation. To gain more insight into the host response to Francisella infection during the acute stage, transcriptomic analysis was performed on lung tissue from mice exposed to virulent (Francisella tularensis ssp tularensis SchuS4). Despite an extensive transcriptional response in the lungs of animals as early as 4 hrs post-exposure, Francisella tularensis was associated with an almost complete lack of induction of immune-related genes during the initial 24 hrs post-exposure. This broad subversion of innate immune responses was particularly evident when compared to the pulmonary inflammatory response induced by other lethal (Yersinia pestis) and non-lethal (Legionella pneumophila, Pseudomonas aeruginosa) pulmonary infections. However, the unique induction of a subset of inflammation-related genes suggests a role for dysregulation of lymphocyte function and anti-inflammatory pathways in the extreme virulence of Francisella. Subsequent activation of a classical inflammatory response 48 hrs post-exposure was associated with altered abundance of Francisella-specific transcripts, including those associated with bacterial surface components. In summary, virulent Francisella induces a unique pulmonary inflammatory response characterized by temporal regulation of innate immune pathways correlating with altered bacterial gene expression patterns. This study represents the first simultaneous measurement of both host and Francisella transcriptome changes that occur during in vivo infection and identifies potential bacterial virulence factors responsible for regulation of host inflammatory pathways.

The Warburg effect, which originally described increased production of lactate in cancer, is associated with diverse cellular processes such as angiogenesis, hypoxia, polarization of macrophages and activation of T cells. This phenomenon is intimately linked to several diseases including neoplasia, sepsis and autoimmune diseases1,2. Lactate, which is converted from pyruvate in tumour cells, is widely known as an energy source and metabolic by-product. However, its non-metabolic functions in physiology and disease remain unknown. Here we show that lactate-derived lactylation of histone lysine residues serves as an epigenetic modification that directly stimulates gene transcription from chromatin. We identify 28 lactylation sites on core histones in human and mouse cells. Hypoxia and bacterial challenges induce the production of lactate by glycolysis, and this acts as a precursor that stimulates histone lactylation. Using M1 macrophages that have been exposed to bacteria as a model system, we show that histone lactylation has different temporal dynamics from acetylation. In the late phase of M1 macrophage polarization, increased histone lactylation induces homeostatic genes that are involved in wound healing, including Arg1. Collectively, our results suggest that an endogenous 'lactate clock' in bacterially challenged M1 macrophages turns on gene expression to promote homeostasis. Histone lactylation thus represents an opportunity to improve our understanding of the functions of lactate and its role in diverse pathophysiological conditions, including infection and cancer.

Riemerella anatipestifer is an important pathogen of waterfowl, which causes septicemia anserum exsudativa in ducks. In this study, an AS87_03730 gene deletion R. anatipestifer mutant Yb2ΔAS87_03730 was constructed to investigate the role of AS87_03730 on R. anatipestifer virulence and gene regulation. By deleting a 708-bp fragment from AS87_03730, the mutant Yb2ΔAS87_03730 showed a significant decreased growth rate in TSB and invasion capacity in Vero cells, compared to wild-type strain Yb2. Moreover, the median lethal dose (LD50) of Yb2ΔAS87_03730 was 1.24 × 10(7) colony forming units (CFU), which is about 80-fold attenuated than that of Yb2 (LD50 = 1.53 × 10(5) CFU). Furthermore, RNA-Seq analysis and Real-time PCR indicated 19 up-regulated and two down-regulated genes in Yb2ΔAS87_03730. Functional analysis revealed that 12 up-regulated genes were related to "Translation, ribosomal structure and biogenesis", two were classified into "Cell envelope biogenesis, outer membrane", one was involved in "Amino acid transport and metabolism", and the other four had unknown functions. Polymerase chain reaction and sequence analysis indicated that the AS87_03730 gene is highly conserved among R. anatipestifer strains, as the percent sequence identity was over 93.5%. This study presents evidence that AS87_03730 gene is involved in bacterial virulence and gene regulation of R. anatipestifer.

Bidirectional interactions between the immune and nervous systems are well established in vertebrates. Insects show similar neuro-immune-behavioral interactions to those seen in vertebrates. Using quantitative real-time PCR, we present evidence that gene expression in the honeybee head is influenced by activation of the immune system 8h after a bacterial challenge with Escherichia coli. Seven genes were selected for quantitative analysis in order to cover both typical functions of the head such as exocrine secretion (mrjp3 and mrjp4) and olfactory processes (obp17) as well as more general processes such as structural functions (mlc2 and paramyosin), stress response (ERp60) and energy housekeeping (enolase). In this way, we show at the molecular level that the immune system functions as a sensory organ in insects -- as it does in vertebrates -- which signals to the head that a bacterial infection is present, and leads to regulation of expression of several genes in the head by a yet unidentified mechanism.

Bacterial chromosomes are organized into polycistronic cotranscribed operons, but the evolutionary pressures maintaining them are unclear. We hypothesized that operons alter gene expression noise characteristics, resulting in selection for or against maintaining operons depending on network architecture. Mathematical models for 6 functional classes of network modules showed that three classes exhibited decreased noise and 3 exhibited increased noise with same-operon cotranscription of interacting proteins. Noise reduction was often associated with a decreased chance of reaching an ultrasensitive threshold. Stochastic simulations of the lac operon demonstrated that the predicted effects of transcriptional coupling hold for a complex network module. We employed bioinformatic analysis to find overrepresentation of noise-minimizing operon organization compared with randomized controls. Among constitutively expressed physically interacting protein pairs, higher coupling frequencies appeared at lower expression levels, where noise effects are expected to be dominant. Our results thereby suggest an important role for gene expression noise, in many cases interacting with an ultrasensitive switch, in maintaining or selecting for operons in bacterial chromosomes.

Chromosomal integration of recombinant genes is desirable compared with expression from plasmids due to increased stability, reduced cell-to-cell variability, and elimination of the need for antibiotics for plasmid maintenance. Here, we present a new approach for tuning pathway gene expression levels via random integration and high-throughput screening. We demonstrate multiplexed gene integration and expression-level optimization for isobutanol production in Escherichia coli The integrated strains could, with far lower expression levels than plasmid-based expression, produce high titers (10.0 ± 0.9 g/liter isobutanol in 48 hours) and yields (69% of the theoretical maximum). Close examination of pathway expression in the top-performing, as well as other isolates, reveals the complexity of cellular metabolism and regulation, underscoring the need for precise optimization while integrating pathway genes into the chromosome. We expect this method for pathway integration and optimization can be readily extended to a wide range of pathways and chassis to create robust and efficient production strains.

In diverse bacterial species, the global regulator Hfq contributes to post-transcriptional networks that control expression of numerous genes. Hfq of the opportunistic pathogen Pseudomonas aeruginosa inhibits translation of target transcripts by forming a regulatory complex with the catabolite repression protein Crc. This repressive complex acts as part of an intricate mechanism of preferred nutrient utilisation. We describe high-resolution cryo-EM structures of the assembly of Hfq and Crc bound to the translation initiation site of a target mRNA. The core of the assembly is formed through interactions of two cognate RNAs, two Hfq hexamers and a Crc pair. Additional Crc protomers are recruited to the core to generate higher-order assemblies with demonstrated regulatory activity in vivo. This study reveals how Hfq cooperates with a partner protein to regulate translation, and provides a structural basis for an RNA code that guides global regulators to interact cooperatively and regulate different RNA targets.

Type I interferons (IFNs) induce hundreds of IFN-stimulated genes (ISGs) in response to viral infection. Induction of these ISGs must be regulated for an efficient and controlled antiviral response, but post-transcriptional controls of these genes have not been well defined. Here, we identify a role for the RNA base modification N6-methyladenosine (m6A) in the regulation of ISGs. Using ribosome profiling and quantitative mass spectrometry, coupled with m6A-immunoprecipitation and sequencing, we identify a subset of ISGs, including IFITM1, whose translation is enhanced by m6A and the m6A methyltransferase proteins METTL3 and METTL14. We further determine that the m6A reader YTHDF1 increases the expression of IFITM1 in an m6A-binding-dependent manner. Importantly, we find that the m6A methyltransferase complex promotes the antiviral activity of type I IFN. Thus, these studies identify m6A as having a role in post-transcriptional control of ISG translation during the type I IFN response for antiviral restriction.

Stimulation of innate immunity by bacterial molecular patterns can induce an enhanced cellular immune response to pathogens that are associated with innate immune memory shaped by epigenetic changes. Immunological memory can be expressed in the acceleration/intensification of inflammation, as well as in the exact opposite-to maintain tolerance and non-response to a repeated stimulus. Tolerance is one of the central concepts of immunity and is ensured by the consistency of all parts of the immune response. The severe consequences of inflammation force researchers to study in detail all stages of the downstream pathways that are activated after exposure to a stimulus, while the formation of non-response to a pro-inflammatory stimulus has not yet received a detailed description. Elucidation of the mechanism of tolerance is an urgent task for the prevention and treatment of inflammatory diseases. The aim of this investigation was to study the dynamic changes in the gene expression of A20 and ATF3, the inflammation suppressors, against the background of the expression of the genes of the innate immunity receptors TLR4 and NOD2 and the pro-inflammatory cytokine TNF-α under the influence of TLR4 and NOD2 agonists, lipopolysaccharide (LPS) and glucosaminylmuramyl dipeptide (GMDP). The mechanism of inflammation regulation by bioregulators of bacterial origin-LPS and GMDP-was evaluated in vitro in human peripheral blood mononuclear cells and in vivo after i.p. administration of LPS and GMDP to mice. Gene expression was assessed by RT-PCR. Innate immune receptors and the pro-inflammatory cytokine TNF-α were found to develop early in response to LPS and GMDP, both in vitro and in vivo. Genes of cytosolic proteins controlling inflammation (A20 and ATF3) were expressed later. Prior exposure of the innate immune system to LPS and muramyl peptides may modulate host defense against acute inflammation. As a result of the study, new data were obtained on dynamic changes in deubiquitinase A20 and the transcription factor ATF3, which are involved in the limitation and suppression of inflammatory reactions caused by fragments of bacterial cell walls-LPS and GMDP. Thus, bioregulators of bacterial origin LPS and GMDP, along with pro-inflammatory factors, activate the expression of genes that suppress inflammation, which should be considered when analyzing data from studies of the pro-inflammatory properties of LPS and GMDP and when developing drugs based on them.

Inducible gene expression based upon Tet repressor (tet regulation) is a broadly applied tool in molecular genetics. In its original environment, Tet repressor (TetR) negatively controls tetracycline (tc) resistance in bacteria. In the presence of tc, TetR is induced and detaches from its cognate DNA sequence tetO, so that a tc antiporter protein is expressed. In this article, we provide a comprehensive overview about tet regulation in bacteria and illustrate the parameters of different regulatory architectures. While some of these set-ups rely on natural tet-control regions like those found on transposon Tn10, highly efficient variations of this system have recently been adapted to different Gram-negative and Gram-positive bacteria. Novel tet-controllable artificial or hybrid promoters were employed for target gene expression. They are controlled by regulators expressed at different levels either in a constitutive or in an autoregulated manner. The resulting tet systems have been used for various purposes. We discuss integrative elements vested with tc-sensitive promoters, as well as tet regulation in Gram-negative and Gram-positive bacteria for analytical purposes and for protein overproduction. Also the use of TetR as an in vivo biosensor for tetracyclines or as a regulatory device in synthetic biology constructs is outlined. Technical specifications underlying different regulatory set-ups are highlighted, and finally recent developments concerning variations of TetR are presented, which may expand the use of prokaryotic tet systems in the future.

The Escherichia coli ClpYQ (HslUV) complex is an ATP-dependent protease, and the clpQ⁺Y⁺ (hslV⁺U⁺) operon encodes two heat shock proteins, ClpQ and ClpY, respectively. The transcriptional (op) or translational (pr) clpQ⁺::lacZ fusion gene was constructed, with the clpQ⁺Y⁺ promoter fused to a lacZ reporter gene. The clpQ⁺::lacZ (op or pr) fusion gene was each crossed into lambda phage. The λlpQ⁺::lacZ⁺ (op), a transcriptional fusion gene, was used to form lysogens in the wild-type, rpoH or/and rpoS mutants. Upon shifting the temperature up from 30 ° C to 42 ° C, the wild-type λclpQ⁺::lacZ⁺ (op) demonstrates an increased β-galactosidase (βGal) activity. However, the βGal activity of clpQ⁺::lacZ⁺ (op) was decreased in the rpoH and rpoH rpoS mutants but not in the rpoS mutant. The levels of clpQ⁺::lacZ⁺ mRNA transcripts correlated well to their βGal activity. Similarly, the expression of the clpQ⁺::lacZ⁺ gene fusion was nearly identical to the clpQ⁺Y⁺ transcript under the in vivo condition. The clpQ(m1)::lacZ⁺, containing a point mutation in the -10 promoter region for RpoH binding, showed decreased βGal activity, independent of activation by RpoH. We conclude that RpoH itself regulates clpQ⁺Y⁺ gene expression. In addition, the clpQ⁺Y⁺ message carries a conserved 71 bp at the 5' untranslated region (5'UTR) that is predicted to form the stem-loop structure by analysis of its RNA secondary structure. The clpQ(m2)Δ40::lacZ⁺, with a 40 bp deletion in the 5'UTR, showed a decreased βGal activity. In addition, from our results, it is suggested that this stem-loop structure is necessary for the stability of the clpQ⁺Y⁺ message.

Riemerella anatipestifer is an important pathogen that causes septicemia anserum exsudativa in ducks. Lipopolysaccharide (LPS) is considered to be a major virulence factor of R. anatipestifer. To identify genes involved in LPS biosynthesis, we screened a library of random Tn4351 transposon mutants using a monoclonal antibody against R. anatipestifer serotype 1 LPS (anti-LPS MAb). A mutant strain RA1067 which lost the reactivity in an indirect ELISA was obtained. Southern blot and sequencing analyses indicated a single Tn4351 was inserted at 116 bp in the M949_RS01915 gene in the RA1067 chromosomal DNA. Silver staining and Western blot analyses indicated that the RA1067 LPS was defected compared to the wild-type strain CH3 LPS. The RA1067 displayed a significant decreased growth rate at the late stage of growth in TSB in comparison with CH3. In addition, RA1067 showed higher susceptibility to complement-dependent killing, more than 360-fold attenuated virulence based on the median lethal dose determination, increased bacterial adhesion and invasion capacities to Vero cells and significantly decreased blood bacterial loads in RA1067 infected ducks, when compared to the CH3. An animal experiment indicated that inactivated RA1067 cells was effective in cross-protecting of the ducks from challenging with R. anatipestifer strains WJ4 (serotype 1), Yb2 (serotype 2) and HXb2 (serotype 10), further confirming the alteration of the RA1067 antigenicity. Moreover, RNA-Seq analysis and real-time PCR verified two up-regulated and three down-regulated genes in RA1067. Our findings demonstrate that the M949_RS01915 gene is associated to bacterial antigenicity, pathogenicity and gene regulation of R. anatipestifer.

The bean bug Riptortus pedestris possesses a specialized symbiotic organ in a posterior region of the midgut, where numerous crypts harbor extracellular betaproteobacterial symbionts of the genus Burkholderia. Second instar nymphs orally acquire the symbiont from the environment, and the symbiont infection benefits the host by facilitating growth and by occasionally conferring insecticide resistance. Here we performed comparative transcriptomic analyses of insect genes expressed in symbiotic and non-symbiotic regions of the midgut dissected from Burkholderia-infected and uninfected R. pedestris. Expression sequence tag analysis of cDNA libraries and quantitative reverse transcription PCR identified a number of insect genes expressed in symbiosis- or aposymbiosis-associated patterns. For example, genes up-regulated in symbiotic relative to aposymbiotic individuals, including many cysteine-rich secreted protein genes and many cathepsin protease genes, are likely to play a role in regulating the symbiosis. Conversely, genes up-regulated in aposymbiotic relative to symbiotic individuals, including a chicken-type lysozyme gene and a defensin-like protein gene, are possibly involved in regulation of non-symbiotic bacterial infections. Our study presents the first transcriptomic data on gut symbiotic organ of a stinkbug, which provides initial clues to understanding of molecular mechanisms underlying the insect-bacterium gut symbiosis and sheds light on several intriguing commonalities between endocellular and extracellular symbiotic associations.

Understanding the mechanisms of lipid accumulation in microorganisms is important for several reasons. In addition to providing insight into assembly of biological membranes, lipid accumulation has important applications in the production of renewable fuels and chemicals. The photosynthetic bacterium Rhodobacter sphaeroides is an attractive organism to study lipid accumulation, as it has the ability to increase membrane production at low O2 tensions. Under these conditions, R. sphaeroides develops invaginations of the cytoplasmic membrane to increase its membrane surface area for housing of the membrane-bound components of its photosynthetic apparatus. Here we use fatty acid levels as a reporter of membrane lipid content. We show that, under low-O2 and anaerobic conditions, the total fatty acid content per cell increases 3-fold. We also find that the increases in the amount of fatty acid and photosynthetic pigment per cell are correlated as O2 tensions or light intensity are changed. To ask if lipid and pigment accumulation were genetically separable, we analyzed strains with mutations in known photosynthetic regulatory pathways. While a strain lacking AppA failed to induce photosynthetic pigment-protein complex accumulation, it increased fatty acid content under low-O2 conditions. We also found that an intact PrrBA pathway is required for low-O2-induced fatty acid accumulation. Our findings suggest a previously unknown role of R. sphaeroides transcriptional regulators in increasing fatty acid and phospholipid accumulation in response to decreased O2 tension.

Abundance and diversity of bacteria and their viral predators, bacteriophages (phages), in the digestive tract are associated with human health. Particularly intriguing is the long-term coexistence of these two antagonistic populations. We performed genome-wide RNA sequencing on a human enteroaggregative Escherichia coli isolate to identify genes differentially expressed between in vitro conditions and in murine intestines. We experimentally demonstrated that four of these differentially expressed genes modified the interactions between E. coli and three virulent phages by either increasing or decreasing its susceptibility/resistance pattern and also by interfering with biofilm formation. Therefore, the regulation of bacterial genes expression during the colonization of the digestive tract influences the coexistence of phages and bacteria, highlighting the intricacy of tripartite relationships between phages, bacteria, and the animal host in intestinal homeostasis.

Understanding gene function and regulation is essential for the interpretation, prediction, and ultimate design of cell responses to changes in the environment. An important step toward meeting the challenge of understanding gene function and regulation is the identification of sets of genes that are always co-expressed. These gene sets, Atomic Regulons (ARs), represent fundamental units of function within a cell and could be used to associate genes of unknown function with cellular processes and to enable rational genetic engineering of cellular systems. Here, we describe an approach for inferring ARs that leverages large-scale expression data sets, gene context, and functional relationships among genes. We computed ARs for Escherichia coli based on 907 gene expression experiments and compared our results with gene clusters produced by two prevalent data-driven methods: Hierarchical clustering and k-means clustering. We compared ARs and purely data-driven gene clusters to the curated set of regulatory interactions for E. coli found in RegulonDB, showing that ARs are more consistent with gold standard regulons than are data-driven gene clusters. We further examined the consistency of ARs and data-driven gene clusters in the context of gene interactions predicted by Context Likelihood of Relatedness (CLR) analysis, finding that the ARs show better agreement with CLR predicted interactions. We determined the impact of increasing amounts of expression data on AR construction and find that while more data improve ARs, it is not necessary to use the full set of gene expression experiments available for E. coli to produce high quality ARs. In order to explore the conservation of co-regulated gene sets across different organisms, we computed ARs for Shewanella oneidensis, Pseudomonas aeruginosa, Thermus thermophilus, and Staphylococcus aureus, each of which represents increasing degrees of phylogenetic distance from E. coli. Comparison of the organism-specific ARs showed that the consistency of AR gene membership correlates with phylogenetic distance, but there is clear variability in the regulatory networks of closely related organisms. As large scale expression data sets become increasingly common for model and non-model organisms, comparative analyses of atomic regulons will provide valuable insights into fundamental regulatory modules used across the bacterial domain.