iModulonMiner and PyModulon: Software for unsupervised mining of gene expression compendia.

Public gene expression databases are a rapidly expanding resource of organism responses to diverse perturbations, presenting both an opportunity and a challenge for bioinformatics workflows to extract actionable knowledge of transcription regulatory network function. Here, we introduce a five-step computational pipeline, called iModulonMiner, to compile, process, curate, analyze, and characterize the totality of RNA-seq data for a given organism or cell type. This workflow is centered around the data-driven computation of co-regulated gene sets using Independent Component Analysis, called iModulons, which have been shown to have broad applications. As a demonstration, we applied this workflow to generate the iModulon structure of Bacillus subtilis using all high-quality, publicly-available RNA-seq data. Using this structure, we predicted regulatory interactions for multiple transcription factors, identified groups of co-expressed genes that are putatively regulated by undiscovered transcription factors, and predicted properties of a recently discovered single-subunit phage RNA polymerase. We also present a Python package, PyModulon, with functions to characterize, visualize, and explore computed iModulons. The pipeline, available at https://github.com/SBRG/iModulonMiner, can be readily applied to diverse organisms to gain a rapid understanding of their transcriptional regulatory network structure and condition-specific activity.

The Escherichia coli Fur pan-regulon has few conserved but many unique regulatory targets.

While global transcription factors (TFs) have been studied extensively in Escherichia coli model strains, conservation and diversity in TF regulation between strains is still unknown. Here we use a combination of ChIP-exo-to define ferric uptake regulator (Fur) binding sites-and differential gene expression-to define the Fur regulon in nine E. coli strains. We then define a pan-regulon consisting of 469 target genes that includes all Fur target genes in all nine strains. The pan-regulon is then divided into the core regulon (target genes found in all the strains, n = 36), the accessory regulon (target found in two to eight strains, n = 158) and the unique regulon (target genes found in one strain, n = 275). Thus, there is a small set of Fur regulated genes common to all nine strains, but a large number of regulatory targets unique to a particular strain. Many of the unique regulatory targets are genes unique to that strain. This first-established pan-regulon reveals a common core of conserved regulatory targets and significant diversity in transcriptional regulation amongst E. coli strains, reflecting diverse niche specification and strain history.

Adaptive laboratory evolution and independent component analysis disentangle complex vancomycin adaptation trajectories.

Human infections with methicillin-resistant Staphylococcus aureus (MRSA) are commonly treated with vancomycin, and strains with decreased susceptibility, designated as vancomycin-intermediate S. aureus (VISA), are associated with treatment failure. Here, we profiled the phenotypic, mutational, and transcriptional landscape of 10 VISA strains adapted by laboratory evolution from one common MRSA ancestor, the USA300 strain JE2. Using functional and independent component analysis, we found that: 1) despite the common genetic background and environmental conditions, the mutational landscape diverged between evolved strains and included mutations previously associated with vancomycin resistance (in vraT, graS, vraFG, walKR, and rpoBCD) as well as novel adaptive mutations (SAUSA300_RS04225, ssaA, pitAR, and sagB); 2) the first wave of mutations affected transcriptional regulators and the second affected genes involved in membrane biosynthesis; 3) expression profiles were predominantly strain-specific except for sceD and lukG, which were the only two genes significantly differentially expressed in all clones; 4) three independent virulence systems (φSa3, SaeR, and T7SS) featured as the most transcriptionally perturbed gene sets across clones; 5) there was a striking variation in oxacillin susceptibility across the evolved lineages (from a 10-fold increase to a 63-fold decrease) that also arose in clinical MRSA isolates exposed to vancomycin and correlated with susceptibility to teichoic acid inhibitors; and 6) constitutive expression of the VraR regulon explained cross-susceptibility, while mutations in walK were associated with cross-resistance. Our results show that adaptation to vancomycin involves a surprising breadth of mutational and transcriptional pathways that affect antibiotic susceptibility and possibly the clinical outcome of infections.

Metabolic and genetic basis for auxotrophies in Gram-negative species.

Auxotrophies constrain the interactions of bacteria with their environment, but are often difficult to identify. Here, we develop an algorithm (AuxoFind) using genome-scale metabolic reconstruction to predict auxotrophies and apply it to a series of available genome sequences of over 1,300 Gram-negative strains. We identify 54 auxotrophs, along with the corresponding metabolic and genetic basis, using a pangenome approach, and highlight auxotrophies conferring a fitness advantage in vivo. We show that the metabolic basis of auxotrophy is species-dependent and varies with 1) pathway structure, 2) enzyme promiscuity, and 3) network redundancy. Various levels of complexity constitute the genetic basis, including 1) deleterious single-nucleotide polymorphisms (SNPs), in-frame indels, and deletions; 2) single/multigene deletion; and 3) movement of mobile genetic elements (including prophages) combined with genomic rearrangements. Fourteen out of 19 predictions agree with experimental evidence, with the remaining cases highlighting shortcomings of sequencing, assembly, annotation, and reconstruction that prevent predictions of auxotrophies. We thus develop a framework to identify the metabolic and genetic basis for auxotrophies in Gram-negatives.

Revealing 29 sets of independently modulated genes in Staphylococcus aureus, their regulators, and role in key physiological response.

The ability of Staphylococcus aureus to infect many different tissue sites is enabled, in part, by its transcriptional regulatory network (TRN) that coordinates its gene expression to respond to different environments. We elucidated the organization and activity of this TRN by applying independent component analysis to a compendium of 108 RNA-sequencing expression profiles from two S. aureus clinical strains (TCH1516 and LAC). ICA decomposed the S. aureus transcriptome into 29 independently modulated sets of genes (i-modulons) that revealed: 1) High confidence associations between 21 i-modulons and known regulators; 2) an association between an i-modulon and σS, whose regulatory role was previously undefined; 3) the regulatory organization of 65 virulence factors in the form of three i-modulons associated with AgrR, SaeR, and Vim-3; 4) the roles of three key transcription factors (CodY, Fur, and CcpA) in coordinating the metabolic and regulatory networks; and 5) a low-dimensional representation, involving the function of few transcription factors of changes in gene expression between two laboratory media (RPMI, cation adjust Mueller Hinton broth) and two physiological media (blood and serum). This representation of the TRN covers 842 genes representing 76% of the variance in gene expression that provides a quantitative reconstruction of transcriptional modules in S. aureus, and a platform enabling its full elucidation.

Adaptive evolution reveals a tradeoff between growth rate and oxidative stress during naphthoquinone-based aerobic respiration.

Evolution fine-tunes biological pathways to achieve a robust cellular physiology. Two and a half billion years ago, rapidly rising levels of oxygen as a byproduct of blooming cyanobacterial photosynthesis resulted in a redox upshift in microbial energetics. The appearance of higher-redox-potential respiratory quinone, ubiquinone (UQ), is believed to be an adaptive response to this environmental transition. However, the majority of bacterial species are still dependent on the ancient respiratory quinone, naphthoquinone (NQ). Gammaproteobacteria can biosynthesize both of these respiratory quinones, where UQ has been associated with aerobic lifestyle and NQ with anaerobic lifestyle. We engineered an obligate NQ-dependent γ-proteobacterium, Escherichia coli ΔubiC, and performed adaptive laboratory evolution to understand the selection against the use of NQ in an oxic environment and also the adaptation required to support the NQ-driven aerobic electron transport chain. A comparative systems-level analysis of pre- and postevolved NQ-dependent strains revealed a clear shift from fermentative to oxidative metabolism enabled by higher periplasmic superoxide defense. This metabolic shift was driven by the concerted activity of 3 transcriptional regulators (PdhR, RpoS, and Fur). Analysis of these findings using a genome-scale model suggested that resource allocation to reactive oxygen species (ROS) mitigation results in lower growth rates. These results provide a direct elucidation of a resource allocation tradeoff between growth rate and ROS mitigation costs associated with NQ usage under oxygen-replete condition.

OxyR Is a Convergent Target for Mutations Acquired during Adaptation to Oxidative Stress-Prone Metabolic States.

Oxidative stress is concomitant with aerobic metabolism. Thus, bacterial genomes encode elaborate mechanisms to achieve redox homeostasis. Here we report that the peroxide-sensing transcription factor, oxyR, is a common mutational target using bacterial species belonging to two genera, Escherichia coli and Vibrio natriegens, in separate growth conditions implemented during laboratory evolution. The mutations clustered in the redox active site, dimer interface, and flexible redox loop of the protein. These mutations favor the oxidized conformation of OxyR that results in constitutive expression of the genes it regulates. Independent component analysis of the transcriptome revealed that the constitutive activity of OxyR reduces DNA damage from reactive oxygen species, as inferred from the activity of the SOS response regulator LexA. This adaptation to peroxide stress came at a cost of lower growth, as revealed by calculations of proteome allocation using genome-scale models of metabolism and macromolecular expression. Further, identification of similar sequence changes in natural isolates of E. coli indicates that adaptation to oxidative stress through genetic changes in oxyR can be a common occurrence.

A computational knowledge-base elucidates the response of Staphylococcus aureus to different media types.

S. aureus is classified as a serious threat pathogen and is a priority that guides the discovery and development of new antibiotics. Despite growing knowledge of S. aureus metabolic capabilities, our understanding of its systems-level responses to different media types remains incomplete. Here, we develop a manually reconstructed genome-scale model (GEM-PRO) of metabolism with 3D protein structures for S. aureus USA300 str. JE2 containing 854 genes, 1,440 reactions, 1,327 metabolites and 673 3-dimensional protein structures. Computations were in 85% agreement with gene essentiality data from random barcode transposon site sequencing (RB-TnSeq) and 68% agreement with experimental physiological data. Comparisons of computational predictions with experimental observations highlight: 1) cases of non-essential biomass precursors; 2) metabolic genes subject to transcriptional regulation involved in Staphyloxanthin biosynthesis; 3) the essentiality of purine and amino acid biosynthesis in synthetic physiological media; and 4) a switch to aerobic fermentation upon exposure to extracellular glucose elucidated as a result of integrating time-course of quantitative exo-metabolomics data. An up-to-date GEM-PRO thus serves as a knowledge-based platform to elucidate S. aureus' metabolic response to its environment.

BOFdat: Generating biomass objective functions for genome-scale metabolic models from experimental data.

Genome-scale metabolic models (GEMs) are mathematically structured knowledge bases of metabolism that provide phenotypic predictions from genomic information. GEM-guided predictions of growth phenotypes rely on the accurate definition of a biomass objective function (BOF) that is designed to include key cellular biomass components such as the major macromolecules (DNA, RNA, proteins), lipids, coenzymes, inorganic ions and species-specific components. Despite its importance, no standardized computational platform is currently available to generate species-specific biomass objective functions in a data-driven, unbiased fashion. To fill this gap in the metabolic modeling software ecosystem, we implemented BOFdat, a Python package for the definition of a Biomass Objective Function from experimental data. BOFdat has a modular implementation that divides the BOF definition process into three independent modules defined here as steps: 1) the coefficients for major macromolecules are calculated, 2) coenzymes and inorganic ions are identified and their stoichiometric coefficients estimated, 3) the remaining species-specific metabolic biomass precursors are algorithmically extracted in an unbiased way from experimental data. We used BOFdat to reconstruct the BOF of the Escherichia coli model iML1515, a gold standard in the field. The BOF generated by BOFdat resulted in the most concordant biomass composition, growth rate, and gene essentiality prediction accuracy when compared to other methods. Installation instructions for BOFdat are available in the documentation and the source code is available on GitHub (https://github.com/jclachance/BOFdat).

Strain-Specific Metabolic Requirements Revealed by a Defined Minimal Medium for Systems Analyses of Staphylococcus aureus.

Staphylococcus aureus is a Gram-positive pathogenic bacterium that colonizes an estimated one-third of the human population and can cause a wide spectrum of disease, ranging from superficial skin infections to life-threatening sepsis. The adaptive mechanisms that contribute to the success of this pathogen remain obscure partially due to a lack of knowledge of its metabolic requirements. Systems biology approaches can be extremely useful in predicting and interpreting metabolic phenotypes; however, such approaches rely on a chemically defined minimal medium as a basis to investigate the requirements of the cell. In this study, a chemically defined minimal medium formulation, termed synthetic minimal medium (SMM), was investigated and validated to support growth of three S. aureus strains: LAC and TCH1516 (USA300 lineage), as well as D592 (USA100 lineage). The formulated SMM was used in an adaptive laboratory evolution experiment to probe the various mutational trajectories of all three strains leading to optimized growth capabilities. The evolved strains were phenotypically characterized for their growth rate and antimicrobial susceptibility. Strains were also resequenced to examine the genetic basis for observed changes in phenotype and to design follow-up metabolite supplementation assays. Our results reveal evolutionary trajectories that arose from strain-specific metabolic requirements. SMM and the evolved strains can also serve as important tools to study antibiotic resistance phenotypes of S. aureusIMPORTANCE As researchers try to understand and combat the development of antibiotic resistance in pathogens, there is a growing need to thoroughly understand the physiology and metabolism of the microbes. Staphylococcus aureus is a threatening pathogen with increased antibiotic resistance and well-studied virulence mechanisms. However, the adaptive mechanisms used by this pathogen to survive environmental stresses remain unclear, mostly due to the lack of information about its metabolic requirements. Defining the minimal metabolic requirements for S. aureus growth is a first step toward unraveling the mechanisms by which it adapts to metabolic stresses. Here, we present the development of a chemically defined minimal medium supporting growth of three S. aureus strains, and we reveal key genetic mutations contributing to improved growth in minimal medium.

Elucidating the CodY regulon in Staphylococcus aureus USA300 substrains TCH1516 and LAC.

CodY is a conserved broad-acting transcription factor that regulates the expression of genes related to amino acid metabolism and virulence in Gram-positive bacteria. Here, we performed the first in vivo determination of CodY target genes using a novel CodY monoclonal antibody in methicillin-resistant Staphylococcus aureus (MRSA) USA300. Our results showed (i) the same 135 CodY promoter binding sites regulating the 165 target genes identified in two closely related virulent S. aureus USA300 TCH1516 and LAC strains; (ii) the differential binding intensity for the same target genes under the same conditions was due to sequence differences in the same CodY-binding site in the two strains; (iii) a CodY regulon comprising 72 target genes that are differentially regulated relative to a CodY deletion strain, representing genes that are mainly involved in amino acid transport and metabolism, inorganic ion transport and metabolism, transcription and translation, and virulence, all based on transcriptomic data; and (iv) CodY systematically regulated central metabolic flux to generate branched-chain amino acids (BCAAs) by mapping the CodY regulon onto a genome-scale metabolic model of S. aureus. Our study performed the first system-level analysis of CodY in two closely related USA300 TCH1516 and LAC strains, revealing new insights into the similarities and differences of CodY regulatory roles between the closely related strains. IMPORTANCE With the increasing availability of whole-genome sequences for many strains within the same pathogenic species, a comparative analysis of key regulators is needed to understand how the different strains uniquely coordinate metabolism and expression of virulence. To successfully infect the human host, Staphylococcus aureus USA300 relies on the transcription factor CodY to reorganize metabolism and express virulence factors. While CodY is a known key transcription factor, its target genes are not characterized on a genome-wide basis. We performed a comparative analysis to describe the transcriptional regulation of CodY between two dominant USA300 strains. This study motivates the characterization of common pathogenic strains and an evaluation of the possibility of developing specialized treatments for major strains circulating in the population.

Phenotypic response of yeast metabolic network to availability of proteinogenic amino acids.

Genome-scale metabolism can best be described as a highly interconnected network of biochemical reactions and metabolites. The flow of metabolites, i.e., flux, throughout these networks can be predicted and analyzed using approaches such as flux balance analysis (FBA). By knowing the network topology and employing only a few simple assumptions, FBA can efficiently predict metabolic functions at the genome scale as well as microbial phenotypes. The network topology is represented in the form of genome-scale metabolic models (GEMs), which provide a direct mapping between network structure and function via the enzyme-coding genes and corresponding metabolic capacity. Recently, the role of protein limitations in shaping metabolic phenotypes have been extensively studied following the reconstruction of enzyme-constrained GEMs. This framework has been shown to significantly improve the accuracy of predicting microbial phenotypes, and it has demonstrated that a global limitation in protein availability can prompt the ubiquitous metabolic strategy of overflow metabolism. Being one of the most abundant and differentially expressed proteome sectors, metabolic proteins constitute a major cellular demand on proteinogenic amino acids. However, little is known about the impact and sensitivity of amino acid availability with regards to genome-scale metabolism. Here, we explore these aspects by extending on the enzyme-constrained GEM framework by also accounting for the usage of amino acids in expressing the metabolic proteome. Including amino acids in an enzyme-constrained GEM of Saccharomyces cerevisiae, we demonstrate that the expanded model is capable of accurately reproducing experimental amino acid levels. We further show that the metabolic proteome exerts variable demands on amino acid supplies in a condition-dependent manner, suggesting that S. cerevisiae must have evolved to efficiently fine-tune the synthesis of amino acids for expressing its metabolic proteins in response to changes in the external environment. Finally, our results demonstrate how the metabolic network of S. cerevisiae is robust towards perturbations of individual amino acids, while simultaneously being highly sensitive when the relative amino acid availability is set to mimic a priori distributions of both yeast and non-yeast origins.

Pan-Genome Analysis of Transcriptional Regulation in Six Salmonella enterica Serovar Typhimurium Strains Reveals Their Different Regulatory Structures.

Establishing transcriptional regulatory networks (TRNs) in bacteria has been limited to well-characterized model strains. Using machine learning methods, we established the transcriptional regulatory networks of six Salmonella enterica serovar Typhimurium strains from their transcriptomes. By decomposing a compendia of RNA sequencing (RNA-seq) data with independent component analysis, we obtained 400 independently modulated sets of genes, called iModulons. We (i) performed pan-genome analysis of the phylogroup structure of S. Typhimurium and analyzed the iModulons against this background, (ii) revealed different genetic signatures in pathogenicity islands that explained phenotypes, (iii) discovered three transport iModulons linked to antibiotic resistance, (iv) described concerted responses to cationic antimicrobial peptides, and (v) uncovered new regulons. Thus, by combining pan-genome and transcriptomic analytics, we revealed variations in TRNs across six strains of serovar Typhimurium. IMPORTANCE Salmonella enterica serovar Typhimurium is a pathogen involved in human nontyphoidal infections. Treating S. Typhimurium infections is difficult due to the species's dynamic adaptation to its environment, which is dictated by a complex transcriptional regulatory network (TRN) that is different across strains. In this study, we describe the use of independent component analysis to characterize the differential TRNs across the S. Typhimurium pan-genome using a compendium of high-quality RNA-seq data. This approach provided unprecedented insights into the differences between regulation of key cellular functions and pathogenicity in the different strains. The study provides an impetus to initiate a large-scale effort to reveal the TRN differences between the major phylogroups of the pathogenic bacteria, which could fundamentally impact personalizing treatments of bacterial pathogens.

Path to improving the life cycle and quality of genome-scale models of metabolism.

Genome-scale models of metabolism (GEMs) are key computational tools for the systems-level study of metabolic networks. Here, we describe the "GEM life cycle," which we subdivide into four stages: inception, maturation, specialization, and amalgamation. We show how different types of GEM reconstruction workflows fit in each stage and proceed to highlight two fundamental bottlenecks for GEM quality improvement: GEM maturation and content removal. We identify common characteristics contributing to increasing quality of maturing GEMs drawing from past independent GEM maturation efforts. We then shed some much-needed light on the latent and unrecognized but pervasive issue of content removal, demonstrating the substantial effects of model pruning on its solution space. Finally, we propose a novel framework for content removal and associated confidence-level assignment which will help guide future GEM development efforts, reduce duplication of effort across groups, potentially aid automated reconstruction platforms, and boost the reproducibility of model development.

Pangenome Analytics Reveal Two-Component Systems as Conserved Targets in ESKAPEE Pathogens.

The two-component system (TCS) helps bacteria sense and respond to environmental stimuli through histidine kinases and response regulators. TCSs are the largest family of multistep signal transduction processes, and they are involved in many important cellular processes such as antibiotic resistance, pathogenicity, quorum sensing, osmotic stress, and biofilms. Here, we perform the first comprehensive study to highlight the role of TCSs as potential drug targets against ESKAPEE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter spp., and Escherichia coli) pathogens through annotation, mapping, pangenomic status, gene orientation, and sequence variation analysis. The distribution of the TCSs is group specific with regard to Gram-positive and Gram-negative bacteria, except for KdpDE. The TCSs among ESKAPEE pathogens form closed pangenomes, except for Pseudomonas aeruginosa Furthermore, their conserved nature due to closed pangenomes might make them good drug targets. Fitness score analysis suggests that any mutation in some TCSs such as BaeSR, ArcBA, EvgSA, and AtoSC, etc., might be lethal to the cell. Taken together, the results of this pangenomic assessment of TCSs reveal a range of strategies deployed by the ESKAPEE pathogens to manifest pathogenicity and antibiotic resistance. This study further suggests that the conserved features of TCSs might make them an attractive group of potential targets with which to address antibiotic resistance.IMPORTANCE The ESKAPEE pathogens are the leading cause of health care-associated infections worldwide. Two-component systems (TCSs) can be used as effective targets against pathogenic bacteria since they are ubiquitous and manage various vital functions such as antibiotic resistance, virulence, biofilms, quorum sensing, and pH balance, among others. This study provides a comprehensive overview of the pangenomic status of the TCSs among ESKAPEE pathogens. The annotation and pangenomic analysis of TCSs show that they are significantly distributed and conserved among the pathogens, as most of them form closed pangenomes. Furthermore, our analysis also reveals that the removal of the TCSs significantly affects the fitness of the cell. Hence, they may be used as promising drug targets against bacteria.

Environmental conditions dictate differential evolution of vancomycin resistance in Staphylococcus aureus.

While microbiological resistance to vancomycin in Staphylococcus aureus is rare, clinical vancomycin treatment failures are common, and methicillin-resistant S. aureus (MRSA) strains isolated from patients after prolonged vancomycin treatment failure remain susceptible. Adaptive laboratory evolution was utilized to uncover mutational mechanisms associated with MRSA vancomycin resistance in a physiological medium as well as a bacteriological medium used in clinical susceptibility testing. Sequencing of resistant clones revealed shared and media-specific mutational outcomes, with an overlap in cell wall regulons (walKRyycHI, vraSRT). Evolved strains displayed similar properties to resistant clinical isolates in their genetic and phenotypic traits. Importantly, resistant phenotypes that developed in physiological media did not translate into resistance in bacteriological media. Further, a bacteriological media-specific mechanism for vancomycin resistance associated with a mutated mprF was confirmed. This study bridges the gap between the understanding of clinical and microbiological vancomycin resistance in S. aureus and expands the number of allelic variants (18 ± 4 mutations for the top 5 mutated genes) that result in vancomycin resistance phenotypes.

Identifying the effect of vancomycin on health care-associated methicillin-resistant Staphylococcus aureus strains using bacteriological and physiological media.

BACKGROUND: The evolving antibiotic-resistant behavior of health care-associated methicillin-resistant Staphylococcus aureus (HA-MRSA) USA100 strains are of major concern. They are resistant to a broad class of antibiotics such as macrolides, aminoglycosides, fluoroquinolones, and many more. FINDINGS: The selection of appropriate antibiotic susceptibility examination media is very important. Thus, we use bacteriological (cation-adjusted Mueller-Hinton broth) as well as physiological (R10LB) media to determine the effect of vancomycin on USA100 strains. The study includes the profiling behavior of HA-MRSA USA100 D592 and D712 strains in the presence of vancomycin through various high-throughput assays. The US100 D592 and D712 strains were characterized at sub-inhibitory concentrations through growth curves, RNA sequencing, bacterial cytological profiling, and exo-metabolomics high throughput experiments. CONCLUSIONS: The study reveals the vancomycin resistance behavior of HA-MRSA USA100 strains in dual media conditions using wide-ranging experiments.

A workflow for generating multi-strain genome-scale metabolic models of prokaryotes.

Genome-scale models (GEMs) of bacterial strains' metabolism have been formulated and used over the past 20 years. Recently, with the number of genome sequences exponentially increasing, multi-strain GEMs have proved valuable to define the properties of a species. Here, through four major stages, we extend the original Protocol used to generate a GEM for a single strain to enable multi-strain GEMs: (i) obtain or generate a high-quality model of a reference strain; (ii) compare the genome sequence between a reference strain and target strains to generate a homology matrix; (iii) generate draft strain-specific models from the homology matrix; and (iv) manually curate draft models. These multi-strain GEMs can be used to study pan-metabolic capabilities and strain-specific differences across a species, thus providing insights into its range of lifestyles. Unlike the original Protocol, this procedure is scalable and can be partly automated with the Supplementary Jupyter notebook Tutorial. This Protocol Extension joins the ranks of other comparable methods for generating models such as CarveMe and KBase. This extension of the original Protocol takes on the order of weeks to multiple months to complete depending on the availability of a suitable reference model.

Elucidation of Regulatory Modes for Five Two-Component Systems in Escherichia coli Reveals Novel Relationships.

Escherichia coli uses two-component systems (TCSs) to respond to environmental signals. TCSs affect gene expression and are parts of E. coli's global transcriptional regulatory network (TRN). Here, we identified the regulons of five TCSs in E. coli MG1655: BaeSR and CpxAR, which were stimulated by ethanol stress; KdpDE and PhoRB, induced by limiting potassium and phosphate, respectively; and ZraSR, stimulated by zinc. We analyzed RNA-seq data using independent component analysis (ICA). ChIP-exo data were used to validate condition-specific target gene binding sites. Based on these data, we do the following: (i) identify the target genes for each TCS; (ii) show how the target genes are transcribed in response to stimulus; and (iii) reveal novel relationships between TCSs, which indicate noncognate inducers for various response regulators, such as BaeR to iron starvation, CpxR to phosphate limitation, and PhoB and ZraR to cell envelope stress. Our understanding of the TRN in E. coli is thus notably expanded.IMPORTANCE E. coli is a common commensal microbe found in the human gut microenvironment; however, some strains cause diseases like diarrhea, urinary tract infections, and meningitis. E. coli's two-component systems (TCSs) modulate target gene expression, especially related to virulence, pathogenesis, and antimicrobial peptides, in response to environmental stimuli. Thus, it is of utmost importance to understand the transcriptional regulation of TCSs to infer bacterial environmental adaptation and disease pathogenicity. Utilizing a combinatorial approach integrating RNA sequencing (RNA-seq), independent component analysis, chromatin immunoprecipitation coupled with exonuclease treatment (ChIP-exo), and data mining, we suggest five different modes of TCS transcriptional regulation. Our data further highlight noncognate inducers of TCSs, which emphasizes the cross-regulatory nature of TCSs in E. coli and suggests that TCSs may have a role beyond their cognate functionalities. In summary, these results can lead to an understanding of the metabolic capabilities of bacteria and correctly predict complex phenotype under diverse conditions, especially when further incorporated with genome-scale metabolic models.

Systems Biology and Pangenome of Salmonella O-Antigens.

O-antigens are glycopolymers in lipopolysaccharides expressed on the cell surface of Gram-negative bacteria. Variability in the O-antigen structure constitutes the basis for the establishment of the serotyping schema. We pursued a two-pronged approach to define the basis for O-antigen structural diversity. First, we developed a bottom-up systems biology approach to O-antigen metabolism by building a reconstruction of Salmonella O-antigen biosynthesis and used it to (i) update 410 existing Salmonella strain-specific metabolic models, (ii) predict a strain's serogroup and its O-antigen glycan synthesis capability (yielding 98% agreement with experimental data), and (iii) extend our workflow to more than 1,400 Gram-negative strains. Second, we used a top-down pangenome analysis to elucidate the genetic basis for intraserogroup O-antigen structural variations. We assembled a database of O-antigen gene islands from over 11,000 sequenced Salmonella strains, revealing (i) that gene duplication, pseudogene formation, gene deletion, and bacteriophage insertion elements occur ubiquitously across serogroups; (ii) novel serotypes in the group O:4 B2 variant, as well as an additional genotype variant for group O:4, and (iii) two novel O-antigen gene islands in understudied subspecies. We thus comprehensively defined the genetic basis for O-antigen diversity.IMPORTANCE Lipopolysaccharides are a major component of the outer membrane in Gram-negative bacteria. They are composed of a conserved lipid structure that is embedded in the outer leaflet of the outer membrane and a polysaccharide known as the O-antigen. O-antigens are highly variable in structure across strains of a species and are crucial to a bacterium's interactions with its environment. They constitute the first line of defense against both the immune system and bacteriophage infections and have been shown to mediate antimicrobial resistance. The significance of our research is in identifying the metabolic and genetic differences within and across O-antigen groups in Salmonella strains. Our effort constitutes a first step toward characterizing the O-antigen metabolic network across Gram-negative organisms and a comprehensive overview of genetic variations in Salmonella.