Transcriptional regulatory networks in bacteria: from input signals to output responses

Regulation of TCA cycle genes by srbA sRNA: Impacts on Pseudomonas aeruginosa virulence and survival

Pseudomonas aeruginosa, an opportunistic bacterial pathogen of public health concern, is known for its metabolic versatility, adaptability in harsh environment, and pathogenic aggressiveness. P. aeruginosa relies on various regulatory networks modulated by small non-coding RNAs, which in turn influence different physiological traits such as metabolism, stress response, and pathogenesis. In this study, srbA sRNA has been shown to play a diverse role in regulating cellular metabolism and the production of different virulence factors in P. aeruginosa. srbA was found to control the TCA cycle, a key regulatory pathway for cellular metabolism and energy production, by regulating three main enzymes: citrate synthase (gltA), isocitrate dehydrogenase (icd), and α-ketoglutarate dehydrogenase E1 subunit (sucA) at both the transcriptional and translational levels. By modulating the TCA cycle, srbA could help the bacteria to adapt nutritional stress by lowering energy consumption. Additionally, srbA has been found to differentially regulate production of various virulence factors such as rhamnolipid, elastase, LasA protease, and pyocyanin under both nutrient-rich and nutrient-limiting conditions. It could also influence motilities in P. aeruginosa, linked to biofilm formation and pathogenicity. Thus, srbA might hold a promise in the research area for identifying virulence pathways and developing novel therapeutic targets to combat the global pathogenic threat of P. aeruginosa.

Exploring the Genome-wide Expression Level of the Bacterial Strain Belonging to Bacillus safensis (MM19) Against Phomopsis viticola

Rhizobacteria have the ability to compete with pathogenic microorganisms and contribute to plant immunity and defense mechanisms. Their growth and survival in the rhizosphere ensure a biological balance in favor of the host plant. The differential gene expression profiles of B. safensis (MM19) revealed significantly increased expression of prominent genes related to thiamine biosynthesis involving various metabolites and enzymes that participate in the suppression of mycelium growth and pathogen inhibition. Correspondingly, the expression of three major genes (HOG1, FUS3, SGI) involved in the virulence of P. viticola was assessed using qPCR analysis. HOG1 was the highest expressed gene in the pathogen when it was co-cultivated with MM19. Based on these findings, we performed molecular docking and dynamics analysis to explore the interaction between HOG1 and thiamine, as well as expression network analysis constructed using Cytoscape. The functional genomic data related to thiamine biosynthesis and the corresponding pathways ensure a priming role in the antagonistic behavior of B. safensis (MM19) against P. viticola as a support for plant immunity.

Functional and structural characterization of an IclR family transcription factor for the development of dicarboxylic acid biosensors

Prokaryotic transcription factors (TFs) regulate gene expression in response to small molecules, thus representing promising candidates as versatile small molecule-detecting biosensors valuable for synthetic biology applications. The engineering of such biosensors requires thorough in vitro and in vivo characterization of TF ligand response as well as detailed molecular structure information. In this work, we functionally and structurally characterize the Pca regulon regulatory protein (PcaR) transcription factor belonging to the IclR transcription factor family. Here, we present in vitro functional analysis of the ligand profile of PcaR and the construction of genetic circuits for the characterization of PcaR as an in vivo biosensor in the model eukaryote Saccharomyces cerevisiae. We report the crystal structures of PcaR in the apo state and in complex with one of its ligands, succinate, which suggests the mechanism of dicarboxylic acid recognition by this transcription factor. This work contributes key structural and functional insights enabling the engineering of PcaR for dicarboxylic acid biosensors, in addition to providing more insights into the IclR family of regulators.

Constitutive Activation of RpoH and the Addition of L-arabinose Influence Antibiotic Sensitivity of PHL628 E. coli

Antibiotics are used to combat the ever-present threat of infectious diseases, but bacteria are continually evolving an assortment of defenses that enable their survival against even the most potent treatments. While the demand for novel antibiotic agents is high, the discovery of a new agent is exceedingly rare. We chose to focus on understanding how different signal transduction pathways in the gram-negative bacterium Escherichia coli (E. coli) influence the sensitivity of the organism to antibiotics from three different classes: tetracycline, chloramphenicol, and levofloxacin. Using the PHL628 strain of E. coli, we exogenously overexpressed two transcription factors, FliA and RpoH.I54N (a constitutively active mutant), to determine their influence on the minimum inhibitory concentration (MIC) and minimum duration of killing (MDK) concentration for each of the studied antibiotics. We hypothesized that activating these pathways, which upregulate genes that respond to specific stressors, could mitigate bacterial response to antibiotic treatment. We also compared the exogenous overexpression of the constitutively active RpoH mutant to thermal heat shock that has feedback loops maintained. While FliA overexpression had no impact on MIC or antibiotic tolerance, RpoH.I54N overexpression reduced the MIC for tetracycline and chloramphenicol but had no independent impact on antibiotic tolerance. Thermal heat shock alone also did not affect MIC or antibiotic tolerance. L-arabinose, the small molecule used to induce expression in our system, unexpectedly independently increased the MICs for tetracycline (>2-fold) and levofloxacin (3-fold). Additionally, the combination of thermal heat shock and arabinose provided a synergistic, 5-fold increase in MIC for chloramphenicol. Arabinose increased the tolerance, as assessed by MDK99, for chloramphenicol (2-fold) and levofloxacin (4-fold). These experiments highlight the potential of the RpoH pathway to modulate antibiotic sensitivity and the emerging implication of arabinose in enhanced MIC and antibiotic tolerance.

Independent component analysis of Corynebacterium glutamicum transcriptomes reveals its transcriptional regulatory network

Gene expression in bacteria is regulated by multiple transcription factors. Clarifying the regulation mechanism of gene expression is necessary to understand bacterial physiological activities. To further understand the structure of the transcriptional regulatory network of Corynebacterium glutamicum, we applied independent component analysis, an unsupervised machine learning algorithm, to the high-quality C. glutamicum gene expression profile which includes 263 samples from 29 independent projects. We obtained 87 robust independent regulatory modules (iModulons). These iModulons explain 76.7% of the variance in the expression profile and constitute the quantitative transcriptional regulatory network of C. glutamicum. By analyzing the constituent genes in iModulons, we identified potential targets for 20 transcription factors. We also captured the changes in iModulon activities under different growth rates and dissolved oxygen concentrations, demonstrating the ability of iModulons to comprehensively interpret transcriptional responses to environmental changes. In summary, this study provides a genome-scale quantitative transcriptional regulatory network for C. glutamicum and informs future research on complex changes in the transcriptome.

Divergent directed evolution of a TetR-type repressor towards aromatic molecules

Reprogramming cellular behaviour is one of the hallmarks of synthetic biology. To this end, prokaryotic allosteric transcription factors (aTF) have been repurposed as versatile tools for processing small molecule signals into cellular responses. Expanding the toolbox of aTFs that recognize new inducer molecules is of considerable interest in many applications. Here, we first establish a resorcinol responsive aTF-based biosensor in Escherichia coli using the TetR-family repressor RolR from Corynebacterium glutamicum. We then perform an iterative walk along the fitness landscape of RolR to identify new inducer specificities, namely catechol, methyl catechol, caffeic acid, protocatechuate, L-DOPA, and the tumour biomarker homovanillic acid. Finally, we demonstrate the versatility of these engineered aTFs by transplanting them into the model eukaryote Saccharomyces cerevisiae. This work provides a framework for efficient aTF engineering to expand ligand specificity towards novel molecules on laboratory timescales, which, more broadly, is invaluable across a wide range of applications such as protein and metabolic engineering, as well as point-of-care diagnostics.

Transcription factor expression levels and environmental signals constrain transcription factor innovation

Evolutionary innovation of transcription factors frequently drives phenotypic diversification and adaptation to environmental change. Transcription factors can gain or lose connections to target genes, resulting in novel regulatory responses and phenotypes. However the frequency of functional adaptation varies between different regulators, even when they are closely related. To identify factors influencing propensity for innovation, we utilise a Pseudomonas fluorescens SBW25 strain rendered incapable of flagellar mediated motility in soft-agar plates via deletion of the flagellar master regulator (fleQ ). This bacterium can evolve to rescue flagellar motility via gene regulatory network rewiring of an alternative transcription factor to rescue activity of FleQ. Previously, we have identified two members (out of 22) of the RpoN-dependent enhancer binding protein (RpoN-EBP) family of transcription factors (NtrC and PFLU1132) that are capable of innovating in this way. These two transcription factors rescue motility repeatably and reliably in a strict hierarchy – with NtrC the only route in a ∆fleQ background, and PFLU1132 the only route in a ∆fleQ ∆ntrC background. However, why other members in the same transcription factor family have not been observed to rescue flagellar activity is unclear. Previous work shows that protein homology cannot explain this pattern within the protein family (RpoN-EBPs), and mutations in strains that rescued motility suggested high levels of transcription factor expression and activation drive innovation. We predict that mutations that increase expression of the transcription factor are vital to unlock evolutionary potential for innovation. Here, we construct titratable expression mutant lines for 11 of the RpoN-EBPs in P. fluorescens . We show that in five additional RpoN-EBPs (FleR, HbcR, GcsR, DctD, AauR and PFLU2209), high expression levels result in different mutations conferring motility rescue, suggesting alternative rewiring pathways. Our results indicate that expression levels (and not protein homology) of RpoN-EBPs are a key constraining factor in determining evolutionary potential for innovation. This suggests that transcription factors that can achieve high expression through few mutational changes, or transcription factors that are active in the selective environment, are more likely to innovate and contribute to adaptive gene regulatory network evolution.

On the evolution of natural product biosynthesis

Natural products are the raw material for drug discovery programmes. Bioactive natural products are used extensively in medicine and agriculture and have found utility as antibiotics, immunosuppressives, anti-cancer drugs and anthelminthics. Remarkably, the natural role and what mechanisms drive evolution of these molecules is relatively poorly understood. The exponential increase in genome and chemical data in recent years, coupled with technical advances in bioinformatics and genetics have enabled progress to be made in understanding the evolution of biosynthetic gene clusters and the products of their enzymatic machinery. Here we discuss the diversity of natural products, incorporating the mechanisms that govern evolution of metabolic pathways and how this can be applied to biosynthetic gene clusters. We build on the nomenclature of natural products in terms of primary, integrated, secondary and specialised metabolism and place this within an ecology-evolutionary-developmental biology framework. This eco-evo-devo framework we believe will help to clarify the nature and use of the term specialised metabolites in the future.

Yeast Mannan-Rich Fraction Modulates Endogenous Reactive Oxygen Species Generation and Antibiotic Sensitivity in Resistant E. coli

Mannan-rich fraction (MRF) isolated from Saccharomyces cerevisiae has been studied for its beneficial impact on animal intestinal health. Herein, we examined how MRF affected the formation of reactive oxygen species (ROS), impacting antibiotic susceptibility in resistant Escherichia coli through the modulation of bacterial metabolism. The role of MRF in effecting proteomic change was examined using a proteomics-based approach. The results showed that MRF, when combined with bactericidal antibiotic treatment, increased ROS production in resistant E. coli by 59.29 ± 4.03% compared to the control (p ≤ 0.05). We further examined the effect of MRF alone and in combination with antibiotic treatment on E. coli growth and explored how MRF potentiates bacterial susceptibility to antibiotics via proteomic changes in key metabolic pathways. Herein we demonstrated that MRF supplementation in the growth media of ampicillin-resistant E. coli had a significant impact on the normal translational control of the central metabolic pathways, including those involved in the glycolysis-TCA cycle (p ≤ 0.05).

Navigating Environmental Transitions: the Role of Phenotypic Variation in Bacterial Responses

The ability of bacteria to respond to changes in their environment is critical to their survival, allowing them to withstand stress, form complex communities, and induce virulence responses during host infection. A remarkable feature of many of these bacterial responses is that they are often variable across individual cells, despite occurring in an isogenic population exposed to a homogeneous environmental change, a phenomenon known as phenotypic heterogeneity. Phenotypic heterogeneity can enable bet-hedging or division of labor strategies that allow bacteria to survive fluctuating conditions. Investigating the significance of phenotypic heterogeneity in environmental transitions requires dynamic, single-cell data. Technical advances in quantitative single-cell measurements, imaging, and microfluidics have led to a surge of publications on this topic. Here, we review recent discoveries on single-cell bacterial responses to environmental transitions of various origins and complexities, from simple diauxic shifts to community behaviors in biofilm formation to virulence regulation during infection. We describe how these studies firmly establish that this form of heterogeneity is prevalent and a conserved mechanism by which bacteria cope with fluctuating conditions. We end with an outline of current challenges and future directions for the field. While it remains challenging to predict how an individual bacterium will respond to a given environmental input, we anticipate that capturing the dynamics of the process will begin to resolve this and facilitate rational perturbation of environmental responses for therapeutic and bioengineering purposes.

The Spatial Organization of Bacterial Transcriptional Regulatory Networks

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

Microbial interactions play an important role in regulating the effects of plant species on soil bacterial diversity

Plant species and microbial interactions have significant impacts on the diversity of bacterial communities. However, few studies have explored interactions among these factors, such the role of microbial interactions in regulating the effects of plant species on soil bacterial diversity. We assumed that plant species not only affect bacterial community diversity directly, but also influence bacterial community diversity indirectly through changing microbial interactions. Specifically, we collected soil samples associated with three different plant species, one evergreen shrub (Rhododendron simsii) and the other two deciduous shrubs (Dasiphora fruticosa and Salix oritrepha). Soil bacterial community composition and diversity were examined by high-throughput sequencing. Moreover, soil bacterial antagonistic interactions and soil edaphic characteristics were evaluated. We used structural equation modeling (SEM) to disentangle and compare the direct effect of different plant species on soil bacterial community diversity, and their indirect effects through influence on soil edaphic characteristics and microbial antagonistic interactions. The results showed that (1) Plant species effects on soil bacterial diversity were significant; (2) Plant species effects on soil microbial antagonistic interactions were significant; and (3) there was not only a significant direct plant species effect on bacterial diversity, but also a significant indirect effect on bacterial diversity through influence on microbial antagonistic interactions. Our study reveals the difference among plant species in their effects on soil microbial antagonistic interactions and highlights the vital role of microbial interactions on shaping soil microbial community diversity.

Targeting the Essential Transcription Factor HP1043 of Helicobacter pylori: A Drug Repositioning Study

Antibiotic-resistant bacterial pathogens are a very challenging problem nowadays. Helicobacter pylori is one of the most widespread and successful human pathogens since it colonizes half of the world population causing chronic and atrophic gastritis, peptic ulcer, mucosa-associated lymphoid tissue-lymphoma, and even gastric adenocarcinoma. Moreover, it displays resistance to numerous antibiotics. One of the H. pylori pivotal transcription factors, HP1043, plays a fundamental role in regulating essential cellular processes. Like other bacterial transcription factors, HP1043 does not display a eukaryote homolog. These characteristics make HP1043 a promising candidate to develop novel antibacterial strategies. Drug repositioning is a relatively recent strategy employed in drug development; testing approved drugs on new targets considerably reduces the time and cost of this process. The combined computational and in vitro approach further reduces the number of compounds to be tested in vivo. Our aim was to identify a subset of known drugs able to prevent HP1043 binding to DNA promoters. This result was reached through evaluation by molecular docking the binding capacity of about 14,350 molecules on the HP1043 dimer in both conformations, bound and unbound to the DNA. Employing an ad hoc pipeline including MMGBSA molecular dynamics, a selection of seven drugs was obtained. These were tested in vitro by electrophoretic mobility shift assay to evaluate the HP1043–DNA interaction. Among these, three returned promising results showing an appreciable reduction of the DNA-binding activity of HP1043. Overall, we applied a computational methodology coupled with experimental validation of the results to screen a large number of known drugs on one of the H. pylori essential transcription factors. This methodology allowed a rapid reduction of the number of drugs to be tested, and the drug repositioning approach considerably reduced the drug design costs. Identified drugs do not belong to the same pharmaceutical category and, by computational analysis, bound different cavities, but all display a reduction of HP1043 binding activity on the DNA.

A Versatile Transcription Factor Biosensor System Responsive to Multiple Aromatic and Indole Inducers

Allosteric transcription factor (aTF) biosensors are valuable tools for engineering microbes toward a multitude of applications in metabolic engineering, biotechnology, and synthetic biology. One of the challenges toward constructing functional and diverse biosensors in engineered microbes is the limited toolbox of identified and characterized aTFs. To overcome this, extensive bioprospecting of aTFs from sequencing databases, as well as aTF ligand-specificity engineering are essential in order to realize their full potential as biosensors for novel applications. In this work, using the TetR-family repressor CmeR from Campylobacter jejuni, we construct aTF genetic circuits that function as salicylate biosensors in the model organisms Escherichia coli and Saccharomyces cerevisiae. In addition to salicylate, we demonstrate the responsiveness of CmeR-regulated promoters to multiple aromatic and indole inducers. This relaxed ligand specificity of CmeR makes it a useful tool for detecting molecules in many metabolic engineering applications, as well as a good target for directed evolution to engineer proteins that are able to detect new and diverse chemistries.

Machine learning from Pseudomonas aeruginosa transcriptomes identifies independently modulated sets of genes associated with known transcriptional regulators

The transcriptional regulatory network (TRN) of Pseudomonas aeruginosa coordinates cellular processes in response to stimuli. We used 364 transcriptomes (281 publicly available + 83 in-house generated) to reconstruct the TRN of P. aeruginosa using independent component analysis. We identified 104 independently modulated sets of genes (iModulons) among which 81 reflect the effects of known transcriptional regulators. We identified iModulons that (i) play an important role in defining the genomic boundaries of biosynthetic gene clusters (BGCs), (ii) show increased expression of the BGCs and associated secretion systems in nutrient conditions that are important in cystic fibrosis, (iii) show the presence of a novel ribosomally synthesized and post-translationally modified peptide (RiPP) BGC which might have a role in P. aeruginosa virulence, (iv) exhibit interplay of amino acid metabolism regulation and central metabolism across different carbon sources and (v) clustered according to their activity changes to define iron and sulfur stimulons. Finally, we compared the identified iModulons of P. aeruginosa with those previously described in Escherichia coli to observe conserved regulons across two Gram-negative species. This comprehensive TRN framework encompasses the majority of the transcriptional regulatory machinery in P. aeruginosa, and thus should prove foundational for future research into its physiological functions.

Interaction with mammalian enteric viruses alters outer membrane vesicle production and content by commensal bacteria

Intestinal commensal bacteria contribute to maintaining gut homeostasis. Disruptions to the commensal flora are linked to the development and persistence of disease. The importance of these organisms is further demonstrated by the widespread ability of enteric viruses to exploit commensal bacteria to enhance viral infection. These viruses interact directly with commensal bacteria, and while the impact of this interaction on viral infection is well described for several viruses, the impact on the commensal bacteria has yet to be explored. In this article, we demonstrate, for the first time, that enteric viruses alter the gene expression and phenotype of individual commensal bacteria. Human and murine norovirus interaction with bacteria resulted in genome-wide differential gene expression and marked changes in the surface architecture of the bacterial cells. Furthermore, the interaction of the virus with bacteria led to increased production of smaller outer membrane vesicles (OMVs). Enhanced production of smaller vesicles was also observed when noroviruses were incubated with other commensal bacteria, indicating a potentially broad impact of norovirus interaction. The vesicle production observed in the in vivo model followed a similar trend where an increased quantity of smaller bacterial vesicles was observed in stool collected from virus-infected mice compared to mock-infected mice. Furthermore, changes in vesicle size were linked to changes in protein content and abundance, indicating that viral binding induced a shift in the mechanism of the OMV biogenesis. Collectively, these data demonstrate that enteric viruses induce specific changes in bacterial gene expression, leading to changes in bacterial extracellular vesicle production that can potentially impact host responses to infection.

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

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

Expression of the ace operon in Escherichia coli is triggered in response to growth rate-dependent flux-signal of ATP

The signal that triggers the expression of the ace operon and, in turn, the transition of central metabolism's architecture from acetogenic to gluconeogenic in Escherichia coli remains elusive despite extensive research both in vivo and in vitro. Here, with the aid of flux analysis together with measurements of the enzymic activity of isocitrate lyase (ICL) and its aceA-messenger ribonucleuc acid (mRNA) transcripts, we provide credible evidence suggesting that the expression of the ace operon in E. coli is triggered in response to growth rate-dependent threshold flux-signal of adenosine triphosphate (ATP). Flux analysis revealed that the shortfall in ATP supply observed as the growth rate ($\mu $) diminishes from µmax to ≤ 0.43h-1 ($ \pm 0.02;n4)\ $is partially redressed by up-regulating flux through succinyl CoA synthetase. Unlike glycerol and glucose, pyruvate cannot feed directly into the two glycolytic ATP-generating reactions catalyzed by phosphoglycerokinase and pyruvate kinase. On the other hand, glycerol, which upon its conversion to D-glyceraldehyde, feeds into the phosphorylation and dephosphorylation parts of glycolysis including the substrate-level phosphorylation-ATP generating reactions, thus preventing ATP flux from dropping to the critical threshold signal required to trigger the acetate-diauxic switch until glycerol is fully consumed. The mRNA transcriptional patterns of key gluconeogenic enzymes, namely, ackA, acetate kinase; pta, phosphotransacetylase; acs, acetyl CoA synthetase and aceA, ICL, suggest that the pyruvate phenotype is better equipped than the glycerol phenotype for the switch from acetogenic to gluconeogenic metabolism.

System-level understanding of gene expression and regulation for engineering secondary metabolite production in Streptomyces

The gram-positive bacterium, Streptomyces, is noticed for its ability to produce a wide array of pharmaceutically active compounds through secondary metabolism. To discover novel bioactive secondary metabolites and increase the production, Streptomyces species have been extensively studied for the past decades. Among the cellular components, RNA molecules play important roles as the messengers for gene expression and diverse regulations taking place at the RNA level. Thus, the analysis of RNA-level regulation is critical to understanding the regulation of Streptomyces’ metabolism and secondary metabolite production. A dramatic advance in Streptomyces research was made recently, by exploiting high-throughput technology to systematically understand RNA levels. In this review, we describe the current status of the system-wide investigation of Streptomyces in terms of RNA, toward expansion of its genetic potential for secondary metabolite synthesis.

Programmatic access to bacterial regulatory networks with regutools

SUMMARY

RegulonDB has collected, harmonized and centralized data from hundreds of experiments for nearly two decades and is considered a point of reference for transcriptional regulation in Escherichia coli K12. Here, we present the regutools R package to facilitate programmatic access to RegulonDB data in computational biology. regutools gives researchers the possibility of writing reproducible workflows with automated queries to RegulonDB. The regutools package serves as a bridge between RegulonDB data and the Bioconductor ecosystem by reusing the data structures and statistical methods powered by other Bioconductor packages. We demonstrate the integration of regutools with Bioconductor by analyzing transcription factor DNA binding sites and transcriptional regulatory networks from RegulonDB. We anticipate that regutools will serve as a useful building block in our progress to further our understanding of gene regulatory networks.

AVAILABILITY AND IMPLEMENTATION

regutools is an R package available through Bioconductor at bioconductor.org/packages/regutools.

Networks Consolidate the Core Concepts of Evolution by Natural Selection

Microbiology has unraveled rich evidence of ongoing reticulate evolutionary processes and complex interactions both within and between cells. These phenomena feature real biological networks, which can logically be analyzed using network-based tools. It is thus not surprising that network sciences, a field independent from evolutionary biology and microbiology, have recently pervasively infused their methods into both fields. Importantly, network tools bring forward observations enhancing the understanding of three core evolutionary concepts: variation, fitness, and heredity. Consequently, our work shows how network sciences can enhance evolutionary theory by explaining the evolution by natural selection of a broad diversity of units of selection, while updating the popular figure of Darwin's tree of life with a comprehensive sketch of the networks of evolution.

Comparative genomics analysis of Nitriliruptoria reveals the genomic differences and salt adaptation strategies

The group Nitriliruptoria, recently classified as a separate class of phylum Actinobacteria, has five members at present, which belong to halophilic or halotolerant Actinobacteria. Here, we sequenced the genomes of Egicoccus halophilus EGI 80432^T and Egibacter rhizosphaerae EGI 80759^T, and performed a comparative genomics approach to analyze the genomic differences and salt adaptation mechanisms in Nitriliruptoria. Phylogenetic analysis suggested that Euzebya tangerina F10^T has a closer phylogenetic relationship to Euzebya rosea DSW09^T, while genomic analysis revealed highest genomic similarity with Nitriliruptor alkaliphilus ANL-iso2^T and E. halophilus EGI 80432^T. Genomic differences of Nitriliruptoria were mainly observed in genome size, gene contents, and the amounts of gene in per functional categories. Furthermore, our analysis also revealed that Nitriliruptoria possess similar synthesis systems of solutes, such as trehalose, glutamine, glutamate, and proline. On the other hand, each member of Nitriliruptoria species possesses specific mechanisms, K⁺ influx and efflux, betaine and ectoine synthesis, and compatible solutes transport to survive in various high-salt environments.

Contrasting effects of isocitrate dehydrogenase deletion on fluxes through enzymes of central metabolism in Escherichia coli

Flux analysis is central to understanding cellular metabolism and successful manipulation of metabolic fluxes in microbial cell-factories. Isocitrate dehydrogenase (ICDH) deletion conferred contrasting effects on fluxes through substrate-level phosphorylation (SLP) reactions. While significantly increasing flux through pyruvate kinase, it diminishes flux through succinyl CoA synthetase and upregulates phosphotransacetylase (PTA) and acetate kinase (AK). In addition to acetate, the ICDH-less strain excretes pyruvate, citrate and isocitrate. While efflux to acetate excretion by the Escherichia coli parental strain and its ICDH-less derivative is a reflection of high throughput of glycolytic intermediates, excretion of pyruvate is a reflection of high throughput via pyruvate kinase. On the other hand, citrate and isocitrate excretion is a reflection of truncating the Krebs cycle at the level of ICDH. Furthermore, another striking finding is the inability of the ICDH-less cultures to utilize acetate as a source of carbon despite the availability of an adequate supply of extracellular glutamate (for biosynthesis) and elevated levels of AK and PTA (for acetate uptake). This striking observation is now explicable in the light of the newly proposed hypothesis that the expression of the ace operon enzymes is controlled in response to a minimum threshold signal (ATP), which could not be achieved in the ICDH-less strain.

Disruption of the NlpD lipoprotein of the plague pathogen Yersinia pestis affects iron acquisition and the activity of the twin-arginine translocation system

We have previously shown that the cell morphogenesis NlpD lipoprotein is essential for virulence of the plague bacteria, Yersinia pestis. To elucidate the role of NlpD in Y. pestis pathogenicity, we conducted a whole-genome comparative transcriptome analysis of the wild-type Y. pestis strain and an nlpD mutant under conditions mimicking early stages of infection. The analysis suggested that NlpD is involved in three phenomena: (i) Envelope stability/integrity evidenced by compensatory up-regulation of the Cpx and Psp membrane stress-response systems in the mutant; (ii) iron acquisition, supported by modulation of iron metabolism genes and by limited growth in iron-deprived medium; (iii) activity of the twin-arginine (Tat) system, which translocates folded proteins across the cytoplasmic membrane. Virulence studies of Y. pestis strains mutated in individual Tat components clearly indicated that the Tat system is central in Y. pestis pathogenicity and substantiated the assumption that NlpD essentiality in iron utilization involves the activity of the Tat system. This study reveals a new role for NlpD in Tat system activity and iron assimilation suggesting a modality by which this lipoprotein is involved in Y. pestis pathogenesis.

We have previously shown that the NlpD lipoprotein, which is involved in the regulation of cell morphogenesis, is essential for virulence of the plague bacteria, Yersinia pestis. To uncover the role of NlpD in Y. pestis pathogenicity, we conducted a whole-genome comparative transcriptome analysis as well as phenotypic and virulence evaluation analyses of the nlpD and related mutants. The study reveals a new role for the Y. pestis NlpD lipoprotein in iron assimilation and Tat system activity.

Multifaceted Stoichiometry Control of Bacterial Operons Revealed by Deep Proteome Quantification

More than half of the protein-coding genes in bacteria are organized in polycistronic operons composed of two or more genes. It remains under debate whether the operon organization maintains the stoichiometric expression of the genes within an operon. In this study, we performed a label-free data-independent acquisition hyper reaction monitoring mass-spectrometry (HRM-MS) experiment to quantify the Escherichia coli proteome in exponential phase and quantified 93.6% of the cytosolic proteins, covering 67.9% and 56.0% of the translating polycistronic operons in BW25113 and MG1655 strains, respectively. We found that the translational regulation contributes largely to the proteome complexity: the shorter operons tend to be more tightly controlled for stoichiometry than longer operons; the operons which mainly code for complexes is more tightly controlled for stoichiometry than the operons which mainly code for metabolic pathways. The gene interval (distance between adjacent genes in one operon) may serve as a regulatory factor for stoichiometry. The catalytic efficiency might be a driving force for differential expression of enzymes encoded in one operon. These results illustrated the multifaceted nature of the operon regulation: the operon unified transcriptional level and gene-specific translational level. This multi-level regulation benefits the host by optimizing the efficiency of the productivity of metabolic pathways and maintenance of different types of protein complexes.

Essential Role of σ Factor RpoF in Flagellar Biosynthesis and Flagella-Mediated Motility of Acidithiobacillus caldus

Acidithiobacillaceae, an important family of acidophilic and chemoautotrophic sulfur or iron oxidizers, participate in geobiochemical circulation of the elements and drive the release of heavy metals in mining associated habitats. Because of their environmental adaptability and energy metabolic systems, Acidithiobacillus spp. have become the dominant bacteria used in bioleaching for heavy metal recovery. Flagella-driven motility is associated with bacterial chemotaxis and bacterial responses to environmental stimuli. However, little is known about how the flagellum of Acidithiobacillus spp. is regulated and how the flagellum affects the growth of these chemoautotrophic bacteria. In this study, we analyzed the flagellar gene clusters in Acidithiobacillus strains and uncovered the close relationship between flagella and the sulfur-oxidizing systems (Sox system). The σ²⁸ gene (rpoF) knockout and overexpression strains of Acidithiobacillus caldus were constructed. Scanning electron microscopy shows that A. caldus ΔrpoF cells lacked flagella, indicating the essential role of RpoF in regulating flagella synthesis in these chemoautotrophic bacteria. Motility analysis suggests that the deletion of rpoF resulted in the reduction of swarming capability, while this capability was enhanced in the rpoF overexpression strain. Both static cultivation and low concentration of energy substrates (elemental sulfur or tetrathionate) led to weak growth of A. caldus ΔrpoF cells. The deletion of rpoF promoted bacterial attachment to the surface of elemental sulfur in static cultivation. The absence of RpoF caused an obvious change in transcription profile, including genes in flagellar cluster and those involved in biofilm formation. These results provide an understanding on the regulation of flagellar hierarchy and the flagellar function in these sulfur or iron oxidizers.

Genome sequence and comparative genomics of Rhizobium sp. Td3, a novel plant growth promoting phosphate solubilizing Cajanus cajan symbiont

Rhizobium sp. Td3 is a Sesbania plant growth promoting, Cajanus cajan nodulating rhizobia. Studying its whole genome was important as it is a potent phosphate solubilizer with constitutive gluconic acid production ability through operation of the periplasmic glucose oxidation pathway even under conditions of catabolite repression. This is in contrast to the other explored phosphate solubilizers. Rhizobial isolates sequenced so far are known to lack components of the direct glucose oxidation pathway and cannot produce gluconic acid on its own. Here, we present the genome sequence of Rhizobium sp. Td3. Genome comprises of a single chromosome of size 5,606,547 bp (5.6 Mb) with no symbiotic plasmid. Rhizobium leguminosarum bv. viciae USDA2370 was the closest whole genome known. 109 genes responsible for diverse plant growth promoting activities like P solubilization, synthesis of acetoin, nitric oxide, indole-3 acetic acid, exopolysaccharide, siderophore and trehalose have been identified. Flagellar proteins, genes encoding antibiotic and metal resistance, enzymes required for combating oxidative stress as well as attachment and colonization in the plant rhizosphere are also present. Availability of genome sequence of such a versatile plant growth promoting agent will help in exploiting all the phyto-beneficial traits of Td3 for its use as a biofertilizer.

LexA Binds to Transcription Regulatory Site of Cell Division Gene ftsZ in Toxic Cyanobacterium Microcystis aeruginosa

Previously, we showed that DNA replication and cell division in toxic cyanobacterium Microcystis aeruginosa are coordinated by transcriptional regulation of cell division gene ftsZ and that an unknown protein specifically bound upstream of ftsZ (BpFz; DNA-binding protein to an upstream site of ftsZ) during successful DNA replication and cell division. Here, we purified BpFz from M. aeruginosa strain NIES-298 using DNA-affinity chromatography and gel-slicing combined with gel electrophoresis mobility shift assay (EMSA). The N-terminal amino acid sequence of BpFz was identified as TNLESLTQ, which was identical to that of transcription repressor LexA from NIES-843. EMSA analysis using mutant probes showed that the sequence GTACTAN₃GTGTTC was important in LexA binding. Comparison of the upstream regions of lexA in the genomes of closely related cyanobacteria suggested that the sequence TASTRNNNNTGTWC could be a putative LexA recognition sequence (LexA box). Searches for TASTRNNNNTGTWC as a transcriptional regulatory site (TRS) in the genome of M. aeruginosa NIES-843 showed that it was present in genes involved in cell division, photosynthesis, and extracellular polysaccharide biosynthesis. Considering that BpFz binds to the TRS of ftsZ during normal cell division, LexA may function as a transcriptional activator of genes related to cell reproduction in M. aeruginosa, including ftsZ. This may be an example of informality in the control of bacterial cell division.

An Orthogonal Permease-Inducer-Repressor Feedback Loop Shows Bistability

Feedback loops in biological networks, among others, enable differentiation and cell cycle progression, and increase robustness in signal transduction. In natural networks, feedback loops are often complex and intertwined, making it challenging to identify which loops are mainly responsible for an observed behavior. However, minimal synthetic replicas could allow for such identification. Here, we engineered a synthetic permease-inducer-repressor system in Saccharomyces cerevisiae to analyze if a transport-mediated positive feedback loop could be a core mechanism for the switch-like behavior in the regulation of metabolic gene networks such as the S. cerevisiae GAL system or the Escherichia coli lac operon. We characterized the synthetic circuit using deterministic and stochastic mathematical models. Similar to its natural counterparts, our synthetic system shows bistable and hysteretic behavior, and the inducer concentration range for bistability as well as the switching rates between the two stable states depend on the repressor concentration. Our results indicate that a generic permease-inducer-repressor circuit with a single feedback loop is sufficient to explain the experimentally observed bistable behavior of the natural systems. We anticipate that the approach of reimplementing natural systems with orthogonal parts to identify crucial network components is applicable to other natural systems such as signaling pathways.

Crystal structure and DNA binding activity of a PadR family transcription regulator from hypervirulent Clostridium difficile R20291

Background

Clostridium difficile is a spore-forming obligate anaerobe that can remain viable for extended periods, even in the presence of antibiotics, which contributes to the persistence of this bacterium as a human pathogen during host-to-host transmission and in hospital environments. We examined the structure and function of a gene product with the locus tag CDR20291_0991 (cdPadR1) as part of our broader goal aimed at elucidating transcription regulatory mechanisms involved in virulence and antibiotic resistance of the recently emergent hypervirulent C. difficile strain R20291. cdPadR1 is genomically positioned near genes that are involved in stress response and virulence. In addition, it was previously reported that cdPadR1 and a homologue from the historical C. difficile strain 630 (CD630_1154) were differentially expressed when exposed to stressors, including antibiotics.

Results

The crystal structure of cdPadR1 was determined to 1.9 Å resolution, which revealed that it belongs to the PadR-s2 subfamily of PadR transcriptional regulators. cdPadR1 binds its own promoter and other promoter regions from within the C. difficile R20291 genome. DNA binding experiments demonstrated that cdPadR1 binds a region comprised of inverted repeats and an AT-rich core with the predicted specific binding motif, GTACTAT(N₂)ATTATA(N)AGTA, within its own promoter that is also present in 200 other regions in the C. difficile R20291 genome. Mutation of the highly conserved W in α4 of the effector binding/oligomerization domain, which is predicted to be involved in multi-drug recognition and dimerization in other PadR-s2 proteins, resulted in alterations of cdPadR1 binding to the predicted binding motif, potentially due to loss of higher order oligomerization.

Conclusions

Our results indicate that cdPadR1 binds a region within its own promoter consisting of the binding motif GTACTAT(N₂)ATTATA(N)AGTA and seems to associate non-specifically with longer DNA fragments in vitro, which may facilitate promoter and motif searching. This suggests that cdPadR1 acts as a transcriptional auto-regulator, binding specific sites within its own promoter, and is part of a broad gene regulatory network involved, in part, with environmental stress response, antibiotic resistance and virulence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0850-0) contains supplementary material, which is available to authorized users.

Catabolite control protein A has an important role in the metabolic regulation of Streptococcus suis type 2 according to iTRAQ-based quantitative proteomic analysis

The catabolite control protein A (ccpA) regulates the carbon metabolism in Streptococcus suis type 2 and has pleiotropic regulatory functions in bacterial virulence and transcription. The present study systematically investigated ccpA activity in Streptococcus suis type 2 using isobaric tag for relative and absolute quantification (iTRAQ) liquid chromatography‑tandem mass spectrometry‑based proteomics. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses demonstrated that ccpA is an important protein for the regulation of metabolism, virulence and immune pathways in Streptococcus suis type 2. The present study therefore expanded the current understanding of the effects of ccpA on virulence, metabolic regulation and transcription in Streptococcus suis type 2 and other important pathogens.

How to train your microbe: methods for dynamically characterizing gene networks

Gene networks regulate biological processes dynamically. However, researchers have largely relied upon static perturbations, such as growth media variations and gene knockouts, to elucidate gene network structure and function. Thus, much of the regulation on the path from DNA to phenotype remains poorly understood. Recent studies have utilized improved genetic tools, hardware, and computational control strategies to generate precise temporal perturbations outside and inside of live cells. These experiments have, in turn, provided new insights into the organizing principles of biology. Here, we introduce the major classes of dynamical perturbations that can be used to study gene networks, and discuss technologies available for creating them in a wide range of microbial pathways.

Mesoscopic model and free energy landscape for protein-DNA binding sites: analysis of cyanobacterial promoters

The identification of protein binding sites in promoter sequences is a key problem to understand and control regulation in biochemistry and biotechnological processes. We use a computational method to analyze promoters from a given genome. Our approach is based on a physical model at the mesoscopic level of protein-DNA interaction based on the influence of DNA local conformation on the dynamics of a general particle along the chain. Following the proposed model, the joined dynamics of the protein particle and the DNA portion of interest, only characterized by its base pair sequence, is simulated. The simulation output is analyzed by generating and analyzing the Free Energy Landscape of the system. In order to prove the capacity of prediction of our computational method we have analyzed nine promoters of Anabaena PCC 7120. We are able to identify the transcription starting site of each of the promoters as the most populated macrostate in the dynamics. The developed procedure allows also to characterize promoter macrostates in terms of thermo-statistical magnitudes (free energy and entropy), with valuable biological implications. Our results agree with independent previous experimental results. Thus, our methods appear as a powerful complementary tool for identifying protein binding sites in promoter sequences.

Binding of specific proteins to particular sites in the DNA sequence is a fundamental issue for gene regulation in molecular biology and genetic engineering. A deep understanding of cell physiology requires the analysis of a plethora of genes involving characterization of their promoter architectures that determine their regulation and gene transcription. In order to locate the promoter elements of a given gene, experimental determination of its transcription start site (TSS) is required. This is an expensive, time-consuming task that, depending on our requirements, could be simplified using computational analysis as a first approach. Nevertheless, most computational methods lack a physical basis on the protein-DNA interaction mechanism. We adopt here this strategy, by using a simple model for protein-DNA interaction to find TSS in a bunch of cyanobacteria promoters. We make use of physical tools to characterize these TSS and to relate them with biological properties as the relative strength of the promoter. Our study shows how a model based on a coarse-grained description of a biomolecule can give valuable insight on its biological function.

Comprehensive mapping of the Escherichia coli flagellar regulatory network

Flagellar synthesis is a highly regulated process in all motile bacteria. In Escherichia coli and related species, the transcription factor FlhDC is the master regulator of a multi-tiered transcription network. FlhDC activates transcription of a number of genes, including some flagellar genes and the gene encoding the alternative Sigma factor FliA. Genes whose expression is required late in flagellar assembly are primarily transcribed by FliA, imparting temporal regulation of transcription and coupling expression to flagellar assembly. In this study, we use ChIP-seq and RNA-seq to comprehensively map the E. coli FlhDC and FliA regulons. We define a surprisingly restricted FlhDC regulon, including two novel regulated targets and two binding sites not associated with detectable regulation of surrounding genes. In contrast, we greatly expand the known FliA regulon. Surprisingly, 30 of the 52 FliA binding sites are located inside genes. Two of these intragenic promoters are associated with detectable noncoding RNAs, while the others either produce highly unstable RNAs or are inactive under these conditions. Together, our data redefine the E. coli flagellar regulatory network, and provide new insight into the temporal orchestration of gene expression that coordinates the flagellar assembly process.

Flagella are surface-associated appendages that propel bacteria and are involved in diverse functions such as chemotaxis, surface attachment, and host cell invasion. Flagella are incredibly complex macromolecular machines that are energetically costly to produce, assemble, and power. Flagellar production is tightly regulated and flagellar components are only synthesized when flagellar motility is advantageous. Regulation also ensures that flagellar components are produced in roughly the same order in which they are needed, increasing efficiency of the assembly process. The transcriptional regulation of flagellar genes has been studied extensively in the model organism Escherichia coli; however, many questions remain. We have used an unbiased, genome-wide approach to comprehensively identify all of the binding sites and regulatory targets for the two key regulators of flagellar synthesis, FlhDC and FliA. Our results redefine the flagellar regulatory network, and suggest that FliA binds many sites that are not associated with productive transcription. This work is important because it suggests possible new functions for FliA outside of the transcription of canonical mRNAs, and it provides new insight into the temporal orchestration of gene expression that coordinates the flagellar assembly process.

The CcpA regulon of Streptococcus suis reveals novel insights into the regulation of the streptococcal central carbon metabolism by binding of CcpA to two distinct binding motifs

Streptococcus suis (S. suis) is a neglected zoonotic streptococcus causing fatal diseases in humans and in pigs. The transcriptional regulator CcpA (catabolite control protein A) is involved in the metabolic adaptation to different carbohydrate sources and virulence of S. suis and other pathogenic streptococci. In this study, we determined the DNA binding characteristics of CcpA and identified the CcpA regulon during growth of S. suis. Electrophoretic mobility shift analyses showed promiscuous DNA binding of CcpA to cognate cre sites in vitro. In contrast, sequencing of immunoprecipitated chromatin revealed two specific consensus motifs, a pseudo-palindromic cre motif (WWGAAARCGYTTTCWW) and a novel cre2 motif (TTTTYHWDHHWWTTTY), within the regulatory elements of the genes directly controlled by CcpA. Via these elements CcpA regulates expression of genes involved in carbohydrate uptake and conversion, and in addition in important metabolic pathways of the central carbon metabolism, like glycolysis, mixed-acid fermentation, and the fragmentary TCA cycle. Furthermore, our analyses provide evidence that CcpA regulates the genes of the central carbon metabolism by binding either the pseudo-palindromic cre motif or the cre2 motif in a HPr(Ser)∼P independent conformation.

Regulation Systems of Bacteria such as Escherichia coli in Response to Nutrient Limitation and Environmental Stresses

An overview was made to understand the regulation system of a bacterial cell such as Escherichia coli in response to nutrient limitation such as carbon, nitrogen, phosphate, sulfur, ion sources, and environmental stresses such as oxidative stress, acid shock, heat shock, and solvent stresses. It is quite important to understand how the cell detects environmental signals, integrate such information, and how the cell system is regulated. As for catabolite regulation, F1,6B P (FDP), PEP, and PYR play important roles in enzyme level regulation together with transcriptional regulation by such transcription factors as Cra, Fis, CsrA, and cAMP-Crp. αKG plays an important role in the coordinated control between carbon (C)- and nitrogen (N)-limitations, where αKG inhibits enzyme I (EI) of phosphotransferase system (PTS), thus regulating the glucose uptake rate in accordance with N level. As such, multiple regulation systems are co-ordinated for the cell synthesis and energy generation against nutrient limitations and environmental stresses. As for oxidative stress, the TCA cycle both generates and scavenges the reactive oxygen species (ROSs), where NADPH produced at ICDH and the oxidative pentose phosphate pathways play an important role in coping with oxidative stress. Solvent resistant mechanism was also considered for the stresses caused by biofuels and biochemicals production in the cell.

Murine vaginal colonization model for investigating asymptomatic mucosal carriage of Streptococcus pyogenes

While many virulence factors promoting Streptococcus pyogenes invasive disease have been described, specific streptococcal factors and host properties influencing asymptomatic mucosal carriage remain uncertain. To address the need for a refined model of prolonged S. pyogenes asymptomatic mucosal colonization, we have adapted a preestrogenized murine vaginal colonization model for S. pyogenes. In this model, derivatives of strains HSC5, SF370, JRS4, NZ131, and MEW123 established a reproducible, asymptomatic colonization of the vaginal mucosa over a period of typically 3 to 4 weeks' duration at a relatively high colonization efficiency. Prior treatment with estradiol prolonged streptococcal colonization and was associated with reduced inflammation in the colonized vaginal epithelium as well as a decreased leukocyte presence in vaginal fluid compared to the levels of inflammation and leukocyte presence in non-estradiol-treated control mice. The utility of our model for investigating S. pyogenes factors contributing to mucosal carriage was verified, as a mutant with a mutation in the transcriptional regulator catabolite control protein A (CcpA) demonstrated significant impairment in vaginal colonization. An assessment of in vivo transcriptional activity in the CcpA(-) strain for several known CcpA-regulated genes identified significantly elevated transcription of lactate oxidase (lctO) correlating with excessive generation of hydrogen peroxide to self-lethal levels. Deletion of lctO did not impair colonization, but deletion of lctO in a CcpA(-) strain prolonged carriage, exceeding even that of the wild-type strain. Thus, while LctO is not essential for vaginal colonization, its dysregulation is deleterious, highlighting the critical role of CcpA in promoting mucosal colonization. The vaginal colonization model should prove effective for future analyses of S. pyogenes mucosal colonization.

Neutral forces acting on intragenomic variability shape the Escherichia coli regulatory network topology

Cis-regulatory networks (CRNs) play a central role in cellular decision making. Like every other biological system, CRNs undergo evolution, which shapes their properties by a combination of adaptive and nonadaptive evolutionary forces. Teasing apart these forces is an important step toward functional analyses of the different components of CRNs, designing regulatory perturbation experiments, and constructing synthetic networks. Although tests of neutrality and selection based on molecular sequence data exist, no such tests are currently available based on CRNs. In this work, we present a unique genotype model of CRNs that is grounded in a genomic context and demonstrate its use in identifying portions of the CRN with properties explainable by neutral evolutionary forces at the system, subsystem, and operon levels. We leverage our model against experimentally derived data from Escherichia coli. The results of this analysis show statistically significant and substantial neutral trends in properties previously identified as adaptive in origin--degree distribution, clustering coefficient, and motifs--within the E. coli CRN. Our model captures the tightly coupled genome-interactome of an organism and enables analyses of how evolutionary events acting at the genome level, such as mutation, and at the population level, such as genetic drift, give rise to neutral patterns that we can quantify in CRNs.

Metabolic Regulation of a Bacterial Cell System with Emphasis on Escherichia coli Metabolism

It is quite important to understand the overall metabolic regulation mechanism of bacterial cells such as Escherichia coli from both science (such as biochemistry) and engineering (such as metabolic engineering) points of view. Here, an attempt was made to clarify the overall metabolic regulation mechanism by focusing on the roles of global regulators which detect the culture or growth condition and manipulate a set of metabolic pathways by modulating the related gene expressions. For this, it was considered how the cell responds to a variety of culture environments such as carbon (catabolite regulation), nitrogen, and phosphate limitations, as well as the effects of oxygen level, pH (acid shock), temperature (heat shock), and nutrient starvation.

Evolutionary consequences of antibiotic use for the resistome, mobilome and microbial pangenome

The widespread use and abuse of antibiotic therapy has evolutionary and ecological consequences, some of which are only just beginning to be examined. One well known consequence is the fixation of mutations and lateral gene transfer (LGT) events that confer antibiotic resistance. Sequential selection events, driven by different classes of antibiotics, have resulted in the assembly of diverse resistance determinants and mobile DNAs into novel genetic elements of ever-growing complexity and flexibility. These novel plasmids, integrons, and genomic islands have now become fixed at high frequency in diverse cell lineages by human antibiotic use. Consequently they can be regarded as xenogenetic pollutants, analogous to xenobiotic compounds, but with the critical distinction that they replicate rather than degrade when released to pollute natural environments. Antibiotics themselves must also be regarded as pollutants, since human production overwhelms natural synthesis, and a major proportion of ingested antibiotic is excreted unchanged into waste streams. Such antibiotic pollutants have non-target effects, raising the general rates of mutation, recombination, and LGT in all the microbiome, and simultaneously providing the selective force to fix such changes. This has the consequence of recruiting more genes into the resistome and mobilome, and of increasing the overlap between these two components of microbial genomes. Thus the human use and environmental release of antibiotics is having second order effects on the microbial world, because these small molecules act as drivers of bacterial evolution. Continued pollution with both xenogenetic elements and the selective agents that fix such elements in populations has potentially adverse consequences for human welfare.

Gene network requirements for regulation of metabolic gene expression to a desired state

Gene circuits that control metabolism should restore metabolic functions upon environmental changes. Whether gene networks are capable of steering metabolism to optimal states is an open question. Here we present a method to identify such optimal gene networks. We show that metabolic network optimisation over a range of environments results in an input-output relationship for the gene network that guarantees optimal metabolic states. Optimal control is possible if the gene network can achieve this input-output relationship. We illustrate our approach with the best-studied regulatory network in yeast, the galactose network. We find that over the entire range of external galactose concentrations, the regulatory network is able to optimally steer galactose metabolism. Only a few gene network parameters affect this optimal regulation. The other parameters can be tuned independently for optimisation of other functions, such as fast and low-noise gene expression. This study highlights gene network plasticity, evolvability, and modular functionality.

How to cluster gene expression dynamics in response to environmental signals

Organisms usually cope with change in the environment by altering the dynamic trajectory of gene expression to adjust the complement of active proteins. The identification of particular sets of genes whose expression is adaptive in response to environmental changes helps to understand the mechanistic base of gene-environment interactions essential for organismic development. We describe a computational framework for clustering the dynamics of gene expression in distinct environments through Gaussian mixture fitting to the expression data measured at a set of discrete time points. We outline a number of quantitative testable hypotheses about the patterns of dynamic gene expression in changing environments and gene-environment interactions causing developmental differentiation. The future directions of gene clustering in terms of incorporations of the latest biological discoveries and statistical innovations are discussed. We provide a set of computational tools that are applicable to modeling and analysis of dynamic gene expression data measured in multiple environments.

Counteractive balancing of transcriptome expression involving CodY and CovRS in Streptococcus pyogenes

Streptococcus pyogenes (group A streptococcus [GAS]) responds to environmental changes in a manner that results in an adaptive regulation of the transcriptome. The objective of the present study was to understand how two global transcriptional regulators, CodY and CovRS, coordinate the transcriptional network in S. pyogenes. Results from expression microarray data and quantitative reverse transcription-PCR (qRT-PCR) showed that the global regulator CodY controls the expression of about 250 genes, or about 17% of the genome of strain NZ131. Additionally, the codY gene was shown to be negatively autoregulated, with its protein binding directly to the promoter region with a CodY binding site. In further studies, the influence of codY, covRS, and codY-covRS mutations on gene expression was analyzed in growth phase-dependent conditions using C medium, reported to mimic nutritional abundance and famine conditions similar to those found during host GAS infection. Additional biological experiments of several virulence phenotypes, including pilin production, biofilm formation, and NAD glycohydrolase activity, demonstrated the role that both CodY and CovRS play in their regulation. Correlation analysis of the overall data revealed that, in exponentially growing cells, CodY and CovRS act in opposite directions, with CodY stimulating and CovRS repressing a substantial fraction of the core genome, including many virulence factors. This is the first report of counteractive balancing of transcriptome expression by global transcription regulators and provides important insight into how GAS modulates gene expression by integrating important extracellular and intracellular information.