1
|
Eble H, Joswig M, Lamberti L, Ludington WB. Master regulators of biological systems in higher dimensions. Proc Natl Acad Sci U S A 2023; 120:e2300634120. [PMID: 38096409 PMCID: PMC10743376 DOI: 10.1073/pnas.2300634120] [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: 01/12/2023] [Accepted: 10/23/2023] [Indexed: 12/18/2023] Open
Abstract
A longstanding goal of biology is to identify the key genes and species that critically impact evolution, ecology, and health. Network analysis has revealed keystone species that regulate ecosystems and master regulators that regulate cellular genetic networks. Yet these studies have focused on pairwise biological interactions, which can be affected by the context of genetic background and other species present, generating higher-order interactions. The important regulators of higher-order interactions are unstudied. To address this, we applied a high-dimensional geometry approach that quantifies epistasis in a fitness landscape to ask how individual genes and species influence the interactions in the rest of the biological network. We then generated and also reanalyzed 5-dimensional datasets (two genetic, two microbiome). We identified key genes (e.g., the rbs locus and pykF) and species (e.g., Lactobacilli) that control the interactions of many other genes and species. These higher-order master regulators can induce or suppress evolutionary and ecological diversification by controlling the topography of the fitness landscape. Thus, we provide a method and mathematical justification for exploration of biological networks in higher dimensions.
Collapse
Affiliation(s)
- Holger Eble
- Chair of Discrete Mathematics/Geometry, Technical University Berlin, Berlin10623, Germany
| | - Michael Joswig
- Chair of Discrete Mathematics/Geometry, Technical University Berlin, Berlin10623, Germany
- Max Planck Institute for Mathematics in the Sciences, Leipzig04103, Germany
| | - Lisa Lamberti
- Department of Biosystems Science and Engineering, Federal Institute of Technology (ETH Zürich), Basel4058, Switzerland
- Swiss Institute of Bioinformatics, Basel4058, Switzerland
| | - William B. Ludington
- Department of Biosphere Sciences and Engineering, Carnegie Institution for Science, Baltimore, MD21218
- Department of Biology, Johns Hopkins University, Baltimore, MD21218
| |
Collapse
|
2
|
Servajean R, Bitbol AF. Impact of population size on early adaptation in rugged fitness landscapes. Philos Trans R Soc Lond B Biol Sci 2023; 378:20220045. [PMID: 37004726 PMCID: PMC10067268 DOI: 10.1098/rstb.2022.0045] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2023] Open
Abstract
Owing to stochastic fluctuations arising from finite population size, known as genetic drift, the ability of a population to explore a rugged fitness landscape depends on its size. In the weak mutation regime, while the mean steady-state fitness increases with population size, we find that the height of the first fitness peak encountered when starting from a random genotype displays various behaviours versus population size, even among small and simple rugged landscapes. We show that the accessibility of the different fitness peaks is key to determining whether this height overall increases or decreases with population size. Furthermore, there is often a finite population size that maximizes the height of the first fitness peak encountered when starting from a random genotype. This holds across various classes of model rugged landscapes with sparse peaks, and in some experimental and experimentally inspired ones. Thus, early adaptation in rugged fitness landscapes can be more efficient and predictable for relatively small population sizes than in the large-size limit. This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.
Collapse
Affiliation(s)
- Richard Servajean
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- SIB Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Anne-Florence Bitbol
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- SIB Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| |
Collapse
|
3
|
Liu W, Tang S, Peng J, Zhu Y, Pan L, Wang J, Peng X, Cheng H, Chen Z, Wang Y, Zhou H. Enhancing lactose recognition of a key enzyme in 2'-fucosyllactose synthesis: α-1,2-fucosyltransferase. JOURNAL OF THE SCIENCE OF FOOD AND AGRICULTURE 2023; 103:1303-1314. [PMID: 36116126 DOI: 10.1002/jsfa.12224] [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: 07/06/2022] [Revised: 09/12/2022] [Accepted: 09/18/2022] [Indexed: 06/15/2023]
Abstract
BACKGROUND 2'-Fucosyllactose, a representative oligosaccharide in human milk, is an emerging and promising food and pharmaceutical ingredient due to its powerful health benefits, such as participating in immune regulation, regulation of intestinal flora, etc. To enable economically viable production of 2'-fucosyllactose, different biosynthesis strategies using precursors and pathway enzymes have been developed. The α-1,2-fucosyltransferases are an essential part involved in these strategies, but their strict substrate selectivity and unsatisfactory substrate tolerance are one of the key roadblocks limiting biosynthesis. RESULTS To tackle this issue, a semi-rational manipulation combining computer-aided designing and screening with biochemical experiments were adopted. The mutant had a 100-fold increase in catalytic efficiency compared to the wild-type. The highest 2'-fucosyllactose yield was up to 0.65 mol mol-1 lactose with a productivity of 2.56 g mL-1 h-1 performed by enzymatic catalysis in vitro. Further analysis revealed that the interactions between the mutant and substrates were reduced. The crucial contributions of wild-type and mutant to substrate recognition ability were closely related to their distinct phylotypes in terms of amino acid preference. CONCLUSION It is envisioned that the engineered α-1,2-fucosyltransferase could be harnessed to relieve constraints imposed on the bioproduction of 2'-fucosyllactose and lay a theoretical foundation for elucidating the substrate recognition mechanisms of fucosyltransferases. © 2022 Society of Chemical Industry.
Collapse
Affiliation(s)
- Wenxian Liu
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
| | - Shizhe Tang
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
| | - Jing Peng
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
| | - Yuling Zhu
- Changsha Yunkang Biotechnology Co Ltd, Changsha, P. R. China
| | - Lina Pan
- Ausnutria Institute Food & Nutrition, Ausnutria Dairy China Co Ltd, Changsha, P. R. China
| | - Jiaqi Wang
- Ausnutria Institute Food & Nutrition, Ausnutria Dairy China Co Ltd, Changsha, P. R. China
| | - Xiaoyu Peng
- Ausnutria Institute Food & Nutrition, Ausnutria Dairy China Co Ltd, Changsha, P. R. China
| | - Haina Cheng
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
- Key Laboratory of Biometallurgy of Ministry of Education, Central South University, Changsha, P. R. China
| | - Zhu Chen
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
- Key Laboratory of Biometallurgy of Ministry of Education, Central South University, Changsha, P. R. China
| | - Yuguang Wang
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
- Key Laboratory of Biometallurgy of Ministry of Education, Central South University, Changsha, P. R. China
| | - Hongbo Zhou
- School of Minerals Processing and Bioengineering, Central South University, Changsha, P. R. China
- Key Laboratory of Biometallurgy of Ministry of Education, Central South University, Changsha, P. R. China
| |
Collapse
|
4
|
Abstract
One core goal of genetics is to systematically understand the mapping between the DNA sequence of an organism (genotype) and its measurable characteristics (phenotype). Understanding this mapping is often challenging because of interactions between mutations, where the result of combining several different mutations can be very different than the sum of their individual effects. Here we provide a statistical framework for modeling complex genetic interactions of this type. The key idea is to ask how fast the effects of mutations change when introducing the same mutation in increasingly distant genetic backgrounds. We then propose a model for phenotypic prediction that takes into account this tendency for the effects of mutations to be more similar in nearby genetic backgrounds. Contemporary high-throughput mutagenesis experiments are providing an increasingly detailed view of the complex patterns of genetic interaction that occur between multiple mutations within a single protein or regulatory element. By simultaneously measuring the effects of thousands of combinations of mutations, these experiments have revealed that the genotype–phenotype relationship typically reflects not only genetic interactions between pairs of sites but also higher-order interactions among larger numbers of sites. However, modeling and understanding these higher-order interactions remains challenging. Here we present a method for reconstructing sequence-to-function mappings from partially observed data that can accommodate all orders of genetic interaction. The main idea is to make predictions for unobserved genotypes that match the type and extent of epistasis found in the observed data. This information on the type and extent of epistasis can be extracted by considering how phenotypic correlations change as a function of mutational distance, which is equivalent to estimating the fraction of phenotypic variance due to each order of genetic interaction (additive, pairwise, three-way, etc.). Using these estimated variance components, we then define an empirical Bayes prior that in expectation matches the observed pattern of epistasis and reconstruct the genotype–phenotype mapping by conducting Gaussian process regression under this prior. To demonstrate the power of this approach, we present an application to the antibody-binding domain GB1 and also provide a detailed exploration of a dataset consisting of high-throughput measurements for the splicing efficiency of human pre-mRNA 5′ splice sites, for which we also validate our model predictions via additional low-throughput experiments.
Collapse
|
5
|
Gabzi T, Pilpel Y, Friedlander T. Fitness landscape analysis of a tRNA gene reveals that the wild type allele is sub-optimal, yet mutationally robust. Mol Biol Evol 2022; 39:6670756. [PMID: 35976926 PMCID: PMC9447856 DOI: 10.1093/molbev/msac178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
Fitness landscape mapping and the prediction of evolutionary trajectories on these landscapes are major tasks in evolutionary biology research. Evolutionary dynamics is tightly linked to the landscape topography, but this relation is not straightforward. Here, we analyze a fitness landscape of a yeast tRNA gene, previously measured under four different conditions. We find that the wild type allele is sub-optimal, and 8–10% of its variants are fitter. We rule out the possibilities that the wild type is fittest on average on these four conditions or located on a local fitness maximum. Notwithstanding, we cannot exclude the possibility that the wild type might be fittest in some of the many conditions in the complex ecology that yeast lives at. Instead, we find that the wild type is mutationally robust (“flat”), while more fit variants are typically mutationally fragile. Similar observations of mutational robustness or flatness have been so far made in very few cases, predominantly in viral genomes.
Collapse
Affiliation(s)
- Tzahi Gabzi
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Yitzhak Pilpel
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Tamar Friedlander
- The Robert H. Smith Institute of Plant Sciences and Genetics in Agriculture Faculty of Agriculture, Hebrew University of Jerusalem, 229 Herzl St., Rehovot 7610001, Israel
| |
Collapse
|
6
|
Yan W, Wang B, Chan E, Mitchell-Olds T. Genetic architecture and adaptation of flowering time among environments. THE NEW PHYTOLOGIST 2021; 230:1214-1227. [PMID: 33484593 PMCID: PMC8193995 DOI: 10.1111/nph.17229] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 01/07/2021] [Indexed: 05/17/2023]
Abstract
The genetic basis of flowering time changes across environments, and pleiotropy may limit adaptive evolution of populations in response to local conditions. However, little information is known about how genetic architecture changes among environments. We used genome-wide association studies (GWAS) in Boechera stricta (Graham) Al-Shehbaz, a relative of Arabidopsis, to examine flowering variation among environments and associations with climate conditions in home environments. Also, we used molecular population genetics to search for evidence of historical natural selection. GWAS found 47 significant quantitative trait loci (QTLs) that influence flowering time in one or more environments, control plastic changes in phenology between experiments, or show associations with climate in sites of origin. Genetic architecture of flowering varied substantially among environments. We found that some pairs of QTLs showed similar patterns of pleiotropy across environments. A large-effect QTL showed molecular signatures of adaptive evolution and is associated with climate in home environments. The derived allele at this locus causes later flowering and predominates in sites with greater water availability. This work shows that GWAS of climate associations and ecologically important traits across diverse environments can be combined with molecular signatures of natural selection to elucidate ecological genetics of adaptive evolution.
Collapse
Affiliation(s)
- Wenjie Yan
- College of Marine Life Sciences, Ocean University of China, Qingdao 266003, China
- Department of Biology, Duke University, Box 90338, Durham, NC 27708, USA
| | - Baosheng Wang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Guangzhou, 510650 China
| | - Emily Chan
- Department of Biology, Duke University, Box 90338, Durham, NC 27708, USA
| | | |
Collapse
|
7
|
The evolutionary scaling of cellular traits imposed by the drift barrier. Proc Natl Acad Sci U S A 2020; 117:10435-10444. [PMID: 32345718 DOI: 10.1073/pnas.2000446117] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Owing to internal homeostatic mechanisms, cellular traits may experience long periods of stable selective pressures, during which the stochastic forces of drift and mutation conspire to generate variation. However, even in the face of invariant selection, the drift barrier defined by the genetic effective population size, which is negatively associated with organism size, can have a substantial influence on the location and dispersion of the long-term steady-state distribution of mean phenotypes. In addition, for multilocus traits, the multiplicity of alternative, functionally equivalent states can draw mean phenotypes away from selective optima, even in the absence of mutation bias. Using a framework for traits with an additive genetic basis, it is shown that 1) optimal phenotypic states may be only rarely achieved; 2) gradients of mean phenotypes with respect to organism size (i.e., allometric relationships) are likely to be molded by differences in the power of random genetic drift across the tree of life; and 3) for any particular set of population-genetic conditions, significant variation in mean phenotypes may exist among lineages exposed to identical selection pressures. These results provide a potentially useful framework for understanding numerous aspects of cellular diversification and illustrate the risks of interpreting such variation in a purely adaptive framework.
Collapse
|