1
|
Zhang R, Drummond AJ, Mendes FK. Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters. Syst Biol 2024; 73:102-124. [PMID: 38085256 PMCID: PMC11129596 DOI: 10.1093/sysbio/syad067] [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: 11/06/2023] [Revised: 03/23/2023] [Accepted: 11/07/2023] [Indexed: 05/28/2024] Open
Abstract
Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.
Collapse
Affiliation(s)
- Rong Zhang
- Programme in Emerging Infectious Diseases, Duke-NUS Medical School 169857, Singapore
| | - Alexei J Drummond
- Centre for Computational Evolution, The University of Auckland, Auckland 1010, New Zealand
- School of Biological Sciences, The University of Auckland, Auckland 1010, New Zealand
| | - Fábio K Mendes
- Department of Biology, Washington University in St. Louis, St. Louis, MO 63130, USA
| |
Collapse
|
2
|
Shao Y, Magee AF, Vasylyeva TI, Suchard MA. Scalable gradients enable Hamiltonian Monte Carlo sampling for phylodynamic inference under episodic birth-death-sampling models. PLoS Comput Biol 2024; 20:e1011640. [PMID: 38551979 PMCID: PMC11006205 DOI: 10.1371/journal.pcbi.1011640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 04/10/2024] [Accepted: 03/10/2024] [Indexed: 04/09/2024] Open
Abstract
Birth-death models play a key role in phylodynamic analysis for their interpretation in terms of key epidemiological parameters. In particular, models with piecewise-constant rates varying at different epochs in time, to which we refer as episodic birth-death-sampling (EBDS) models, are valuable for their reflection of changing transmission dynamics over time. A challenge, however, that persists with current time-varying model inference procedures is their lack of computational efficiency. This limitation hinders the full utilization of these models in large-scale phylodynamic analyses, especially when dealing with high-dimensional parameter vectors that exhibit strong correlations. We present here a linear-time algorithm to compute the gradient of the birth-death model sampling density with respect to all time-varying parameters, and we implement this algorithm within a gradient-based Hamiltonian Monte Carlo (HMC) sampler to alleviate the computational burden of conducting inference under a wide variety of structures of, as well as priors for, EBDS processes. We assess this approach using three different real world data examples, including the HIV epidemic in Odesa, Ukraine, seasonal influenza A/H3N2 virus dynamics in New York state, America, and Ebola outbreak in West Africa. HMC sampling exhibits a substantial efficiency boost, delivering a 10- to 200-fold increase in minimum effective sample size per unit-time, in comparison to a Metropolis-Hastings-based approach. Additionally, we show the robustness of our implementation in both allowing for flexible prior choices and in modeling the transmission dynamics of various pathogens by accurately capturing the changing trend of viral effective reproductive number.
Collapse
Affiliation(s)
- Yucai Shao
- Department of Biostatistics, University of California, Los Angeles, California, United States of America
| | - Andrew F. Magee
- Department of Biomathematics, University of California, Los Angeles, California, United States of America
| | - Tetyana I. Vasylyeva
- Department of Medicine, University of California San Diego, La Jolla, California, United States of America
- Department of Population Health and Disease Prevention, University of California Irvine, Irvine, California, United States of America
| | - Marc A. Suchard
- Department of Biostatistics, University of California, Los Angeles, California, United States of America
- Department of Biomathematics, University of California, Los Angeles, California, United States of America
- Department of Human Genetics, Universtiy of California, Los Angeles, California, United States of America
| |
Collapse
|
3
|
Lambert S, Voznica J, Morlon H. Deep Learning from Phylogenies for Diversification Analyses. Syst Biol 2023; 72:1262-1279. [PMID: 37556735 DOI: 10.1093/sysbio/syad044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 06/20/2023] [Accepted: 08/08/2023] [Indexed: 08/11/2023] Open
Abstract
Birth-death (BD) models are widely used in combination with species phylogenies to study past diversification dynamics. Current inference approaches typically rely on likelihood-based methods. These methods are not generalizable, as a new likelihood formula must be established each time a new model is proposed; for some models, such a formula is not even tractable. Deep learning can bring solutions in such situations, as deep neural networks can be trained to learn the relation between simulations and parameter values as a regression problem. In this paper, we adapt a recently developed deep learning method from pathogen phylodynamics to the case of diversification inference, and we extend its applicability to the case of the inference of state-dependent diversification models from phylogenies associated with trait data. We demonstrate the accuracy and time efficiency of the approach for the time-constant homogeneous BD model and the Binary-State Speciation and Extinction model. Finally, we illustrate the use of the proposed inference machinery by reanalyzing a phylogeny of primates and their associated ecological role as seed dispersers. Deep learning inference provides at least the same accuracy as likelihood-based inference while being faster by several orders of magnitude, offering a promising new inference approach for the deployment of future models in the field.
Collapse
Affiliation(s)
- Sophia Lambert
- Institut de Biologie de l'École Normale Supérieure, École Normale Supérieure, CNRS, INSERM, Université Paris Sciences et Lettres, 46 Rue d'Ulm, 75005 Paris, France
- Institute of Ecology and Evolution, Department of Biology, 5289 University of Oregon, Eugene, OR 97403, USA
| | - Jakub Voznica
- Institut Pasteur, Université Paris Cité, Unité Bioinformatique Evolutive, 25-28 Rue du Dr Roux, 75015 Paris, France
- Unité de Biologie Computationnelle, USR 3756 CNRS, 25-28 Rue du Dr Roux, 75015 Paris, France
| | - Hélène Morlon
- Institut de Biologie de l'École Normale Supérieure, École Normale Supérieure, CNRS, INSERM, Université Paris Sciences et Lettres, 46 Rue d'Ulm, 75005 Paris, France
| |
Collapse
|
4
|
Shao Y, Magee AF, Vasylyeva TI, Suchard MA. Scalable gradients enable Hamiltonian Monte Carlo sampling for phylodynamic inference under episodic birth-death-sampling models. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.31.564882. [PMID: 37961423 PMCID: PMC10634968 DOI: 10.1101/2023.10.31.564882] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Birth-death models play a key role in phylodynamic analysis for their interpretation in terms of key epidemiological parameters. In particular, models with piecewise-constant rates varying at different epochs in time, to which we refer as episodic birth-death-sampling (EBDS) models, are valuable for their reflection of changing transmission dynamics over time. A challenge, however, that persists with current time-varying model inference procedures is their lack of computational efficiency. This limitation hinders the full utilization of these models in large-scale phylodynamic analyses, especially when dealing with high-dimensional parameter vectors that exhibit strong correlations. We present here a linear-time algorithm to compute the gradient of the birth-death model sampling density with respect to all time-varying parameters, and we implement this algorithm within a gradient-based Hamiltonian Monte Carlo (HMC) sampler to alleviate the computational burden of conducting inference under a wide variety of structures of, as well as priors for, EBDS processes. We assess this approach using three different real world data examples, including the HIV epidemic in Odesa, Ukraine, seasonal influenza A/H3N2 virus dynamics in New York state, America, and Ebola outbreak in West Africa. HMC sampling exhibits a substantial efficiency boost, delivering a 10- to 200-fold increase in minimum effective sample size per unit-time, in comparison to a Metropolis-Hastings-based approach. Additionally, we show the robustness of our implementation in both allowing for flexible prior choices and in modeling the transmission dynamics of various pathogens by accurately capturing the changing trend of viral effective reproductive number.
Collapse
Affiliation(s)
- Yucai Shao
- Department of Biostatistics, Jonathan and Karin Fielding School of Public Health, University of California, Los Angeles, United States
| | - Andrew F. Magee
- Department of Biomathematics, David Geffen School of Medicine at UCLA, University of California, Los Angeles, United States
| | - Tetyana I. Vasylyeva
- Department of Medicine, University of California San Diego, La Jolla, United States
- Department of Population Health and Disease Prevention, University of California Irvine, Irvine, United States
| | - Marc A. Suchard
- Department of Biostatistics, Jonathan and Karin Fielding School of Public Health, University of California, Los Angeles, United States
- Department of Biomathematics, David Geffen School of Medicine at UCLA, University of California, Los Angeles, United States
- Department of Human Genetics, David Geffen School of Medicine at UCLA, Universtiy of California, Los Angeles, United States
| |
Collapse
|
5
|
Villaverde T, Larridon I, Shah T, Fowler RM, Chau JH, Olmstead RG, Sanmartín I. Phylogenomics sheds new light on the drivers behind a long-lasting systematic riddle: the figwort family Scrophulariaceae. THE NEW PHYTOLOGIST 2023; 240:1601-1615. [PMID: 36869601 DOI: 10.1111/nph.18845] [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: 10/13/2022] [Accepted: 02/23/2023] [Indexed: 06/18/2023]
Abstract
The figwort family, Scrophulariaceae, comprises c. 2000 species whose evolutionary relationships at the tribal level have proven difficult to resolve, hindering our ability to understand their origin and diversification. We designed a specific probe kit for Scrophulariaceae, targeting 849 nuclear loci and obtaining plastid regions as by-products. We sampled c. 87% of the genera described in the family and use the nuclear dataset to estimate evolutionary relationships, timing of diversification, and biogeographic patterns. Ten tribes, including two new tribes, Androyeae and Camptolomeae, are supported, and the phylogenetic positions of Androya, Camptoloma, and Phygelius are unveiled. Our study reveals a major diversification at c. 60 million yr ago in some Gondwanan landmasses, where two different lineages diversified, one of which gave rise to nearly 81% of extant species. A Southern African origin is estimated for most modern-day tribes, with two exceptions, the American Leucophylleae, and the mainly Australian Myoporeae. The rapid mid-Eocene diversification is aligned with geographic expansion within southern Africa in most tribes, followed by range expansion to tropical Africa and multiple dispersals out of Africa. Our robust phylogeny provides a framework for future studies aimed at understanding the role of macroevolutionary patterns and processes that generated Scrophulariaceae diversity.
Collapse
Affiliation(s)
- Tamara Villaverde
- Real Jardín Botánico (CSIC), Plaza de Murillo, 2, Madrid, 28014, Spain
| | - Isabel Larridon
- Royal Botanic Gardens, Kew, Richmond, TW9 3AE, UK
- Systematic and Evolutionary Botany Lab, Department of Biology, Ghent University, K.L. Ledeganckstraat 35, 9000, Ghent, Belgium
| | - Toral Shah
- Royal Botanic Gardens, Kew, Richmond, TW9 3AE, UK
- Department of Life Sciences, Imperial College London, Silwood Park Campus, Ascot, SL5 7PY, UK
| | - Rachael M Fowler
- School of BioSciences, The University of Melbourne, Parkville, Vic., 3010, Australia
| | - John H Chau
- Department of Zoology, Centre for Ecological Genomics and Wildlife Conservation, University of Johannesburg, Auckland Park, 2006, South Africa
| | - Richard G Olmstead
- Department of Biology and Burke Museum, University of Washington, Seattle, WA, 98155, USA
| | - Isabel Sanmartín
- Real Jardín Botánico (CSIC), Plaza de Murillo, 2, Madrid, 28014, Spain
| |
Collapse
|
6
|
Kopperud BT, Magee AF, Höhna S. Rapidly changing speciation and extinction rates can be inferred in spite of nonidentifiability. Proc Natl Acad Sci U S A 2023; 120:e2208851120. [PMID: 36757894 PMCID: PMC9963352 DOI: 10.1073/pnas.2208851120] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Accepted: 12/01/2022] [Indexed: 02/10/2023] Open
Abstract
The birth-death model is commonly used to infer speciation and extinction rates by fitting the model to phylogenetic trees with exclusively extant taxa. Recently, it was demonstrated that speciation and extinction rates are not identifiable if the rates are allowed to vary freely over time. The group of birth-death models that have the same likelihood is called a congruence class, and there is no statistical evidence to favor one model over the other. This issue has led researchers to question if and what patterns can reliably be inferred from phylogenies of only extant taxa and whether time-variable birth-death models should be fitted at all. We explore the congruence class in the context of several empirical phylogenies as well as hypothetical scenarios. For these empirical phylogenies, we assume that we inferred the true congruence class. Thus, our conclusions apply to any empirical phylogeny for which we robustly inferred the true congruence class. When we summarize shared patterns in the congruence class, we show that strong directional trends in speciation and extinction rates are shared among most models. Therefore, we conclude that the inference of strong directional trends is robust. Conversely, estimates of constant rates or gentle slopes are not robust and must be treated with caution. Interestingly, the space of valid speciation rates is narrower and more limited in contrast to extinction rates, which are less constrained. These results provide further evidence and insights that speciation rates can be estimated more reliably than extinction rates.
Collapse
Affiliation(s)
- Bjørn T. Kopperud
- GeoBio-Center LMU, Ludwig-Maximilians-Universität München, Munich80333, Germany
- Department of Earth and Environmental Sciences, Palaeontology & Geobiology, Ludwig-Maximilians-Universität München, Munich80333, Germany
| | - Andrew F. Magee
- Department of Biostatistics, University of California, Los Angeles, CA90095
| | - Sebastian Höhna
- GeoBio-Center LMU, Ludwig-Maximilians-Universität München, Munich80333, Germany
- Department of Earth and Environmental Sciences, Palaeontology & Geobiology, Ludwig-Maximilians-Universität München, Munich80333, Germany
| |
Collapse
|
7
|
Höhna S, Kopperud BT, Magee AF. CRABS: Congruent rate analyses in birth–death scenarios. Methods Ecol Evol 2022. [DOI: 10.1111/2041-210x.13997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Sebastian Höhna
- GeoBio‐Center LMU Ludwig‐Maximilians‐Universität München Munich Germany
- Department of Earth and Environmental Sciences, Paleontology & Geobiology Ludwig‐Maximilians‐Universität München Munich Germany
| | - Bjørn T. Kopperud
- GeoBio‐Center LMU Ludwig‐Maximilians‐Universität München Munich Germany
- Department of Earth and Environmental Sciences, Paleontology & Geobiology Ludwig‐Maximilians‐Universität München Munich Germany
| | - Andrew F. Magee
- Department of Human Genetics University of California Los Angeles California USA
| |
Collapse
|
8
|
Roycroft E, Moritz C, Rowe KC, Moussalli A, Eldridge MDB, Portela Miguez R, Piggott MP, Potter S. Sequence Capture From Historical Museum Specimens: Maximizing Value for Population and Phylogenomic Studies. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.931644] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
The application of high-throughput, short-read sequencing to degraded DNA has greatly increased the feasibility of generating genomic data from historical museum specimens. While many published studies report successful sequencing results from historical specimens; in reality, success and quality of sequence data can be highly variable. To examine predictors of sequencing quality, and methodological approaches to improving data accuracy, we generated and analyzed genomic sequence data from 115 historically collected museum specimens up to 180 years old. Data span both population genomic and phylogenomic scales, including historically collected specimens from 34 specimens of four species of Australian rock-wallabies (genus Petrogale) and 92 samples from 79 specimens of Australo-Papuan murine rodents (subfamily Murinae). For historical rodent specimens, where the focus was sampling for phylogenomics, we found that regardless of specimen age, DNA sequence libraries prepared from toe pad or bone subsamples performed significantly better than those taken from the skin (in terms of proportion of reads on target, number of loci captured, and data accuracy). In total, 93% of DNA libraries from toe pad or bone subsamples resulted in reliable data for phylogenetic inference, compared to 63% of skin subsamples. For skin subsamples, proportion of reads on target weakly correlated with collection year. Then using population genomic data from rock-wallaby skins as a test case, we found substantial improvement in final data quality by mapping to a high-quality “closest sister” de novo assembly from fresh tissues, compared to mapping to a sample-specific historical de novo assembly. Choice of mapping approach also affected final estimates of the number of segregating sites and Watterson's θ, both important parameters for population genomic inference. The incorporation of accurate and reliable sequence data from historical specimens has important outcomes for evolutionary studies at both population and phylogenomic scales. By assessing the outcomes of different approaches to specimen subsampling, library preparation and bioinformatic processing, our results provide a framework for increasing sequencing success for irreplaceable historical specimens.
Collapse
|
9
|
González-Miguéns R, Soler-Zamora C, Useros F, Nogal-Prata S, Berney C, Blanco-Rotea A, Carrasco-Braganza MI, de Salvador-Velasco D, Guillén-Oterino A, Tenorio-Rodríguez D, Velázquez D, Heger TJ, Sanmartín I, Lara E. Cyphoderia ampulla (Cyphoderiidae: Rhizaria), a tale of freshwater sailors. The causes and consequences of ecological transitions through the salinity barrier in a family of benthic protists. Mol Ecol 2022; 31:2644-2663. [PMID: 35262986 PMCID: PMC9311665 DOI: 10.1111/mec.16424] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 02/17/2022] [Accepted: 03/02/2022] [Indexed: 11/29/2022]
Abstract
The salinity barrier that separates marine and freshwater biomes is probably the most important division in biodiversity on Earth. Those organisms that successfully performed this transition had access to new ecosystems while undergoing changes in selective pressure, which often led to major shifts in diversification rates. While these transitions have been extensively investigated in animals, the tempo, mode, and outcome of crossing the salinity barrier have been scarcely studied in other eukaryotes. Here, we reconstructed the evolutionary history of the species complex Cyphoderia ampulla (Euglyphida: Cercozoa: Rhizaria) based on DNA sequences from the nuclear SSU rRNA gene and the mitochondrial cytochrome oxidase subunit I gene, obtained from publicly available environmental DNA data (GeneBank, EukBank) and isolated organisms. A tree calibrated with euglyphid fossils showed that four independent transitions towards freshwater systems occurred from the Mid Miocene onwards, coincident with important fluctuations in sea level. Ancestral trait reconstructions indicated that the whole family Cyphoderiidae had a marine origin and suggest that ancestors of the freshwater forms were euryhaline and lived in environments with fluctuating salinity. Diversification rates did not show any obvious increase concomitant with ecological transitions, but morphometric analyses indicated that species increased in size and homogenized their morphology after colonizing the new environments. This suggests adaptation to changes in selective pressure exerted by life in freshwater sediments.
Collapse
Affiliation(s)
| | - Carmen Soler-Zamora
- Real Jardín Botánico de Madrid (RJB-CSIC), Plaza Murillo 2, 28014, Madrid, Spain
| | - Fernando Useros
- Real Jardín Botánico de Madrid (RJB-CSIC), Plaza Murillo 2, 28014, Madrid, Spain
| | - Sandra Nogal-Prata
- Real Jardín Botánico de Madrid (RJB-CSIC), Plaza Murillo 2, 28014, Madrid, Spain
| | - Cédric Berney
- Université de la Sorbonne CNRS, Station Biologique de Roscoff, UMR 7144, ECOMAP, 29680, Roscoff, France.,Research Federation for the study of Global Ocean Systems Ecology and Evolution, FR2022/Tara GOSEE, 10, Paris, France
| | - Andrés Blanco-Rotea
- Estación Biológica Internacional Duero-Douro, (EUROPARQUES-EBI), Buque hidrográfico Helios-Cousteau en el Lago de Sanabria, 49632, Ribadelago, Castilla y León, Spain
| | - María Isabel Carrasco-Braganza
- Estación Biológica Internacional Duero-Douro, (EUROPARQUES-EBI), Buque hidrográfico Helios-Cousteau en el Lago de Sanabria, 49632, Ribadelago, Castilla y León, Spain
| | - David de Salvador-Velasco
- Estación Biológica Internacional Duero-Douro, (EUROPARQUES-EBI), Buque hidrográfico Helios-Cousteau en el Lago de Sanabria, 49632, Ribadelago, Castilla y León, Spain
| | - Antonio Guillén-Oterino
- Estación Biológica Internacional Duero-Douro, (EUROPARQUES-EBI), Buque hidrográfico Helios-Cousteau en el Lago de Sanabria, 49632, Ribadelago, Castilla y León, Spain
| | - Daniel Tenorio-Rodríguez
- Estación Biológica Internacional Duero-Douro, (EUROPARQUES-EBI), Buque hidrográfico Helios-Cousteau en el Lago de Sanabria, 49632, Ribadelago, Castilla y León, Spain
| | - David Velázquez
- Dpt. of Biology, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Thierry J Heger
- Soil Science and Environment Group, CHANGINS, University of Applied Sciences and Arts Western Switzerland, Route de Duillier 50, 1260, Nyon, Switzerland
| | - Isabel Sanmartín
- Real Jardín Botánico de Madrid (RJB-CSIC), Plaza Murillo 2, 28014, Madrid, Spain
| | - Enrique Lara
- Real Jardín Botánico de Madrid (RJB-CSIC), Plaza Murillo 2, 28014, Madrid, Spain
| |
Collapse
|
10
|
Palazzesi L, Hidalgo O, Barreda VD, Forest F, Höhna S. The rise of grasslands is linked to atmospheric CO 2 decline in the late Palaeogene. Nat Commun 2022; 13:293. [PMID: 35022396 PMCID: PMC8755714 DOI: 10.1038/s41467-021-27897-y] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2017] [Accepted: 12/13/2021] [Indexed: 01/25/2023] Open
Abstract
Grasslands are predicted to experience a major biodiversity change by the year 2100. A better understanding of how grasslands have responded to past environmental changes will help predict the outcome of current and future environmental changes. Here, we explore the relationship between past atmospheric CO2 and temperature fluctuations and the shifts in diversification rate of Poaceae (grasses) and Asteraceae (daisies), two exceptionally species-rich grassland families (~11,000 and ~23,000 species, respectively). To this end, we develop a Bayesian approach that simultaneously estimates diversification rates through time from time-calibrated phylogenies and correlations between environmental variables and diversification rates. Additionally, we present a statistical approach that incorporates the information of the distribution of missing species in the phylogeny. We find strong evidence supporting a simultaneous increase in diversification rates for grasses and daisies after the most significant reduction of atmospheric CO2 in the Cenozoic (~34 Mya). The fluctuations of paleo-temperatures, however, appear not to have had a significant relationship with the diversification of these grassland families. Overall, our results shed new light on our understanding of the origin of grasslands in the context of past environmental changes. A better understanding of how grasslands have responded to past environmental changes will help predict the outcomes of future changes. This study explores past climatic fluctuations and shifts in the diversification rate of grasses and daisies, finding strong evidence for a simultaneous increase in their diversification rates following a reduction of atmospheric CO2 in the Cenozoic.
Collapse
Affiliation(s)
- Luis Palazzesi
- Museo Argentino de Ciencias Naturales & Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, C1405DJR, Argentina. .,Jodrell Laboratory, Royal Botanic Gardens, Kew, Richmond, Surrey, TW9 3DS, UK.
| | - Oriane Hidalgo
- Jodrell Laboratory, Royal Botanic Gardens, Kew, Richmond, Surrey, TW9 3DS, UK.,Institut Botánic de Barcelona (IBB, CSIC-Ajuntament de Barcelona), Catalonia, Spain
| | - Viviana D Barreda
- Museo Argentino de Ciencias Naturales & Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, C1405DJR, Argentina
| | - Félix Forest
- Jodrell Laboratory, Royal Botanic Gardens, Kew, Richmond, Surrey, TW9 3DS, UK
| | - Sebastian Höhna
- GeoBio-Center, Ludwig-Maximilians-Universität München, Richard-Wagner-Str. 10, 80333, Munich, Germany. .,Department of Earth and Environmental Sciences, Paleontology & Geobiology, Ludwig-Maximilians-Universität München, Richard-Wagner-Str. 10, 80333, Munich, Germany.
| |
Collapse
|
11
|
Höhna S, Landis MJ, Huelsenbeck JP. Parallel power posterior analyses for fast computation of marginal likelihoods in phylogenetics. PeerJ 2021; 9:e12438. [PMID: 34760401 PMCID: PMC8570164 DOI: 10.7717/peerj.12438] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 10/15/2021] [Indexed: 11/30/2022] Open
Abstract
In Bayesian phylogenetic inference, marginal likelihoods can be estimated using several different methods, including the path-sampling or stepping-stone-sampling algorithms. Both algorithms are computationally demanding because they require a series of power posterior Markov chain Monte Carlo (MCMC) simulations. Here we introduce a general parallelization strategy that distributes the power posterior MCMC simulations and the likelihood computations over available CPUs. Our parallelization strategy can easily be applied to any statistical model despite our primary focus on molecular substitution models in this study. Using two phylogenetic example datasets, we demonstrate that the runtime of the marginal likelihood estimation can be reduced significantly even if only two CPUs are available (an average performance increase of 1.96x). The performance increase is nearly linear with the number of available CPUs. We record a performance increase of 13.3x for cluster nodes with 16 CPUs, representing a substantial reduction to the runtime of marginal likelihood estimations. Hence, our parallelization strategy enables the estimation of marginal likelihoods to complete in a feasible amount of time which previously needed days, weeks or even months. The methods described here are implemented in our open-source software RevBayes which is available from http://www.RevBayes.com.
Collapse
Affiliation(s)
- Sebastian Höhna
- GeoBio-Center, Ludwig-Maximilians-Universität München, Munich, Germany.,Department of Earth and Environmental Sciences, Paleontology & Geobiology, Ludwig-Maximilians- Universität München, Munich, Germany
| | - Michael J Landis
- Department of Biology, Washington University in St. Louis, St. Louis, United States of America
| | - John P Huelsenbeck
- Department of Integrative Biology, University of California,, Berkeley, United States of America
| |
Collapse
|
12
|
Harmon LJ, Pennell MW, Henao-Diaz LF, Rolland J, Sipley BN, Uyeda JC. Causes and Consequences of Apparent Timescaling Across All Estimated Evolutionary Rates. ANNUAL REVIEW OF ECOLOGY, EVOLUTION, AND SYSTEMATICS 2021. [DOI: 10.1146/annurev-ecolsys-011921-023644] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Evolutionary rates play a central role in connecting micro- and macroevolution. All evolutionary rate estimates, including rates of molecular evolution, trait evolution, and lineage diversification, share a similar scaling pattern with time: The highest rates are those measured over the shortest time interval. This creates a disconnect between micro- and macroevolution, although the pattern is the opposite of what some might expect: Patterns of change over short timescales predict that evolution has tremendous potential to create variation and that potential is barely tapped by macroevolution. In this review, we discuss this shared scaling pattern across evolutionary rates. We break down possible explanations for scaling into two categories, estimation error and model misspecification, and discuss how both apply to each type of rate. We also discuss the consequences of this ubiquitous pattern, which can lead to unexpected results when comparing ratesover different timescales. Finally, after addressing purely statistical concerns, we explore a few possibilities for a shared unifying explanation across the three types of rates that results from a failure to fully understand and account for how biological processes scale over time.
Collapse
Affiliation(s)
- Luke J. Harmon
- Institute for Bioinformatics and Evolutionary Studies (IBEST) and Department of Biological Sciences, University of Idaho, Moscow, Idaho 83844, USA
| | - Matthew W. Pennell
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - L. Francisco Henao-Diaz
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - Jonathan Rolland
- Laboratoire Evolution et Diversité Biologique, CNRS, UMR5174, Université Toulouse III–Paul Sabatier, 31062 Toulouse, France
| | - Breanna N. Sipley
- Program for Bioinformatics and Computational Biology, University of Idaho, Moscow, Idaho 83844, USA
| | - Josef C. Uyeda
- Department of Biological Sciences, Virginia Tech University, Blacksburg, Virginia 24061, USA
| |
Collapse
|
13
|
Carruthers T, Scotland RW. The implications of interrelated assumptions on estimates of divergence times and rates of diversification. Syst Biol 2021; 70:1181-1199. [PMID: 33760070 DOI: 10.1093/sysbio/syab021] [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: 04/29/2020] [Revised: 03/16/2021] [Accepted: 03/22/2021] [Indexed: 11/15/2022] Open
Abstract
Phylogenies are increasingly being used as a basis to provide insight into macroevolutionary history. Here, we use simulation experiments and empirical analyses to evaluate methods that use phylogenies as a basis to make estimates of divergence times and rates of diversification. This is the first study to present a comprehensive assessment of the key variables that underpin analyses in this field - including substitution rates, speciation rates, and extinction, plus character sampling and taxon sampling. We show that in unrealistically simplistic cases (where substitution rates and speciation rates are constant, and where there is no extinction), increased character and taxon sampling lead to more accurate and precise parameter estimates. By contrast, in more complex but realistic cases (where substitution rates, speciation rates, and extinction rates vary), gains in accuracy and precision from increased character and taxon sampling are far more limited. The lack of accuracy and precision even occurs when using methods that are designed to account for more complex cases, such as relaxed clocks, fossil calibrations, and models that allow speciation rates and extinction rates to vary. The problem also persists when analysing genomic scale datasets. These results suggest two interrelated problems that occur when the processes that generated the data are more complex. First, methodological assumptions are more likely to be violated. Second, limitations in the information content of the data become more important.
Collapse
Affiliation(s)
- Tom Carruthers
- Royal Botanic Gardens Kew, Richmond, London, TW9 3AE, United Kingdom
| | - Robert W Scotland
- Department of Plant Sciences, University of Oxford, South Parks Road, Oxford, OX1 3RB, United Kingdom
| |
Collapse
|
14
|
Hay EM, Poulin R, Jorge F. Macroevolutionary dynamics of parasite diversification: A reality check. J Evol Biol 2020; 33:1758-1769. [PMID: 33047407 DOI: 10.1111/jeb.13714] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2020] [Revised: 09/10/2020] [Accepted: 09/21/2020] [Indexed: 12/21/2022]
Abstract
Parasitism is often invoked as a factor explaining the variation in diversification rates across the tree of life, while also representing up to half of Earth's diversity. Yet, patterns and processes of parasite diversification remain mostly unknown. In this study, we assess the patterns of parasite diversification and specifically determine the role of life-history traits (i.e. life cycle complexity and host range) and major coevolutionary events in driving diversification across eight phylogenetic datasets spanning taxonomically different parasite groups. Aware of the degree of incomplete sampling among all parasite phylogenies, we also tested the impact of sampling bias on estimates of diversification. We show that the patterns and rates of parasite diversification differ among taxa according to life cycle complexity and to some extent major host transitions. Only directly transmitted parasites were found to be influenced by an effect of major host transitions on diversification rates. Although parasitism may be a main factor responsible for heterogeneity in diversification among the tree of life, the high degree of incomplete parasite phylogenies remains an obstacle when modelling diversification dynamics. Nevertheless, we provide the first comparative test of parasite diversification, revealing some consistent patterns and insight into the processes that shape it.
Collapse
Affiliation(s)
- Eleanor M Hay
- School of Biological Sciences, Monash University, Clayton, VIC, Australia
| | - Robert Poulin
- Department of Zoology, University of Otago, Dunedin, New Zealand
| | - Fátima Jorge
- Department of Zoology, University of Otago, Dunedin, New Zealand
| |
Collapse
|
15
|
Magee AF, Höhna S, Vasylyeva TI, Leaché AD, Minin VN. Locally adaptive Bayesian birth-death model successfully detects slow and rapid rate shifts. PLoS Comput Biol 2020; 16:e1007999. [PMID: 33112848 PMCID: PMC7652323 DOI: 10.1371/journal.pcbi.1007999] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 11/09/2020] [Accepted: 05/28/2020] [Indexed: 11/18/2022] Open
Abstract
Birth-death processes have given biologists a model-based framework to answer questions about changes in the birth and death rates of lineages in a phylogenetic tree. Therefore birth-death models are central to macroevolutionary as well as phylodynamic analyses. Early approaches to studying temporal variation in birth and death rates using birth-death models faced difficulties due to the restrictive choices of birth and death rate curves through time. Sufficiently flexible time-varying birth-death models are still lacking. We use a piecewise-constant birth-death model, combined with both Gaussian Markov random field (GMRF) and horseshoe Markov random field (HSMRF) prior distributions, to approximate arbitrary changes in birth rate through time. We implement these models in the widely used statistical phylogenetic software platform RevBayes, allowing us to jointly estimate birth-death process parameters, phylogeny, and nuisance parameters in a Bayesian framework. We test both GMRF-based and HSMRF-based models on a variety of simulated diversification scenarios, and then apply them to both a macroevolutionary and an epidemiological dataset. We find that both models are capable of inferring variable birth rates and correctly rejecting variable models in favor of effectively constant models. In general the HSMRF-based model has higher precision than its GMRF counterpart, with little to no loss of accuracy. Applied to a macroevolutionary dataset of the Australian gecko family Pygopodidae (where birth rates are interpretable as speciation rates), the GMRF-based model detects a slow decrease whereas the HSMRF-based model detects a rapid speciation-rate decrease in the last 12 million years. Applied to an infectious disease phylodynamic dataset of sequences from HIV subtype A in Russia and Ukraine (where birth rates are interpretable as the rate of accumulation of new infections), our models detect a strongly elevated rate of infection in the 1990s.
Collapse
Affiliation(s)
- Andrew F. Magee
- Department of Biology, University of Washington, Seattle, WA, 98195, USA
| | - Sebastian Höhna
- GeoBio-Center, Ludwig-Maximilians-Universität München, 80333 Munich, Germany
- Department of Earth and Environmental Sciences, Paleontology & Geobiology, Ludwig-Maximilians-Universität München, 80333 Munich, Germany
| | | | - Adam D. Leaché
- Department of Biology, University of Washington, Seattle, WA, 98195, USA
| | - Vladimir N. Minin
- Department of Statistics, University of California, Irvine, CA, 92697, USA
| |
Collapse
|
16
|
Sun M, Folk RA, Gitzendanner MA, Soltis PS, Chen Z, Soltis DE, Guralnick RP. Estimating rates and patterns of diversification with incomplete sampling: a case study in the rosids. AMERICAN JOURNAL OF BOTANY 2020; 107:895-909. [PMID: 32519354 PMCID: PMC7384126 DOI: 10.1002/ajb2.1479] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 03/03/2020] [Indexed: 05/03/2023]
Abstract
PREMISE Recent advances in generating large-scale phylogenies enable broad-scale estimation of species diversification. These now common approaches typically are characterized by (1) incomplete species coverage without explicit sampling methodologies and/or (2) sparse backbone representation, and usually rely on presumed phylogenetic placements to account for species without molecular data. We used empirical examples to examine the effects of incomplete sampling on diversification estimation and provide constructive suggestions to ecologists and evolutionary biologists based on those results. METHODS We used a supermatrix for rosids and one well-sampled subclade (Cucurbitaceae) as empirical case studies. We compared results using these large phylogenies with those based on a previously inferred, smaller supermatrix and on a synthetic tree resource with complete taxonomic coverage. Finally, we simulated random and representative taxon sampling and explored the impact of sampling on three commonly used methods, both parametric (RPANDA and BAMM) and semiparametric (DR). RESULTS We found that the impact of sampling on diversification estimates was idiosyncratic and often strong. Compared to full empirical sampling, representative and random sampling schemes either depressed or inflated speciation rates, depending on methods and sampling schemes. No method was entirely robust to poor sampling, but BAMM was least sensitive to moderate levels of missing taxa. CONCLUSIONS We suggest caution against uncritical modeling of missing taxa using taxonomic data for poorly sampled trees and in the use of summary backbone trees and other data sets with high representative bias, and we stress the importance of explicit sampling methodologies in macroevolutionary studies.
Collapse
Affiliation(s)
- Miao Sun
- Florida Museum of Natural HistoryUniversity of FloridaGainesvilleFlorida32611USA
- State Key Laboratory of Systematic and Evolutionary BotanyInstitute of BotanyChinese Academy of SciencesBeijing100093China
- Department of BioscienceAarhus UniversityAarhus8000Denmark
| | - Ryan A. Folk
- Department of Biological SciencesMississippi State UniversityMississippi StateMississippi39762USA
| | - Matthew A. Gitzendanner
- Department of BiologyUniversity of FloridaGainesvilleFlorida32611USA
- Biodiversity InstituteUniversity of FloridaGainesvilleFlorida32611USA
| | - Pamela S. Soltis
- Florida Museum of Natural HistoryUniversity of FloridaGainesvilleFlorida32611USA
- Biodiversity InstituteUniversity of FloridaGainesvilleFlorida32611USA
- Genetics InstituteUniversity of FloridaGainesvilleFlorida32608USA
| | - Zhiduan Chen
- State Key Laboratory of Systematic and Evolutionary BotanyInstitute of BotanyChinese Academy of SciencesBeijing100093China
| | - Douglas E. Soltis
- Florida Museum of Natural HistoryUniversity of FloridaGainesvilleFlorida32611USA
- Department of BiologyUniversity of FloridaGainesvilleFlorida32611USA
- Biodiversity InstituteUniversity of FloridaGainesvilleFlorida32611USA
- Genetics InstituteUniversity of FloridaGainesvilleFlorida32608USA
| | - Robert P. Guralnick
- Florida Museum of Natural HistoryUniversity of FloridaGainesvilleFlorida32611USA
- Biodiversity InstituteUniversity of FloridaGainesvilleFlorida32611USA
| |
Collapse
|
17
|
Renner MAM, Foster CSP, Miller JT, Murphy DJ. Increased diversification rates are coupled with higher rates of climate space exploration in Australian Acacia (Caesalpinioideae). THE NEW PHYTOLOGIST 2020; 226:609-622. [PMID: 31792997 DOI: 10.1111/nph.16349] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Accepted: 11/17/2019] [Indexed: 06/10/2023]
Abstract
Australia is an excellent setting to explore relationships between climate change and diversification dynamics. Aridification since the Eocene has resulted in spectacular radiations within one or more Australian biomes. Acacia is the largest plant genus on the Australian continent, with around 1000 species, and is present in all biomes. We investigated the macroevolutionary dynamics of Acacia within climate space. We analysed phylogenetic and climatic data for 503 Acacia species to estimate a time-calibrated phylogeny and central climatic tendencies for BioClim layers from 132 000 herbarium specimens. Diversification rate heterogeneity and rates of climate space exploration were tested. We inferred two diversification rate increases, both associated with significantly higher rates of climate space exploration. Observed spikes in climate disparity within the Pleistocene correspond with onset of Pleistocene glacial-interglacial cycling. Positive time dependency in environmental disparity applies in the basal grade of Acacia, though climate space exploration rates were lower. Incongruence between rates of climate space exploration and disparity suggests different Acacia lineages have experienced different macroevolutionary processes. The second diversification rate increase is associated with a south-east Australian mesic lineage, suggesting adaptations to progressively aridifying environments and ability to transition into mesic environments contributed to Acacia's dominance across Australia.
Collapse
Affiliation(s)
- Matt A M Renner
- Royal Botanic Garden and Domain Trust, Sydney, NSW, 2000, Australia
| | - Charles S P Foster
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, 2006, Australia
| | - Joseph T Miller
- Global Biodiversity Information Facility, DK-2100, Copenhagen, Denmark
| | - Daniel J Murphy
- Royal Botanic Gardens Victoria, Melbourne, 3004, VIC, Australia
| |
Collapse
|
18
|
Lozano-Fernandez J, Tanner AR, Puttick MN, Vinther J, Edgecombe GD, Pisani D. A Cambrian-Ordovician Terrestrialization of Arachnids. Front Genet 2020; 11:182. [PMID: 32218802 PMCID: PMC7078165 DOI: 10.3389/fgene.2020.00182] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Accepted: 02/14/2020] [Indexed: 12/13/2022] Open
Abstract
Understanding the temporal context of terrestrialization in chelicerates depends on whether terrestrial groups, the traditional Arachnida, have a single origin and whether or not horseshoe crabs are primitively or secondarily marine. Molecular dating on a phylogenomic tree that recovers arachnid monophyly, constrained by 27 rigorously vetted fossil calibrations, estimates that Arachnida originated during the Cambrian or Ordovician. After the common ancestor colonized the land, the main lineages appear to have rapidly radiated in the Cambrian-Ordovician boundary interval, coinciding with high rates of molecular evolution. The highest rates of arachnid diversification are detected between the Permian and Early Cretaceous. A pattern of ancient divergence estimates for terrestrial arthropod groups in the Cambrian while the oldest fossils are Silurian (seen in both myriapods and arachnids) is mirrored in the molecular and fossil records of land plants. We suggest the discrepancy between molecular and fossil evidence for terrestrialization is likely driven by the extreme sparseness of terrestrial sediments in the rock record before the late Silurian.
Collapse
Affiliation(s)
- Jesus Lozano-Fernandez
- School of Biological Sciences, University of Bristol, Bristol, United Kingdom
- School of Earth Sciences, University of Bristol, Bristol, United Kingdom
| | - Alastair R. Tanner
- School of Biological Sciences, University of Bristol, Bristol, United Kingdom
| | - Mark N. Puttick
- Department of Biology and Biochemistry, Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Jakob Vinther
- School of Biological Sciences, University of Bristol, Bristol, United Kingdom
- School of Earth Sciences, University of Bristol, Bristol, United Kingdom
| | | | - Davide Pisani
- School of Biological Sciences, University of Bristol, Bristol, United Kingdom
- School of Earth Sciences, University of Bristol, Bristol, United Kingdom
| |
Collapse
|
19
|
Matschiner M. Selective Sampling of Species and Fossils Influences Age Estimates Under the Fossilized Birth-Death Model. Front Genet 2019; 10:1064. [PMID: 31737047 PMCID: PMC6836569 DOI: 10.3389/fgene.2019.01064] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 10/03/2019] [Indexed: 01/24/2023] Open
Abstract
The fossilized birth-death (FBD) model allows the estimation of species divergence times from molecular and fossil information in a coherent framework of diversification and fossil sampling. Some assumptions of the FBD model, however, are difficult to meet in phylogenetic analyses of highly diverse groups. Here, I use simulations to assess the impact of extreme model violations, including diversified sampling of species and the exclusive use of the oldest fossils per clade, on divergence times estimated with the FBD model. My results demonstrate that selective sampling of fossils can produce dramatically overestimated divergence times when the FBD model is used for inference, due to an interplay of underestimates for the model parameters net diversification rate, turnover, and fossil-sampling proportion. In contrast, divergence times estimated with CladeAge, a method that uses information about the oldest fossils per clade together with estimates of sampling and diversification rates, are accurate under these conditions. Practitioners of Bayesian divergence-time estimation should therefore ensure that the dataset conforms to the expectations of the FBD model, or estimates of sampling and diversification rates should be obtained a priori so that CladeAge can be used for the inference.
Collapse
Affiliation(s)
- Michael Matschiner
- Department of Palaentology and Museum, University of Zurich, Zurich, Switzerland
- Centre of Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| |
Collapse
|
20
|
Culshaw V, Stadler T, Sanmartín I. Exploring the power of Bayesian birth-death skyline models to detect mass extinction events from phylogenies with only extant taxa. Evolution 2019; 73:1133-1150. [PMID: 31017656 PMCID: PMC6767073 DOI: 10.1111/evo.13753] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2017] [Revised: 03/30/2019] [Accepted: 04/14/2019] [Indexed: 01/03/2023]
Abstract
Mass extinction events (MEEs), defined as significant losses of species diversity in significantly short time periods, have attracted the attention of biologists because of their link to major environmental change. MEEs have traditionally been studied through the fossil record, but the development of birth-death models has made it possible to detect their signature based on extant-taxa phylogenies. Most birth-death models consider MEEs as instantaneous events where a high proportion of species are simultaneously removed from the tree ("single pulse" approach), in contrast to the paleontological record, where MEEs have a time duration. Here, we explore the power of a Bayesian Birth-Death Skyline (BDSKY) model to detect the signature of MEEs through changes in extinction rates under a "time-slice" approach. In this approach, MEEs are time intervals where the extinction rate is greater than the speciation rate. Results showed BDSKY can detect and locate MEEs but that precision and accuracy depend on the phylogeny's size and MEE intensity. Comparisons of BDSKY with the single-pulse Bayesian model, CoMET, showed a similar frequency of Type II error and neither model exhibited Type I error. However, while CoMET performed better in detecting and locating MEEs for smaller phylogenies, BDSKY showed higher accuracy in estimating extinction and speciation rates.
Collapse
Affiliation(s)
- Victoria Culshaw
- Real Jardín Botánico (RJB)CSICPlaza de Murillo 228014MadridSpain
| | - Tanja Stadler
- Department of Biosystems Science and EngineeringEidgenössische Technische Hochschule Zürich4058BaselSwitzerland
| | - Isabel Sanmartín
- Real Jardín Botánico (RJB)CSICPlaza de Murillo 228014MadridSpain
| |
Collapse
|
21
|
Abstract
For centuries, biologists have been captivated by the vast disparity in species richness between different groups of organisms. Variation in diversity is widely attributed to differences between groups in how fast they speciate or go extinct. Such macroevolutionary rates have been estimated for thousands of groups and have been correlated with an incredible variety of organismal traits. Here we analyze a large collection of phylogenetic trees and fossil time series and describe a hidden generality among these seemingly idiosyncratic results: speciation and extinction rates follow a scaling law in which both depend on the age of the group in which they are measured, with the fastest rates in the youngest clades. Using a series of simulations and sensitivity analyses, we demonstrate that the time dependency is unlikely to be a result of simple statistical artifacts. As such, this time scaling is likely a genuine feature of the tree of life, hinting that the dynamics of biodiversity over deep time may be driven in part by surprisingly simple and general principles.
Collapse
|
22
|
Sarver BA, Pennell MW, Brown JW, Keeble S, Hardwick KM, Sullivan J, Harmon LJ. The choice of tree prior and molecular clock does not substantially affect phylogenetic inferences of diversification rates. PeerJ 2019; 7:e6334. [PMID: 30886768 PMCID: PMC6421065 DOI: 10.7717/peerj.6334] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2018] [Accepted: 12/23/2018] [Indexed: 11/20/2022] Open
Abstract
Comparative methods allow researchers to make inferences about evolutionary processes and patterns from phylogenetic trees. In Bayesian phylogenetics, estimating a phylogeny requires specifying priors on parameters characterizing the branching process and rates of substitution among lineages, in addition to others. Accordingly, characterizing the effect of prior selection on phylogenies is an active area of research. The choice of priors may systematically bias phylogenetic reconstruction and, subsequently, affect conclusions drawn from the resulting phylogeny. Here, we focus on the impact of priors in Bayesian phylogenetic inference and evaluate how they affect the estimation of parameters in macroevolutionary models of lineage diversification. Specifically, we simulate trees under combinations of tree priors and molecular clocks, simulate sequence data, estimate trees, and estimate diversification parameters (e.g., speciation and extinction rates) from these trees. When substitution rate heterogeneity is large, diversification rate estimates deviate substantially from those estimated under the simulation conditions when not captured by an appropriate choice of relaxed molecular clock. However, in general, we find that the choice of tree prior and molecular clock has relatively little impact on the estimation of diversification rates insofar as the sequence data are sufficiently informative and substitution rate heterogeneity among lineages is low-to-moderate.
Collapse
Affiliation(s)
- Brice A.J. Sarver
- Department of Biological Sciences and Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID, USA
| | - Matthew W. Pennell
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Joseph W. Brown
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | - Sara Keeble
- Department of Molecular and Computational Biology, University of Southern California, Los Angeles, CA, USA
| | - Kayla M. Hardwick
- Department of Biological Sciences and Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID, USA
| | - Jack Sullivan
- Department of Biological Sciences and Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID, USA
| | - Luke J. Harmon
- Department of Biological Sciences and Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID, USA
| |
Collapse
|
23
|
Ribeiro E, Davis AM, Rivero-Vega RA, Ortí G, Betancur-R R. Post-Cretaceous bursts of evolution along the benthic-pelagic axis in marine fishes. Proc Biol Sci 2018; 285:20182010. [PMID: 30963906 PMCID: PMC6304066 DOI: 10.1098/rspb.2018.2010] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Accepted: 11/21/2018] [Indexed: 01/25/2023] Open
Abstract
Ecological opportunity arising in the aftermath of mass extinction events is thought to be a powerful driver of evolutionary radiations. Here, we assessed how the wake of the Cretaceous-Palaeogene (K-Pg) mass extinction shaped diversification dynamics in a clade of mostly marine fishes (Carangaria), which comprises a disparate array of benthic and pelagic dwellers including some of the most astonishing fish forms (e.g. flatfishes, billfishes, remoras, archerfishes). Analyses of lineage diversification show time-heterogeneous rates of lineage diversification in carangarians, with highest rates reached during the Palaeocene. Likewise, a remarkable proportion of Carangaria's morphological variation originated early in the history of the group and in tandem with a marked incidence of habitat shifts. Taken together, these results suggest that all major lineages and body plans in Carangaria originated in an early burst shortly after the K-Pg mass extinction, which ultimately allowed the occupation of newly released niches along the benthic-pelagic habitat axis.
Collapse
Affiliation(s)
- Emanuell Ribeiro
- Department of Biology, University of Puerto Rico, Rio Piedras, PO Box 23360, San Juan, Puerto Rico 00931, USA
- Department of Biology, The University of Oklahoma, 730 Van Vleet Oval, Room 314, Norman, OK 73019, USA
| | - Aaron M. Davis
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, PO Box 37012, MRC 159, Washington, DC 20013-7012, USA
- Centre for Tropical Water and Aquatic Ecosystem Research (TropWATER), and School of Marine and Tropical Biology, James Cook University, Townsville, Queensland 4811, Australia
| | - Rafael A. Rivero-Vega
- Department of Biology, University of Puerto Rico, Rio Piedras, PO Box 23360, San Juan, Puerto Rico 00931, USA
| | - Guillermo Ortí
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, PO Box 37012, MRC 159, Washington, DC 20013-7012, USA
- Department of Biological Sciences, The George Washington University, 2023 G Street NW, Washington, DC 20052, USA
| | - Ricardo Betancur-R
- Department of Biology, University of Puerto Rico, Rio Piedras, PO Box 23360, San Juan, Puerto Rico 00931, USA
- Department of Biology, The University of Oklahoma, 730 Van Vleet Oval, Room 314, Norman, OK 73019, USA
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, PO Box 37012, MRC 159, Washington, DC 20013-7012, USA
| |
Collapse
|
24
|
Freyman WA, Höhna S. Stochastic Character Mapping of State-Dependent Diversification Reveals the Tempo of Evolutionary Decline in Self-Compatible Onagraceae Lineages. Syst Biol 2018; 68:505-519. [DOI: 10.1093/sysbio/syy078] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2018] [Revised: 11/05/2018] [Accepted: 11/13/2018] [Indexed: 11/13/2022] Open
Affiliation(s)
- William A Freyman
- Department of Integrative Biology, University of California, Berkeley, 3040 Valley Life Sciences Building #3140, CA 94720, USA
| | - Sebastian Höhna
- Division of Evolutionary Biology, Ludwig-Maximilians-Universität München, Geschwister-Scholl-Platz 1, 80539 Munich, Germany
| |
Collapse
|
25
|
May MR, Höhna S, Moore BR. A Bayesian approach for detecting the impact of mass‐extinction events on molecular phylogenies when rates of lineage diversification may vary. Methods Ecol Evol 2016. [DOI: 10.1111/2041-210x.12563] [Citation(s) in RCA: 77] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Michael R. May
- Department of Evolution and Ecology University of California Davis CA95616USA
| | - Sebastian Höhna
- Department of Evolution and Ecology University of California Davis CA95616USA
- Department of Integrative Biology University of California Berkeley CA94720USA
- Department of Statistics University of California Berkeley CA94720USA
- Department of Mathematics Stockholm University Stockholm SE‐106 91Sweden
| | - Brian R. Moore
- Department of Evolution and Ecology University of California Davis CA95616USA
| |
Collapse
|
26
|
May MR, Moore BR. How Well Can We Detect Lineage-Specific Diversification-Rate Shifts? A Simulation Study of Sequential AIC Methods. Syst Biol 2016; 65:1076-1084. [PMID: 27037081 PMCID: PMC5066061 DOI: 10.1093/sysbio/syw026] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2014] [Accepted: 03/25/2016] [Indexed: 11/26/2022] Open
Abstract
Evolutionary biologists have long been fascinated by the extreme differences in species numbers across branches of the Tree of Life. This has motivated the development of statistical methods for detecting shifts in the rate of lineage diversification across the branches of phylogenic trees. One of the most frequently used methods, MEDUSA, explores a set of diversification-rate models, where each model assigns branches of the phylogeny to a set of diversification-rate categories. Each model is first fit to the data, and the Akaike information criterion (AIC) is then used to identify the optimal diversification model. Surprisingly, the statistical behavior of this popular method is uncharacterized, which is a concern in light of: (1) the poor performance of the AIC as a means of choosing among models in other phylogenetic contexts; (2) the ad hoc algorithm used to visit diversification models, and; (3) errors that we reveal in the likelihood function used to fit diversification models to the phylogenetic data. Here, we perform an extensive simulation study demonstrating that MEDUSA (1) has a high false-discovery rate (on average, spurious diversification-rate shifts are identified ≈30% of the time), and (2) provides biased estimates of diversification-rate parameters. Understanding the statistical behavior of MEDUSA is critical both to empirical researchers—in order to clarify whether these methods can make reliable inferences from empirical datasets—and to theoretical biologists—in order to clarify the specific problems that need to be solved in order to develop more reliable approaches for detecting shifts in the rate of lineage diversification. [Akaike information criterion; extinction; lineage-specific diversification rates; phylogenetic model selection; speciation.]
Collapse
Affiliation(s)
- Michael R May
- Department of Evolution and Ecology, University of California, Davis, Davis, CA 95616, USA
| | - Brian R Moore
- Department of Evolution and Ecology, University of California, Davis, Davis, CA 95616, USA
| |
Collapse
|
27
|
Zhao J, Teufel AI, Liberles DA, Liu L. A generalized birth and death process for modeling the fates of gene duplication. BMC Evol Biol 2015; 15:275. [PMID: 26643106 PMCID: PMC4672517 DOI: 10.1186/s12862-015-0539-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2015] [Accepted: 11/10/2015] [Indexed: 01/15/2023] Open
Abstract
Background Accurately estimating the timing and mode of gene duplications along the evolutionary history of species can provide invaluable information about underlying mechanisms by which the genomes of organisms evolved and the genes with novel functions arose. Mechanistic models have previously been introduced that allow for probabilistic inference of the evolutionary mechanism for duplicate gene retention based upon the average rate of loss over time of the duplicate. However, there is currently no probabilistic model embedded in a birth-death modeling framework that can take into account the effects of different evolutionary mechanisms of gene retention when analyzing gene family data. Results In this study, we describe a generalized birth-death process for modeling the fates of gene duplication. Use of mechanistic models in a phylogenetic framework requires an age-dependent birth-death process. Starting with a single population corresponding to the lineage of a phylogenetic tree and with an assumption of a clock that starts ticking for each duplicate at its birth, an age-dependent birth-death process is developed by extending the results from the time-dependent birth-death process. The implementation of such models in a full phylogenetic framework is expected to enable large scale probabilistic analysis of duplicates in comparative genomic studies. Conclusions We develop an age-dependent birth-death model for understanding the mechanisms of gene retention, which allows a gene loss rate dependent on each duplication event. Simulation results indicate that different mechanisms of gene retentions produce distinct likelihood functions, which can be used with genomic data to quantitatively distinguish those mechanisms.
Collapse
Affiliation(s)
- Jing Zhao
- Department of Statistics, University of Georgia, 101 Cedar Street, Athens, GA, 30602, USA.
| | - Ashley I Teufel
- Department of Molecular Biology, University of Wyoming, Laramie, WY, 82071, USA. .,Center for Computational Genetics and Genomics and Department of Biology, Temple University, Philadelphia, PA, 19122, USA.
| | - David A Liberles
- Department of Molecular Biology, University of Wyoming, Laramie, WY, 82071, USA. .,Center for Computational Genetics and Genomics and Department of Biology, Temple University, Philadelphia, PA, 19122, USA.
| | - Liang Liu
- Department of Statistics, University of Georgia, 101 Cedar Street, Athens, GA, 30602, USA. .,Institute of Bioinformatics, University of Georgia, Athens, GA, 30602, USA.
| |
Collapse
|
28
|
Höhna S, May MR, Moore BR. TESS: an R package for efficiently simulating phylogenetic trees and performing Bayesian inference of lineage diversification rates. Bioinformatics 2015; 32:789-91. [DOI: 10.1093/bioinformatics/btv651] [Citation(s) in RCA: 95] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2015] [Accepted: 10/30/2015] [Indexed: 11/12/2022] Open
|
29
|
Höhna S. The time-dependent reconstructed evolutionary process with a key-role for mass-extinction events. J Theor Biol 2015; 380:321-31. [DOI: 10.1016/j.jtbi.2015.06.005] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2014] [Revised: 03/04/2015] [Accepted: 06/02/2015] [Indexed: 10/23/2022]
|
30
|
Janzen T, Höhna S, Etienne RS. Approximate Bayesian Computation of diversification rates from molecular phylogenies: introducing a new efficient summary statistic, the
nLTT. Methods Ecol Evol 2015. [DOI: 10.1111/2041-210x.12350] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Thijs Janzen
- Groningen Institute for Evolutionary Life Sciences University of Groningen 9700 CC Groningen the Netherlands
| | - Sebastian Höhna
- Department of Evolution and Ecology University of California Davis, Storer Hall, One Shields Avenue Davis CA 95616 USA
| | - Rampal S. Etienne
- Groningen Institute for Evolutionary Life Sciences University of Groningen 9700 CC Groningen the Netherlands
| |
Collapse
|
31
|
Höhna S, Heath TA, Boussau B, Landis MJ, Ronquist F, Huelsenbeck JP. Probabilistic graphical model representation in phylogenetics. Syst Biol 2014; 63:753-71. [PMID: 24951559 PMCID: PMC4184382 DOI: 10.1093/sysbio/syu039] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Recent years have seen a rapid expansion of the model space explored in statistical phylogenetics, emphasizing the need for new approaches to statistical model representation and software development. Clear communication and representation of the chosen model is crucial for: (i) reproducibility of an analysis, (ii) model development, and (iii) software design. Moreover, a unified, clear and understandable framework for model representation lowers the barrier for beginners and nonspecialists to grasp complex phylogenetic models, including their assumptions and parameter/variable dependencies. Graphical modeling is a unifying framework that has gained in popularity in the statistical literature in recent years. The core idea is to break complex models into conditionally independent distributions. The strength lies in the comprehensibility, flexibility, and adaptability of this formalism, and the large body of computational work based on it. Graphical models are well-suited to teach statistical models, to facilitate communication among phylogeneticists and in the development of generic software for simulation and statistical inference. Here, we provide an introduction to graphical models for phylogeneticists and extend the standard graphical model representation to the realm of phylogenetics. We introduce a new graphical model component, tree plates, to capture the changing structure of the subgraph corresponding to a phylogenetic tree. We describe a range of phylogenetic models using the graphical model framework and introduce modules to simplify the representation of standard components in large and complex models. Phylogenetic model graphs can be readily used in simulation, maximum likelihood inference, and Bayesian inference using, for example, Metropolis–Hastings or Gibbs sampling of the posterior distribution. [Computation; graphical models; inference; modularization; statistical phylogenetics; tree plate.]
Collapse
Affiliation(s)
- Sebastian Höhna
- Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia;Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia;
| | - Tracy A Heath
- Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia;Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Bastien Boussau
- Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia;Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Michael J Landis
- Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Fredrik Ronquist
- Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - John P Huelsenbeck
- Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia;Department of Mathematics, Stockholm University, Stockholm, SE-106 91 Stockholm, Sweden; Department of Evolution and Ecology, University of California, Davis, Storer Hall, One Shields Avenue, Davis, CA 95616, USA; Department of Integrative Biology, University of California, Berkeley, CA 94720, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA; Bioinformatics and Evolutionary Genomics, Université de Lyon, Villeurbanne, France; Department of Bioinformatics and Genetics, Swedish Museum of Natural History, SE-10405 Stockholm, Sweden; and Department of Biological Science, King Abdulaziz University, Jeddah, Saudi Arabia
| |
Collapse
|