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.

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.

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.

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.

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.

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.

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.

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.

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.

The Spider Exposure Pathway and the Potential Risk to Arachnivorous Birds.

There is growing concern over the health of North American birds, with evidence suggesting substantial population declines. Spiders are prominent dietary items for many bird species and mediate the transfer of contaminants to arachnivorous birds that consume them. Few studies have investigated the potential risk the spider exposure pathway poses to these birds because most studies have focused on piscivores. In the present study, we developed new chronic and acute As, Cd, Cu, Pb, Ni, Se, Zn, and MeHg spider-based avian wildlife values (SBAWVs) for multiple adult and nestling birds (primarily passerines) and then used the newly generated SBAWVs to characterize the risk to birds across 2 study areas: 1) 5 reaches in the southern Appalachian Mountains, an area with substantial mercury deposition but minimal anthropogenic impact, and 2) 4 reaches adjacent to the Emory River, an area impacted by the largest fly coal-ash spill in US history. We identified MeHg and Cu, Pb, Se, and Zn as contaminants of potential concern (COPC) at the Appalachian Mountain and Emory River study areas, respectively, based on dietary exposure of aquatic contaminants via riparian spiders. The identification of COPC at both study areas due to dietary spider exposure is notable not only because the spider exposure pathway has largely been uninvestigated at these sites but also because the aquatic systems in both areas have been studied extensively. Significant differences in MeHg concentrations were detected among spider taxa and suggest that the selection of spider taxa can impact risk characterization. These results indicate that the spider exposure pathway is important to consider when assessing potential risk, particularly for passerine birds. Environ Toxicol Chem 2020;39:2314-2324. © 2020 SETAC.

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.

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.

Profiling the effect of nafcillin on HA-MRSA D712 using bacteriological and physiological media.

Staphylococcus aureus strains have been continuously evolving resistance to numerous classes of antibiotics including methicillin, vancomycin, daptomycin and linezolid, compounding the enormous healthcare and economic burden of the pathogen. Cation-adjusted Mueller-Hinton broth (CA-MHB) is the standard bacteriological media for measuring antibiotic susceptibility in the clinical lab, but the use of media that more closely mimic the physiological state of the patient, e.g. mammalian tissue culture media, can in certain circumstances reveal antibiotic activities that may be more predictive of effectiveness in vivo. In the current study, we use both types of media to explore antibiotic resistance phenomena in hospital-acquired USA100 lineage methicillin-resistant, vancomycin-intermediate Staphylococcus aureus (MRSA/VISA) strain D712 via multidimensional high throughput analysis of growth rates, bacterial cytological profiling, RNA sequencing, and exo-metabolomics (HPLC and LC-MS). Here, we share data generated from these assays to shed light on the antibiotic resistance behavior of MRSA/VISA D712 in both bacteriological and physiological media.

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.

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.

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

Characterization of CA-MRSA TCH1516 exposed to nafcillin in bacteriological and physiological media.

Cation adjusted-Mueller Hinton Broth (CA-MHB) is the standard bacteriological medium utilized in the clinic for the determination of antibiotic susceptibility. However, a growing number of literature has demonstrated that media conditions can cause a substantial difference in the efficacy of antibiotics and antimicrobials. Recent studies have also shown that minimum inhibitory concentration (MIC) tests performed in standard cell culture media (e.g. RPMI and DMEM) are more indicative of in vivo antibiotic efficacy, presumably because they are a better proxy for the human host's physiological conditions. The basis for the bacterial media dependent susceptibility to antibiotics remains undefined. To address this question, we characterized the physiological response of methicillin-resistant Staphylococcus aureus (MRSA) during exposure to sub-inhibitory concentrations of the beta-lactam antibiotic nafcillin in either CA-MHB or RPMI + 10% LB (R10LB). Here, we present high quality transcriptomic, exo-metabolomic and morphological data paired with growth and susceptibility results for MRSA cultured in either standard bacteriologic or more physiologic relevant medium.