1
|
Sanchez-Gorostiaga A, Bajić D, Osborne ML, Poyatos JF, Sanchez A. High-order interactions distort the functional landscape of microbial consortia. PLoS Biol 2019; 17:e3000550. [PMID: 31830028 PMCID: PMC6932822 DOI: 10.1371/journal.pbio.3000550] [Citation(s) in RCA: 76] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2019] [Revised: 12/26/2019] [Accepted: 11/15/2019] [Indexed: 12/11/2022] Open
Abstract
Understanding the link between community composition and function is a major challenge in microbial population biology, with implications for the management of natural microbiomes and the design of synthetic consortia. Specifically, it is poorly understood whether community functions can be quantitatively predicted from traits of species in monoculture. Inspired by the study of complex genetic interactions, we have examined how the amylolytic rate of combinatorial assemblages of six starch-degrading soil bacteria depend on the separate functional contributions from each species and their interactions. Filtering our results through the theory of biochemical kinetics, we show that this simple function is additive in the absence of interactions among community members. For about half of the combinatorially assembled consortia, the amylolytic function is dominated by pairwise and higher-order interactions. For the other half, the function is additive despite the presence of strong competitive interactions. We explain the mechanistic basis of these findings and propose a quantitative framework that allows us to separate the effect of behavioral and population dynamics interactions. Our results suggest that the functional robustness of a consortium to pairwise and higher-order interactions critically affects our ability to predict and bottom-up engineer ecosystem function in complex communities.
Collapse
Affiliation(s)
- Alicia Sanchez-Gorostiaga
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Microbial Sciences Institute, Yale University, West Haven, Connecticut, United States of America
| | - Djordje Bajić
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Microbial Sciences Institute, Yale University, West Haven, Connecticut, United States of America
| | - Melisa L. Osborne
- The Rowland Institute at Harvard, Harvard University, Cambridge, Massachusetts, United States of America
- Biological Design Center, Boston University, Boston, Massachusetts, United States of America
| | - Juan F. Poyatos
- The Rowland Institute at Harvard, Harvard University, Cambridge, Massachusetts, United States of America
- Logic of Genomic Systems Laboratory, Spanish National Biotechnology Centre (CNB-CSIC), Madrid, Spain
| | - Alvaro Sanchez
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Microbial Sciences Institute, Yale University, West Haven, Connecticut, United States of America
- The Rowland Institute at Harvard, Harvard University, Cambridge, Massachusetts, United States of America
| |
Collapse
|
2
|
Sanchez-Gorostiaga A, Bajić D, Osborne ML, Poyatos JF, Sanchez A. High-order interactions distort the functional landscape of microbial consortia. PLoS Biol 2019; 17:e3000550. [PMID: 31830028 DOI: 10.1101/333534] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2019] [Revised: 12/26/2019] [Accepted: 11/15/2019] [Indexed: 05/23/2023] Open
Abstract
Understanding the link between community composition and function is a major challenge in microbial population biology, with implications for the management of natural microbiomes and the design of synthetic consortia. Specifically, it is poorly understood whether community functions can be quantitatively predicted from traits of species in monoculture. Inspired by the study of complex genetic interactions, we have examined how the amylolytic rate of combinatorial assemblages of six starch-degrading soil bacteria depend on the separate functional contributions from each species and their interactions. Filtering our results through the theory of biochemical kinetics, we show that this simple function is additive in the absence of interactions among community members. For about half of the combinatorially assembled consortia, the amylolytic function is dominated by pairwise and higher-order interactions. For the other half, the function is additive despite the presence of strong competitive interactions. We explain the mechanistic basis of these findings and propose a quantitative framework that allows us to separate the effect of behavioral and population dynamics interactions. Our results suggest that the functional robustness of a consortium to pairwise and higher-order interactions critically affects our ability to predict and bottom-up engineer ecosystem function in complex communities.
Collapse
Affiliation(s)
- Alicia Sanchez-Gorostiaga
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Microbial Sciences Institute, Yale University, West Haven, Connecticut, United States of America
| | - Djordje Bajić
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Microbial Sciences Institute, Yale University, West Haven, Connecticut, United States of America
| | - Melisa L Osborne
- The Rowland Institute at Harvard, Harvard University, Cambridge, Massachusetts, United States of America
- Biological Design Center, Boston University, Boston, Massachusetts, United States of America
| | - Juan F Poyatos
- The Rowland Institute at Harvard, Harvard University, Cambridge, Massachusetts, United States of America
- Logic of Genomic Systems Laboratory, Spanish National Biotechnology Centre (CNB-CSIC), Madrid, Spain
| | - Alvaro Sanchez
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Microbial Sciences Institute, Yale University, West Haven, Connecticut, United States of America
- The Rowland Institute at Harvard, Harvard University, Cambridge, Massachusetts, United States of America
| |
Collapse
|
3
|
Gan J, Cao Y, Jiang L, Wu R. Mapping Covariation Quantitative Trait Loci That Control Organ Growth and Whole-Plant Biomass. FRONTIERS IN PLANT SCIENCE 2019; 10:719. [PMID: 31214231 PMCID: PMC6558071 DOI: 10.3389/fpls.2019.00719] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/09/2018] [Accepted: 05/16/2019] [Indexed: 06/09/2023]
Abstract
Covariation between organ growth and biomass accumulation plays an important role in plants. Plant to capture optimal fitness in nature, which depend coordinate and interact for distinct organs such as leaves, stems, and roots. Although many studies have focused on plant growth or biomass allocation, detailed information on the genetic mechanism of coordinated variation is lacking. Here, we expand a new mapping model based on functional mapping to detect covariation quantitative trait loci (QTLs) that govern development of plant organs and whole biomass, which, via a series of hypothesis tests, allows quantification of how QTLs regulate covariation between organ growth and biomass accumulation. The model was implemented to analyze leaf number data and the whole dry weight of recombinant inbred lines (RILs) of Arabidopsis. Two key QTLs related to growth and biomass allocation that reside within biologically meaningful genes, CRA1 and HIPP25, are characterized. These two genes may control covariation between two traits. The new model will enable the elucidation of the genetic architecture underlying growth and biomass accumulation, which may enhance our understanding of fitness development in plants.
Collapse
Affiliation(s)
- Jingwen Gan
- Center for Computational Biology, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Yige Cao
- Center for Computational Biology, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Libo Jiang
- Center for Computational Biology, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Rongling Wu
- Center for Computational Biology, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- Center for Statistical Genetics, Departments of Public Health Sciences and Statistics, The Pennsylvania State University, Hershey, PA, United States
| |
Collapse
|
4
|
Pedruzzi G, Barlukova A, Rouzine IM. Evolutionary footprint of epistasis. PLoS Comput Biol 2018; 14:e1006426. [PMID: 30222748 PMCID: PMC6177197 DOI: 10.1371/journal.pcbi.1006426] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2018] [Revised: 10/09/2018] [Accepted: 08/09/2018] [Indexed: 11/18/2022] Open
Abstract
Variation of an inherited trait across a population cannot be explained by additive contributions of relevant genes, due to epigenetic effects and biochemical interactions (epistasis). Detecting epistasis in genomic data still represents a significant challenge that requires a better understanding of epistasis from the mechanistic point of view. Using a standard Wright-Fisher model of bi-allelic asexual population, we study how compensatory epistasis affects the process of adaptation. The main result is a universal relationship between four haplotype frequencies of a single site pair in a genome, which depends only on the epistasis strength of the pair defined regarding Darwinian fitness. We demonstrate the existence, at any time point, of a quasi-equilibrium between epistasis and disorder (entropy) caused by random genetic drift and mutation. We verify the accuracy of these analytic results by Monte-Carlo simulation over a broad range of parameters, including the topology of the interacting network. Thus, epistasis assists the evolutionary transit through evolutionary hurdles leaving marks at the level of haplotype disequilibrium. The method allows determining selection coefficient for each site and the epistasis strength of each pair from a sequence set. The resulting ability to detect clusters of deleterious mutation close to full compensation is essential for biomedical applications. These findings help to understand the role of epistasis in multiple compensatory mutations in viral resistance to antivirals and immune response.
Collapse
Affiliation(s)
- Gabriele Pedruzzi
- Sorbonne Université, Institute de Biologie Paris-Seine, Laboratoire de Biologie Computationelle et Quantitative, Paris, France
| | - Ayuna Barlukova
- Sorbonne Université, Institute de Biologie Paris-Seine, Laboratoire de Biologie Computationelle et Quantitative, Paris, France
| | - Igor M. Rouzine
- Sorbonne Université, Institute de Biologie Paris-Seine, Laboratoire de Biologie Computationelle et Quantitative, Paris, France
- * E-mail:
| |
Collapse
|
5
|
Gosik K, Sun L, Chinchilli VM, Wu R. An Ultrahigh-Dimensional Mapping Model of High-order Epistatic Networks for Complex Traits. Curr Genomics 2018; 19:384-394. [PMID: 30065614 PMCID: PMC6030858 DOI: 10.2174/1389202919666171218162210] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2017] [Revised: 03/28/2017] [Accepted: 05/04/2017] [Indexed: 11/23/2022] Open
Abstract
BACKGROUND Genetic interactions involving more than two loci have been thought to affect quantitatively inherited traits and diseases more pervasively than previously appreciated. However, the detection of such high-order interactions to chart a complete portrait of genetic architecture has not been well explored. METHODS We present an ultrahigh-dimensional model to systematically characterize genetic main effects and interaction effects of various orders among all possible markers in a genetic mapping or association study. The model was built on the extension of a variable selection procedure, called iFORM, derived from forward selection. The model shows its unique power to estimate the magnitudes and signs of high-order epistatic effects, in addition to those of main effects and pairwise epistatic effects. RESULTS The statistical properties of the model were tested and validated through simulation studies. By analyzing a real data for shoot growth in a mapping population of woody plant, mei (Prunus mume), we demonstrated the usefulness and utility of the model in practical genetic studies. The model has identified important high-order interactions that contribute to shoot growth for mei. CONCLUSION The model provides a tool to precisely construct genotype-phenotype maps for quantitative traits by identifying any possible high-order epistasis which is often ignored in the current genetic literature.
Collapse
Affiliation(s)
- Kirk Gosik
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA17033, USA
| | - Lidan Sun
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA17033, USA
| | - Vernon M. Chinchilli
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA17033, USA
| | - Rongling Wu
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA17033, USA
| |
Collapse
|
6
|
Kulwal PL. Trait Mapping Approaches Through Linkage Mapping in Plants. PLANT GENETICS AND MOLECULAR BIOLOGY 2018; 164:53-82. [DOI: 10.1007/10_2017_49] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
7
|
Gosik K, Kong L, Chinchilli VM, Wu R. iFORM/eQTL: an ultrahigh-dimensional platform for inferring the global genetic architecture of gene transcripts. Brief Bioinform 2017; 18:250-259. [PMID: 26944084 DOI: 10.1093/bib/bbw014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2015] [Indexed: 01/03/2023] Open
Abstract
Knowledge about how changes in gene expression are encoded by expression quantitative trait loci (eQTLs) is a key to construct the genotype-phenotype map for complex traits or diseases. Traditional eQTL mapping is to associate one transcript with a single marker at a time, thereby limiting our inference about a complete picture of the genetic architecture of gene expression. Here, we implemented an ultrahigh-dimensional variable selection model to build a computing platform that can systematically scan main effects and interaction effects among all possible loci and identify a set of significant eQTLs modulating differentiation and function of gene expression. This platform, named iFORM/eQTL, was assembled by forward-selection-based procedures to tackle complex covariance structures of gene-gene interactions. iFORM/eQTL can particularly discern the role of cis-QTLs, trans-QTLs and their epistatic interactions in gene expression. Results from the reanalysis of a published genetic and genomic data set through iFORM/eQTL gain new discoveries on the genetic origin of gene expression differentiation in Caenorhabditis elegans, which could not be detected by a traditional one-locus/one-transcript analysis approach.
Collapse
Affiliation(s)
- Kirk Gosik
- Department of Statistics, The Pennsylvania State University, University Park, PA, USA
| | - Lan Kong
- Department of Public Health Sciences Pennsylvania State University College of Medicine, Hershey, PA, USA
| | - Vernon M Chinchilli
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA, USA
| | - Rongling Wu
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA, USA
| |
Collapse
|
8
|
Jiang L, Zhang M, Sang M, Ye M, Wu R. Evo-Devo-EpiR: a genome-wide search platform for epistatic control on the evolution of development. Brief Bioinform 2017; 18:754-760. [PMID: 27473062 DOI: 10.1093/bib/bbw062] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2016] [Indexed: 11/14/2022] Open
Abstract
Evo-devo is a theory proposed to study how phenotypes evolve by comparing the developmental processes of different organisms or the same organism experiencing changing environments. It has been recognized that nonallelic interactions at different genes or quantitative trait loci, known as epistasis, may play a pivotal role in the evolution of development, but it has proven difficult to quantify and elucidate this role into a coherent picture. We implement a high-dimensional genome-wide association study model into the evo-devo paradigm and pack it into the R-based Evo-Devo-EpiR, aimed at facilitating the genome-wide landscaping of epistasis for the diversification of phenotypic development. By analyzing a high-throughput assay of DNA markers and their pairs simultaneously, Evo-Devo-EpiR is equipped with a capacity to systematically characterize various epistatic interactions that impact on the pattern and timing of development and its evolution. Enabling a global search for all possible genetic interactions for developmental processes throughout the whole genome, Evo-Devo-EpiR provides a computational tool to illustrate a precise genotype-phenotype map at interface between epistasis, development and evolution.
Collapse
|
9
|
Bessonov K, Gusareva ES, Van Steen K. A cautionary note on the impact of protocol changes for genome-wide association SNP × SNP interaction studies: an example on ankylosing spondylitis. Hum Genet 2015; 134:761-73. [DOI: 10.1007/s00439-015-1560-7] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2015] [Accepted: 04/26/2015] [Indexed: 12/11/2022]
|
10
|
Abstract
Genome-wide association studies (GWASs) have become the focus of the statistical analysis of complex traits in humans, successfully shedding light on several aspects of genetic architecture and biological aetiology. Single-nucleotide polymorphisms (SNPs) are usually modelled as having additive, cumulative and independent effects on the phenotype. Although evidently a useful approach, it is often argued that this is not a realistic biological model and that epistasis (that is, the statistical interaction between SNPs) should be included. The purpose of this Review is to summarize recent directions in methodology for detecting epistasis and to discuss evidence of the role of epistasis in human complex trait variation. We also discuss the relevance of epistasis in the context of GWASs and potential hazards in the interpretation of statistical interaction terms.
Collapse
|
11
|
Tong C, Shen L, Lv Y, Wang Z, Wang X, Feng S, Li X, Sui Y, Pang X, Wu R. Structural mapping: how to study the genetic architecture of a phenotypic trait through its formation mechanism. Brief Bioinform 2012; 15:43-53. [PMID: 23104859 DOI: 10.1093/bib/bbs067] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Traditional approaches for genetic mapping are to simply associate the genotypes of a quantitative trait locus (QTL) with the phenotypic variation of a complex trait. A more mechanistic strategy has emerged to dissect the trait phenotype into its structural components and map specific QTLs that control the mechanistic and structural formation of a complex trait. We describe and assess such a strategy, called structural mapping, by integrating the internal structural basis of trait formation into a QTL mapping framework. Electrical impedance spectroscopy (EIS) has been instrumental for describing the structural components of a phenotypic trait and their interactions. By building robust mathematical models on circuit EIS data and embedding these models within a mixture model-based likelihood for QTL mapping, structural mapping implements the EM algorithm to obtain maximum likelihood estimates of QTL genotype-specific EIS parameters. The uniqueness of structural mapping is to make it possible to test a number of hypotheses about the pattern of the genetic control of structural components. We validated structural mapping by analyzing an EIS data collected for QTL mapping of frost hardiness in a controlled cross of jujube trees. The statistical properties of parameter estimates were examined by simulation studies. Structural mapping can be a powerful alternative for genetic mapping of complex traits by taking account into the biological and physical mechanisms underlying their formation.
Collapse
Affiliation(s)
- Chunfa Tong
- *These authors contributed equally to this work
| | | | | | | | | | | | | | | | | | | |
Collapse
|