Genomic surveillance as a scalable framework for precision phage therapy against antibiotic-resistant pathogens.

Phage therapy is gaining increasing interest in the fight against critically antibiotic-resistant nosocomial pathogens. However, the narrow host range of bacteriophages hampers the development of broadly effective phage therapeutics and demands precision approaches. Here, we combine large-scale phylogeographic analysis with high-throughput phage typing to guide the development of precision phage cocktails targeting carbapenem-resistant Acinetobacter baumannii, a top-priority pathogen. Our analysis reveals that a few strain types dominate infections in each world region, with their geographical distribution remaining stable within 6 years. As we demonstrate in Eastern Europe, this spatiotemporal distribution enables preemptive preparation of region-specific phage collections that target most local infections. Finally, we showcase the efficacy of phage cocktails against prevalent strain types using in vitro and animal infection models. Ultimately, genomic surveillance identifies patients benefiting from the same phages across geographical scales, thus providing a scalable framework for precision phage therapy.

Cell-cell metabolite exchange creates a pro-survival metabolic environment that extends lifespan.

Metabolism is deeply intertwined with aging. Effects of metabolic interventions on aging have been explained with intracellular metabolism, growth control, and signaling. Studying chronological aging in yeast, we reveal a so far overlooked metabolic property that influences aging via the exchange of metabolites. We observed that metabolites exported by young cells are re-imported by chronologically aging cells, resulting in cross-generational metabolic interactions. Then, we used self-establishing metabolically cooperating communities (SeMeCo) as a tool to increase metabolite exchange and observed significant lifespan extensions. The longevity of the SeMeCo was attributable to metabolic reconfigurations in methionine consumer cells. These obtained a more glycolytic metabolism and increased the export of protective metabolites that in turn extended the lifespan of cells that supplied them with methionine. Our results establish metabolite exchange interactions as a determinant of cellular aging and show that metabolically cooperating cells can shape the metabolic environment to extend their lifespan.

Widespread natural selection on metabolite levels in humans.

Natural selection acts ubiquitously on complex human traits, predominantly constraining the occurrence of extreme phenotypes (stabilizing selection). These constraints propagate to DNA sequence variants associated with traits under selection. The genetic signatures of such evolutionary events can thus be detected via combining effect size estimates from genetic association studies and the corresponding allele frequencies. Although this approach has been successfully applied to high-level traits, the prevalence and mode of selection acting on molecular traits remain poorly understood. Here, we estimate the action of natural selection on genetic variants associated with metabolite levels, an important layer of molecular traits. By leveraging summary statistics of published genome-wide association studies with large sample sizes, we find strong evidence of stabilizing selection for 15 out of 97 plasma metabolites, with nonessential amino acids displaying especially strong selection signatures. Mendelian randomization analysis reveals that metabolites under stronger stabilizing selection display larger effects on a range of clinically relevant complex traits, suggesting that maintaining a disease-free profile may be an important source of selective constraints on the metabolome. Metabolites under strong stabilizing selection in humans are also more conserved in their concentrations among diverse mammalian species, suggesting shared selective forces across micro- and macroevolutionary timescales. Overall, this study demonstrates that variation in metabolite levels among humans is frequently shaped by natural selection and this may act through their causal impact on disease susceptibility.

The metabolic domestication syndrome of budding yeast.

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

Principles of metabolome conservation in animals.

Metabolite levels shape cellular physiology and disease susceptibility, yet the general principles governing metabolome evolution are largely unknown. Here, we introduce a measure of conservation of individual metabolite levels among related species. By analyzing multispecies tissue metabolome datasets in phylogenetically diverse mammals and fruit flies, we show that conservation varies extensively across metabolites. Three major functional properties, metabolite abundance, essentiality, and association with human diseases predict conservation, highlighting a striking parallel between the evolutionary forces driving metabolome and protein sequence conservation. Metabolic network simulations recapitulated these general patterns and revealed that abundant metabolites are highly conserved due to their strong coupling to key metabolic fluxes in the network. Finally, we show that biomarkers of metabolic diseases can be distinguished from other metabolites simply based on evolutionary conservation, without requiring any prior clinical knowledge. Overall, this study uncovers simple rules that govern metabolic evolution in animals and implies that most tissue metabolome differences between species are permitted, rather than favored by natural selection. More broadly, our work paves the way toward using evolutionary information to identify biomarkers, as well as to detect pathogenic metabolome alterations in individual patients.

Phase-separated ribosome-nascent chain complexes in genotoxic stress response.

Assemblysomes are EDTA- and RNase-resistant ribonucleoprotein (RNP) complexes of paused ribosomes with protruding nascent polypeptide chains. They have been described in yeast and human cells for the proteasome subunit Rpt1, and the disordered amino-terminal part of the nascent chain was found to be indispensable for the accumulation of the Rpt1-RNP into assemblysomes. Motivated by this, to find other assemblysome-associated RNPs we used bioinformatics to rank subunits of Saccharomyces cerevisiae protein complexes according to their amino-terminal disorder propensity. The results revealed that gene products involved in DNA repair are enriched among the top candidates. The Sgs1 DNA helicase was chosen for experimental validation. We found that indeed nascent chains of Sgs1 form EDTA-resistant RNP condensates, assemblysomes by definition. Moreover, upon exposure to UV, SGS1 mRNA shifted from assemblysomes to polysomes, suggesting that external stimuli are regulators of assemblysome dynamics. We extended our studies to human cell lines. The BLM helicase, ortholog of yeast Sgs1, was identified upon sequencing assemblysome-associated RNAs from the MCF7 human breast cancer cell line, and mRNAs encoding DNA repair proteins were overall enriched. Using the radiation-resistant A549 cell line, we observed by transmission electron microscopy that 1,6-hexanediol, an agent known to disrupt phase-separated condensates, depletes ring ribosome structures compatible with assemblysomes from the cytoplasm of cells and makes the cells more sensitive to X-ray treatment. Taken together, these findings suggest that assemblysomes may be a component of the DNA damage response from yeast to human.

Chloroplastic ascorbate modifies plant metabolism and may act as a metabolite signal regardless of oxidative stress.

Ascorbate (Asc) is a major plant metabolite that plays crucial roles in various processes, from reactive oxygen scavenging to epigenetic regulation. However, to what extent and how Asc modulates metabolism is largely unknown. We investigated the consequences of chloroplastic and total cellular Asc deficiencies by studying chloroplastic Asc transporter mutant lines lacking PHOSPHATE TRANSPORTER 4; 4 and the Asc-deficient vtc2-4 mutant of Arabidopsis (Arabidopsis thaliana). Under regular growth conditions, both Asc deficiencies caused minor alterations in photosynthesis, with no apparent signs of oxidative damage. In contrast, metabolomics analysis revealed global and largely overlapping alterations in the metabolome profiles of both Asc-deficient mutants, suggesting that chloroplastic Asc modulates plant metabolism. We observed significant alterations in amino acid metabolism, particularly in arginine metabolism, activation of nucleotide salvage pathways, and changes in secondary metabolism. In addition, proteome-wide analysis of thermostability revealed that Asc may interact with enzymes involved in arginine metabolism, the Calvin-Benson cycle, and several photosynthetic electron transport components. Overall, our results suggest that, independent of oxidative stress, chloroplastic Asc modulates the activity of diverse metabolic pathways in vascular plants and may act as an internal metabolite signal.

mulea: An R package for enrichment analysis using multiple ontologies and empirical false discovery rate.

Traditional gene set enrichment analyses are typically limited to a few ontologies and do not account for the interdependence of gene sets or terms, resulting in overcorrected p-values. To address these challenges, we introduce mulea, an R package offering comprehensive overrepresentation and functional enrichment analysis. mulea employs a progressive empirical false discovery rate (eFDR) method, specifically designed for interconnected biological data, to accurately identify significant terms within diverse ontologies. mulea expands beyond traditional tools by incorporating a wide range of ontologies, encompassing Gene Ontology, pathways, regulatory elements, genomic locations, and protein domains. This flexibility enables researchers to tailor enrichment analysis to their specific questions, such as identifying enriched transcriptional regulators in gene expression data or overrepresented protein domains in protein sets. To facilitate seamless analysis, mulea provides gene sets (in standardised GMT format) for 27 model organisms, covering 22 ontology types from 16 databases and various identifiers resulting in almost 900 files. Additionally, the muleaData ExperimentData Bioconductor package simplifies access to these pre-defined ontologies. Finally, mulea's architecture allows for easy integration of user-defined ontologies, or GMT files from external sources (e.g., MSigDB or Enrichr), expanding its applicability across diverse research areas. mulea is distributed as a CRAN R package downloadable from https://cran.r-project.org/web/packages/mulea/ and https://github.com/ELTEbioinformatics/mulea . It offers researchers a powerful and flexible toolkit for functional enrichment analysis, addressing limitations of traditional tools with its progressive eFDR and by supporting a variety of ontologies. Overall, mulea fosters the exploration of diverse biological questions across various model organisms.

Plasticity and Stereotypic Rewiring of the Transcriptome Upon Bacterial Evolution of Antibiotic Resistance.

Bacterial evolution of antibiotic resistance frequently has deleterious side effects on microbial growth, virulence, and susceptibility to other antimicrobial agents. However, it is unclear how these trade-offs could be utilized for manipulating antibiotic resistance in the clinic, not least because the underlying molecular mechanisms are poorly understood. Using laboratory evolution, we demonstrate that clinically relevant resistance mutations in Escherichia coli constitutively rewire a large fraction of the transcriptome in a repeatable and stereotypic manner. Strikingly, lineages adapted to functionally distinct antibiotics and having no resistance mutations in common show a wide range of parallel gene expression changes that alter oxidative stress response, iron homeostasis, and the composition of the bacterial outer membrane and cell surface. These common physiological alterations are associated with changes in cell morphology and enhanced sensitivity to antimicrobial peptides. Finally, the constitutive transcriptomic changes induced by resistance mutations are largely distinct from those induced by antibiotic stresses in the wild type. This indicates a limited role for genetic assimilation of the induced antibiotic stress response during resistance evolution. Our work suggests that diverse resistance mutations converge on similar global transcriptomic states that shape genetic susceptibility to antimicrobial compounds.

Underground metabolism as a rich reservoir for pathway engineering.

MOTIVATION: Bioproduction of value-added compounds is frequently achieved by utilizing enzymes from other species. However, expression of such heterologous enzymes can be detrimental due to unexpected interactions within the host cell. Recently, an alternative strategy emerged, which relies on recruiting side activities of host enzymes to establish new biosynthetic pathways. Although such low-level 'underground' enzyme activities are prevalent, it remains poorly explored whether they may serve as an important reservoir for pathway engineering. RESULTS: Here, we use genome-scale modeling to estimate the theoretical potential of underground reactions for engineering novel biosynthetic pathways in Escherichia coli. We found that biochemical reactions contributed by underground enzyme activities often enhance the in silico production of compounds with industrial importance, including several cases where underground activities are indispensable for production. Most of these new capabilities can be achieved by the addition of one or two underground reactions to the native network, suggesting that only a few side activities need to be enhanced during implementation. Remarkably, we find that the contribution of underground reactions to the production of value-added compounds is comparable to that of heterologous reactions, underscoring their biotechnological potential. Taken together, our genome-wide study demonstrates that exploiting underground enzyme activities could be a promising addition to the toolbox of industrial strain development. AVAILABILITY AND IMPLEMENTATION: The data and scripts underlying this article are available on GitHub at https://github.com/pappb/Kovacs-et-al-Underground-metabolism. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Suboptimal Global Transcriptional Response Increases the Harmful Effects of Loss-of-Function Mutations.

The fitness impact of loss-of-function mutations is generally assumed to reflect the loss of specific molecular functions associated with the perturbed gene. Here, we propose that rewiring of the transcriptome upon deleterious gene inactivation is frequently nonspecific and mimics stereotypic responses to external environmental change. Consequently, transcriptional response to gene deletion could be suboptimal and incur an extra fitness cost. Analysis of the transcriptomes of ∼1,500 single-gene deletion Saccharomyces cerevisiae strains supported this scenario. First, most transcriptomic changes are not specific to the deleted gene but are rather triggered by perturbations in functionally diverse genes. Second, gene deletions that alter the expression of dosage-sensitive genes are especially harmful. Third, by elevating the expression level of downregulated genes, we could experimentally mitigate the fitness defect of gene deletions. Our work shows that rewiring of genomic expression upon gene inactivation shapes the harmful effects of mutations.

Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations.

Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology.

Unmatched Level of Molecular Convergence among Deeply Divergent Complex Multicellular Fungi.

Convergent evolution is pervasive in nature, but it is poorly understood how various constraints and natural selection limit the diversity of evolvable phenotypes. Here, we analyze the transcriptome across fruiting body development to understand the independent evolution of complex multicellularity in the two largest clades of fungi-the Agarico- and Pezizomycotina. Despite >650 My of divergence between these clades, we find that very similar sets of genes have convergently been co-opted for complex multicellularity, followed by expansions of their gene families by duplications. Over 82% of shared multicellularity-related gene families were expanding in both clades, indicating a high prevalence of convergence also at the gene family level. This convergence is coupled with a rich inferred repertoire of multicellularity-related genes in the most recent common ancestor of the Agarico- and Pezizomycotina, consistent with the hypothesis that the coding capacity of ancestral fungal genomes might have promoted the repeated evolution of complex multicellularity. We interpret this repertoire as an indication of evolutionary predisposition of fungal ancestors for evolving complex multicellular fruiting bodies. Our work suggests that evolutionary convergence may happen not only when organisms are closely related or are under similar selection pressures, but also when ancestral genomic repertoires render certain evolutionary trajectories more likely than others, even across large phylogenetic distances.

Directed evolution of multiple genomic loci allows the prediction of antibiotic resistance.

Antibiotic development is frequently plagued by the rapid emergence of drug resistance. However, assessing the risk of resistance development in the preclinical stage is difficult. Standard laboratory evolution approaches explore only a small fraction of the sequence space and fail to identify exceedingly rare resistance mutations and combinations thereof. Therefore, new rapid and exhaustive methods are needed to accurately assess the potential of resistance evolution and uncover the underlying mutational mechanisms. Here, we introduce directed evolution with random genomic mutations (DIvERGE), a method that allows an up to million-fold increase in mutation rate along the full lengths of multiple predefined loci in a range of bacterial species. In a single day, DIvERGE generated specific mutation combinations, yielding clinically significant resistance against trimethoprim and ciprofloxacin. Many of these mutations have remained previously undetected or provide resistance in a species-specific manner. These results indicate pathogen-specific resistance mechanisms and the necessity of future narrow-spectrum antibacterial treatments. In contrast to prior claims, we detected the rapid emergence of resistance against gepotidacin, a novel antibiotic currently in clinical trials. Based on these properties, DIvERGE could be applicable to identify less resistance-prone antibiotics at an early stage of drug development. Finally, we discuss potential future applications of DIvERGE in synthetic and evolutionary biology.

Limited Evolutionary Conservation of the Phenotypic Effects of Antibiotic Resistance Mutations.

Multidrug-resistant clinical isolates are common in certain pathogens, but rare in others. This pattern may be due to the fact that mutations shaping resistance have species-specific effects. To investigate this issue, we transferred a range of resistance-conferring mutations and a full resistance gene into Escherichia coli and closely related bacteria. We found that resistance mutations in one bacterial species frequently provide no resistance, in fact even yielding drug hypersensitivity in close relatives. In depth analysis of a key gene involved in aminoglycoside resistance (trkH) indicated that preexisting mutations in other genes-intergenic epistasis-underlie such extreme differences in mutational effects between species. Finally, reconstruction of adaptive landscapes under multiple antibiotic stresses revealed that mutations frequently provide multidrug resistance or elevated drug susceptibility (i.e., collateral sensitivity) only with certain combinations of other resistance mutations. We conclude that resistance and collateral sensitivity are contingent upon the genetic makeup of the bacterial population, and such contingency could shape the long-term fate of resistant bacteria. These results underlie the importance of species-specific treatment strategies.

Enzyme promiscuity shapes adaptation to novel growth substrates.

Evidence suggests that novel enzyme functions evolved from low-level promiscuous activities in ancestral enzymes. Yet, the evolutionary dynamics and physiological mechanisms of how such side activities contribute to systems-level adaptations are not well characterized. Furthermore, it remains untested whether knowledge of an organism's promiscuous reaction set, or underground metabolism, can aid in forecasting the genetic basis of metabolic adaptations. Here, we employ a computational model of underground metabolism and laboratory evolution experiments to examine the role of enzyme promiscuity in the acquisition and optimization of growth on predicted non-native substrates in Escherichia coli K-12 MG1655. After as few as approximately 20 generations, evolved populations repeatedly acquired the capacity to grow on five predicted non-native substrates-D-lyxose, D-2-deoxyribose, D-arabinose, m-tartrate, and monomethyl succinate. Altered promiscuous activities were shown to be directly involved in establishing high-efficiency pathways. Structural mutations shifted enzyme substrate turnover rates toward the new substrate while retaining a preference for the primary substrate. Finally, genes underlying the phenotypic innovations were accurately predicted by genome-scale model simulations of metabolism with enzyme promiscuity.

Correction: Phenotypic heterogeneity promotes adaptive evolution.

[This corrects the article DOI: 10.1371/journal.pbio.2000644.].

Phenotypic heterogeneity promotes adaptive evolution.

Genetically identical cells frequently display substantial heterogeneity in gene expression, cellular morphology and physiology. It has been suggested that by rapidly generating a subpopulation with novel phenotypic traits, phenotypic heterogeneity (or plasticity) accelerates the rate of adaptive evolution in populations facing extreme environmental challenges. This issue is important as cell-to-cell phenotypic heterogeneity may initiate key steps in microbial evolution of drug resistance and cancer progression. Here, we study how stochastic transitions between cellular states influence evolutionary adaptation to a stressful environment in yeast Saccharomyces cerevisiae. We developed inducible synthetic gene circuits that generate varying degrees of expression stochasticity of an antifungal resistance gene. We initiated laboratory evolutionary experiments with genotypes carrying different versions of the genetic circuit by exposing the corresponding populations to gradually increasing antifungal stress. Phenotypic heterogeneity altered the evolutionary dynamics by transforming the adaptive landscape that relates genotype to fitness. Specifically, it enhanced the adaptive value of beneficial mutations through synergism between cell-to-cell variability and genetic variation. Our work demonstrates that phenotypic heterogeneity is an evolving trait when populations face a chronic selection pressure. It shapes evolutionary trajectories at the genomic level and facilitates evolutionary rescue from a deteriorating environmental stress.

The dawn of evolutionary genome engineering.

Genome engineering strategies--such as genome editing, reduction and shuffling, and de novo genome synthesis--enable the modification of specific genomic locations in a directed and combinatorial manner. These approaches offer an unprecedented opportunity to study central evolutionary issues in which natural genetic variation is limited or biased, which sheds light on the evolutionary forces driving complex and extremely slowly evolving traits; the selective constraints on genome architecture; and the reconstruction of ancestral states of cellular structures and networks.