1
|
Overcast I, Noguerales V, Meramveliotakis E, Andújar C, Arribas P, Creedy TJ, Emerson BC, Vogler AP, Papadopoulou A, Morlon H. Inferring the ecological and evolutionary determinants of community genetic diversity. Mol Ecol 2023; 32:6093-6109. [PMID: 37221561 DOI: 10.1111/mec.16958] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Revised: 04/06/2023] [Accepted: 04/12/2023] [Indexed: 05/25/2023]
Abstract
Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.
Collapse
Affiliation(s)
- Isaac Overcast
- Institut de Biologie de l'ENS (IBENS), Ecole Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
- Department of Vertebrate Zoology, American Museum of Natural History, New York, New York, USA
| | - Víctor Noguerales
- Instituto de Productos Naturales y Agrobiología (IPNA-CSIC), San Cristóbal de La Laguna, Spain
- Department of Biological Sciences, University of Cyprus, Nicosia, Cyprus
| | | | - Carmelo Andújar
- Instituto de Productos Naturales y Agrobiología (IPNA-CSIC), San Cristóbal de La Laguna, Spain
| | - Paula Arribas
- Instituto de Productos Naturales y Agrobiología (IPNA-CSIC), San Cristóbal de La Laguna, Spain
| | - Thomas J Creedy
- Department of Life Sciences, Natural History Museum, London, UK
| | - Brent C Emerson
- Instituto de Productos Naturales y Agrobiología (IPNA-CSIC), San Cristóbal de La Laguna, Spain
| | - Alfried P Vogler
- Department of Life Sciences, Natural History Museum, London, UK
- Department of Life Sciences, Imperial College London, Ascot, UK
| | - Anna Papadopoulou
- Department of Biological Sciences, University of Cyprus, Nicosia, Cyprus
| | - Hélène Morlon
- Institut de Biologie de l'ENS (IBENS), Ecole Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
| |
Collapse
|
2
|
Danesh G, Saulnier E, Gascuel O, Choisy M, Alizon S. TiPS
: Rapidly simulating trajectories and phylogenies from compartmental models. Methods Ecol Evol 2022. [DOI: 10.1111/2041-210x.14038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Gonché Danesh
- MIVEGEC, CNRS, IRD Université de Montpellier Montpellie France
| | - Emma Saulnier
- MIVEGEC, CNRS, IRD Université de Montpellier Montpellie France
| | | | - Marc Choisy
- Centre for Tropical Medicine and Global Health Nuffield Department of Medicine, University of Oxford Oxford UK
- Oxford University Clinical Research Unit Ho Chi Minh City Vietnam
| | - Samuel Alizon
- MIVEGEC, CNRS, IRD Université de Montpellier Montpellie France
- Center for Interdisciplinary Research in Biology (CIRB) College de France, CNRS, INSERM, Université PSL Paris France
| |
Collapse
|
3
|
Lamarins A, Fririon V, Folio D, Vernier C, Daupagne L, Labonne J, Buoro M, Lefèvre F, Piou C, Oddou‐Muratorio S. Importance of interindividual interactions in eco-evolutionary population dynamics: The rise of demo-genetic agent-based models. Evol Appl 2022; 15:1988-2001. [PMID: 36540635 PMCID: PMC9753837 DOI: 10.1111/eva.13508] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 10/26/2022] [Accepted: 10/28/2022] [Indexed: 11/29/2022] Open
Abstract
The study of eco-evolutionary dynamics, that is of the intertwinning between ecological and evolutionary processes when they occur at comparable time scales, is of growing interest in the current context of global change. However, many eco-evolutionary studies overlook the role of interindividual interactions, which are hard to predict and yet central to selective values. Here, we aimed at putting forward models that simulate interindividual interactions in an eco-evolutionary framework: the demo-genetic agent-based models (DG-ABMs). Being demo-genetic, DG-ABMs consider the feedback loop between ecological and evolutionary processes. Being agent-based, DG-ABMs follow populations of interacting individuals with sets of traits that vary among the individuals. We argue that the ability of DG-ABMs to take into account the genetic heterogeneity-that affects individual decisions/traits related to local and instantaneous conditions-differentiates them from analytical models, another type of model largely used by evolutionary biologists to investigate eco-evolutionary feedback loops. Based on the review of studies employing DG-ABMs and explicitly or implicitly accounting for competitive, cooperative or reproductive interactions, we illustrate that DG-ABMs are particularly relevant for the exploration of fundamental, yet pressing, questions in evolutionary ecology across various levels of organization. By jointly modelling the effects of management practices and other eco-evolutionary processes on interindividual interactions and population dynamics, DG-ABMs are also effective prospective and decision support tools to evaluate the short- and long-term evolutionary costs and benefits of management strategies and to assess potential trade-offs. Finally, we provide a list of the recent practical advances of the ABM community that should facilitate the development of DG-ABMs.
Collapse
Affiliation(s)
- Amaïa Lamarins
- E2S UPPA, INRAE, ECOBIOPUniversité de Pau et des Pays de l'AdourSaint‐Pée‐sur‐NivelleFrance,Management of Diadromous Fish in their Environment, OFB, INRAE, Institut AgroUniv Pau & Pays Adour/E2S UPPARennesFrance
| | - Victor Fririon
- INRAE, UR 629 Ecologie des Forêts Méditerranéennes, URFMAvignonFrance
| | - Dorinda Folio
- E2S UPPA, INRAE, ECOBIOPUniversité de Pau et des Pays de l'AdourSaint‐Pée‐sur‐NivelleFrance
| | - Camille Vernier
- CIRAD, UMR CBGP, INRAE, IRD, Montpellier SupAgroUniv. MontpellierMontpellierFrance
| | - Léa Daupagne
- E2S UPPA, INRAE, ECOBIOPUniversité de Pau et des Pays de l'AdourSaint‐Pée‐sur‐NivelleFrance
| | - Jacques Labonne
- E2S UPPA, INRAE, ECOBIOPUniversité de Pau et des Pays de l'AdourSaint‐Pée‐sur‐NivelleFrance
| | - Mathieu Buoro
- E2S UPPA, INRAE, ECOBIOPUniversité de Pau et des Pays de l'AdourSaint‐Pée‐sur‐NivelleFrance
| | - François Lefèvre
- INRAE, UR 629 Ecologie des Forêts Méditerranéennes, URFMAvignonFrance
| | - Cyril Piou
- CIRAD, UMR CBGP, INRAE, IRD, Montpellier SupAgroUniv. MontpellierMontpellierFrance
| | - Sylvie Oddou‐Muratorio
- E2S UPPA, INRAE, ECOBIOPUniversité de Pau et des Pays de l'AdourSaint‐Pée‐sur‐NivelleFrance
| |
Collapse
|
4
|
Feeding ecology has shaped the evolution of modern sharks. Curr Biol 2021; 31:5138-5148.e4. [PMID: 34614390 DOI: 10.1016/j.cub.2021.09.028] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Revised: 06/05/2021] [Accepted: 09/09/2021] [Indexed: 11/20/2022]
Abstract
Sharks are iconic predators in today's oceans, yet their modern diversity has ancient origins. In particular, present hypotheses suggest that a combination of mass extinction, global climate change, and competition has regulated the community structure of dominant mackerel (Lamniformes) and ground (Carcharhiniformes) sharks over the last 66 million years. However, while these scenarios advocate an interplay of major abiotic and biotic events, the precise drivers remain obscure. Here, we focus on the role of feeding ecology using a geometric morphometric analysis of 3,837 fossil and extant shark teeth. Our results reveal that morphological segregation rather than competition has characterized lamniform and carcharhiniform evolution. Moreover, although lamniforms suffered a long-term disparity decline potentially linked to dietary "specialization," their recent disparity rivals that of "generalist" carcharhiniforms. We further confirm that low eustatic sea levels impacted lamniform disparity across the end-Cretaceous mass extinction. Adaptations to changing prey availability and the proliferation of coral reef habitats during the Paleogene also likely facilitated carcharhiniform dispersals and cladogenesis, underpinning their current taxonomic dominance. Ultimately, we posit that trophic partitioning and resource utilization shaped past shark ecology and represent critical determinants for their future species survivorship.
Collapse
|
5
|
Stevenson DS, Wallace R. Biogeographical Modeling of Alien Worlds. ASTROBIOLOGY 2021; 21:831-844. [PMID: 33904766 DOI: 10.1089/ast.2020.2304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In this article, we partially quantify the biological potential of an exoplanet. We employ a variety of biogeographical analyses, placing biological evolution in the context of the geological evolution of the planet as a whole. Terrestrial (as in Earthly) biodiversity is tightly constrained in terms of species richness by its environment. An organism's habitable environment may be considered its niche space or hypervolume in terms of the physical characteristics in which that organism can survive and reproduce. This fundamental niche forms the broader space in which the organism realizes its true niche in terms of its interactions with other species. Many of the physical characteristics can be determined from astrophysical constraints and are thus amenable for dissection. However, the geographical space that organisms occupy is driven by the geological evolution of a sizable telluric planet. In turn, this is driven by the progressive differentiation of its interior to produce increasingly felsic crust. Using a variety of available models, we can then constrain the available space that species can inhabit using species-area relationships. By considering a combination of astrophysical constraints and geographical space, we partially quantify the numbers of species that can inhabit the landscape that geology provides. Finally, we also identify a correlation between geomorphological scale and speciation, which, if validated, will allow further dissection of species diversity on alien worlds.
Collapse
Affiliation(s)
- David S Stevenson
- Science Department, Carlton le Willows Academy, Nottingham, United Kingdom
| | - Rodrick Wallace
- Division of Epidemiology, The New York State Psychiatric Institute at Columbia University, New York Psychiatric Institute, New York, USA
| |
Collapse
|
6
|
Abstract
Genealogical tree modeling is essential for estimating evolutionary parameters in population genetics and phylogenetics. Recent mathematical results concerning ranked genealogies without leaf labels unlock opportunities in the analysis of evolutionary trees. In particular, comparisons between ranked genealogies facilitate the study of evolutionary processes of different organisms sampled at multiple time periods. We propose metrics on ranked tree shapes and ranked genealogies for lineages isochronously and heterochronously sampled. Our proposed tree metrics make it possible to conduct statistical analyses of ranked tree shapes and timed ranked tree shapes or ranked genealogies. Such analyses allow us to assess differences in tree distributions, quantify estimation uncertainty, and summarize tree distributions. We show the utility of our metrics via simulations and an application in infectious diseases.
Collapse
Affiliation(s)
- Jaehee Kim
- Department of Biology, Stanford University, Stanford, CA 94305
| | | | - Julia A Palacios
- Department of Statistics, Stanford University, Stanford, CA 94305;
- Department of Biomedical Data Science, Stanford School of Medicine, Stanford, CA 94305
| |
Collapse
|
7
|
Fraser D, Soul LC, Tóth AB, Balk MA, Eronen JT, Pineda-Munoz S, Shupinski AB, Villaseñor A, Barr WA, Behrensmeyer AK, Du A, Faith JT, Gotelli NJ, Graves GR, Jukar AM, Looy CV, Miller JH, Potts R, Lyons SK. Investigating Biotic Interactions in Deep Time. Trends Ecol Evol 2020; 36:61-75. [PMID: 33067015 DOI: 10.1016/j.tree.2020.09.001] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2020] [Revised: 08/31/2020] [Accepted: 09/02/2020] [Indexed: 11/16/2022]
Abstract
Recent renewed interest in using fossil data to understand how biotic interactions have shaped the evolution of life is challenging the widely held assumption that long-term climate changes are the primary drivers of biodiversity change. New approaches go beyond traditional richness and co-occurrence studies to explicitly model biotic interactions using data on fossil and modern biodiversity. Important developments in three primary areas of research include analysis of (i) macroevolutionary rates, (ii) the impacts of and recovery from extinction events, and (iii) how humans (Homo sapiens) affected interactions among non-human species. We present multiple lines of evidence for an important and measurable role of biotic interactions in shaping the evolution of communities and lineages on long timescales.
Collapse
Affiliation(s)
- Danielle Fraser
- Palaeobiology, Canadian Museum of Nature, Ottawa, ON, Canada; Biology and Earth Sciences, Carleton University, Ottawa, ON, Canada; Department of Paleobiology and Evolution of Terrestrial Ecosystems Program, Smithsonian Institution, National Museum of Natural History, Washington, DC , USA.
| | - Laura C Soul
- Department of Paleobiology and Evolution of Terrestrial Ecosystems Program, Smithsonian Institution, National Museum of Natural History, Washington, DC , USA
| | - Anikó B Tóth
- Centre for Ecosystem Science, School of Biological, Earth and Environmental Sciences, UNSW, Sydney, NSW, Australia
| | - Meghan A Balk
- Department of Paleobiology and Evolution of Terrestrial Ecosystems Program, Smithsonian Institution, National Museum of Natural History, Washington, DC , USA; BIO5 Institute, University of Arizona, Tucson, AZ, USA
| | - Jussi T Eronen
- Ecosystems and Environment Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland; Helsinki Institute of Sustainability Science, Faculty of Biological and Environmental Sciences, Ecosystems and Environment Research Programme, University of Helsinki, Helsinki, Finland; BIOS research Unit, Helsinki, Finland
| | - Silvia Pineda-Munoz
- Department of Earth and Atmospheric Sciences, Georgia Institute of Technology, Atlanta, GA, USA
| | | | - Amelia Villaseñor
- Department of Anthropology, University of Arkansas, Fayetteville, AR, USA
| | - W Andrew Barr
- Department of Paleobiology and Evolution of Terrestrial Ecosystems Program, Smithsonian Institution, National Museum of Natural History, Washington, DC , USA; Center for the Advanced Study of Human Paleobiology, Department of Anthropology, The George Washington University, Washington, DC, USA
| | - Anna K Behrensmeyer
- Department of Paleobiology and Evolution of Terrestrial Ecosystems Program, Smithsonian Institution, National Museum of Natural History, Washington, DC , USA
| | - Andrew Du
- Department of Anthropology and Geography, Colorado State University, Fort Collins, CO, USA
| | - J Tyler Faith
- Natural History Museum of Utah, University of Utah, Salt Lake City, UT,USA; Department of Anthropology, University of Utah, Salt Lake City, UT, USA
| | | | - Gary R Graves
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, USA; Center for Macroecology, Evolution and Climate, GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
| | - Advait M Jukar
- Department of Paleobiology and Evolution of Terrestrial Ecosystems Program, Smithsonian Institution, National Museum of Natural History, Washington, DC , USA
| | - Cindy V Looy
- Department of Integrative Biology, Museum of Paleontology, University and Jepson Herbaria, University of California-Berkeley, Berkeley, CA , USA
| | - Joshua H Miller
- Department of Geology, University of Cincinnati, Cincinnati, OH, USA
| | - Richard Potts
- Human Origins Program, National Museum of Natural History, Smithsonian Institution, Washington, DC , USA
| | - S Kathleen Lyons
- School of Biological Sciences, University of Nebraska-Lincoln, Lincoln, NE, USA
| |
Collapse
|
8
|
Princepe D, De Aguiar MAM. Modeling Mito-nuclear Compatibility and Its Role in Species Identification. Syst Biol 2020; 70:133-144. [PMID: 32497198 DOI: 10.1093/sysbio/syaa044] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Revised: 05/26/2020] [Accepted: 05/27/2020] [Indexed: 01/27/2023] Open
Abstract
Mitochondrial genetic material (mtDNA) is widely used for phylogenetic reconstruction and as a barcode for species identification. The utility of mtDNA in these contexts derives from its particular molecular properties, including its high evolutionary rate, uniparental inheritance, and small size. But mtDNA may also play a fundamental role in speciation-as suggested by recent observations of coevolution with the nuclear DNA, along with the fact that respiration depends on coordination of genes from both sources. Here, we study how mito-nuclear interactions affect the accuracy of species identification by mtDNA, as well as the speciation process itself. We simulate the evolution of a population of individuals who carry a recombining nuclear genome and a mitochondrial genome inherited maternally. We compare a null model fitness landscape that lacks any mito-nuclear interaction against a scenario in which interactions influence fitness. Fitness is assigned to individuals according to their mito-nuclear compatibility, which drives the coevolution of the nuclear and mitochondrial genomes. Depending on the model parameters, the population breaks into distinct species and the model output then allows us to analyze the accuracy of mtDNA barcode for species identification. Remarkably, we find that species identification by mtDNA is equally accurate in the presence or absence of mito-nuclear coupling and that the success of the DNA barcode derives mainly from population geographical isolation during speciation. Nevertheless, selection imposed by mito-nuclear compatibility influences the diversification process and leaves signatures in the genetic content and spatial distribution of the populations, in three ways. First, speciation is delayed and the resulting phylogenetic trees are more balanced. Second, clades in the resulting phylogenetic tree correlate more strongly with the spatial distribution of species and clusters of more similar mtDNA's. Third, there is a substantial increase in the intraspecies mtDNA similarity, decreasing the number of alleles substitutions per locus and promoting the conservation of genetic information. We compare the evolutionary patterns observed in our model to empirical data from copepods (Tigriopus californicus). We find good qualitative agreement in the geographic patterns and the topology of the phylogenetic tree, provided the model includes selection based on mito-nuclear interactions. These results highlight the role of mito-nuclear compatibility in the speciation process and its reconstruction from genetic data.[Mito-nuclear coevolution; mtDNA barcode; parapatry; phylogeny.].
Collapse
Affiliation(s)
| | - Marcus A M De Aguiar
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas - 13083-859, Campinas, SP, Brazil
| |
Collapse
|
9
|
Verboom GA, Boucher FC, Ackerly DD, Wootton LM, Freyman WA. Species Selection Regime and Phylogenetic Tree Shape. Syst Biol 2020; 69:774-794. [DOI: 10.1093/sysbio/syz076] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Revised: 11/08/2019] [Accepted: 11/11/2019] [Indexed: 11/12/2022] Open
Abstract
Abstract
Species selection, the effect of heritable traits in generating between-lineage diversification rate differences, provides a valuable conceptual framework for understanding the relationship between traits, diversification, and phylogenetic tree shape. An important challenge, however, is that the nature of real diversification landscapes—curves or surfaces which describe the propensity of species-level lineages to diversify as a function of one or more traits—remains poorly understood. Here, we present a novel, time-stratified extension of the QuaSSE model in which speciation/extinction rate is specified as a static or temporally shifting Gaussian or skewed-Gaussian function of the diversification trait. We then use simulations to show that the generally imbalanced nature of real phylogenetic trees, as well as their generally greater than expected frequency of deep branching events, are typical outcomes when diversification is treated as a dynamic, trait-dependent process. Focusing on four basic models (Gaussian-speciation with and without background extinction; skewed-speciation; Gaussian-extinction), we also show that particular features of the species selection regime produce distinct tree shape signatures and that, consequently, a combination of tree shape metrics has the potential to reveal the species selection regime under which a particular lineage diversified. We evaluate this idea empirically by comparing the phylogenetic trees of plant lineages diversifying within climatically and geologically stable environments of the Greater Cape Floristic Region, with those of lineages diversifying in environments that have experienced major change through the Late Miocene-Pliocene. Consistent with our expectations, the trees of lineages diversifying in a dynamic context are less balanced, show a greater concentration of branching events close to the present, and display stronger diversification rate-trait correlations. We suggest that species selection plays an important role in shaping phylogenetic trees but recognize the need for an explicit probabilistic framework within which to assess the likelihoods of alternative diversification scenarios as explanations of a particular tree shape. [Cape flora; diversification landscape; environmental change; gamma statistic; species selection; time-stratified QuaSSE model; trait-dependent diversification; tree imbalance.]
Collapse
Affiliation(s)
- G Anthony Verboom
- Bolus Herbarium and Department of Biological Sciences, University of Cape Town, Private Bag, Rondebosch 7700, South Africa
| | - Florian C Boucher
- Department of Botany and Zoology, University of Stellenbosch, Private Bag X1, Matieland 7602, South Africa
- Université Grenoble Alpes, CNRS, Laboratoire d’Ecologie Alpine (LECA), 2233 Rue de la Piscine, FR-38000 Grenoble, France
| | - David D Ackerly
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
- Department of Environmental Sciences, Policy, and Management, University of California, Berkeley, CA 94720, USA
| | - Lara M Wootton
- Bolus Herbarium and Department of Biological Sciences, University of Cape Town, Private Bag, Rondebosch 7700, South Africa
| | - William A Freyman
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
- Department of Ecology, Evolution and Behavior, University of Minnesota, Saint Paul, MN 55108, USA
| |
Collapse
|
10
|
Chang JM, Floden EW, Herrero J, Gascuel O, Di Tommaso P, Notredame C. Incorporating alignment uncertainty into Felsenstein's phylogenetic bootstrap to improve its reliability. Bioinformatics 2019; 37:1506-1514. [PMID: 30726875 PMCID: PMC8275982 DOI: 10.1093/bioinformatics/btz082] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Revised: 12/12/2018] [Accepted: 02/05/2019] [Indexed: 12/30/2022] Open
Abstract
Motivation Most evolutionary analyses are based on pre-estimated multiple sequence alignment. Wong et al. established the existence of an uncertainty induced by multiple sequence alignment when reconstructing phylogenies. They were able to show that in many cases different aligners produce different phylogenies, with no simple objective criterion sufficient to distinguish among these alternatives. Results We demonstrate that incorporating MSA induced uncertainty into bootstrap sampling can significantly increase correlation between clade correctness and its corresponding bootstrap value. Our procedure involves concatenating several alternative multiple sequence alignments of the same sequences, produced using different commonly used aligners. We then draw bootstrap replicates while favoring columns of the more unique aligner among the concatenated aligners. We named this concatenation and bootstrapping method, Weighted Partial Super Bootstrap (wpSBOOT). We show on three simulated datasets of 16, 32 and 64 tips that our method improves the predictive power of bootstrap values. We also used as a benchmark an empirical collection of 853 1-to-1 orthologous genes from seven yeast species and found wpSBOOT to significantly improve discrimination capacity between topologically correct and incorrect trees. Bootstrap values of wpSBOOT are comparable to similar readouts estimated using a single method. However, for reduced trees by 50% and 95% bootstrap thresholds, wpSBOOT comes out the lowest Type I error (less FP). Availability The automated generation of replicates has been implemented in the T-Coffee package, which is available as open source freeware available from www.tcoffee.org. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Jia-Ming Chang
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SD, United Kingdom
| | - Evan W Floden
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Javier Herrero
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SD, United Kingdom
| | - Olivier Gascuel
- Unité Bioinformatique Evolutive, Centre de Bioinformatique, Biostatistique et Biologie Intégrative (C3BI)-USR 3756 CNRS and Institut Pasteur, Paris, France
| | - Paolo Di Tommaso
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Cedric Notredame
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
11
|
Costa CLN, Lemos-Costa P, Marquitti FMD, Fernandes LD, Ramos MF, Schneider DM, Martins AB, de Aguiar MAM. Signatures of Microevolutionary Processes in Phylogenetic Patterns. Syst Biol 2018; 68:131-144. [PMID: 29939352 DOI: 10.1093/sysbio/syy049] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2017] [Accepted: 06/13/2018] [Indexed: 11/13/2022] Open
Abstract
Phylogenetic trees are representations of evolutionary relationships among species and contain signatures of the processes responsible for the speciation events they display. Inferring processes from tree properties, however, is challenging. To address this problem, we analyzed a spatially-explicit model of speciation where genome size and mating range can be controlled. We simulated parapatric and sympatric (narrow and wide mating range, respectively) radiations and constructed their phylogenetic trees, computing structural properties such as tree balance and speed of diversification. We showed that parapatric and sympatric speciation are well separated by these structural tree properties. Balanced trees with constant rates of diversification only originate in sympatry and genome size affected both the balance and the speed of diversification of the simulated trees. Comparison with empirical data showed that most of the evolutionary radiations considered to have developed in parapatry or sympatry are in good agreement with model predictions. Even though additional forces other than spatial restriction of gene flow, genome size, and genetic incompatibilities, do play a role in the evolution of species formation, the microevolutionary processes modeled here capture signatures of the diversification pattern of evolutionary radiations, regarding the symmetry and speed of diversification of lineages.
Collapse
Affiliation(s)
- Carolina L N Costa
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil.,Instituto de Biologia, Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil
| | - Paula Lemos-Costa
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil.,Instituto de Biologia, Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil
| | - Flavia M D Marquitti
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil
| | - Lucas D Fernandes
- Departamento de Entomologia e Acarologia, Escola Superior de Agricultura Luiz de Queiroz, Universidade de São Paulo, 13418-900 Piracicaba, SP, Brazil
| | - Marlon F Ramos
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil
| | - David M Schneider
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil
| | - Ayana B Martins
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil.,Department of Fish Ecology & Evolution, Centre of Ecology, Evolution and Biogeochemistry, Eawag Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
| | - Marcus A M de Aguiar
- Instituto de Física 'Gleb Wataghin', Universidade Estadual de Campinas, Unicamp, 13083-970 Campinas, SP, Brazil
| |
Collapse
|
12
|
The Species Problem from the Modeler's Point of View. Bull Math Biol 2018; 81:878-898. [PMID: 30535845 DOI: 10.1007/s11538-018-00536-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2018] [Accepted: 11/27/2018] [Indexed: 01/12/2023]
Abstract
How to define a partition of individuals into species is a long-standing question called the species problem in systematics. Here, we focus on this problem in the thought experiment where individuals reproduce clonally and both the differentiation process and the population genealogies are explicitly known. We specify three desirable properties of species partitions: (A) Heterotypy between species, (B) Homotypy within species and (M) Genealogical monophyly of each species. We then ask: How and when is it possible to delineate species in a way satisfying these properties? We point out that the three desirable properties cannot in general be satisfied simultaneously, but that any two of them can. We mathematically prove the existence of the finest partition satisfying (A) and (M) and the coarsest partition satisfying (B) and (M). For each of them, we propose a simple algorithm to build the associated phylogeny out of the genealogy. The ways we propose to phrase the species problem shed new light on the interaction between the genealogical and phylogenetic scales in modeling work. The two definitions centered on the monophyly property can readily be used at a higher taxonomic level as well, e.g., to cluster species into monophyletic genera.
Collapse
|
13
|
Herrera-Alsina L, Pigot AL, Hildenbrandt H, Etienne RS. The influence of ecological and geographic limits on the evolution of species distributions and diversity. Evolution 2018; 72:1978-1991. [PMID: 30055007 PMCID: PMC6220796 DOI: 10.1111/evo.13563] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2017] [Revised: 07/17/2018] [Accepted: 07/18/2018] [Indexed: 01/11/2023]
Abstract
The role of ecological limits in regulating the distribution and diversification of species remains controversial. Although such limits must ultimately arise from constraints on local species coexistence, this spatial context is missing from most macroevolutionary models. Here, we develop a stochastic, spatially explicit model of species diversification to explore the phylogenetic and biogeographic patterns expected when local diversity is bounded. We show how local ecological limits, by regulating opportunities for range expansion and thus rates of speciation and extinction, lead to temporal slowdowns in diversification and predictable differences in equilibrium diversity between regions. However, our models also show that even when regions have identical diversity limits, the dynamics of diversification and total number of species supported at equilibrium can vary dramatically depending on the relative size of geographic and local ecological niche space. Our model predicts that small regions with higher local ecological limits support a higher standing diversity and more balanced phylogenetic trees than large geographic areas with more stringent constraints on local coexistence. Our findings highlight how considering the spatial context of diversification can provide new insights into the role of ecological limits in driving variation in biodiversity across space, time, and clades.
Collapse
Affiliation(s)
- Leonel Herrera-Alsina
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, 9700 CC, The Netherlands
| | - Alex L Pigot
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, 9700 CC, The Netherlands.,Centre for Biodiversity and Environment Research, Department of Genetics, Evolution and Environment, University College London, London, WC1E 6BT, United Kingdom
| | - Hanno Hildenbrandt
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, 9700 CC, The Netherlands
| | - Rampal S Etienne
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, 9700 CC, The Netherlands
| |
Collapse
|
14
|
Clade diversification dynamics and the biotic and abiotic controls of speciation and extinction rates. Nat Commun 2018; 9:3013. [PMID: 30068945 PMCID: PMC6070539 DOI: 10.1038/s41467-018-05419-7] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Accepted: 06/21/2018] [Indexed: 11/23/2022] Open
Abstract
How ecological interactions, genetic processes, and environmental variability jointly shape the evolution of species diversity remains a challenging problem in biology. We developed an individual-based model of clade diversification to predict macroevolutionary dynamics when resource competition, genetic differentiation, and landscape fluctuations interact. Diversification begins with a phase of geographic adaptive radiation. Extinction rates rise sharply at the onset of the next phase. In this phase of niche self-structuring, speciation and extinction processes, albeit driven by biotic mechanisms (competition and hybridization), have essentially constant rates, determined primarily by the abiotic pace of landscape dynamics. The final phase of diversification begins when intense competition prevents dispersing individuals from establishing new populations. Species’ ranges shrink, causing negative diversity-dependence of speciation rates. These results show how ecological and microevolutionary processes shape macroevolutionary dynamics and rates; they caution against the notion of ecological limits to diversity, and suggest new directions for the phylogenetic analysis of diversification. The history and patterns of species diversity are shaped by a variety of ecological and evolutionary factors. Here, the authors develop a computational model to predict clade diversification dynamics and rates of speciation and extinction under the influences of resource competition, genetic differentiation, and random landscape fluctuation.
Collapse
|
15
|
Hoppe A, Türpitz S, Steel M. Species notions that combine phylogenetic trees and phenotypic partitions. J Math Biol 2018. [PMID: 30043242 DOI: 10.1101/075580] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
A recent paper (Manceau and Lambert in bioRxiv, 2017. https://doi.org/10.1101/075580 ) developed a novel approach for describing two well-defined notions of 'species' based on a phylogenetic tree and a phenotypic partition. In this paper, we explore some further combinatorial properties of this approach and describe an extension that allows an arbitrary number of phenotypic partitions to be combined with a phylogenetic tree for these two species notions.
Collapse
Affiliation(s)
- Anica Hoppe
- Institute of Mathematics and Computer Science, Ernst-Moritz-Arndt University, Greifswald, Germany
| | - Sonja Türpitz
- Institute of Mathematics and Computer Science, Ernst-Moritz-Arndt University, Greifswald, Germany
| | - Mike Steel
- Biomathematics Research Centre, University of Canterbury, Christchurch, New Zealand.
| |
Collapse
|
16
|
Olivares I, Karger DN, Kessler M. Assessing species saturation: conceptual and methodological challenges. Biol Rev Camb Philos Soc 2018; 93:1874-1890. [PMID: 29733121 DOI: 10.1111/brv.12424] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2017] [Revised: 04/04/2018] [Accepted: 04/10/2018] [Indexed: 11/29/2022]
Abstract
Is there a maximum number of species that can coexist? Intuitively, we assume an upper limit to the number of species in a given assemblage, or that a lineage can produce, but defining and testing this limit has proven problematic. Herein, we first outline seven general challenges of studies on species saturation, most of which are independent of the actual method used to assess saturation. Among these are the challenge of defining saturation conceptually and operationally, the importance of setting an appropriate referential system, and the need to discriminate among patterns, processes and mechanisms. Second, we list and discuss the methodological approaches that have been used to study species saturation. These approaches vary in time and spatial scales, and in the variables and assumptions needed to assess saturation. We argue that assessing species saturation is possible, but that many studies conducted to date have conceptual and methodological flaws that prevent us from currently attaining a good idea of the occurrence of species saturation.
Collapse
Affiliation(s)
- Ingrid Olivares
- Department of Systematic and Evolutionary Botany, University of Zurich, Zollikerstrasse 107, 8008 Zurich, Switzerland
| | - Dirk N Karger
- Department of Systematic and Evolutionary Botany, University of Zurich, Zollikerstrasse 107, 8008 Zurich, Switzerland.,Swiss Federal Research Institute WSL, Zürcherstrasse 111, 8903 Birmensdorf, Switzerland
| | - Michael Kessler
- Department of Systematic and Evolutionary Botany, University of Zurich, Zollikerstrasse 107, 8008 Zurich, Switzerland
| |
Collapse
|
17
|
Maliet O, Gascuel F, Lambert A. Ranked Tree Shapes, Nonrandom Extinctions, and the Loss of Phylogenetic Diversity. Syst Biol 2018; 67:1025-1040. [DOI: 10.1093/sysbio/syy030] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2016] [Accepted: 04/08/2018] [Indexed: 11/13/2022] Open
Affiliation(s)
- Odile Maliet
- Institut de Biologie de l’École Normale Supérieure (IBENS), École Normale Supérieure, CNRS, INSERM, PSL Research University, Paris, France
- ED 227, Sorbonne Universités, Paris, France
| | - Fanny Gascuel
- Institut de Biologie de l’École Normale Supérieure (IBENS), École Normale Supérieure, CNRS, INSERM, PSL Research University, Paris, France
- ED 227, Sorbonne Universités, Paris, France
- Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, PSL Research University, Paris, France
| | - Amaury Lambert
- Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, PSL Research University, Paris, France
- Laboratoire Probabilités, Statistique et Modélisation (LPSM), Sorbonne Université, CNRS, Paris, France
| |
Collapse
|
18
|
Chalmandrier L, Albouy C, Descombes P, Sandel B, Faurby S, Svenning JC, Zimmermann NE, Pellissier L. Comparing spatial diversification and meta-population models in the Indo-Australian Archipelago. ROYAL SOCIETY OPEN SCIENCE 2018; 5:171366. [PMID: 29657753 PMCID: PMC5882677 DOI: 10.1098/rsos.171366] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/12/2017] [Accepted: 01/31/2018] [Indexed: 06/08/2023]
Abstract
Reconstructing the processes that have shaped the emergence of biodiversity gradients is critical to understand the dynamics of diversification of life on Earth. Islands have traditionally been used as model systems to unravel the processes shaping biological diversity. MacArthur and Wilson's island biogeographic model predicts diversity to be based on dynamic interactions between colonization and extinction rates, while treating islands themselves as geologically static entities. The current spatial configuration of islands should influence meta-population dynamics, but long-term geological changes within archipelagos are also expected to have shaped island biodiversity, in part by driving diversification. Here, we compare two mechanistic models providing inferences on species richness at a biogeographic scale: a mechanistic spatial-temporal model of species diversification and a spatial meta-population model. While the meta-population model operates over a static landscape, the diversification model is driven by changes in the size and spatial configuration of islands through time. We compare the inferences of both models to floristic diversity patterns among land patches of the Indo-Australian Archipelago. Simulation results from the diversification model better matched observed diversity than a meta-population model constrained only by the contemporary landscape. The diversification model suggests that the dynamic re-positioning of islands promoting land disconnection and reconnection induced an accumulation of particularly high species diversity on Borneo, which is central within the island network. By contrast, the meta-population model predicts a higher diversity on the mainlands, which is less compatible with empirical data. Our analyses highlight that, by comparing models with contrasting assumptions, we can pinpoint the processes that are most compatible with extant biodiversity patterns.
Collapse
Affiliation(s)
- Loïc Chalmandrier
- Landscape Ecology, Institute of Terrestrial Ecosystems, ETH Zürich, Zurich, Switzerland
- Swiss Federal Research Institute WSL, 8903 Birmensdorf, Switzerland
| | - Camille Albouy
- Landscape Ecology, Institute of Terrestrial Ecosystems, ETH Zürich, Zurich, Switzerland
- Swiss Federal Research Institute WSL, 8903 Birmensdorf, Switzerland
| | - Patrice Descombes
- Landscape Ecology, Institute of Terrestrial Ecosystems, ETH Zürich, Zurich, Switzerland
- Swiss Federal Research Institute WSL, 8903 Birmensdorf, Switzerland
| | - Brody Sandel
- Department of Biology, Santa Clara University, 500 El Camino Real, Santa Clara, CA 95053, USA
| | - Soren Faurby
- Department of Biogeography and Global Change, Museo Nacional de Ciencias Naturales, CSIC, Madrid, Spain
- Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, SE 405 30 Gothenburg, Sweden
| | - Jens-Christian Svenning
- Section for Ecoinformatics and Biodiversity, Department of Bioscience, Aarhus University, 8000 Aarhus C, Denmark
- Center for Biodiversity Dynamics in a Changing World (BIOCHANGE), Aarhus University, Ny Munkegade 114, Aarhus, Denmark
| | | | - Loïc Pellissier
- Landscape Ecology, Institute of Terrestrial Ecosystems, ETH Zürich, Zurich, Switzerland
- Swiss Federal Research Institute WSL, 8903 Birmensdorf, Switzerland
| |
Collapse
|
19
|
Albert JS, Schoolmaster DR, Tagliacollo V, Duke-Sylvester SM. Barrier Displacement on a Neutral Landscape: Toward a Theory of Continental Biogeography. Syst Biol 2018; 66:167-182. [PMID: 27590192 DOI: 10.1093/sysbio/syw080] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2016] [Accepted: 08/23/2016] [Indexed: 01/07/2023] Open
Abstract
Macroevolutionary theory posits three processes leading to lineage diversification and the formation of regional biotas: dispersal (species geographic range expansion), speciation (species lineage splitting), and extinction (species lineage termination). The Theory of Island Biogeography (TIB) predicts species richness values using just two of these processes; dispersal and extinction. Yet most species on Earth live on continents or continental shelves, and the dynamics of evolutionary diversification at regional and continental scales are qualitatively different from those that govern the formation of species richness on biogeographic islands. Certain geomorphological processes operating perennially on continental platforms displace barriers to gene flow and organismal dispersal, and affect all three terms of macroevolutionary diversification. For example, uplift of a dissected landscape and river capture both merge and separate portions of adjacent areas, allowing dispersal and larger geographic ranges, vicariant speciation and smaller geographic ranges, and extinction when range sizes are subdivided below a minimum persistence threshold. The TIB also does not predict many biogeographic and phylogenetic patterns widely observed in continentally distributed taxa, including: (i) power function-like species-area relationships; (ii) log-normal distribution of species geographic range sizes, in which most species have restricted ranges (are endemic) and few species have broad ranges (are cosmopolitan); (iii) mid-domain effects with more species toward the geographic center, and more early-branching, species-poor clades toward the geographic periphery; (iv) exponential rates of net diversification with log-linear accumulation of lineages through geological time; and (v) power function-like relationships between species-richness and clade diversity, in which most clades are species-poor and few clades are species-rich. Current theory does not provide a robust mechanistic framework to connect these seemingly disparate patterns. Here we present SEAMLESS (Spatially Explicit Area Model of Landscape Evolution by SimulationS) that generates clade diversification by moving geographic barriers on a continuous, neutral landscape. SEAMLESS is a neutral Landscape Evolution Model (LEM) that treats species and barriers as functionally equivalent with respect to model parameters. SEAMLESS differs from other model-based biogeographic methods (e.g., Lagrange, GeoSSE, BayArea, and BioGeoBEARS) by modeling properties of dispersal barriers rather than areas, and by modeling the evolution of species lineages on a continuous landscape, rather than the evolution of geographic ranges along branches of a phylogeny. SEAMLESS shows how dispersal is required to maintain species richness and avoid clade-wide extinction, demonstrates that ancestral range size does not predict species richness, and provides a unified explanation for the suite of commonly observed biogeographic and phylogenetic patterns listed above. SEAMLESS explains how a simple barrier-displacement mechanism affects lineage diversification under neutral conditions, and is advanced here toward the formulation of a general theory of continental biogeography. [Diversification, extinction, geodispersal, macroevolution, river capture, vicariance.].
Collapse
Affiliation(s)
- James S Albert
- Department of Biology, University of Louisiana at Lafayette, 104 E. University Circle, Lafayette, LA 70503, USA
| | | | - Victor Tagliacollo
- Universidade Federal do Tocantins Avenida NS 15, 109 Norte Palmas, Tocantins 77001-090, Brazil
| | - Scott M Duke-Sylvester
- Department of Biology, University of Louisiana at Lafayette, 104 E. University Circle, Lafayette, LA 70503, USA
| |
Collapse
|
20
|
Biodiversity Dynamics on Islands: Explicitly Accounting for Causality in Mechanistic Models. DIVERSITY-BASEL 2017. [DOI: 10.3390/d9030030] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
21
|
Inferring epidemiological parameters from phylogenies using regression-ABC: A comparative study. PLoS Comput Biol 2017; 13:e1005416. [PMID: 28263987 PMCID: PMC5358897 DOI: 10.1371/journal.pcbi.1005416] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2016] [Revised: 03/20/2017] [Accepted: 02/16/2017] [Indexed: 02/06/2023] Open
Abstract
Inferring epidemiological parameters such as the R0 from time-scaled phylogenies is a timely challenge. Most current approaches rely on likelihood functions, which raise specific issues that range from computing these functions to finding their maxima numerically. Here, we present a new regression-based Approximate Bayesian Computation (ABC) approach, which we base on a large variety of summary statistics intended to capture the information contained in the phylogeny and its corresponding lineage-through-time plot. The regression step involves the Least Absolute Shrinkage and Selection Operator (LASSO) method, which is a robust machine learning technique. It allows us to readily deal with the large number of summary statistics, while avoiding resorting to Markov Chain Monte Carlo (MCMC) techniques. To compare our approach to existing ones, we simulated target trees under a variety of epidemiological models and settings, and inferred parameters of interest using the same priors. We found that, for large phylogenies, the accuracy of our regression-ABC is comparable to that of likelihood-based approaches involving birth-death processes implemented in BEAST2. Our approach even outperformed these when inferring the host population size with a Susceptible-Infected-Removed epidemiological model. It also clearly outperformed a recent kernel-ABC approach when assuming a Susceptible-Infected epidemiological model with two host types. Lastly, by re-analyzing data from the early stages of the recent Ebola epidemic in Sierra Leone, we showed that regression-ABC provides more realistic estimates for the duration parameters (latency and infectiousness) than the likelihood-based method. Overall, ABC based on a large variety of summary statistics and a regression method able to perform variable selection and avoid overfitting is a promising approach to analyze large phylogenies. Given the rapid evolution of many pathogens, analysing their genomes by means of phylogenies can inform us about how they spread. This is the focus of the field known as “phylodynamics”. Most existing methods inferring epidemiological parameters from virus phylogenies are limited by the difficulty of handling complex likelihood functions, which commonly incorporate latent variables. Here, we use an alternative method known as regression-based Approximate Bayesian Computation (ABC), which circumvents this problem by using simulations and dataset comparisons. Since phylogenies are difficult to compare to one another, we introduce many summary statistics to describe them and take advantage of current machine learning techniques able to perform variable selection. We show that the accuracy we reach is comparable to that of existing methods. This accuracy increases with phylogeny size and can even be higher than that of existing methods for some parameters. Overall, regression-based ABC opens new perspectives to infer epidemiological parameters from large phylogenies.
Collapse
|
22
|
Weber MG, Wagner CE, Best RJ, Harmon LJ, Matthews B. Evolution in a Community Context: On Integrating Ecological Interactions and Macroevolution. Trends Ecol Evol 2017; 32:291-304. [PMID: 28215448 DOI: 10.1016/j.tree.2017.01.003] [Citation(s) in RCA: 105] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2016] [Revised: 01/10/2017] [Accepted: 01/11/2017] [Indexed: 11/28/2022]
Abstract
Despite a conceptual understanding that evolution and species interactions are inextricably linked, it remains challenging to study ecological and evolutionary dynamics together over long temporal scales. In this review, we argue that, despite inherent challenges associated with reconstructing historical processes, the interplay of ecology and evolution is central to our understanding of macroevolution and community coexistence, and cannot be safely ignored in community and comparative phylogenetic studies. We highlight new research avenues that foster greater consideration of both ecological and evolutionary dynamics as processes that occur along branches of phylogenetic trees. By promoting new ways forward using this perspective, we hope to inspire further integration that creatively co-utilizes phylogenies and ecological data to study eco-evolutionary dynamics over time and space.
Collapse
Affiliation(s)
- Marjorie G Weber
- Department of Plant Biology, Michigan State University, East Lansing, MI 48823, USA.
| | - Catherine E Wagner
- Biodiversity Institute and Department of Botany, University of Wyoming, 1000 E. University Avenue, Laramie, WY 82071, USA
| | - Rebecca J Best
- Eawag, Department of Aquatic Ecology, Center for Ecology, Evolution and Biogeochemistry, Kastanienbaum, Switzerland; School of Earth Sciences and Environmental Sustainability, Northern Arizona University, 525 S. Beaver Street, Flagstaff, AZ 86011, USA
| | - Luke J Harmon
- Eawag, Department of Aquatic Ecology, Center for Ecology, Evolution and Biogeochemistry, Kastanienbaum, Switzerland; Department of Biological Sciences, University of Idaho, Moscow, ID, USA
| | - Blake Matthews
- Eawag, Department of Aquatic Ecology, Center for Ecology, Evolution and Biogeochemistry, Kastanienbaum, Switzerland
| |
Collapse
|
23
|
Doebeli M, Ispolatov Y, Simon B. Towards a mechanistic foundation of evolutionary theory. eLife 2017; 6. [PMID: 28198700 PMCID: PMC5333952 DOI: 10.7554/elife.23804] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2016] [Accepted: 02/13/2017] [Indexed: 11/17/2022] Open
Abstract
Most evolutionary thinking is based on the notion of fitness and related ideas such as fitness landscapes and evolutionary optima. Nevertheless, it is often unclear what fitness actually is, and its meaning often depends on the context. Here we argue that fitness should not be a basal ingredient in verbal or mathematical descriptions of evolution. Instead, we propose that evolutionary birth-death processes, in which individuals give birth and die at ever-changing rates, should be the basis of evolutionary theory, because such processes capture the fundamental events that generate evolutionary dynamics. In evolutionary birth-death processes, fitness is at best a derived quantity, and owing to the potential complexity of such processes, there is no guarantee that there is a simple scalar, such as fitness, that would describe long-term evolutionary outcomes. We discuss how evolutionary birth-death processes can provide useful perspectives on a number of central issues in evolution. DOI:http://dx.doi.org/10.7554/eLife.23804.001
Collapse
Affiliation(s)
- Michael Doebeli
- Department of Zoology and Department of Mathematics, University of British Columbia, Vancouver, Canada
| | - Yaroslav Ispolatov
- Departamento de Fisica, Universidad de Santiago de Chile, Santiago, Chile
| | - Burt Simon
- Department of Mathematical and Statistical Sciences, University of Colorado, Denver, United States
| |
Collapse
|
24
|
Missa O, Dytham C, Morlon H. Understanding how biodiversity unfolds through time under neutral theory. Philos Trans R Soc Lond B Biol Sci 2016; 371:20150226. [PMID: 26977066 DOI: 10.1098/rstb.2015.0226] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Theoretical predictions for biodiversity patterns are typically derived under the assumption that ecological systems have reached a dynamic equilibrium. Yet, there is increasing evidence that various aspects of ecological systems, including (but not limited to) species richness, are not at equilibrium. Here, we use simulations to analyse how biodiversity patterns unfold through time. In particular, we focus on the relative time required for various biodiversity patterns (macroecological or phylogenetic) to reach equilibrium. We simulate spatially explicit metacommunities according to the Neutral Theory of Biodiversity (NTB) under three modes of speciation, which differ in how evenly a parent species is split between its two daughter species. We find that species richness stabilizes first, followed by species area relationships (SAR) and finally species abundance distributions (SAD). The difference in timing of equilibrium between these different macroecological patterns is the largest when the split of individuals between sibling species at speciation is the most uneven. Phylogenetic patterns of biodiversity take even longer to stabilize (tens to hundreds of times longer than species richness) so that equilibrium predictions from neutral theory for these patterns are unlikely to be relevant. Our results suggest that it may be unwise to assume that biodiversity patterns are at equilibrium and provide a first step in studying how these patterns unfold through time.
Collapse
Affiliation(s)
- Olivier Missa
- Institute of Biology, Ecole Normale Supérieure, 46 rue d'Ulm, Paris 75005, France
| | - Calvin Dytham
- Biology Department, University of York, Wentworth Way, York YO10 5DD, UK
| | - Hélène Morlon
- Institute of Biology, Ecole Normale Supérieure, 46 rue d'Ulm, Paris 75005, France
| |
Collapse
|
25
|
Doebeli M, Ispolatov I. Diversity and Coevolutionary Dynamics in High-Dimensional Phenotype Spaces. Am Nat 2016; 189:105-120. [PMID: 28107053 DOI: 10.1086/689891] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
Abstract
We study macroevolutionary dynamics by extending microevolutionary competition models to long timescales. It has been shown that for a general class of competition models, gradual evolutionary change in continuous phenotypes (evolutionary dynamics) can be nonstationary and even chaotic when the dimension of the phenotype space in which the evolutionary dynamics unfold is high. It has also been shown that evolutionary diversification can occur along nonequilibrium trajectories in phenotype space. We combine these lines of thinking by studying long-term coevolutionary dynamics of emerging lineages in multidimensional phenotype spaces. We use a statistical approach to investigate the evolutionary dynamics of many different systems. We find (1) that, for a given dimension of phenotype space, the coevolutionary dynamics tend to be fast and nonstationary for an intermediate number of coexisting lineages but tend to stabilize as the evolving communities reach a saturation level of diversity and (2) that the amount of diversity at the saturation level increases rapidly (exponentially) with the dimension of phenotype space. These results have implications for theoretical perspectives on major macroevolutionary patterns such as adaptive radiation, long-term temporal patterns of phenotypic changes, and the evolution of diversity.
Collapse
|
26
|
Integrating selection, niche, and diversification into a hierarchical conceptual framework. ORG DIVERS EVOL 2016. [DOI: 10.1007/s13127-016-0299-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|