1
|
Ardell S, Martsul A, Johnson MS, Kryazhimskiy S. Environment-independent distribution of mutational effects emerges from microscopic epistasis. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.11.18.567655. [PMID: 38014325 PMCID: PMC10680819 DOI: 10.1101/2023.11.18.567655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Predicting how new mutations alter phenotypes is difficult because mutational effects vary across genotypes and environments. Recently discovered global epistasis, where the fitness effects of mutations scale with the fitness of the background genotype, can improve predictions, but how the environment modulates this scaling is unknown. We measured the fitness effects of ~100 insertion mutations in 42 strains of Saccharomyces cerevisiae in six laboratory environments and found that the global-epistasis scaling is nearly invariant across environments. Instead, the environment tunes one global parameter, the background fitness at which most mutations switch sign. As a consequence, the distribution of mutational effects is predictable across genotypes and environments. Our results suggest that the effective dimensionality of genotype-to-phenotype maps across environments is surprisingly low.
Collapse
Affiliation(s)
- Sarah Ardell
- Department of Ecology, Behavior and Evolution, University of California San Diego, La Jolla, CA 92093
| | - Alena Martsul
- Department of Ecology, Behavior and Evolution, University of California San Diego, La Jolla, CA 92093
| | - Milo S. Johnson
- Department of Integrative Biology, University of California Berkeley, Berkeley, CA 94720
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Sergey Kryazhimskiy
- Department of Ecology, Behavior and Evolution, University of California San Diego, La Jolla, CA 92093
| |
Collapse
|
2
|
Diaz-Colunga J, Skwara A, Vila JCC, Bajic D, Sanchez A. Global epistasis and the emergence of function in microbial consortia. Cell 2024; 187:3108-3119.e30. [PMID: 38776921 DOI: 10.1016/j.cell.2024.04.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Revised: 12/06/2023] [Accepted: 04/16/2024] [Indexed: 05/25/2024]
Abstract
The many functions of microbial communities emerge from a complex web of interactions between organisms and their environment. This poses a significant obstacle to engineering microbial consortia, hindering our ability to harness the potential of microorganisms for biotechnological applications. In this study, we demonstrate that the collective effect of ecological interactions between microbes in a community can be captured by simple statistical models that predict how adding a new species to a community will affect its function. These predictive models mirror the patterns of global epistasis reported in genetics, and they can be quantitatively interpreted in terms of pairwise interactions between community members. Our results illuminate an unexplored path to quantitatively predicting the function of microbial consortia from their composition, paving the way to optimizing desirable community properties and bringing the tasks of predicting biological function at the genetic, organismal, and ecological scales under the same quantitative formalism.
Collapse
Affiliation(s)
- Juan Diaz-Colunga
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT 06511, USA; Microbial Sciences Institute, Yale University, New Haven, CT 06511, USA; Department of Microbial Biotechnology, National Center for Biotechnology CNB-CSIC, 28049 Madrid, Spain; Institute of Functional Biology and Genomics IBFG-CSIC, University of Salamanca, 37007 Salamanca, Spain.
| | - Abigail Skwara
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT 06511, USA; Microbial Sciences Institute, Yale University, New Haven, CT 06511, USA
| | - Jean C C Vila
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT 06511, USA; Microbial Sciences Institute, Yale University, New Haven, CT 06511, USA; Department of Biology, Stanford University, Stanford, CA 94305, USA
| | - Djordje Bajic
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT 06511, USA; Microbial Sciences Institute, Yale University, New Haven, CT 06511, USA; Department of Biotechnology, Delft University of Technology, Delft 2628 CD, the Netherlands.
| | - Alvaro Sanchez
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT 06511, USA; Microbial Sciences Institute, Yale University, New Haven, CT 06511, USA; Department of Microbial Biotechnology, National Center for Biotechnology CNB-CSIC, 28049 Madrid, Spain; Institute of Functional Biology and Genomics IBFG-CSIC, University of Salamanca, 37007 Salamanca, Spain.
| |
Collapse
|
3
|
Martin NS, Schaper S, Camargo CQ, Louis AA. Non-Poissonian Bursts in the Arrival of Phenotypic Variation Can Strongly Affect the Dynamics of Adaptation. Mol Biol Evol 2024; 41:msae085. [PMID: 38693911 PMCID: PMC11156200 DOI: 10.1093/molbev/msae085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Revised: 03/01/2024] [Accepted: 04/17/2024] [Indexed: 05/03/2024] Open
Abstract
Modeling the rate at which adaptive phenotypes appear in a population is a key to predicting evolutionary processes. Given random mutations, should this rate be modeled by a simple Poisson process, or is a more complex dynamics needed? Here we use analytic calculations and simulations of evolving populations on explicit genotype-phenotype maps to show that the introduction of novel phenotypes can be "bursty" or overdispersed. In other words, a novel phenotype either appears multiple times in quick succession or not at all for many generations. These bursts are fundamentally caused by statistical fluctuations and other structure in the map from genotypes to phenotypes. Their strength depends on population parameters, being highest for "monomorphic" populations with low mutation rates. They can also be enhanced by additional inhomogeneities in the mapping from genotypes to phenotypes. We mainly investigate the effect of bursts using the well-studied genotype-phenotype map for RNA secondary structure, but find similar behavior in a lattice protein model and in Richard Dawkins's biomorphs model of morphological development. Bursts can profoundly affect adaptive dynamics. Most notably, they imply that fitness differences play a smaller role in determining which phenotype fixes than would be the case for a Poisson process without bursts.
Collapse
Affiliation(s)
- Nora S Martin
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Steffen Schaper
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Chico Q Camargo
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, UK
- Faculty of Environment, Science and Economy, University of Exeter, Exeter EX4 4QF, UK
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, UK
| |
Collapse
|
4
|
Gonzales J, Kim I, Hwang W, Cho JH. Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.24.595776. [PMID: 38826371 PMCID: PMC11142230 DOI: 10.1101/2024.05.24.595776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85β. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85β-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.
Collapse
|
5
|
Hinz A, Amado A, Kassen R, Bank C, Wong A. Unpredictability of the Fitness Effects of Antimicrobial Resistance Mutations Across Environments in Escherichia coli. Mol Biol Evol 2024; 41:msae086. [PMID: 38709811 PMCID: PMC11110942 DOI: 10.1093/molbev/msae086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Revised: 04/11/2024] [Accepted: 04/30/2024] [Indexed: 05/08/2024] Open
Abstract
The evolution of antimicrobial resistance (AMR) in bacteria is a major public health concern, and antibiotic restriction is often implemented to reduce the spread of resistance. These measures rely on the existence of deleterious fitness effects (i.e. costs) imposed by AMR mutations during growth in the absence of antibiotics. According to this assumption, resistant strains will be outcompeted by susceptible strains that do not pay the cost during the period of restriction. The fitness effects of AMR mutations are generally studied in laboratory reference strains grown in standard growth environments; however, the genetic and environmental context can influence the magnitude and direction of a mutation's fitness effects. In this study, we measure how three sources of variation impact the fitness effects of Escherichia coli AMR mutations: the type of resistance mutation, the genetic background of the host, and the growth environment. We demonstrate that while AMR mutations are generally costly in antibiotic-free environments, their fitness effects vary widely and depend on complex interactions between the mutation, genetic background, and environment. We test the ability of the Rough Mount Fuji fitness landscape model to reproduce the empirical data in simulation. We identify model parameters that reasonably capture the variation in fitness effects due to genetic variation. However, the model fails to accommodate the observed variation when considering multiple growth environments. Overall, this study reveals a wealth of variation in the fitness effects of resistance mutations owing to genetic background and environmental conditions, which will ultimately impact their persistence in natural populations.
Collapse
Affiliation(s)
- Aaron Hinz
- Department of Biology, Carleton University, Ottawa, ON K1S 5B6, Canada
- Department of Biology, University of Ottawa, Ottawa, ON K1N 6N5, Canada
- Department of Biology, McGill University, Montreal, QC H3A 1B1, Canada
| | - André Amado
- Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
- Division of Theoretical Ecology and Evolution, Swiss Institute of Bioinformatics, Lausanne, Switzerland
- Evolutionary Dynamics Group, Gulbenkian Science Institute, Oeiras, Portugal
| | - Rees Kassen
- Department of Biology, University of Ottawa, Ottawa, ON K1N 6N5, Canada
- Department of Biology, McGill University, Montreal, QC H3A 1B1, Canada
| | - Claudia Bank
- Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
- Division of Theoretical Ecology and Evolution, Swiss Institute of Bioinformatics, Lausanne, Switzerland
- Evolutionary Dynamics Group, Gulbenkian Science Institute, Oeiras, Portugal
| | - Alex Wong
- Department of Biology, Carleton University, Ottawa, ON K1S 5B6, Canada
| |
Collapse
|
6
|
Baverstock K. Responses to commentaries on "The gene: An appraisal". PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2024; 188:31-42. [PMID: 38360273 DOI: 10.1016/j.pbiomolbio.2024.02.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2023] [Accepted: 02/02/2024] [Indexed: 02/17/2024]
Abstract
The central conclusions of "The Gene: An Appraisal" are that genetic variance does not underpin biological evolution, and, therefore, that genes are not Mendel's units of inheritance. In this response, I will address the criticisms I have received via commentaries on that paper by defending the following statements: 1. Epistasis does not explain the power-law fitness profile of the Long-Term Evolution Experiment (LTEE). The data from the evolution of natural systems displays the power-law form ubiquitously. Epistasis plays no role in evolution. 2. The common characteristics of living things (natural systems) are described by the principle of least action in de Maupertuis's original form, which is synonymous with the 2nd law of thermodynamics and Newton's 2nd law of motion in its complete form, i.e., F = dp/dt. Organisms strive to achieve free energy balance with their environments. 3. Based on an appraisal of the scientific environment between 1880 and 1911, I conclude that Johannsen's genotype conception was perhaps, the only option available to him. 4. The power-law fitness profile of the LTEE falsifies Fisher's Genetical Theory of Natural Selection, Johannsen's genotype conception, and the idea that 'Darwinian evolution' is an exception to the generic thermodynamic process of evolution in natural systems. 5. The use of the technique of genome-wide association to identify the causes and the likelihoods of inherited common diseases and behavioural traits is a 'wild goose chase' because genes are not the units of inheritance.
Collapse
Affiliation(s)
- Keith Baverstock
- Department of Environmental and Biological Sciences, University of Eastern Finland, Kuopio Campus, Kuopio, Finland.
| |
Collapse
|
7
|
Venkataraman P, Nagendra P, Ahlawat N, Brajesh RG, Saini S. Convergent genetic adaptation of Escherichia coli in minimal media leads to pleiotropic divergence. Front Mol Biosci 2024; 11:1286824. [PMID: 38660375 PMCID: PMC11039892 DOI: 10.3389/fmolb.2024.1286824] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 02/15/2024] [Indexed: 04/26/2024] Open
Abstract
Adaptation in an environment can either be beneficial, neutral or disadvantageous in another. To test the genetic basis of pleiotropic behaviour, we evolved six lines of E. coli independently in environments where glucose and galactose were the sole carbon sources, for 300 generations. All six lines in each environment exhibit convergent adaptation in the environment in which they were evolved. However, pleiotropic behaviour was observed in several environmental contexts, including other carbon environments. Genome sequencing reveals that mutations in global regulators rpoB and rpoC cause this pleiotropy. We report three new alleles of the rpoB gene, and one new allele of the rpoC gene. The novel rpoB alleles confer resistance to Rifampicin, and alter motility. Our results show how single nucleotide changes in the process of adaptation in minimal media can lead to wide-scale pleiotropy, resulting in changes in traits that are not under direct selection.
Collapse
Affiliation(s)
| | | | | | | | - Supreet Saini
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
8
|
Li Y, Barton JP. Correlated Allele Frequency Changes Reveal Clonal Structure and Selection in Temporal Genetic Data. Mol Biol Evol 2024; 41:msae060. [PMID: 38507665 PMCID: PMC10986812 DOI: 10.1093/molbev/msae060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 02/02/2024] [Accepted: 03/15/2024] [Indexed: 03/22/2024] Open
Abstract
In evolving populations where the rate of beneficial mutations is large, subpopulations of individuals with competing beneficial mutations can be maintained over long times. Evolution with this kind of clonal structure is commonly observed in a wide range of microbial and viral populations. However, it can be difficult to completely resolve clonal dynamics in data. This is due to limited read lengths in high-throughput sequencing methods, which are often insufficient to directly measure linkage disequilibrium or determine clonal structure. Here, we develop a method to infer clonal structure using correlated allele frequency changes in time-series sequence data. Simulations show that our method recovers true, underlying clonal structures when they are known and accurately estimate linkage disequilibrium. This information can then be combined with other inference methods to improve estimates of the fitness effects of individual mutations. Applications to data suggest novel clonal structures in an E. coli long-term evolution experiment, and yield improved predictions of the effects of mutations on bacterial fitness and antibiotic resistance. Moreover, our method is computationally efficient, requiring orders of magnitude less run time for large data sets than existing methods. Overall, our method provides a powerful tool to infer clonal structures from data sets where only allele frequencies are available, which can also improve downstream analyses.
Collapse
Affiliation(s)
- Yunxiao Li
- Department of Physics and Astronomy, University of California, Riverside, CA 92521, USA
| | - John P Barton
- Department of Physics and Astronomy, University of California, Riverside, CA 92521, USA
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15260, USA
| |
Collapse
|
9
|
DelaFuente J, Diaz-Colunga J, Sanchez A, San Millan A. Global epistasis in plasmid-mediated antimicrobial resistance. Mol Syst Biol 2024; 20:311-320. [PMID: 38409539 PMCID: PMC10987494 DOI: 10.1038/s44320-024-00012-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Revised: 01/03/2024] [Accepted: 01/04/2024] [Indexed: 02/28/2024] Open
Abstract
Antimicrobial resistance (AMR) in bacteria is a major public health threat and conjugative plasmids play a key role in the dissemination of AMR genes among bacterial pathogens. Interestingly, the association between AMR plasmids and pathogens is not random and certain associations spread successfully at a global scale. The burst of genome sequencing has increased the resolution of epidemiological programs, broadening our understanding of plasmid distribution in bacterial populations. Despite the immense value of these studies, our ability to predict future plasmid-bacteria associations remains limited. Numerous empirical studies have recently reported systematic patterns in genetic interactions that enable predictability, in a phenomenon known as global epistasis. In this perspective, we argue that global epistasis patterns hold the potential to predict interactions between plasmids and bacterial genomes, thereby facilitating the prediction of future successful associations. To assess the validity of this idea, we use previously published data to identify global epistasis patterns in clinically relevant plasmid-bacteria associations. Furthermore, using simple mechanistic models of antibiotic resistance, we illustrate how global epistasis patterns may allow us to generate new hypotheses on the mechanisms associated with successful plasmid-bacteria associations. Collectively, we aim at illustrating the relevance of exploring global epistasis in the context of plasmid biology.
Collapse
Affiliation(s)
| | - Juan Diaz-Colunga
- Centro Nacional de Biotecnología (CNB-CSIC), Madrid, Spain
- Institute of Functional Biology & Genomics, IBFG - CSIC, Universidad de Salamanca, Salamanca, Spain
| | - Alvaro Sanchez
- Centro Nacional de Biotecnología (CNB-CSIC), Madrid, Spain.
- Institute of Functional Biology & Genomics, IBFG - CSIC, Universidad de Salamanca, Salamanca, Spain.
| | - Alvaro San Millan
- Centro Nacional de Biotecnología (CNB-CSIC), Madrid, Spain.
- Centro de Investigación Biológica en Red de Epidemiología y Salud Pública (CIBERESP), Instituto de Salud Carlos III, Madrid, Spain.
| |
Collapse
|
10
|
Mehra P, Hintze A. Reducing Epistasis and Pleiotropy Can Avoid the Survival of the Flattest Tragedy. BIOLOGY 2024; 13:193. [PMID: 38534462 DOI: 10.3390/biology13030193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2024] [Revised: 03/10/2024] [Accepted: 03/15/2024] [Indexed: 03/28/2024]
Abstract
This study investigates whether reducing epistasis and pleiotropy enhances mutational robustness in evolutionary adaptation, utilizing an indirect encoded model within the "survival of the flattest" (SoF) fitness landscape. By simulating genetic variations and their phenotypic consequences, we explore organisms' adaptive mechanisms to maintain positions on higher, narrower evolutionary peaks amidst environmental and genetic pressures. Our results reveal that organisms can indeed sustain their advantageous positions by minimizing the complexity of genetic interactions-specifically, by reducing the levels of epistasis and pleiotropy. This finding suggests a counterintuitive strategy for evolutionary stability: simpler genetic architectures, characterized by fewer gene interactions and multifunctional genes, confer a survival advantage by enhancing mutational robustness. This study contributes to our understanding of the genetic underpinnings of adaptability and robustness, challenging traditional views that equate complexity with fitness in dynamic environments.
Collapse
Affiliation(s)
- Priyanka Mehra
- Department for MicroData Analytics, Dalarna University, 791 88 Falun, Sweden
| | - Arend Hintze
- Department for MicroData Analytics, Dalarna University, 791 88 Falun, Sweden
- BEACON Center for the Study of Evolution in Action, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
11
|
Park Y, Metzger BP, Thornton JW. The simplicity of protein sequence-function relationships. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.09.02.556057. [PMID: 37732229 PMCID: PMC10508729 DOI: 10.1101/2023.09.02.556057] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/22/2023]
Abstract
How complicated is the genetic architecture of proteins - the set of causal effects by which sequence determines function? High-order epistatic interactions among residues are thought to be pervasive, making a protein's function difficult to predict or understand from its sequence. Most studies, however, used methods that overestimate epistasis, because they analyze genetic architecture relative to a designated reference sequence - causing measurement noise and small local idiosyncrasies to propagate into pervasive high-order interactions - or have not effectively accounted for global nonlinearity in the sequence-function relationship. Here we present a new reference-free method that jointly estimates global nonlinearity and specific epistatic interactions across a protein's entire genotype-phenotype map. This method yields a maximally efficient explanation of a protein's genetic architecture and is more robust than existing methods to measurement noise, partial sampling, and model misspecification. We reanalyze 20 combinatorial mutagenesis experiments from a diverse set of proteins and find that additive and pairwise effects, along with a simple nonlinearity to account for limited dynamic range, explain a median of 96% of total variance in measured phenotypes (and >92% in every case). Only a tiny fraction of genotypes are strongly affected by third- or higher-order epistasis. Genetic architecture is also sparse: the number of terms required to explain the vast majority of variance is smaller than the number of genotypes by many orders of magnitude. The sequence-function relationship in most proteins is therefore far simpler than previously thought, opening the way for new and tractable approaches to characterize it.
Collapse
Affiliation(s)
- Yeonwoo Park
- Committee on Genetics, Genomics, and Systems Biology, University of Chicago, Chicago, IL 60637
- Current affiliation: Center for RNA Research, Institute for Basic Science, Seoul, Republic of Korea 08826
| | - Brian P.H. Metzger
- Department of Ecology and Evolution, University of Chicago, Chicago, IL 60637
- Current affiliation: Department of Biological Sciences, Purdue University, West Lafayette, IN 47907
| | - Joseph W. Thornton
- Department of Ecology and Evolution, University of Chicago, Chicago, IL 60637
- Department of Human Genetics, University of Chicago, Chicago, IL 60637
| |
Collapse
|
12
|
Duan B, Qiu C, Lockless SW, Sze SH, Kaplan CD. Higher-order epistasis within Pol II trigger loop haplotypes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.20.576280. [PMID: 38293233 PMCID: PMC10827151 DOI: 10.1101/2024.01.20.576280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
RNA polymerase II (Pol II) has a highly conserved domain, the trigger loop (TL), that controls transcription fidelity and speed. We previously probed pairwise genetic interactions between residues within and surrounding the TL and identified widespread incompatibility between TLs of different species when placed in the Saccharomyces cerevisiae Pol II context, indicating epistasis between the TL and its surrounding context. We sought to understand the nature of this incompatibility and probe higher order epistasis internal to the TL. We have employed deep mutational scanning with selected natural TL variants ("haplotypes"), and all possible intermediate substitution combinations between them and the yeast Pol II TL. We identified both positive and negative higher-order residue interactions within example TL haplotypes. Intricate higher-order epistasis formed by TL residues was sometimes only apparent from analysis of intermediate genotypes, emphasizing complexity of epistatic interactions. Furthermore, we distinguished TL substitutions with distinct classes of epistatic patterns, suggesting specific TL residues that potentially influence TL evolution. Our examples of complex residue interactions suggest possible pathways for epistasis to facilitate Pol II evolution.
Collapse
Affiliation(s)
- Bingbing Duan
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, PA 15260
| | - Chenxi Qiu
- Department of Genetics, Harvard Medical School, Boston, MA 02215
| | - Steve W Lockless
- Department of Biology, Texas A&M University, College Station, TX 77843
| | - Sing-Hoi Sze
- Department of Computer Science & Engineering, Texas A&M University, College Station, TX 77843
- Department of Biochemistry & Biophysics, Texas A&M University, College Station, TX 77843
| | - Craig D Kaplan
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, PA 15260
| |
Collapse
|
13
|
Diaz-Colunga J, Sanchez A, Ogbunugafor CB. Environmental modulation of global epistasis in a drug resistance fitness landscape. Nat Commun 2023; 14:8055. [PMID: 38052815 DOI: 10.1038/s41467-023-43806-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2022] [Accepted: 11/21/2023] [Indexed: 12/07/2023] Open
Abstract
Interactions between mutations (epistasis) can add substantial complexity to genotype-phenotype maps, hampering our ability to predict evolution. Yet, recent studies have shown that the fitness effect of a mutation can often be predicted from the fitness of its genetic background using simple, linear relationships. This phenomenon, termed global epistasis, has been leveraged to reconstruct fitness landscapes and infer adaptive trajectories in a wide variety of contexts. However, little attention has been paid to how patterns of global epistasis may be affected by environmental variation, despite this variation frequently being a major driver of evolution. This is particularly relevant for the evolution of drug resistance, where antimicrobial drugs may change the environment faced by pathogens and shape their adaptive trajectories in ways that can be difficult to predict. By analyzing a fitness landscape of four mutations in a gene encoding an essential enzyme of P. falciparum (a parasite cause of malaria), here we show that patterns of global epistasis can be strongly modulated by the concentration of a drug in the environment. Expanding on previous theoretical results, we demonstrate that this modulation can be quantitatively explained by how specific gene-by-gene interactions are modified by drug dose. Importantly, our results highlight the need to incorporate potential environmental variation into the global epistasis framework in order to predict adaptation in dynamic environments.
Collapse
Affiliation(s)
- Juan Diaz-Colunga
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT, 06511, USA.
- Department of Microbial Biotechnology, Spanish National Center for Biotechnology CNB-CSIC, 28049, Madrid, Spain.
- Institute of Functional Biology and Genomics IBFG-CSIC, University of Salamanca, 37007, Salamanca, Spain.
| | - Alvaro Sanchez
- Department of Microbial Biotechnology, Spanish National Center for Biotechnology CNB-CSIC, 28049, Madrid, Spain.
- Institute of Functional Biology and Genomics IBFG-CSIC, University of Salamanca, 37007, Salamanca, Spain.
| | - C Brandon Ogbunugafor
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, CT, 06511, USA.
- Santa Fe Institute, Santa Fe, NM, 87501, USA.
| |
Collapse
|
14
|
Wang Y, Sang M, Feng L, Gragnoli C, Griffin C, Wu R. A pleiotropic-epistatic entangelement model of drug response. Drug Discov Today 2023; 28:103790. [PMID: 37758020 DOI: 10.1016/j.drudis.2023.103790] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2023] [Revised: 09/10/2023] [Accepted: 09/20/2023] [Indexed: 10/03/2023]
Abstract
Because drug response is multifactorial, graph models are uniquely powerful for comprehending its genetic architecture. We deconstruct drug response into many different and interdependent sub-traits, with each sub-trait controlled by multiple genes that act and interact in a complicated manner. The outcome of drug response is the consequence of multileveled intertwined interactions between pleiotropic effects and epistatic effects. Here, we propose a general statistical physics framework to chart the 3D geometric network that codes how epistasis pleiotropically influences a complete set of sub-traits to shape body-drug interactions. This model can dissect the topological architecture of epistatically induced pleiotropic networks (EiPN) and pleiotropically influenced epistatic networks (PiEN). We analyze and interpret the practical implications of the pleiotropic-epistatic entanglement model for pharmacogenomic studies.
Collapse
Affiliation(s)
- Yu Wang
- Center for Computational Biology, Beijing Forestry University, Beijing 100083, China
| | - Mengmeng Sang
- Institute of Reproductive Medicine, Medical School of Nantong University, Nantong, Jiangsu 226019, China
| | - Li Feng
- Fisheries Engineering Institute, Chinese Academy of Fishery Sciences, Beijing 1000141, China
| | - Claudia Gragnoli
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA 17033, USA; Department of Medicine, Creighton University School of Medicine, Omaha, NE 68124, USA; Molecular Biology Laboratory, Bios Biotech Multi-Diagnostic Health Center, Rome 00197, Italy
| | - Christopher Griffin
- Applied Research Laboratory, The Pennsylvania State University, University Park, PA 16802, USA
| | - Rongling Wu
- Center for Computational Biology, Beijing Forestry University, Beijing 100083, China; Beijing Yanqi Lake Institute of Mathematical Sciences and Applications, Beijing 101408, China; Yau Mathematical Sciences Center, Tsinghua University, Beijing 100084, China.
| |
Collapse
|