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.

Bacterial fitness landscapes stratify based on proteome allocation associated with discrete aero-types.

The fitness landscape is a concept commonly used to describe evolution towards optimal phenotypes. It can be reduced to mechanistic detail using genome-scale models (GEMs) from systems biology. We use recently developed GEMs of Metabolism and protein Expression (ME-models) to study the distribution of Escherichia coli phenotypes on the rate-yield plane. We found that the measured phenotypes distribute non-uniformly to form a highly stratified fitness landscape. Systems analysis of the ME-model simulations suggest that this stratification results from discrete ATP generation strategies. Accordingly, we define "aero-types", a phenotypic trait that characterizes how a balanced proteome can achieve a given growth rate by modulating 1) the relative utilization of oxidative phosphorylation, glycolysis, and fermentation pathways; and 2) the differential employment of electron-transport-chain enzymes. This global, quantitative, and mechanistic systems biology interpretation of fitness landscape formed upon proteome allocation offers a fundamental understanding of bacterial physiology and evolution dynamics.

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.

Community formation in wealth-mediated thermodynamic strategy evolution.

We study a dynamical system defined by a repeated game on a 1D lattice, in which the players keep track of their gross payoffs over time in a bank. Strategy updates are governed by a Boltzmann distribution, which depends on the neighborhood bank values associated with each strategy, relative to a temperature scale, which defines the random fluctuations. Players with higher bank values are, thus, less likely to change strategy than players with a lower bank value. For a parameterized rock-paper-scissors game, we derive a condition under which communities of a given strategy form with either fixed or drifting boundaries. We show the effect of a temperature increase on the underlying system and identify surprising properties of this model through numerical simulations.

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.

Adaptive laboratory evolution of Escherichia coli under acid stress.

The ability of Escherichia coli to tolerate acid stress is important for its survival and colonization in the human digestive tract. Here, we performed adaptive laboratory evolution of the laboratory strain E. coli K-12 MG1655 at pH 5.5 in glucose minimal medium. After 800 generations, six independent populations under evolution had reached 18.0 % higher growth rates than their starting strain at pH 5.5, while maintaining comparable growth rates to the starting strain at pH 7. We characterized the evolved strains and found that: (1) whole genome sequencing of isolated clones from each evolved population revealed mutations in rpoC appearing in five of six sequenced clones; and (2) gene expression profiles revealed different strategies to mitigate acid stress, which are related to amino acid metabolism and energy production and conversion. Thus, a combination of adaptive laboratory evolution, genome resequencing and expression profiling revealed, on a genome scale, the strategies that E. coli uses to mitigate acid stress.

Predicting proteome allocation, overflow metabolism, and metal requirements in a model acetogen.

The unique capability of acetogens to ferment a broad range of substrates renders them ideal candidates for the biotechnological production of commodity chemicals. In particular the ability to grow with H2:CO2 or syngas (a mixture of H2/CO/CO2) makes these microorganisms ideal chassis for sustainable bioproduction. However, advanced design strategies for acetogens are currently hampered by incomplete knowledge about their physiology and our inability to accurately predict phenotypes. Here we describe the reconstruction of a novel genome-scale model of metabolism and macromolecular synthesis (ME-model) to gain new insights into the biology of the model acetogen Clostridium ljungdahlii. The model represents the first ME-model of a Gram-positive bacterium and captures all major central metabolic, amino acid, nucleotide, lipid, major cofactors, and vitamin synthesis pathways as well as pathways to synthesis RNA and protein molecules necessary to catalyze these reactions, thus significantly broadens the scope and predictability. Use of the model revealed how protein allocation and media composition influence metabolic pathways and energy conservation in acetogens and accurately predicted secretion of multiple fermentation products. Predicting overflow metabolism is of particular interest since it enables new design strategies, e.g. the formation of glycerol, a novel product for C. ljungdahlii, thus broadening the metabolic capability for this model microbe. Furthermore, prediction and experimental validation of changing secretion rates based on different metal availability opens the window into fermentation optimization and provides new knowledge about the proteome utilization and carbon flux in acetogens.

Laboratory evolution reveals a two-dimensional rate-yield tradeoff in microbial metabolism.

Growth rate and yield are fundamental features of microbial growth. However, we lack a mechanistic and quantitative understanding of the rate-yield relationship. Studies pairing computational predictions with experiments have shown the importance of maintenance energy and proteome allocation in explaining rate-yield tradeoffs and overflow metabolism. Recently, adaptive evolution experiments of Escherichia coli reveal a phenotypic diversity beyond what has been explained using simple models of growth rate versus yield. Here, we identify a two-dimensional rate-yield tradeoff in adapted E. coli strains where the dimensions are (A) a tradeoff between growth rate and yield and (B) a tradeoff between substrate (glucose) uptake rate and growth yield. We employ a multi-scale modeling approach, combining a previously reported coarse-grained small-scale proteome allocation model with a fine-grained genome-scale model of metabolism and gene expression (ME-model), to develop a quantitative description of the full rate-yield relationship for E. coli K-12 MG1655. The multi-scale analysis resolves the complexity of ME-model which hindered its practical use in proteome complexity analysis, and provides a mechanistic explanation of the two-dimensional tradeoff. Further, the analysis identifies modifications to the P/O ratio and the flux allocation between glycolysis and pentose phosphate pathway (PPP) as potential mechanisms that enable the tradeoff between glucose uptake rate and growth yield. Thus, the rate-yield tradeoffs that govern microbial adaptation to new environments are more complex than previously reported, and they can be understood in mechanistic detail using a multi-scale modeling approach.

The genetic basis for adaptation of model-designed syntrophic co-cultures.

Understanding the fundamental characteristics of microbial communities could have far reaching implications for human health and applied biotechnology. Despite this, much is still unknown regarding the genetic basis and evolutionary strategies underlying the formation of viable synthetic communities. By pairing auxotrophic mutants in co-culture, it has been demonstrated that viable nascent E. coli communities can be established where the mutant strains are metabolically coupled. A novel algorithm, OptAux, was constructed to design 61 unique multi-knockout E. coli auxotrophic strains that require significant metabolite uptake to grow. These predicted knockouts included a diverse set of novel non-specific auxotrophs that result from inhibition of major biosynthetic subsystems. Three OptAux predicted non-specific auxotrophic strains-with diverse metabolic deficiencies-were co-cultured with an L-histidine auxotroph and optimized via adaptive laboratory evolution (ALE). Time-course sequencing revealed the genetic changes employed by each strain to achieve higher community growth rates and provided insight into mechanisms for adapting to the syntrophic niche. A community model of metabolism and gene expression was utilized to predict the relative community composition and fundamental characteristics of the evolved communities. This work presents new insight into the genetic strategies underlying viable nascent community formation and a cutting-edge computational method to elucidate metabolic changes that empower the creation of cooperative communities.

Identifying contaminants of potential concern in remote headwater streams of Tennessee's Appalachian Mountains.

The susceptibility of Tennessee's Appalachian Mountains to anthropogenic stressors has remained largely uninvestigated likely due to a lack of known point source contamination. However, a growing body of scientific evidence suggests that depositional inputs can lead to concerning levels of contamination, even in remote areas. To investigate potential concerns, water quality parameters, contaminants in water (nitrogen, TSS, and metals), and contaminants in eastern brook trout (mercury, polychlorinated biphenyls [PCBs], organochlorine [OC] pesticides, dioxins, furans, and phthalates) were measured in four Appalachian Mountain streams from 2015 to 2017. Concentrations were compared to literature and/or model-derived (e.g., biotic ligand model) threshold values to determine whether levels exceeded those acceptable for stream health. Dioxins and furans were detectable in fish tissue at all sites with an average 2,3,7,8-tetrachlorodinbenzodioxin toxicity equivalence (TEQ) of 0.0015 ng/kg. Concentrations of PCBs, phthalates, and organochlorine pesticides were never above analytical quantitation limits, although several OC pesticides (e.g., alpha-chlordane) were detectable in fish. Aluminum concentrations in water were found at levels shown previously to cause mortality in brook trout during acidic rain events. The average whole-body methylmercury concentrations in fish among sites were 0.037 ± 0.003 µg/kg and were on average 75 ± 2% of total mercury.

The Role of Sexual Dimorphism and Tissue Selection in Ecotoxicological Studies Using the Riparian Spider Tetragnatha elongata.

Tetragnathid spiders (Tetragnatha spp.) found in riparian habitats have recently been used as bioindicators of sediment contamination and insect-mediated contaminant flux. We investigated whether sexual dimorphism (size and behavior) influenced the female:male ratio in composite samples, stable isotope ratios (carbon [δ13C], nitrogen [δ15N]), and Hg concentrations in the southern United States. Additionally, we explored whether biomass for contaminant analysis could be preserved by using the legs of tetragnathids as a surrogate for whole-body δ13C and δ15N signatures. We found that female tetragnathids were significantly larger than male spiders and represented a larger proportion of spiders collected at all sites. However, despite the difference in size between sexes, no differences in growth dynamics, isotopic signatures (δ13C and δ15N), or mercury concentrations were observed. It was determined that the leg of a tetragnathid can accurately represent the stable isotope signature of an entire spider.

Mercury Accumulation in Millipedes (Narceus spp.) Living Adjacent to a Southern Appalachian Mountain Stream (USA).

Millipedes are among the most important processors of leaf litter in temperate forests. Through consumption of leaf litter, millipedes may be exposed to mercury that accumulates in leaf tissues prior to senescence. To investigate mercury uptake in millipedes, Narceus spp. were collected from a remote site in the southern Appalachian Mountains, an area known to receive high mercury deposition. Additionally, aquatic primary consumers (larval caddisflies and stoneflies), brook trout (Salvelinus fontinalis) and rainbow trout (Oncorhynchus mykiss) were collected from the same site for comparisons of mercury concentrations and percent methylmercury. Bioaccumulation factors for millipedes were 18.5 and 20.2 for total and methylmercury, respectively. At this site, the mean THg concentration in millipedes was ~ 10 × greater than both brook trout and rainbow trout and ~ 200 × greater than that of aquatic primary consumers. Millipede THg concentrations ranged from 222 to 1620 ng/g ww in an area where EPA fish consumption criteria (300 ng/g MeHg in fish tissue, ww) were not exceeded. The mean percent methylmercury in millipedes was 1.4%, suggesting these animals were accumulating large quantities of inorganic mercury.

Advancing the scale of synthetic biology via cross-species transfer of cellular functions enabled by iModulon engraftment.

Machine learning applied to large compendia of transcriptomic data has enabled the decomposition of bacterial transcriptomes to identify independently modulated sets of genes, such iModulons represent specific cellular functions. The identification of iModulons enables accurate identification of genes necessary and sufficient for cross-species transfer of cellular functions. We demonstrate cross-species transfer of: 1) the biotransformation of vanillate to protocatechuate, 2) a malonate catabolic pathway, 3) a catabolic pathway for 2,3-butanediol, and 4) an antimicrobial resistance to ampicillin found in multiple Pseudomonas species to Escherichia coli. iModulon-based engineering is a transformative strategy as it includes all genes comprising the transferred cellular function, including genes without functional annotation. Adaptive laboratory evolution was deployed to optimize the cellular function transferred, revealing mutations in the host. Combining big data analytics and laboratory evolution thus enhances the level of understanding of systems biology, and synthetic biology for strain design and development.

Ecotoxicological Studies Indicate That Sublethal and Lethal Processes Limit Insect-Mediated Contaminant Flux.

Merolimnic insects can accumulate and transport considerable amounts of aquatic contaminants to terrestrial systems. The rate of contaminant biotransport, termed insect-mediated contaminant flux (IMCF), depends on emergent insect biomass and contaminant accumulation, both functions of environmental concentration. We developed a mathematical model of IMCF and apply it to three ecotoxicological studies obtained through the US Environmental Protection Agency's ECOTOX database to determine at which concentration maximum IMCF occurs. Model results demonstrate that the maximum IMCF depends on competing rates of biomass loss and contaminant accumulation and does not necessarily occur at the highest insect or environmental contaminant concentration. In addition, modeling results suggest that sublethal contaminant effects (e.g., decreased growth) on insect biomass can be an important and potentially underappreciated control on IMCF. Environ Toxicol Chem 2023;42:1982-1992. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Laboratory evolution, transcriptomics, and modeling reveal mechanisms of paraquat tolerance.

Relationships between the genome, transcriptome, and metabolome underlie all evolved phenotypes. However, it has proved difficult to elucidate these relationships because of the high number of variables measured. A recently developed data analytic method for characterizing the transcriptome can simplify interpretation by grouping genes into independently modulated sets (iModulons). Here, we demonstrate how iModulons reveal deep understanding of the effects of causal mutations and metabolic rewiring. We use adaptive laboratory evolution to generate E. coli strains that tolerate high levels of the redox cycling compound paraquat, which produces reactive oxygen species (ROS). We combine resequencing, iModulons, and metabolic models to elucidate six interacting stress-tolerance mechanisms: (1) modification of transport, (2) activation of ROS stress responses, (3) use of ROS-sensitive iron regulation, (4) motility, (5) broad transcriptional reallocation toward growth, and (6) metabolic rewiring to decrease NADH production. This work thus demonstrates the power of iModulon knowledge mapping for evolution analysis.

Laboratory evolution of synthetic electron transport system variants reveals a larger metabolic respiratory system and its plasticity.

The bacterial respiratory electron transport system (ETS) is branched to allow condition-specific modulation of energy metabolism. There is a detailed understanding of the structural and biochemical features of respiratory enzymes; however, a holistic examination of the system and its plasticity is lacking. Here we generate four strains of Escherichia coli harboring unbranched ETS that pump 1, 2, 3, or 4 proton(s) per electron and characterized them using a combination of synergistic methods (adaptive laboratory evolution, multi-omic analyses, and computation of proteome allocation). We report that: (a) all four ETS variants evolve to a similar optimized growth rate, and (b) the laboratory evolutions generate specific rewiring of major energy-generating pathways, coupled to the ETS, to optimize ATP production capability. We thus define an Aero-Type System (ATS), which is a generalization of the aerobic bioenergetics and is a metabolic systems biology description of respiration and its inherent plasticity.

Use of Riparian Spiders as Sentinels of Persistent and Bioavailable Chemical Contaminants in Aquatic Ecosystems: A Review.

Aquatic ecosystems around the world are contaminated with a wide range of anthropogenic chemicals, including metals and organic pollutants, that originate from point and nonpoint sources. Many of these chemical contaminants have complex environmental cycles, are persistent and bioavailable, can be incorporated into aquatic food webs, and pose a threat to the health of wildlife and humans. Identifying appropriate sentinels that reflect bioavailability is critical to assessing and managing aquatic ecosystems impacted by contaminants. The objective of the present study is to review research on riparian spiders as sentinels of persistent and bioavailable chemical contaminants in aquatic ecosystems. Our review of the literature on riparian spiders as sentinels suggests that significant progress has been made during the last two decades of research. We identified 55 published studies conducted around the world in which riparian spiders (primarily of the families Tetragnathidae, Araneidae, Lycosidae, and Pisauridae) were used as sentinels of chemical contamination of lotic, lentic, and estuarine systems. For several contaminants, such as polychlorinated biphenyls (PCBs), Hg, and Se, it is now clear that riparian spiders are appropriate sentinels. However, many contaminants and factors that could impact chemical concentrations in riparian spiders have not been well characterized. Further study of riparian spiders and their potential role as sentinels is critical because it would allow for development of national-scale programs that utilize riparian spiders as sentinels to monitor chemical contaminants in aquatic ecosystems. A riparian spider sentinel program in the United States would be complementary to existing national sentinel programs, including those for fish and immature dragonflies. Environ Toxicol Chem 2022;41:499-514. © 2021 SETAC.

Restoration of fitness lost due to dysregulation of the pyruvate dehydrogenase complex is triggered by ribosomal binding site modifications.

Pyruvate dehydrogenase complex (PDC) functions as the main determinant of the respiro-fermentative balance because it converts pyruvate to acetyl-coenzyme A (CoA), which then enters the TCA (tricarboxylic acid cycle). PDC is repressed by the pyruvate dehydrogenase complex regulator (PdhR) in Escherichia coli. The deletion of the pdhR gene compromises fitness in aerobic environments. We evolve the E. coli pdhR deletion strain to examine its achievable growth rate and the underlying adaptive strategies. We find that (1) optimal proteome allocation to PDC is critical in achieving optimal growth rate; (2) expression of PDC in evolved strains is reduced through mutations in the Shine-Dalgarno sequence; (3) rewiring of the TCA flux and increased reactive oxygen species (ROS) defense occur in the evolved strains; and (4) the evolved strains adapt to an efficient biomass yield. Together, these results show how adaptation can find alternative regulatory mechanisms for a key cellular process if the primary regulatory mode fails.

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.