51
|
Huang J, Bennett J, Flouri T, Leaché AD, Yang Z. Phase Resolution of Heterozygous Sites in Diploid Genomes is Important to Phylogenomic Analysis under the Multispecies Coalescent Model. Syst Biol 2021; 71:334-352. [PMID: 34143216 PMCID: PMC8977997 DOI: 10.1093/sysbio/syab047] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Revised: 06/03/2021] [Accepted: 06/21/2021] [Indexed: 01/01/2023] Open
Abstract
Genome sequencing projects routinely generate haploid consensus sequences from diploid
genomes, which are effectively chimeric sequences with the phase at heterozygous sites
resolved at random. The impact of phasing errors on phylogenomic analyses under the
multispecies coalescent (MSC) model is largely unknown. Here, we conduct a computer
simulation to evaluate the performance of four phase-resolution strategies (the true phase
resolution, the diploid analytical integration algorithm which averages over all phase
resolutions, computational phase resolution using the program PHASE, and random
resolution) on estimation of the species tree and evolutionary parameters in analysis of
multilocus genomic data under the MSC model. We found that species tree estimation is
robust to phasing errors when species divergences were much older than average coalescent
times but may be affected by phasing errors when the species tree is shallow. Estimation
of parameters under the MSC model with and without introgression is affected by phasing
errors. In particular, random phase resolution causes serious overestimation of population
sizes for modern species and biased estimation of cross-species introgression probability.
In general, the impact of phasing errors is greater when the mutation rate is higher, the
data include more samples per species, and the species tree is shallower with recent
divergences. Use of phased sequences inferred by the PHASE program produced small biases
in parameter estimates. We analyze two real data sets, one of East Asian brown frogs and
another of Rocky Mountains chipmunks, to demonstrate that heterozygote phase-resolution
strategies have similar impacts on practical data analyses. We suggest that genome
sequencing projects should produce unphased diploid genotype sequences if fully phased
data are too challenging to generate, and avoid haploid consensus sequences, which have
heterozygous sites phased at random. In case the analytical integration algorithm is
computationally unfeasible, computational phasing prior to population genomic analyses is
an acceptable alternative. [BPP; introgression; multispecies coalescent; phase; species
tree.]
Collapse
Affiliation(s)
- Jun Huang
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK.,Department of Mathematics, Beijing Jiaotong University, Beijing, 100044, P.R. China
| | - Jeremy Bennett
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK.,Department of Ecology and Evolutionary Biology, University of Connecticut, 75 N. Eagleville Road, Unit 3043, Storrs, CT 06269-3043, USA
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK
| | - Adam D Leaché
- Department of Biology & Burke Museum of Natural History and Culture, University of Washington, Seattle, WA 98195-1800, USA
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
52
|
Cai R, Ané C. Assessing the fit of the multi-species network coalescent to multi-locus data. Bioinformatics 2021; 37:634-641. [PMID: 33027508 DOI: 10.1093/bioinformatics/btaa863] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Revised: 09/14/2020] [Accepted: 09/22/2020] [Indexed: 01/25/2023] Open
Abstract
MOTIVATION With growing genome-wide molecular datasets from next-generation sequencing, phylogenetic networks can be estimated using a variety of approaches. These phylogenetic networks include events like hybridization, gene flow or horizontal gene transfer explicitly. However, the most accurate network inference methods are computationally heavy. Methods that scale to larger datasets do not calculate a full likelihood, such that traditional likelihood-based tools for model selection are not applicable to decide how many past hybridization events best fit the data. We propose here a goodness-of-fit test to quantify the fit between data observed from genome-wide multi-locus data, and patterns expected under the multi-species coalescent model on a candidate phylogenetic network. RESULTS We identified weaknesses in the previously proposed TICR test, and proposed corrections. The performance of our new test was validated by simulations on real-world phylogenetic networks. Our test provides one of the first rigorous tools for model selection, to select the adequate network complexity for the data at hand. The test can also work for identifying poorly inferred areas on a network. AVAILABILITY AND IMPLEMENTATION Software for the goodness-of-fit test is available as a Julia package at https://github.com/cecileane/QuartetNetworkGoodnessFit.jl. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Ruoyi Cai
- Department of Statistics, University of Wisconsin - Madison, Madison, WI 53706, USA
| | - Cécile Ané
- Department of Statistics, University of Wisconsin - Madison, Madison, WI 53706, USA.,Department of Botany, University of Wisconsin - Madison, Madison, WI 53706, USA
| |
Collapse
|
53
|
Introgression is widespread in the radiation of carnivorous Nepenthes pitcher plants. Mol Phylogenet Evol 2021; 163:107214. [PMID: 34052438 DOI: 10.1016/j.ympev.2021.107214] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2021] [Revised: 05/14/2021] [Accepted: 05/25/2021] [Indexed: 11/23/2022]
Abstract
Introgression and hybridization are important processes in plant evolution, but they are difficult to study from a phylogenetic perspective, because they conflict with the bifurcating evolutionary history typically depicted in phylogenetic models. The role of hybridization in plant evolution is best documented in the form of allo-polyploidizations. In contrast, homoploid hybridization and introgression are less explored, although they may be crucial in adaptive radiations. Here we employ genome-wide data (ddRAD-seq, transcriptomes) to investigate the evolutionary history of Nepenthes, a radiation of c. 160 species of iconic carnivorous plants mainly from tropical Asia. Our data indicates that the main radiation is only c. 5 million years old, and confirms previous bifurcating phylogenies. However, due to a greatly expanded number of loci, we were able test for the first time the long-standing hypotheses of introgression and historical hybridization. The genus presents one very clear case of organellar capture between two distantly related but sympatric groups. Furthermore, all Nepenthes species show introgression signals in their nuclear genomes, as uncovered by a general survey of ABBA-BABA-like statistics. The ancestor of the rapid main radiation shows ancestry from two deeply diverged lineages, as indicated by phylogenetic network analyses. All major clades of the main radiation show further introgression both within and between each other, as suggested by admixture graphs. Our study supports the hypothesis that rapid adaptive radiations are hotspots of introgression in the tree of life, and highlights the need to consider non-treelike processes in evolutionary studies of Nepenthes in particular.
Collapse
|
54
|
Wang Y, Ogilvie HA, Nakhleh L. Practical Speedup of Bayesian Inference of Species Phylogenies by Restricting the Space of Gene Trees. Mol Biol Evol 2021; 37:1809-1818. [PMID: 32077947 PMCID: PMC7253205 DOI: 10.1093/molbev/msaa045] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Species tree inference from multilocus data has emerged as a powerful paradigm in the postgenomic era, both in terms of the accuracy of the species tree it produces as well as in terms of elucidating the processes that shaped the evolutionary history. Bayesian methods for species tree inference are desirable in this area as they have been shown not only to yield accurate estimates, but also to naturally provide measures of confidence in those estimates. However, the heavy computational requirements of Bayesian inference have limited the applicability of such methods to very small data sets. In this article, we show that the computational efficiency of Bayesian inference under the multispecies coalescent can be improved in practice by restricting the space of the gene trees explored during the random walk, without sacrificing accuracy as measured by various metrics. The idea is to first infer constraints on the trees of the individual loci in the form of unresolved gene trees, and then to restrict the sampler to consider only resolutions of the constrained trees. We demonstrate the improvements gained by such an approach on both simulated and biological data.
Collapse
Affiliation(s)
- Yaxuan Wang
- Computer Science Department, Rice University, Houston, TX
| | - Huw A Ogilvie
- Computer Science Department, Rice University, Houston, TX
| | - Luay Nakhleh
- Computer Science Department, Rice University, Houston, TX
| |
Collapse
|
55
|
Rothfels CJ. Polyploid phylogenetics. THE NEW PHYTOLOGIST 2021; 230:66-72. [PMID: 33491778 DOI: 10.1111/nph.17105] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2020] [Accepted: 10/27/2020] [Indexed: 05/23/2023]
Abstract
Polyploidy is a dominant feature of extant plant diversity. However, major research questions, including whether polyploidy is important to long-term evolution or is just 'evolutionary noise', remain unresolved due to difficulties associated with the generation and analysis of data from polyploid lineages. Many of these difficulties have been recently overcome, such that it is now often relatively straightforward to infer the full and often reticulate phylogenetic history of groups with recently formed polyploids. This nascent field of 'polyploid phylogenetics' allows researchers to tackle long-standing questions of polyploid macroevolution, supplies the foundation for mechanistic models of ploidy change, and provides the opportunity to include a more complete and representative sample of plant taxa in our analyses in general.
Collapse
Affiliation(s)
- Carl J Rothfels
- Department of Integrative Biology, University Herbarium, University of California, Berkeley, CA, 94702, USA
| |
Collapse
|
56
|
Flouri T, Jiao X, Rannala B, Yang Z. A Bayesian Implementation of the Multispecies Coalescent Model with Introgression for Phylogenomic Analysis. Mol Biol Evol 2021; 37:1211-1223. [PMID: 31825513 PMCID: PMC7086182 DOI: 10.1093/molbev/msz296] [Citation(s) in RCA: 58] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Recent analyses suggest that cross-species gene flow or introgression is common in nature, especially during species divergences. Genomic sequence data can be used to infer introgression events and to estimate the timing and intensity of introgression, providing an important means to advance our understanding of the role of gene flow in speciation. Here, we implement the multispecies-coalescent-with-introgression model, an extension of the multispecies-coalescent model to incorporate introgression, in our Bayesian Markov chain Monte Carlo program Bpp. The multispecies-coalescent-with-introgression model accommodates deep coalescence (or incomplete lineage sorting) and introgression and provides a natural framework for inference using genomic sequence data. Computer simulation confirms the good statistical properties of the method, although hundreds or thousands of loci are typically needed to estimate introgression probabilities reliably. Reanalysis of data sets from the purple cone spruce confirms the hypothesis of homoploid hybrid speciation. We estimated the introgression probability using the genomic sequence data from six mosquito species in the Anopheles gambiae species complex, which varies considerably across the genome, likely driven by differential selection against introgressed alleles.
Collapse
Affiliation(s)
- Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| | - Xiyun Jiao
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| | - Bruce Rannala
- Department of Evolution and Ecology, University of California, Davis, Davis, CA
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| |
Collapse
|
57
|
A molecularphylogeny offorktail damselflies(genus Ischnura)revealsa dynamic macroevolutionary history of female colour polymorphisms. Mol Phylogenet Evol 2021; 160:107134. [PMID: 33677008 DOI: 10.1016/j.ympev.2021.107134] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Revised: 02/11/2021] [Accepted: 02/26/2021] [Indexed: 12/18/2022]
Abstract
Colour polymorphisms are popular study systems among biologists interested in evolutionary dynamics, genomics, sexual selection and sexual conflict. In many damselfly groups, such as in the globally distributed genus Ischnura (forktails), sex-limited female colour polymorphisms occur in multiple species. Female-polymorphic species contain two or three female morphs, one of which phenotypically matches the male (androchrome or male mimic) and the other(s) which are phenotypically distinct from the male (heterochrome). These female colour polymorphisms are thought to be maintained by frequency-dependent sexual conflict, but their macroevolutionary histories are unknown, due to the lack of a robust molecular phylogeny. Here, we present the first time-calibrated phylogeny of Ischnura, using a multispecies coalescent approach (StarBEAST2) and incorporating both molecular and fossil data for 41 extant species (55% of the genus). We estimate the age of Ischnura to be between 13.8 and 23.4 millions of years, i.e. Miocene. We infer the ancestral state of this genus as female monomorphism with heterochrome females, with multiple gains and losses of female polymorphisms, evidence of trans-species female polymorphisms and a significant positive relationship between female polymorphism incidence and current geographic range size. Our study provides a robust phylogenetic framework for future research on the dynamic macroevolutionary history of this clade with its extraordinary diversity of sex-limited female polymorphisms.
Collapse
|
58
|
Westbury MV, Le Duc D, Duchêne DA, Krishnan A, Prost S, Rutschmann S, Grau JH, Dalen L, Weyrich A, Norén K, Werdelin L, Dalerum F, Schöneberg T, Hofreiter M. Ecological Specialisation and Evolutionary Reticulation in Extant Hyaenidae. Mol Biol Evol 2021; 38:3884-3897. [PMID: 34426844 PMCID: PMC8382907 DOI: 10.1093/molbev/msab055] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
During the Miocene, Hyaenidae was a highly diverse family of Carnivora that has since been severely reduced to four species: the bone-cracking spotted, striped, and brown hyenas, and the specialized insectivorous aardwolf. Previous studies investigated the evolutionary histories of the spotted and brown hyenas, but little is known about the remaining two species. Moreover, the genomic underpinnings of scavenging and insectivory, defining traits of the extant species, remain elusive. Here, we generated an aardwolf genome and analyzed it together with the remaining three species to reveal their evolutionary relationships, genomic underpinnings of their scavenging and insectivorous lifestyles, and their respective genetic diversities and demographic histories. High levels of phylogenetic discordance suggest gene flow between the aardwolf lineage and the ancestral brown/striped hyena lineage. Genes related to immunity and digestion in the bone-cracking hyenas and craniofacial development in the aardwolf showed the strongest signals of selection, suggesting putative key adaptations to carrion and termite feeding, respectively. A family-wide expansion in olfactory receptor genes suggests that an acute sense of smell was a key early adaptation. Finally, we report very low levels of genetic diversity within the brown and striped hyenas despite no signs of inbreeding, putatively linked to their similarly slow decline in effective population size over the last ∼2 million years. High levels of genetic diversity and more stable population sizes through time are seen in the spotted hyena and aardwolf. Taken together, our findings highlight how ecological specialization can impact the evolutionary history, demographics, and adaptive genetic changes of an evolutionary lineage.
Collapse
Affiliation(s)
- M V Westbury
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, Potsdam, 14476, Germany.,Section for Evolutionary Genomics, The GLOBE Institute, University of Copenhagen, Øster Voldgade 5-7, Copenhagen, Denmark
| | - Diana Le Duc
- Institute of Human Genetics, University Medical Center Leipzig, Leipzig, 04103, Germany.,Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, 04103, Germany
| | - David A Duchêne
- Section for Evolutionary Genomics, The GLOBE Institute, University of Copenhagen, Øster Voldgade 5-7, Copenhagen, Denmark.,Research School of Biology, Australian National University, Canberra, ACT 2601, Australia
| | - Arunkumar Krishnan
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA.,Department of Biological Sciences, Indian Institute of Science Education and Research (IISER) Berhampur, Odisha, 760010, India
| | - Stefan Prost
- LOEWE-Center for Translational Biodiversity Genomics, Senckenberg, 60325, Germany. Frankfurt.,South African National Biodiversity Institute, National Zoological Garden, Pretoria, 0184, South Africa
| | - Sereina Rutschmann
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, Potsdam, 14476, Germany
| | - Jose H Grau
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, Potsdam, 14476, Germany.,amedes Genetics, amedes Medizinische Dienstleistungen, Berlin, Germany
| | - Love Dalen
- Centre for Palaeogenetics, Svante Arrhenius väg 20C, Stockholm, 10691, Sweden.,Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Box 50007, Stockholm, 10405, Sweden
| | - Alexandra Weyrich
- Department of Evolutionary Genetics, Leibniz Institute for Zoo and Wildlife Research (IZW), Berlin, 10315, Germany
| | - Karin Norén
- Department of Zoology, Stockholm University, Stockholm, 106 91, Sweden
| | - Lars Werdelin
- Department of Palaeobiology, Swedish Museum of Natural History, Box 50007, Stockholm, SE-10405, Sweden
| | - Fredrik Dalerum
- Department of Zoology, Stockholm University, Stockholm, 106 91, Sweden.,Research Unit of Biodiversity (UO-CSIC-PA), Mieres Campus, University of Oviedo, Mieres, Asturias, 33600, Spain.,Mammal Research Institute, Department of Zoology and Entomology, University of Pretoria, South Africa
| | - Torsten Schöneberg
- Rudolf Schönheimer Institute of Biochemistry, Molecular Biochemistry, Medical Faculty, Johannisallee 30, Leipzig, 04103, Germany
| | - Michael Hofreiter
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, Potsdam, 14476, Germany
| |
Collapse
|
59
|
Allman ES, Mitchell JD, Rhodes JA. Gene tree discord, simplex plots, and statistical tests under the coalescent. Syst Biol 2021; 71:929-942. [PMID: 33560348 DOI: 10.1093/sysbio/syab008] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 01/31/2021] [Accepted: 02/03/2021] [Indexed: 02/06/2023] Open
Abstract
A simple graphical device, the simplex plot of quartet concordance factors, is introduced to aid in the exploration of a collection of gene trees on a common set of taxa. A single plot summarizes all gene tree discord, and allows for visual comparison to the expected discord from the multispecies coalescent model (MSC) of incomplete lineage sorting on a species tree. A formal statistical procedure is described that can quantify the deviation from expectation for each subset of four taxa, suggesting when the data is not in accord with the MSC, and thus that either gene tree inference error is substantial or a more complex model such as that on a network may be required. If the collection of gene trees is in accord with the MSC, the plots reveal when substantial incomplete lineage sorting is present. Applications to both simulated and empirical multilocus data sets illustrate the insights provided.
Collapse
Affiliation(s)
- Elizabeth S Allman
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK 99709, USA
| | - Jonathan D Mitchell
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK 99709, USA.,Unité Bioinformatique Evolutive, C3BI USR 3756, Institut Pasteur & CNRS, Paris, France
| | - John A Rhodes
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK 99709, USA
| |
Collapse
|
60
|
Lopes F, Oliveira LR, Kessler A, Beux Y, Crespo E, Cárdenas-Alayza S, Majluf P, Sepúlveda M, Brownell RL, Franco-Trecu V, Páez-Rosas D, Chaves J, Loch C, Robertson BC, Acevedo-Whitehouse K, Elorriaga-Verplancken FR, Kirkman SP, Peart CR, Wolf JBW, Bonatto SL. Phylogenomic Discordance in the Eared Seals is best explained by Incomplete Lineage Sorting following Explosive Radiation in the Southern Hemisphere. Syst Biol 2020; 70:786-802. [PMID: 33367817 DOI: 10.1093/sysbio/syaa099] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 12/02/2020] [Accepted: 12/08/2020] [Indexed: 12/25/2022] Open
Abstract
The phylogeny and systematics of fur seals and sea lions (Otariidae) have long been studied with diverse data types, including an increasing amount of molecular data. However, only a few phylogenetic relationships have reached acceptance because of strong gene-tree species tree discordance. Divergence times estimates in the group also vary largely between studies. These uncertainties impeded the understanding of the biogeographical history of the group, such as when and how trans-equatorial dispersal and subsequent speciation events occurred. Here, we used high-coverage genome-wide sequencing for 14 of the 15 species of Otariidae to elucidate the phylogeny of the family and its bearing on the taxonomy and biogeographical history. Despite extreme topological discordance among gene trees, we found a fully supported species tree that agrees with the few well-accepted relationships and establishes monophyly of the genus Arctocephalus. Our data support a relatively recent trans-hemispheric dispersal at the base of a southern clade, which rapidly diversified into six major lineages between 3 and 2.5 Ma. Otaria diverged first, followed by Phocarctos and then four major lineages within Arctocephalus. However, we found Zalophus to be nonmonophyletic, with California (Zalophus californianus) and Steller sea lions (Eumetopias jubatus) grouping closer than the Galapagos sea lion (Zalophus wollebaeki) with evidence for introgression between the two genera. Overall, the high degree of genealogical discordance was best explained by incomplete lineage sorting resulting from quasi-simultaneous speciation within the southern clade with introgresssion playing a subordinate role in explaining the incongruence among and within prior phylogenetic studies of the family. [Hybridization; ILS; phylogenomics; Pleistocene; Pliocene; monophyly.].
Collapse
Affiliation(s)
- Fernando Lopes
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, 90619-900 Porto Alegre, RS, Brazil.,Laboratório de Ecologia de Mamíferos, Universidade do Vale do Rio dos Sinos, São Leopoldo, RS, Brazil
| | - Larissa R Oliveira
- Laboratório de Ecologia de Mamíferos, Universidade do Vale do Rio dos Sinos, São Leopoldo, RS, Brazil.,GEMARS, Grupo de Estudos de Mamíferos Aquáticos do Rio Grande do Sul, 95560-000 Torres, RS, Brazil
| | - Amanda Kessler
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, 90619-900 Porto Alegre, RS, Brazil
| | - Yago Beux
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, 90619-900 Porto Alegre, RS, Brazil
| | - Enrique Crespo
- Centro Nacional Patagónico - CENPAT, CONICET, Puerto Madryn, Argentina
| | - Susana Cárdenas-Alayza
- Centro para la Sostenibilidad Ambiental, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Patricia Majluf
- Centro para la Sostenibilidad Ambiental, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Maritza Sepúlveda
- Centro de Investigación y Gestión de Recursos Naturales (CIGREN), Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
| | - Robert L Brownell
- Southwest Fisheries Science Center, National Oceanic and Atmospheric Administration, NOAA, La Jolla, USA
| | - Valentina Franco-Trecu
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Diego Páez-Rosas
- Colegio de Ciencias Biológicas y Ambientales, COCIBA, Universidad San Francisco de Quito, Quito, Ecuador
| | - Jaime Chaves
- Colegio de Ciencias Biológicas y Ambientales, COCIBA, Universidad San Francisco de Quito, Quito, Ecuador.,Department of Biology, San Francisco State University, 1800 Holloway Ave, San Francisco, CA, USA
| | - Carolina Loch
- Sir John Walsh Research Institute, Faculty of Dentistry, University of Otago, Dunedin, New Zealand
| | | | - Karina Acevedo-Whitehouse
- Unit for Basic and Applied Microbiology, School of Natural Sciences, Universidad Autónoma de Querétaro, Querétaro, Mexico
| | | | - Stephen P Kirkman
- Department of Environmental Affairs, Oceans and Coasts, Cape Town, South Africa
| | - Claire R Peart
- Department Biologie II, Division of Evolutionary Biology, Ludwig-Maximilians-Universität München, Münich, Germany
| | - Jochen B W Wolf
- Department Biologie II, Division of Evolutionary Biology, Ludwig-Maximilians-Universität München, Münich, Germany
| | - Sandro L Bonatto
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, 90619-900 Porto Alegre, RS, Brazil
| |
Collapse
|
61
|
Hibbins MS, Gibson MJS, Hahn MW. Determining the probability of hemiplasy in the presence of incomplete lineage sorting and introgression. eLife 2020; 9:e63753. [PMID: 33345772 PMCID: PMC7800383 DOI: 10.7554/elife.63753] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 12/18/2020] [Indexed: 12/11/2022] Open
Abstract
The incongruence of character states with phylogenetic relationships is often interpreted as evidence of convergent evolution. However, trait evolution along discordant gene trees can also generate these incongruences - a phenomenon known as hemiplasy. Classic comparative methods do not account for discordance, resulting in incorrect inferences about the number, timing, and direction of trait transitions. Biological sources of discordance include incomplete lineage sorting (ILS) and introgression, but only ILS has received theoretical consideration in the context of hemiplasy. Here, we present a model that shows introgression makes hemiplasy more likely, such that methods that account for ILS alone will be conservative. We also present a method and software (HeIST) for making statistical inferences about the probability of hemiplasy and homoplasy in large datasets that contain both ILS and introgression. We apply our methods to two empirical datasets, finding that hemiplasy is likely to contribute to the observed trait incongruences in both.
Collapse
Affiliation(s)
- Mark S Hibbins
- Department of Biology, Indiana UniversityBloomingtonUnited States
| | | | - Matthew W Hahn
- Department of Biology, Indiana UniversityBloomingtonUnited States
- Department of Computer Science, Indiana UniversityBloomingtonUnited States
| |
Collapse
|
62
|
Koch H, DeGiorgio M. Maximum Likelihood Estimation of Species Trees from Gene Trees in the Presence of Ancestral Population Structure. Genome Biol Evol 2020; 12:3977-3995. [PMID: 32022857 PMCID: PMC7061232 DOI: 10.1093/gbe/evaa022] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/23/2020] [Indexed: 11/12/2022] Open
Abstract
Though large multilocus genomic data sets have led to overall improvements in phylogenetic inference, they have posed the new challenge of addressing conflicting signals across the genome. In particular, ancestral population structure, which has been uncovered in a number of diverse species, can skew gene tree frequencies, thereby hindering the performance of species tree estimators. Here we develop a novel maximum likelihood method, termed TASTI (Taxa with Ancestral structure Species Tree Inference), that can infer phylogenies under such scenarios, and find that it has increasing accuracy with increasing numbers of input gene trees, contrasting with the relatively poor performances of methods not tailored for ancestral structure. Moreover, we propose a supertree approach that allows TASTI to scale computationally with increasing numbers of input taxa. We use genetic simulations to assess TASTI's performance in the three- and four-taxon settings and demonstrate the application of TASTI on a six-species Afrotropical mosquito data set. Finally, we have implemented TASTI in an open-source software package for ease of use by the scientific community.
Collapse
Affiliation(s)
- Hillary Koch
- Department of Statistics, Pennsylvania State University
| | - Michael DeGiorgio
- Department of Computer and Electrical Engineering and Computer Science, Florida Atlantic University
| |
Collapse
|
63
|
Blair C, Ané C. Phylogenetic Trees and Networks Can Serve as Powerful and Complementary Approaches for Analysis of Genomic Data. Syst Biol 2020; 69:593-601. [PMID: 31432090 DOI: 10.1093/sysbio/syz056] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2019] [Accepted: 08/15/2019] [Indexed: 11/14/2022] Open
Abstract
Genomic data have had a profound impact on nearly every biological discipline. In systematics and phylogenetics, the thousands of loci that are now being sequenced can be analyzed under the multispecies coalescent model (MSC) to explicitly account for gene tree discordance due to incomplete lineage sorting (ILS). However, the MSC assumes no gene flow post divergence, calling for additional methods that can accommodate this limitation. Explicit phylogenetic network methods have emerged, which can simultaneously account for ILS and gene flow by representing evolutionary history as a directed acyclic graph. In this point of view, we highlight some of the strengths and limitations of phylogenetic networks and argue that tree-based inference should not be blindly abandoned in favor of networks simply because they represent more parameter rich models. Attention should be given to model selection of reticulation complexity, and the most robust conclusions regarding evolutionary history are likely obtained when combining tree- and network-based inference.
Collapse
Affiliation(s)
- Christopher Blair
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY 11201, USA
- Biology PhD Program, CUNY Graduate Center, 365 5th Ave., New York, NY 10016, USA
| | - Cécile Ané
- Department of Botany, University of Wisconsin - Madison, 1300 University Ave, Madison, WI 53706, USA
- Department of Statistics, University of Wisconsin - Madison, 1300 University Ave, Madison, WI 53706, USA
| |
Collapse
|
64
|
Swimming through the sands of the Sahara and Arabian deserts: Phylogeny of sandfish skinks (Scincidae, Scincus) reveals a recent and rapid diversification. Mol Phylogenet Evol 2020; 155:107012. [PMID: 33217580 DOI: 10.1016/j.ympev.2020.107012] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Revised: 11/05/2020] [Accepted: 11/11/2020] [Indexed: 11/18/2022]
Abstract
Large parts of the Sahara Desert and Arabia are covered by sand seas and sand dunes, which are inhabited by specialized animal communities. For example, many lizards have developed adaptations to life in loose sand, including sand-swimming behavior. The best-known sand swimmers of the Saharo-Arabia are the sandfish skinks (genus Scincus). Although there are currently only four Scincus species recognized, their phylogenetic relationships have not yet been addressed in detail. We use eight genetic markers (three mitochondrial, five nuclear) and a complete sampling of species to infer the relationships within the genus. We employ multiple phylogenetic approaches to reconstruct the evolutionary history of these skinks and to assess the level of reticulation at the onset of their radiation. Our results indicate the presence of five strongly supported species-level lineages, four represented by the currently recognized species and the fifth by S. scincus conirostris, which does not form a clade with S. scincus. Based on these results we elevate the Iranian and northern Arabian S. conirostris to the species level. The two Saharan species, S. albifasciatus and S. scincus, are sister in all analyses. Deeper relationships within the genus, however, remained largely unresolved despite the extensive genetic data set. This basal polytomy, together with the fact that we detected no sign of hybridization in the history of the genus, indicates that the diversification of the five Scincus species was rapid, burst-like, and not followed by secondary hybridization events. Divergence time estimations show a Middle Pliocene crown radiation of the genus (3.3 Mya). We hypothesize that the aridification of the Saharo-Arabia that began in the Late Miocene triggered the initial diversification of Scincus, and that the subsequent expansion of sand deserts enabled their dispersal over the large Saharan and Arabian range. We discuss the evolution of body form in sand swimming lizards and ponder how Scincus retained their fully limbed morphology despite being sand swimmers that are typically limbless.
Collapse
|
65
|
Magalhães RF, Lemes P, Santos MTT, Mol RM, Ramos EKS, Oswald CB, Pezzuti TL, Santos FR, Brandão RA, Garcia PCA. Evidence of introgression in endemic frogs from the campo rupestre contradicts the reduced hybridization hypothesis. Biol J Linn Soc Lond 2020. [DOI: 10.1093/biolinnean/blaa142] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Abstract
The campo rupestre ecosystem is considered an old, climatically buffered, infertile landscape. As a consequence, long-term isolation is thought to have played an important role in the diversification of its biota. Here, we tested for hybridization between two endemic leaf frogs from the campo rupestre. We used sequence markers and coalescent models to verify haplotype sharing between the species, to test the existence and direction of gene flow, and to reconstruct the spatiotemporal dynamics of gene flow. Additionally, ecological niche modelling (ENM) was used to assess for potential co-occurrence by overlapping the climatic niche of these species since the middle Pleistocene. We found haplotype sharing and/or lack of differentiation in four nuclear fragments, one of them associated with introgression. The coalescent models support introgressive hybridization unidirectionally from Pithecopus megacephalus to P. ayeaye, occurring ~300 kya. ENM corroborates this scenario, revealing areas of potential environmental niche overlap for the species at about 787 kya. These results contradict the expectation of reduced hybridization, while ENM suggests climatic fluctuation rather than stability for the two species. The reduced hybridization hypothesis needs to be further investigated because our results suggest that it may have unrealistic premises at least for animals.
Collapse
Affiliation(s)
- Rafael F Magalhães
- Departamento de Ciências Naturais, Universidade Federal de São João del-Rei, São João del-Rei, Minas Gerais, Brazil
- Programa de Pós-Graduação em Zoologia, Departamento de Zoologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
- Laboratório de Biodiversidade e Evolução Molecular, Departamento de Genética, Ecologia e Evolução, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
| | - Priscila Lemes
- Laboratorio de Ecologia e Conservação, Departamento de Botânica e Ecologia, Instituto de Biociências, Universidade Federal do Mato Grosso, Cuiabá, Mato Grosso, Brazil
| | - Marcus Thadeu T Santos
- Programa de Pós-Graduação em Zoologia, Instituto de Biociências, Universidade Estadual Paulista (UNESP), Rio Claro, São Paulo, Brazil
| | - Rafael M Mol
- Programa de Pós-Graduação em Zoologia, Departamento de Zoologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
| | - Elisa K S Ramos
- Programa de Pós-Graduação em Genética e Biologia Molecular, Departamento de Genética, Evolução, Microbiologia e Imunologia, Universidade Estadual de Campinas, Campinas, São Paulo, Brazil
| | - Caroline B Oswald
- Programa de Pós-Graduação em Zoologia, Departamento de Zoologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
| | - Tiago L Pezzuti
- Programa de Pós-Graduação em Zoologia, Departamento de Zoologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
| | - Fabrício R Santos
- Programa de Pós-Graduação em Zoologia, Departamento de Zoologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
- Laboratório de Biodiversidade e Evolução Molecular, Departamento de Genética, Ecologia e Evolução, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
| | - Reuber A Brandão
- Laboratório de Fauna e Unidades de Conservação, Departamento de Engenharia Florestal, Universidade de Brasília, Brasília, Distrito Federal, Brazil
| | - Paulo C A Garcia
- Programa de Pós-Graduação em Zoologia, Departamento de Zoologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brazil
| |
Collapse
|
66
|
Rengifo‐Correa L, Abad‐Franch F, Martínez‐Hernández F, Salazar‐Schettino PM, Téllez‐Rendón JL, Villalobos G, Morrone JJ. A biogeographic–ecological approach to disentangle reticulate evolution in the
Triatoma phyllosoma
species group (Heteroptera: Triatominae), vectors of Chagas disease. J ZOOL SYST EVOL RES 2020. [DOI: 10.1111/jzs.12409] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Laura Rengifo‐Correa
- Departamento de Biología Evolutiva Facultad de Ciencias Museo de Zoología ‘Alfonso L. Herrera’Universidad Nacional Autónoma de México Mexico City Mexico
| | - Fernando Abad‐Franch
- Programa de Pós‐graduação em Medicina Tropical Núcleo de Medicina Tropical Facultade Medicina Universidade de Brasília Brasília Brazil
| | | | - Paz M. Salazar‐Schettino
- Laboratorio de Biología de Parásitos Departamento de Microbiología y Parasitología Facultad de Medicina Universidad Nacional Autónoma de México Mexico City Mexico
| | | | - Guiehdani Villalobos
- Departamento de Ecología de Agentes Patógenos Hospital General Dr. Manuel Gea González Mexico City Mexico
| | - Juan J. Morrone
- Departamento de Biología Evolutiva Facultad de Ciencias Museo de Zoología ‘Alfonso L. Herrera’Universidad Nacional Autónoma de México Mexico City Mexico
| |
Collapse
|
67
|
Huang J, Flouri T, Yang Z. A Simulation Study to Examine the Information Content in Phylogenomic Data Sets under the Multispecies Coalescent Model. Mol Biol Evol 2020; 37:3211-3224. [DOI: 10.1093/molbev/msaa166] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
AbstractWe use computer simulation to examine the information content in multilocus data sets for inference under the multispecies coalescent model. Inference problems considered include estimation of evolutionary parameters (such as species divergence times, population sizes, and cross-species introgression probabilities), species tree estimation, and species delimitation based on Bayesian comparison of delimitation models. We found that the number of loci is the most influential factor for almost all inference problems examined. Although the number of sequences per species does not appear to be important to species tree estimation, it is very influential to species delimitation. Increasing the number of sites and the per-site mutation rate both increase the mutation rate for the whole locus and these have the same effect on estimation of parameters, but the sequence length has a greater effect than the per-site mutation rate for species tree estimation. We discuss the computational costs when the data size increases and provide guidelines concerning the subsampling of genomic data to enable the application of full-likelihood methods of inference.
Collapse
Affiliation(s)
- Jun Huang
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
- Department of Mathematics, Beijing Jiaotong University, Beijing, P.R. China
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| |
Collapse
|
68
|
Jiao X, Yang Z. Defining Species When There is Gene Flow. Syst Biol 2020; 70:108-119. [PMID: 32617579 DOI: 10.1093/sysbio/syaa052] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Revised: 06/23/2020] [Accepted: 06/23/2020] [Indexed: 12/20/2022] Open
Abstract
Whatever one's definition of species, it is generally expected that individuals of the same species should be genetically more similar to each other than they are to individuals of another species. Here, we show that in the presence of cross-species gene flow, this expectation may be incorrect. We use the multispecies coalescent model with continuous-time migration or episodic introgression to study the impact of gene flow on genetic differences within and between species and highlight a surprising but plausible scenario in which different population sizes and asymmetrical migration rates cause a genetic sequence to be on average more closely related to a sequence from another species than to a sequence from the same species. Our results highlight the extraordinary impact that even a small amount of gene flow may have on the genetic history of the species. We suggest that contrasting long-term migration rate and short-term hybridization rate, both of which can be estimated using genetic data, may be a powerful approach to detecting the presence of reproductive barriers and to define species boundaries.[Gene flow; introgression; migration; multispecies coalescent; species concept; species delimitation.].
Collapse
Affiliation(s)
- Xiyun Jiao
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
69
|
Zhu J, Liu X, Ogilvie HA, Nakhleh LK. A divide-and-conquer method for scalable phylogenetic network inference from multilocus data. Bioinformatics 2020; 35:i370-i378. [PMID: 31510688 PMCID: PMC6612858 DOI: 10.1093/bioinformatics/btz359] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Motivation Reticulate evolutionary histories, such as those arising in the presence of hybridization, are best modeled as phylogenetic networks. Recently developed methods allow for statistical inference of phylogenetic networks while also accounting for other processes, such as incomplete lineage sorting. However, these methods can only handle a small number of loci from a handful of genomes. Results In this article, we introduce a novel two-step method for scalable inference of phylogenetic networks from the sequence alignments of multiple, unlinked loci. The method infers networks on subproblems and then merges them into a network on the full set of taxa. To reduce the number of trinets to infer, we formulate a Hitting Set version of the problem of finding a small number of subsets, and implement a simple heuristic to solve it. We studied their performance, in terms of both running time and accuracy, on simulated as well as on biological datasets. The two-step method accurately infers phylogenetic networks at a scale that is infeasible with existing methods. The results are a significant and promising step towards accurate, large-scale phylogenetic network inference. Availability and implementation We implemented the algorithms in the publicly available software package PhyloNet (https://bioinfocs.rice.edu/PhyloNet). Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Jiafan Zhu
- Department of Computer Science, Rice University, Houston, TX, USA
| | - Xinhao Liu
- Department of Computer Science, Rice University, Houston, TX, USA
| | - Huw A Ogilvie
- Department of Computer Science, Rice University, Houston, TX, USA
| | - Luay K Nakhleh
- Department of Computer Science, Rice University, Houston, TX, USA.,Department of BioSciences, Rice University, Houston, TX, USA
| |
Collapse
|
70
|
Abstract
Knowing phylogenetic relationships among species is fundamental for many studies in biology. An accurate phylogenetic tree underpins our understanding of the major transitions in evolution, such as the emergence of new body plans or metabolism, and is key to inferring the origin of new genes, detecting molecular adaptation, understanding morphological character evolution and reconstructing demographic changes in recently diverged species. Although data are ever more plentiful and powerful analysis methods are available, there remain many challenges to reliable tree building. Here, we discuss the major steps of phylogenetic analysis, including identification of orthologous genes or proteins, multiple sequence alignment, and choice of substitution models and inference methodologies. Understanding the different sources of errors and the strategies to mitigate them is essential for assembling an accurate tree of life.
Collapse
|
71
|
Crowl AA, Manos PS, McVay JD, Lemmon AR, Lemmon EM, Hipp AL. Uncovering the genomic signature of ancient introgression between white oak lineages (Quercus). THE NEW PHYTOLOGIST 2020; 226:1158-1170. [PMID: 30963585 DOI: 10.1111/nph.15842] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Accepted: 03/29/2019] [Indexed: 05/10/2023]
Abstract
Botanists have long recognised interspecific gene flow as a common occurrence within white oaks (Quercus section Quercus). Historical allele exchange, however, has not been fully characterised and the complex genomic signals resulting from the combination of vertical and horizontal gene transmission may confound phylogenetic inference and obscure our ability to accurately infer the deep evolutionary history of oaks. Using anchored enrichment, we obtained a phylogenomic dataset consisting of hundreds of single-copy nuclear loci. Concatenation, species-tree and network analyses were carried out in an attempt to uncover the genomic signal of ancient introgression and infer the divergent phylogenetic topology for the white oak clade. Locus and site-level likelihood comparisons were then conducted to further explore the introgressed signal within our dataset. Historical, intersectional gene flow is suggested to have occurred between an ancestor of the Eurasian Roburoid lineage and Quercus pontica and North American Dumosae and Prinoideae lineages. Despite extensive time past, our approach proved successful in detecting the genomic signature of ancient introgression. Our results, however, highlight the importance of sampling and the use of a plurality of analytical tools and methods to sufficiently explore genomic datasets, uncover this signal, and accurately infer evolutionary history.
Collapse
Affiliation(s)
- Andrew A Crowl
- Department of Biology, Duke University, Durham, NC, 27708, USA
| | - Paul S Manos
- Department of Biology, Duke University, Durham, NC, 27708, USA
| | - John D McVay
- Department of Biology, Duke University, Durham, NC, 27708, USA
| | - Alan R Lemmon
- Department of Scientific Computing, Florida State University, Dirac Science Library, Tallahassee, FL, 32317, USA
| | - Emily Moriarty Lemmon
- Department of Biological Science, Florida State University, 89 Chieftan Way, Tallahassee, FL, 32317, USA
| | - Andrew L Hipp
- The Morton Arboretum, 4100 Illinois Route 53, Lisle, IL, 60532, USA
- The Field Museum, 1400 S Lake Shore Drive, Chicago, IL, 60605, USA
| |
Collapse
|
72
|
Hipp AL, Manos PS, Hahn M, Avishai M, Bodénès C, Cavender-Bares J, Crowl AA, Deng M, Denk T, Fitz-Gibbon S, Gailing O, González-Elizondo MS, González-Rodríguez A, Grimm GW, Jiang XL, Kremer A, Lesur I, McVay JD, Plomion C, Rodríguez-Correa H, Schulze ED, Simeone MC, Sork VL, Valencia-Avalos S. Genomic landscape of the global oak phylogeny. THE NEW PHYTOLOGIST 2020; 226:1198-1212. [PMID: 31609470 DOI: 10.1111/nph.16162] [Citation(s) in RCA: 97] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Accepted: 06/05/2019] [Indexed: 05/10/2023]
Abstract
The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny. We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny. Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity. The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.
Collapse
Affiliation(s)
- Andrew L Hipp
- The Morton Arboretum, Lisle, IL, 60532-1293, USA
- The Field Museum, Chicago, IL, 60605, USA
| | | | - Marlene Hahn
- The Morton Arboretum, Lisle, IL, 60532-1293, USA
| | - Michael Avishai
- Previously of, The Hebrew University of Jerusalem, Botanical Garden, Zalman Shne'ur St. 1, Jerusalem, Israel
| | | | | | | | - Min Deng
- Shanghai Chenshan Plant Science Research Center, Chinese Academy of Sciences, Shanghai, 201602, China
| | - Thomas Denk
- Swedish Museum of Natural History, Stockholm, 10405, Sweden
| | | | - Oliver Gailing
- Büsgen-Institute, Georg-August-University Göttingen, Göttingen, D-37077, Germany
| | | | - Antonio González-Rodríguez
- Escuela Nacional de Estudios Superiores Unidad Morelia, Universidad Nacional Autónoma de México, Morelia, 58190, México
| | | | - Xiao-Long Jiang
- Shanghai Chenshan Plant Science Research Center, Chinese Academy of Sciences, Shanghai, 201602, China
| | | | | | | | | | - Hernando Rodríguez-Correa
- Escuela Nacional de Estudios Superiores Unidad Morelia, Universidad Nacional Autónoma de México, Morelia, 58190, México
| | - Ernst-Detlef Schulze
- Max Planck Institute for Biogeochemistry, Hans-Knoell-Str. 10, Jena, 07745, Germany
| | | | - Victoria L Sork
- University of California, Los Angeles, Los Angeles, CA, 90095, USA
| | - Susana Valencia-Avalos
- Herbario de la Facultad de Ciencias, Departamento de Biología Comparada, Universidad Nacional Autónoma de México, Circuito Exterior, s.n., Ciudad Universitaria, Coyoacán, México City, CP 04510, México
| |
Collapse
|
73
|
Bernhardt N, Brassac J, Dong X, Willing EM, Poskar CH, Kilian B, Blattner FR. Genome-wide sequence information reveals recurrent hybridization among diploid wheat wild relatives. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2020; 102:493-506. [PMID: 31821649 DOI: 10.1111/tpj.14641] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Revised: 11/13/2019] [Accepted: 11/28/2019] [Indexed: 05/07/2023]
Abstract
Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome-wide sequence information, independent data were obtained from genotyping-by-sequencing and a target-enrichment experiment that returned 244 low-copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent-based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat-group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat-group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.
Collapse
Affiliation(s)
- Nadine Bernhardt
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), 06466, Gatersleben, Germany
| | - Jonathan Brassac
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), 06466, Gatersleben, Germany
| | - Xue Dong
- Max Planck Institute for Plant Breeding Research, 50829, Cologne, Germany
- Plant Germplasm and Genomics Centre, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, 650201, Kunming, Yunnan, China
| | - Eva-Maria Willing
- Max Planck Institute for Plant Breeding Research, 50829, Cologne, Germany
| | - C Hart Poskar
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), 06466, Gatersleben, Germany
| | - Benjamin Kilian
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), 06466, Gatersleben, Germany
- Global Crop Diversity Trust, 53113, Bonn, Germany
| | - Frank R Blattner
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), 06466, Gatersleben, Germany
- German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, 04103, Leipzig, Germany
| |
Collapse
|
74
|
Olave M, Meyer A. Implementing Large Genomic Single Nucleotide Polymorphism Data Sets in Phylogenetic Network Reconstructions: A Case Study of Particularly Rapid Radiations of Cichlid Fish. Syst Biol 2020; 69:848-862. [DOI: 10.1093/sysbio/syaa005] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Revised: 01/09/2020] [Accepted: 01/23/2020] [Indexed: 12/23/2022] Open
Abstract
AbstractThe Midas cichlids of the Amphilophus citrinellus spp. species complex from Nicaragua (13 species) are an extraordinary example of adaptive and rapid radiation ($<$24,000 years old). These cichlids are a very challenging group to infer its evolutionary history in phylogenetic analyses, due to the apparent prevalence of incomplete lineage sorting (ILS), as well as past and current gene flow. Assuming solely a vertical transfer of genetic material from an ancestral lineage to new lineages is not appropriate in many cases of genes transferred horizontally in nature. Recently developed methods to infer phylogenetic networks under such circumstances might be able to circumvent these problems. These models accommodate not just ILS, but also gene flow, under the multispecies network coalescent (MSNC) model, processes that are at work in young, hybridizing, and/or rapidly diversifying lineages. There are currently only a few programs available that implement MSNC for estimating phylogenetic networks. Here, we present a novel way to incorporate single nucleotide polymorphism (SNP) data into the currently available PhyloNetworks program. Based on simulations, we demonstrate that SNPs can provide enough power to recover the true phylogenetic network. We also show that it can accurately infer the true network more often than other similar SNP-based programs (PhyloNet and HyDe). Moreover, our approach results in a faster algorithm compared to the original pipeline in PhyloNetworks, without losing power. We also applied our new approach to infer the phylogenetic network of Midas cichlid radiation. We implemented the most comprehensive genomic data set to date (RADseq data set of 679 individuals and $>$37K SNPs from 19 ingroup lineages) and present estimated phylogenetic networks for this extremely young and fast-evolving radiation of cichlid fish. We demonstrate that the MSNC is more appropriate than the multispecies coalescent alone for the analysis of this rapid radiation. [Genomics; multispecies network coalescent; phylogenetic networks; phylogenomics; RADseq; SNPs.]
Collapse
Affiliation(s)
- Melisa Olave
- Department of Biology, University of Konstanz, 78457 Konstanz, Germany
| | - Axel Meyer
- Department of Biology, University of Konstanz, 78457 Konstanz, Germany
| |
Collapse
|
75
|
Jiao X, Flouri T, Rannala B, Yang Z. The Impact of Cross-Species Gene Flow on Species Tree Estimation. Syst Biol 2020; 69:830-847. [DOI: 10.1093/sysbio/syaa001] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2018] [Revised: 11/12/2019] [Accepted: 01/15/2020] [Indexed: 12/26/2022] Open
Abstract
Abstract
Recent analyses of genomic sequence data suggest cross-species gene flow is common in both plants and animals, posing challenges to species tree estimation. We examine the levels of gene flow needed to mislead species tree estimation with three species and either episodic introgressive hybridization or continuous migration between an outgroup and one ingroup species. Several species tree estimation methods are examined, including the majority-vote method based on the most common gene tree topology (with either the true or reconstructed gene trees used), the UPGMA method based on the average sequence distances (or average coalescent times) between species, and the full-likelihood method based on multilocus sequence data. Our results suggest that the majority-vote method based on gene tree topologies is more robust to gene flow than the UPGMA method based on coalescent times and both are more robust than likelihood assuming a multispecies coalescent (MSC) model with no cross-species gene flow. Comparison of the continuous migration model with the episodic introgression model suggests that a small amount of gene flow per generation can cause drastic changes to the genetic history of the species and mislead species tree methods, especially if the species diverged through radiative speciation events. Estimates of parameters under the MSC with gene flow suggest that African mosquito species in the Anopheles gambiae species complex constitute such an example of extreme impact of gene flow on species phylogeny. [IM; introgression; migration; MSci; multispecies coalescent; species tree.]
Collapse
Affiliation(s)
- Xiyun Jiao
- Department of Genetics, University College London, Gower Street, London WC1E 6BT, UK
| | - Tomáš Flouri
- Department of Genetics, University College London, Gower Street, London WC1E 6BT, UK
| | - Bruce Rannala
- Department of Evolution and Ecology, University of California, Davis, CA 95616, USA
| | - Ziheng Yang
- Department of Genetics, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
76
|
Duckett DJ, Pelletier TA, Carstens BC. Identifying model violations under the multispecies coalescent model using P2C2M.SNAPP. PeerJ 2020; 8:e8271. [PMID: 31949994 PMCID: PMC6956792 DOI: 10.7717/peerj.8271] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2019] [Accepted: 11/22/2019] [Indexed: 11/20/2022] Open
Abstract
Phylogenetic estimation under the multispecies coalescent model (MSCM) assumes all incongruence among loci is caused by incomplete lineage sorting. Therefore, applying the MSCM to datasets that contain incongruence that is caused by other processes, such as gene flow, can lead to biased phylogeny estimates. To identify possible bias when using the MSCM, we present P2C2M.SNAPP. P2C2M.SNAPP is an R package that identifies model violations using posterior predictive simulation. P2C2M.SNAPP uses the posterior distribution of species trees output by the software package SNAPP to simulate posterior predictive datasets under the MSCM, and then uses summary statistics to compare either the empirical data or the posterior distribution to the posterior predictive distribution to identify model violations. In simulation testing, P2C2M.SNAPP correctly classified up to 83% of datasets (depending on the summary statistic used) as to whether or not they violated the MSCM model. P2C2M.SNAPP represents a user-friendly way for researchers to perform posterior predictive model checks when using the popular SNAPP phylogenetic estimation program. It is freely available as an R package, along with additional program details and tutorials.
Collapse
Affiliation(s)
- Drew J Duckett
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, USA
| | | | - Bryan C Carstens
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
77
|
Granados-Aguilar X, Granados Mendoza C, Cervantes CR, Montes JR, Arias S. Unraveling Reticulate Evolution in Opuntia (Cactaceae) From Southern Mexico. FRONTIERS IN PLANT SCIENCE 2020; 11:606809. [PMID: 33519858 PMCID: PMC7838128 DOI: 10.3389/fpls.2020.606809] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 12/10/2020] [Indexed: 05/20/2023]
Abstract
The process of hybridization occurs in approximately 40% of vascular plants, and this exchange of genetic material between non-conspecific individuals occurs unequally among plant lineages, being more frequent in certain groups such as Opuntia (Cactaceae). This genus is known for multiple taxonomic controversies due to widespread polyploidy and probable hybrid origin of several of its species. Southern Mexico species of this genus have been poorly studied despite their great diversity in regions such as the Tehuacán-Cuicatlán Valley which contains around 12% of recognized Mexico's native Opuntia species. In this work, we focus on testing the hybrid status of two putative hybrids from this region, Opuntia tehuacana and Opuntia pilifera, and estimate if hybridization occurs among sampled southern opuntias using two newly identified nuclear intron markers to construct phylogenetic networks with HyDe and Dsuite and perform invariant analysis under the coalescent model with HyDe and Dsuite. For the test of hybrid origin in O. tehuacana, our results could not recover hybridization as proposed in the literature, but we found introgression into O. tehuacana individuals involving O. decumbens and O. huajuapensis. Regarding O. pilifera, we identified O. decumbens as probable parental species, supported by our analysis, which sustains the previous hybridization hypothesis between Nopalea and Basilares clades. Finally, we suggest new hybridization and introgression cases among southern Mexican species involving O. tehuantepecana and O. depressa as parental species of O. velutina and O. decumbens.
Collapse
Affiliation(s)
- Xochitl Granados-Aguilar
- Posgrado en Ciencias Biológicas, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Jardín Botánico, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
- *Correspondence: Xochitl Granados-Aguilar,
| | - Carolina Granados Mendoza
- Departamento de Botánica, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Cristian Rafael Cervantes
- Posgrado en Ciencias Biológicas, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Jardín Botánico, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - José Rubén Montes
- Posgrado en Ciencias Biológicas, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Departamento de Botánica, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Salvador Arias
- Jardín Botánico, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Salvador Arias,
| |
Collapse
|
78
|
Allman ES, Baños H, Rhodes JA. NANUQ: a method for inferring species networks from gene trees under the coalescent model. Algorithms Mol Biol 2019; 14:24. [PMID: 31827592 PMCID: PMC6896299 DOI: 10.1186/s13015-019-0159-2] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Accepted: 11/07/2019] [Indexed: 01/07/2023] Open
Abstract
Species networks generalize the notion of species trees to allow for hybridization or other lateral gene transfer. Under the network multispecies coalescent model, individual gene trees arising from a network can have any topology, but arise with frequencies dependent on the network structure and numerical parameters. We propose a new algorithm for statistical inference of a level-1 species network under this model, from data consisting of gene tree topologies, and provide the theoretical justification for it. The algorithm is based on an analysis of quartets displayed on gene trees, combining several statistical hypothesis tests with combinatorial ideas such as a quartet-based intertaxon distance appropriate to networks, the NeighborNet algorithm for circular split systems, and the Circular Network algorithm for constructing a splits graph.
Collapse
|
79
|
Karimi N, Grover CE, Gallagher JP, Wendel JF, Ané C, Baum DA. Reticulate Evolution Helps Explain Apparent Homoplasy in Floral Biology and Pollination in Baobabs (Adansonia; Bombacoideae; Malvaceae). Syst Biol 2019; 69:462-478. [DOI: 10.1093/sysbio/syz073] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2019] [Revised: 10/24/2019] [Accepted: 10/26/2019] [Indexed: 12/17/2022] Open
Abstract
Abstract
Baobabs (Adansonia) are a cohesive group of tropical trees with a disjunct distribution in Australia, Madagascar, and continental Africa, and diverse flowers associated with two pollination modes. We used custom-targeted sequence capture in conjunction with new and existing phylogenetic comparative methods to explore the evolution of floral traits and pollination systems while allowing for reticulate evolution. Our analyses suggest that relationships in Adansonia are confounded by reticulation, with network inference methods supporting at least one reticulation event. The best supported hypothesis involves introgression between Adansonia rubrostipa and core Longitubae, both of which are hawkmoth pollinated with yellow/red flowers, but there is also some support for introgression between the African lineage and Malagasy Brevitubae, which are both mammal-pollinated with white flowers. New comparative methods for phylogenetic networks were developed that allow maximum-likelihood inference of ancestral states and were applied to study the apparent homoplasy in floral biology and pollination mode seen in Adansonia. This analysis supports a role for introgressive hybridization in morphological evolution even in a clade with highly divergent and geographically widespread species. Our new comparative methods for discrete traits on species networks are implemented in the software PhyloNetworks. [Comparative methods; Hyb-Seq; introgression; network inference; population trees; reticulate evolution; species tree inference; targeted sequence capture.]
Collapse
Affiliation(s)
- Nisa Karimi
- Department of Botany, University of Wisconsin – Madison, 430 Lincoln Drive, Madison, WI 53706, USA
| | - Corrinne E Grover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, 2200 Osborn Drive, Ames, IA 50011, USA
| | - Joseph P Gallagher
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, 2200 Osborn Drive, Ames, IA 50011, USA
- Department of Biology, University of Massachusetts, 611 North Pleasant Street, Amherst, MA 01003, USA
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, 2200 Osborn Drive, Ames, IA 50011, USA
| | - Cécile Ané
- Department of Botany, University of Wisconsin – Madison, 430 Lincoln Drive, Madison, WI 53706, USA
- Department of Statistics, University of Wisconsin – Madison, 1300 University Ave, WI, 53706, USA
| | - David A Baum
- Department of Botany, University of Wisconsin – Madison, 430 Lincoln Drive, Madison, WI 53706, USA
- Wisconsin Institute for Discovery, 330 N Orchard Street, Madison, 430 Lincoln Drive, Madison, WI 53706, USA
| |
Collapse
|
80
|
Dong S, Xiao Y, Kong H, Feng C, Harris A, Yan Y, Kang M. Nuclear loci developed from multiple transcriptomes yield high resolution in phylogeny of scaly tree ferns (Cyatheaceae) from China and Vietnam. Mol Phylogenet Evol 2019; 139:106567. [DOI: 10.1016/j.ympev.2019.106567] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2019] [Revised: 07/10/2019] [Accepted: 07/18/2019] [Indexed: 11/27/2022]
|
81
|
Blanco-Pastor JL, Bertrand YJK, Liberal IM, Wei Y, Brummer EC, Pfeil BE. Evolutionary networks from RADseq loci point to hybrid origins of Medicago carstiensis and Medicago cretacea. AMERICAN JOURNAL OF BOTANY 2019; 106:1219-1228. [PMID: 31535720 DOI: 10.1002/ajb2.1352] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2019] [Accepted: 07/12/2019] [Indexed: 06/10/2023]
Abstract
PREMISE Although hybridization has played an important role in the evolution of many plant species, phylogenetic reconstructions that include hybridizing lineages have been historically constrained by the available models and data. Restriction-site-associated DNA sequencing (RADseq) has been a popular sequencing technique for the reconstruction of hybridization in the next-generation sequencing era. However, the utility of RADseq for the reconstruction of complex evolutionary networks has not been thoroughly investigated. Conflicting phylogenetic relationships in the genus Medicago have been mainly attributed to hybridization, but the specific hybrid origins of taxa have not been yet clarified. METHODS We obtained new molecular data from diploid species of Medicago section Medicago using single-digest RADseq to reconstruct evolutionary networks from gene trees, an approach that is computationally tractable with data sets that include several species and complex hybridization patterns. RESULTS Our analyses revealed that assembly filters to exclusively select a small set of loci with high phylogenetic information led to the most-divergent network topologies. Conversely, alternative clustering thresholds or filters on the number of samples per locus had a lower impact on networks. A strong hybridization signal was detected for M. carstiensis and M. cretacea, while signals were less clear for M. rugosa, M. rhodopea, M. suffruticosa, M. marina, M. scutellata, and M. sativa. CONCLUSIONS Complex network reconstructions from RADseq gene trees were not robust under variations of the assembly parameters and filters. But when the most-divergent networks were discarded, all remaining analyses consistently supported a hybrid origin for M. carstiensis and M. cretacea.
Collapse
Affiliation(s)
- José Luis Blanco-Pastor
- Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, 40530, Göteborg, Sweden
- INRA, Centre Nouvelle-Aquitaine-Poitiers, UR4 (URP3F), 86600, Lusignan, France
| | - Yann J K Bertrand
- Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, 40530, Göteborg, Sweden
- Institute of Botany, Czech Academy of Sciences, Zámek 1, 25243, Průhonice, Czech Republic
| | | | - Yanling Wei
- Plant Breeding Center, Department of Plant Sciences, University of California, Davis, Davis, CA, USA
| | - E Charles Brummer
- Plant Breeding Center, Department of Plant Sciences, University of California, Davis, Davis, CA, USA
| | - Bernard E Pfeil
- Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, 40530, Göteborg, Sweden
| |
Collapse
|
82
|
Vandendorpe J, van Baak CGC, Stelbrink B, Delicado D, Albrecht C, Wilke T. Historical faunal exchange between the Pontocaspian Basin and North America. Ecol Evol 2019; 9:10816-10827. [PMID: 31632651 PMCID: PMC6787871 DOI: 10.1002/ece3.5602] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Revised: 08/05/2019] [Accepted: 08/06/2019] [Indexed: 12/05/2022] Open
Abstract
Ecrobia is a genus of small brackish-water mud snails with an amphi-Atlantic distribution. Interestingly, the species occurring in the northwestern Atlantic, Ecrobia truncata, is more closely related to the Pontocaspian taxa, Ecrobia grimmi and Ecrobia maritima, than to the species occurring in the northeastern Atlantic and Mediterranean Sea. At least three colonization scenarios may account for this peculiar biogeographical pattern: (1) a recent human-mediated dispersal, (2) a historical transatlantic interchange, and (3) a historical transpolar interchange. To test these three scenarios, we used five operational criteria-time of species divergence, first appearance in the fossil record, dispersal limitation as well as environmental filtering and biotic interactions along the potential migration routes. Specifically, we inferred a time-calibrated molecular phylogeny for Ecrobia and reconstructed a paleogeographical map of the Arctic Ocean at 2.5 million years ago (Mya). Based on the five operational criteria, scenarios 1 and 2 can likely be rejected. In contrast, all criteria support scenario 3 (historical transpolar interchange). It is therefore suggested that a bird-mediated and/or ocean current-mediated faunal interchange via the Arctic Ocean occurred during the Late Pliocene or Early Pleistocene. This dispersal was likely facilitated by reduced distances between the Eurasian and North American/Greenland landmasses, marine introgressions, and/or a stepping-stone system of brackish-water habitats in northern Siberia, as well as a lack of competition along the migration route. As for the direction of dispersal, the scientific data presented are not conclusive. However, there is clearly more support for the scenario of dispersal from the Pontocaspian Basin to North America than vice versa. This is the first study providing evidence for a natural faunal exchange between the Pontocaspian Basin and North America via the Arctic Ocean.
Collapse
Affiliation(s)
- Justine Vandendorpe
- Department of Animal Ecology and SystematicsJustus Liebig University GiessenGiessenGermany
| | | | - Björn Stelbrink
- Department of Animal Ecology and SystematicsJustus Liebig University GiessenGiessenGermany
| | - Diana Delicado
- Department of Animal Ecology and SystematicsJustus Liebig University GiessenGiessenGermany
| | - Christian Albrecht
- Department of Animal Ecology and SystematicsJustus Liebig University GiessenGiessenGermany
| | - Thomas Wilke
- Department of Animal Ecology and SystematicsJustus Liebig University GiessenGiessenGermany
| |
Collapse
|
83
|
Zhu J, Nakhleh L. Inference of species phylogenies from bi-allelic markers using pseudo-likelihood. Bioinformatics 2019; 34:i376-i385. [PMID: 29950004 PMCID: PMC6022577 DOI: 10.1093/bioinformatics/bty295] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Motivation Phylogenetic networks represent reticulate evolutionary histories. Statistical methods for their inference under the multispecies coalescent have recently been developed. A particularly powerful approach uses data that consist of bi-allelic markers (e.g. single nucleotide polymorphism data) and allows for exact likelihood computations of phylogenetic networks while numerically integrating over all possible gene trees per marker. While the approach has good accuracy in terms of estimating the network and its parameters, likelihood computations remain a major computational bottleneck and limit the method’s applicability. Results In this article, we first demonstrate why likelihood computations of networks take orders of magnitude more time when compared to trees. We then propose an approach for inference of phylogenetic networks based on pseudo-likelihood using bi-allelic markers. We demonstrate the scalability and accuracy of phylogenetic network inference via pseudo-likelihood computations on simulated data. Furthermore, we demonstrate aspects of robustness of the method to violations in the underlying assumptions of the employed statistical model. Finally, we demonstrate the application of the method to biological data. The proposed method allows for analyzing larger datasets in terms of the numbers of taxa and reticulation events. While pseudo-likelihood had been proposed before for data consisting of gene trees, the work here uses sequence data directly, offering several advantages as we discuss. Availability and implementation The methods have been implemented in PhyloNet (http://bioinfocs.rice.edu/phylonet).
Collapse
Affiliation(s)
- Jiafan Zhu
- Department of Computer Science, Rice University, Houston, TX, USA
| | - Luay Nakhleh
- Department of Computer Science, Rice University, Houston, TX, USA.,Department of BioSciences, Rice University, Houston, TX, USA
| |
Collapse
|
84
|
Wollenberg Valero KC, Marshall JC, Bastiaans E, Caccone A, Camargo A, Morando M, Niemiller ML, Pabijan M, Russello MA, Sinervo B, Werneck FP, Sites JW, Wiens JJ, Steinfartz S. Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians. Genes (Basel) 2019; 10:genes10090646. [PMID: 31455040 PMCID: PMC6769790 DOI: 10.3390/genes10090646] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Revised: 07/21/2019] [Accepted: 08/05/2019] [Indexed: 12/22/2022] Open
Abstract
In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.
Collapse
Affiliation(s)
| | - Jonathon C Marshall
- Department of Zoology, Weber State University, 1415 Edvalson Street, Dept. 2505, Ogden, UT 84401, USA
| | - Elizabeth Bastiaans
- Department of Biology, State University of New York, College at Oneonta, Oneonta, NY 13820, USA
| | - Adalgisa Caccone
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 06520, USA
| | - Arley Camargo
- Centro Universitario de Rivera, Universidad de la República, Ituzaingó 667, Rivera 40000, Uruguay
| | - Mariana Morando
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales (IPEEC, CENPAT-CONICET) Bv. Brown 2915, Puerto Madryn U9120ACD, Argentina
| | - Matthew L Niemiller
- Department of Biological Sciences, The University of Alabama in Huntsville, Huntsville, AL 35899, USA
| | - Maciej Pabijan
- Department of Comparative Anatomy, Institute of Zoology and Biomedical Research, Jagiellonian University, ul. Gronostajowa 9, 30-387 Kraków, Poland
| | - Michael A Russello
- Department of Biology, University of British Columbia, Okanagan Campus, 3247 University Way, Kelowna, BC V1V 1V7, Canada
| | - Barry Sinervo
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Coastal Biology Building, 130 McAllister Way, Santa Cruz, CA 95060, USA
| | - Fernanda P Werneck
- Programa de Coleções Científicas Biológicas, Coordenação de Biodiversidade, Instituto Nacional de Pesquisas da Amazônia, Manaus 69060-000, Brazil
| | - Jack W Sites
- Department of Biological and Marine Sciences, University of Hull, Cottingham Road, Hull HU6 7RX, UK
| | - John J Wiens
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | - Sebastian Steinfartz
- Molecular Evolution and Systematics of Animals, Institute of Biology, University of Leipzig, Talstrasse 33, 04103 Leipzig, Germany
| |
Collapse
|
85
|
Tan M, Long H, Liao B, Cao Z, Yuan D, Tian G, Zhuang J, Yang J. QS-Net: Reconstructing Phylogenetic Networks Based on Quartet and Sextet. Front Genet 2019; 10:607. [PMID: 31396256 PMCID: PMC6667645 DOI: 10.3389/fgene.2019.00607] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2018] [Accepted: 06/11/2019] [Indexed: 01/27/2023] Open
Abstract
Phylogenetic networks are used to estimate evolutionary relationships among biological entities or taxa involving reticulate events such as horizontal gene transfer, hybridization, recombination, and reassortment. In the past decade, many phylogenetic tree and network reconstruction methods have been proposed. Despite that they are highly accurate in reconstructing simple to moderate complex reticulate events, the performance decreases when several reticulate events are present simultaneously. In this paper, we proposed QS-Net, a phylogenetic network reconstruction method taking advantage of information on the relationship among six taxa. To evaluate the performance of QS-Net, we conducted experiments on three artificial sequence data simulated from an evolutionary tree, an evolutionary network involving three reticulate events, and a complex evolutionary network involving five reticulate events. Comparison with popular phylogenetic methods including Neighbor-Joining, Split-Decomposition, Neighbor-Net, and Quartet-Net suggests that QS-Net is comparable with other methods in reconstructing tree-like evolutionary histories, while it outperforms them in reconstructing reticulate events. In addition, we also applied QS-Net in real data including a bacterial taxonomy data consisting of 36 bacterial species and the whole genome sequences of 22 H7N9 influenza A viruses. The results indicate that QS-Net is capable of inferring commonly believed bacterial taxonomy and influenza evolution as well as identifying novel reticulate events. The software QS-Net is publically available at https://github.com/Tmyiri/QS-Net.
Collapse
Affiliation(s)
- Ming Tan
- College of Computer Science and Electronic Engineering, Hunan University, Changsha, China
| | - Haixia Long
- School of Information Science and Technology , Hainan Normal University, Haikou, China
| | - Bo Liao
- College of Computer Science and Electronic Engineering, Hunan University, Changsha, China.,School of Information Science and Technology , Hainan Normal University, Haikou, China
| | - Zhi Cao
- College of Computer Science and Electronic Engineering, Hunan University, Changsha, China
| | - Dawei Yuan
- Geneis (Beijing) Co. Ltd., Beijing, China
| | - Geng Tian
- Geneis (Beijing) Co. Ltd., Beijing, China
| | - Jujuan Zhuang
- Department of Mathematics, Dalian Martine University, Dalian, China
| | - Jialiang Yang
- School of Information Science and Technology , Hainan Normal University, Haikou, China.,Icahn Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, NY, United States
| |
Collapse
|
86
|
Huynh S, Marcussen T, Felber F, Parisod C. Hybridization preceded radiation in diploid wheats. Mol Phylogenet Evol 2019; 139:106554. [PMID: 31288105 DOI: 10.1016/j.ympev.2019.106554] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2019] [Revised: 07/03/2019] [Accepted: 07/03/2019] [Indexed: 01/06/2023]
Abstract
Evolutionary relationships among the Aegilops-Triticum relatives of cultivated wheats have been difficult to resolve owing to incomplete lineage sorting and reticulate evolution. Recent studies have suggested that the wheat D-genome lineage (progenitor of Ae. tauschii) originated through homoploid hybridization between the A-genome lineage (progenitor of Triticum s.str.) and the B-genome lineage (progenitor of Ae. speltoides). This scenario of reticulation has been debated, calling for adequate phylogenetic analyses based on comprehensive sampling. To reconstruct the evolution of Aegilops-Triticum diploids, we here combined high-throughput sequencing of 38 nuclear low-copy loci of multiple accessions of all 13 species with inferences of the species phylogeny using the full-parameterized MCMC_SEQ method. Phylogenies recovered a monophyletic Aegilops-Triticum lineage that began diversifying ~6.6 Ma ago and gave rise to four sublineages, i.e. the A- (2 species), B- (1 species), D- (9 species) and T- (Ae. mutica) genome lineage. Full-parameterized phylogenies as well as patterns of tree dilation and tree compression supported a hybrid origin of the D-genome lineage from A and B ~3.0-4.0 Ma ago, and did not indicate additional hybridization events. Conflicting ABBA-BABA tests suggestive of further reticulation were shown here to result from ancestral population structure rather than hybridization. This comprehensive and dated phylogeny of wheat relatives indicates that the origin of the hybrid D-genome was followed by intense diversification into the majority of extant diploid as well as allopolyploid wild wheats.
Collapse
Affiliation(s)
- Stella Huynh
- Institute of Biology, University of Neuchâtel, Switzerland
| | - Thomas Marcussen
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Norway
| | - François Felber
- Institute of Biology, University of Neuchâtel, Switzerland; Musée et Jardins botaniques cantonaux de Lausanne et Pont-de-Nant, Switzerland
| | | |
Collapse
|
87
|
Fungal species boundaries in the genomics era. Fungal Genet Biol 2019; 131:103249. [PMID: 31279976 DOI: 10.1016/j.fgb.2019.103249] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Revised: 05/21/2019] [Accepted: 06/28/2019] [Indexed: 12/30/2022]
Abstract
Genomic data has opened new possibilities to understand how organisms change over time, and could enable the discovery of previously undescribed species. Although taxonomy used to be based on phenotypes, molecular data has frequently revealed that morphological traits are insufficient to describe biodiversity. Genomics holds the promise of revealing even more genetic discontinuities, but the parameters on how to describe species from genomic data remain unclear. Fungi have been a successful case in which the use of molecular markers has uncovered the existence of genetic boundaries where no crosses are possible. In this minireview, we highlight recent advances, propose a set of standards to use genomic sequences to uncover species boundaries, point out potential pitfalls, and present possible future research directions.
Collapse
|
88
|
Thawornwattana Y, Dalquen D, Yang Z. Coalescent Analysis of Phylogenomic Data Confidently Resolves the Species Relationships in the Anopheles gambiae Species Complex. Mol Biol Evol 2019; 35:2512-2527. [PMID: 30102363 PMCID: PMC6188554 DOI: 10.1093/molbev/msy158] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Deep coalescence and introgression make it challenging to infer phylogenetic relationships among closely related species that arose through radiative speciation events. Despite numerous phylogenetic analyses and the availability of whole genomes, the phylogeny in the Anopheles gambiae species complex has not been confidently resolved. Here we extract over 80, 000 coding and noncoding short segments (called loci) from the genomes of six members of the species complex and use a Bayesian method under the multispecies coalescent model to infer the species tree, which takes into account genealogical heterogeneity across the genome and uncertainty in the gene trees. We obtained a robust estimate of the species tree from the distal region of the X chromosome: (A. merus, ((A. melas, (A. arabiensis, A. quadriannulatus)), (A. gambiae, A. coluzzii))), with A. merus to be the earliest branching species. This species tree agrees with the chromosome inversion phylogeny and provides a parsimonious interpretation of inversion and introgression events. Simulation informed by the real data suggest that the coalescent approach is reliable while the sliding-window analysis used in a previous phylogenomic study generates artifactual species trees. Likelihood ratio test of gene flow revealed strong evidence of autosomal introgression from A. arabiensis into A. gambiae (at the average rate of ∼0.2 migrants per generation), but not in the opposite direction, and introgression of the 3 L chromosomal region from A. merus into A. quadriannulatus. Our results highlight the importance of accommodating incomplete lineage sorting and introgression in phylogenomic analyses of species that arose through recent radiative speciation events.
Collapse
Affiliation(s)
- Yuttapong Thawornwattana
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom.,Department of Microbiology, Faculty of Science, Mahidol University, Bangkok, Thailand
| | - Daniel Dalquen
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom.,Radcliffe Institute for Advanced Studies, Harvard University, Cambridge, MA
| |
Collapse
|
89
|
Bryja J, Colangelo P, Lavrenchenko LA, Meheretu Y, Šumbera R, Bryjová A, Verheyen E, Leirs H, Castiglia R. Diversity and evolution of African Grass Rats (Muridae:
Arvicanthis
)—From radiation in East Africa to repeated colonization of northwestern and southeastern savannas. J ZOOL SYST EVOL RES 2019. [DOI: 10.1111/jzs.12290] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Josef Bryja
- Institute of Vertebrate Biology of the Czech Academy of Sciences Brno Czech Republic
- Department of Botany and Zoology, Faculty of Science Masaryk University Brno Czech Republic
| | - Paolo Colangelo
- National Research Council Institute of Agro‐environmental and Forest Biology (CNR‐IBAF) Rome Italy
| | - Leonid A. Lavrenchenko
- A.N.Severtsov Institute of Ecology and Evolution of the Russian Academy of Sciences Moscow Russia
| | - Yonas Meheretu
- Department of Biology, Institute of Mountain Research and Development Mekelle University Mekelle Ethiopia
| | - Radim Šumbera
- Department of Zoology, Faculty of Science University of South Bohemia České Budějovice Czech Republic
| | - Anna Bryjová
- Institute of Vertebrate Biology of the Czech Academy of Sciences Brno Czech Republic
| | - Erik Verheyen
- Operational Direction Taxonomy and Phylogeny Royal Belgian Institute for Natural Sciences Brussels Belgium
- Evolutionary Ecology Group, Biology Department University of Antwerp Antwerp Belgium
| | - Herwig Leirs
- Evolutionary Ecology Group, Biology Department University of Antwerp Antwerp Belgium
| | - Riccardo Castiglia
- Department of Biology and Biotechnology “Charles Darwin” Sapienza University of Rome Rome Italy
| |
Collapse
|
90
|
Bouckaert R, Vaughan TG, Barido-Sottani J, Duchêne S, Fourment M, Gavryushkina A, Heled J, Jones G, Kühnert D, De Maio N, Matschiner M, Mendes FK, Müller NF, Ogilvie HA, du Plessis L, Popinga A, Rambaut A, Rasmussen D, Siveroni I, Suchard MA, Wu CH, Xie D, Zhang C, Stadler T, Drummond AJ. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLoS Comput Biol 2019; 15:e1006650. [PMID: 30958812 PMCID: PMC6472827 DOI: 10.1371/journal.pcbi.1006650] [Citation(s) in RCA: 1610] [Impact Index Per Article: 322.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2018] [Revised: 04/18/2019] [Accepted: 02/04/2019] [Indexed: 11/18/2022] Open
Abstract
Elaboration of Bayesian phylogenetic inference methods has continued at pace in recent years with major new advances in nearly all aspects of the joint modelling of evolutionary data. It is increasingly appreciated that some evolutionary questions can only be adequately answered by combining evidence from multiple independent sources of data, including genome sequences, sampling dates, phenotypic data, radiocarbon dates, fossil occurrences, and biogeographic range information among others. Including all relevant data into a single joint model is very challenging both conceptually and computationally. Advanced computational software packages that allow robust development of compatible (sub-)models which can be composed into a full model hierarchy have played a key role in these developments. Developing such software frameworks is increasingly a major scientific activity in its own right, and comes with specific challenges, from practical software design, development and engineering challenges to statistical and conceptual modelling challenges. BEAST 2 is one such computational software platform, and was first announced over 4 years ago. Here we describe a series of major new developments in the BEAST 2 core platform and model hierarchy that have occurred since the first release of the software, culminating in the recent 2.5 release.
Collapse
Affiliation(s)
- Remco Bouckaert
- Centre of Computational Evolution, University of Auckland, Auckland, New Zealand
- Max Planck Institute for the Science of Human History, Jena, Germany
| | - Timothy G. Vaughan
- ETH Zürich, Department of Biosystems Science and Engineering, 4058 Basel, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Joëlle Barido-Sottani
- ETH Zürich, Department of Biosystems Science and Engineering, 4058 Basel, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Sebastián Duchêne
- Department of Biochemistry and Molecular Biology, University of Melbourne, Melbourne, Victoria, Australia
| | - Mathieu Fourment
- ithree institute, University of Technology Sydney, Sydney, Australia
| | | | | | - Graham Jones
- Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, SE 405 30 Göteborg, Sweden
| | - Denise Kühnert
- Max Planck Institute for the Science of Human History, Jena, Germany
| | - Nicola De Maio
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Cambridgeshire, UK
| | - Michael Matschiner
- Department of Environmental Sciences, University of Basel, 4051 Basel, Switzerland
| | - Fábio K. Mendes
- Centre of Computational Evolution, University of Auckland, Auckland, New Zealand
| | - Nicola F. Müller
- ETH Zürich, Department of Biosystems Science and Engineering, 4058 Basel, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Huw A. Ogilvie
- Department of Computer Science, Rice University, Houston, TX 77005-1892, USA
| | - Louis du Plessis
- Department of Zoology, University of Oxford, Oxford, OX1 3PS, UK
| | - Alex Popinga
- Centre of Computational Evolution, University of Auckland, Auckland, New Zealand
| | - Andrew Rambaut
- Institute of Evolutionary Biology, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 3FL UK
| | - David Rasmussen
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC 27695, USA
| | - Igor Siveroni
- Department of Infectious Disease Epidemiology, Imperial College London, Norfolk Place, W2 1PG, UK
| | - Marc A. Suchard
- Department of Biomathematics, David Geffen School of Medicine, University of California, Los Angeles, CA, USA
| | - Chieh-Hsi Wu
- Department of Statistics, University of Oxford, OX1 3LB, UK
| | - Dong Xie
- Centre of Computational Evolution, University of Auckland, Auckland, New Zealand
| | - Chi Zhang
- Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, Beijing, China
| | - Tanja Stadler
- ETH Zürich, Department of Biosystems Science and Engineering, 4058 Basel, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Alexei J. Drummond
- Centre of Computational Evolution, University of Auckland, Auckland, New Zealand
| |
Collapse
|
91
|
Bravo GA, Antonelli A, Bacon CD, Bartoszek K, Blom MPK, Huynh S, Jones G, Knowles LL, Lamichhaney S, Marcussen T, Morlon H, Nakhleh LK, Oxelman B, Pfeil B, Schliep A, Wahlberg N, Werneck FP, Wiedenhoeft J, Willows-Munro S, Edwards SV. Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics. PeerJ 2019; 7:e6399. [PMID: 30783571 PMCID: PMC6378093 DOI: 10.7717/peerj.6399] [Citation(s) in RCA: 67] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2018] [Accepted: 01/07/2019] [Indexed: 12/23/2022] Open
Abstract
Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.
Collapse
Affiliation(s)
- Gustavo A. Bravo
- Department of Organismic and Evolutionary Biology, Museum of Comparative Zoology, Harvard University, Cambridge, MA, USA
| | - Alexandre Antonelli
- Department of Organismic and Evolutionary Biology, Museum of Comparative Zoology, Harvard University, Cambridge, MA, USA
- Gothenburg Global Biodiversity Centre, Göteborg, Sweden
- Department of Biological and Environmental Sciences, University of Gothenburg, Göteborg, Sweden
- Gothenburg Botanical Garden, Göteborg, Sweden
| | - Christine D. Bacon
- Gothenburg Global Biodiversity Centre, Göteborg, Sweden
- Department of Biological and Environmental Sciences, University of Gothenburg, Göteborg, Sweden
| | - Krzysztof Bartoszek
- Department of Computer and Information Science, Linköping University, Linköping, Sweden
| | - Mozes P. K. Blom
- Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Stockholm, Sweden
| | - Stella Huynh
- Institut de Biologie, Université de Neuchâtel, Neuchâtel, Switzerland
| | - Graham Jones
- Department of Biological and Environmental Sciences, University of Gothenburg, Göteborg, Sweden
| | - L. Lacey Knowles
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, USA
| | - Sangeet Lamichhaney
- Department of Organismic and Evolutionary Biology, Museum of Comparative Zoology, Harvard University, Cambridge, MA, USA
| | - Thomas Marcussen
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Oslo, Norway
| | - Hélène Morlon
- Institut de Biologie, Ecole Normale Supérieure de Paris, Paris, France
| | - Luay K. Nakhleh
- Department of Computer Science, Rice University, Houston, TX, USA
| | - Bengt Oxelman
- Gothenburg Global Biodiversity Centre, Göteborg, Sweden
- Department of Biological and Environmental Sciences, University of Gothenburg, Göteborg, Sweden
| | - Bernard Pfeil
- Department of Biological and Environmental Sciences, University of Gothenburg, Göteborg, Sweden
| | - Alexander Schliep
- Department of Computer Science and Engineering, Chalmers University of Technology and University of Gothenburg, Göteborg, Sweden
| | | | - Fernanda P. Werneck
- Coordenação de Biodiversidade, Programa de Coleções Científicas Biológicas, Instituto Nacional de Pesquisa da Amazônia, Manaus, AM, Brazil
| | - John Wiedenhoeft
- Department of Computer Science and Engineering, Chalmers University of Technology and University of Gothenburg, Göteborg, Sweden
- Department of Computer Science, Rutgers University, Piscataway, NJ, USA
| | - Sandi Willows-Munro
- School of Life Sciences, University of Kwazulu-Natal, Pietermaritzburg, South Africa
| | - Scott V. Edwards
- Department of Organismic and Evolutionary Biology, Museum of Comparative Zoology, Harvard University, Cambridge, MA, USA
- Gothenburg Centre for Advanced Studies in Science and Technology, Chalmers University of Technology and University of Gothenburg, Göteborg, Sweden
| |
Collapse
|
92
|
Baños H. Identifying Species Network Features from Gene Tree Quartets Under the Coalescent Model. Bull Math Biol 2019; 81:494-534. [PMID: 30094772 PMCID: PMC6344282 DOI: 10.1007/s11538-018-0485-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2017] [Accepted: 07/30/2018] [Indexed: 10/28/2022]
Abstract
We show that many topological features of level-1 species networks are identifiable from the distribution of the gene tree quartets under the network multi-species coalescent model. In particular, every cycle of size at least 4 and every hybrid node in a cycle of size at least 5 are identifiable. This is a step toward justifying the inference of such networks which was recently implemented by Solís-Lemus and Ané. We show additionally how to compute quartet concordance factors for a network in terms of simpler networks, and explore some circumstances in which cycles of size 3 and hybrid nodes in 4-cycles can be detected.
Collapse
Affiliation(s)
- Hector Baños
- University of Alaska Fairbanks, P.O. Box 756660, Fairbanks, AK, 99775-6660, USA.
| |
Collapse
|
93
|
The Timing and Direction of Introgression Under the Multispecies Network Coalescent. Genetics 2019; 211:1059-1073. [PMID: 30670542 DOI: 10.1534/genetics.118.301831] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 01/21/2019] [Indexed: 12/26/2022] Open
Abstract
Introgression is a pervasive biological process, and many statistical methods have been developed to infer its presence from genomic data. However, many of the consequences and genomic signatures of introgression remain unexplored from a methodological standpoint. Here, we develop a model for the timing and direction of introgression based on the multispecies network coalescent, and from it suggest new approaches for testing introgression hypotheses. We suggest two new statistics, D 1 and D 2, which can be used in conjunction with other information to test hypotheses relating to the timing and direction of introgression, respectively. D 1 may find use in evaluating cases of homoploid hybrid speciation (HHS), while D 2 provides a four-taxon test for polarizing introgression. Although analytical expectations for our statistics require a number of assumptions to be met, we show how simulations can be used to test hypotheses about introgression when these assumptions are violated. We apply the D 1 statistic to genomic data from the wild yeast Saccharomyces paradoxus-a proposed example of HHS-demonstrating its use as a test of this model. These methods provide new and powerful ways to address questions relating to the timing and direction of introgression.
Collapse
|
94
|
Advances in Computational Methods for Phylogenetic Networks in the Presence of Hybridization. BIOINFORMATICS AND PHYLOGENETICS 2019. [DOI: 10.1007/978-3-030-10837-3_13] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
95
|
Jones GR. Divergence Estimation in the Presence of Incomplete Lineage Sorting and Migration. Syst Biol 2018; 68:19-31. [PMID: 29800341 DOI: 10.1093/sysbio/syy041] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2017] [Accepted: 05/16/2018] [Indexed: 11/14/2022] Open
Abstract
This article focuses on the problem of estimating a species tree from multilocus data in the presence of incomplete lineage sorting and migration. I develop a mathematical model similar to IMa2 (Hey 2010) for the relevant evolutionary processes which allows both the population size parameters and the migration rates between pairs of species tree branches to be integrated out. I then describe a BEAST2 package DENIM (Divergence estimation notwithstanding ILS and migration) which is based on this model and which uses an approximation to sample from the posterior. The approximation is based on the assumption that migrations are rare, and it only samples from certain regions of the posterior which seem likely given this assumption. The method breaks down if there is a lot of migration. Using simulations, Leaché et al. (2014) showed that using the standard multispecies coalescent model to infer a species tree can result in poor accuracy if migration is present. I reanalyze this simulated data to explore DENIM's performance and demonstrate substantial improvements in accuracy over *BEAST. I also reanalyze an empirical data set.
Collapse
Affiliation(s)
- Graham R Jones
- Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, SE 405 30 Göteborg, Sweden
| |
Collapse
|
96
|
Hey J, Chung Y, Sethuraman A, Lachance J, Tishkoff S, Sousa VC, Wang Y. Phylogeny Estimation by Integration over Isolation with Migration Models. Mol Biol Evol 2018; 35:2805-2818. [PMID: 30137463 PMCID: PMC6231491 DOI: 10.1093/molbev/msy162] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Phylogeny estimation is difficult for closely related populations and species, especially if they have been exchanging genes. We present a hierarchical Bayesian, Markov-chain Monte Carlo method with a state space that includes all possible phylogenies in a full Isolation-with-Migration model framework. The method is based on a new type of genealogy augmentation called a "hidden genealogy" that enables efficient updating of the phylogeny. This is the first likelihood-based method to fully incorporate directional gene flow and genetic drift for estimation of a species or population phylogeny. Application to human hunter-gatherer populations from Africa revealed a clear phylogenetic history, with strong support for gene exchange with an unsampled ghost population, and relatively ancient divergence between a ghost population and modern human populations, consistent with human/archaic divergence. In contrast, a study of five chimpanzee populations reveals a clear phylogeny with several pairs of populations having exchanged DNA, but does not support a history with an unsampled ghost population.
Collapse
Affiliation(s)
- Jody Hey
- Department of Biology, Center for Computational Genetics and Genomics, Temple University, Philadelphia, PA
| | - Yujin Chung
- Department of Biology, Center for Computational Genetics and Genomics, Temple University, Philadelphia, PA
- The Department of Applied Statistics, Kyonggi University, Suwon, South Korea
| | - Arun Sethuraman
- Department of Biology, Center for Computational Genetics and Genomics, Temple University, Philadelphia, PA
- Department of Biological Sciences, California State University San Marcos, San Marcos, CA
| | - Joseph Lachance
- Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA
- Georgia Institute of Technology, Atlanta, GA
| | - Sarah Tishkoff
- Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA
| | - Vitor C Sousa
- Department of Genetics, Rutgers, the State University of New Jersey, Piscataway, NJ
- University of Lisbon, Lisboa, Portugal
| | - Yong Wang
- Department of Genetics, Rutgers, the State University of New Jersey, Piscataway, NJ
- Ancestry, San Francisco, CA
| |
Collapse
|
97
|
Abstract
PhyloNet was released in 2008 as a software package for representing and analyzing phylogenetic networks. At the time of its release, the main functionalities in PhyloNet consisted of measures for comparing network topologies and a single heuristic for reconciling gene trees with a species tree. Since then, PhyloNet has grown significantly. The software package now includes a wide array of methods for inferring phylogenetic networks from data sets of unlinked loci while accounting for both reticulation (e.g., hybridization) and incomplete lineage sorting. In particular, PhyloNet now allows for maximum parsimony, maximum likelihood, and Bayesian inference of phylogenetic networks from gene tree estimates. Furthermore, Bayesian inference directly from sequence data (sequence alignments or biallelic markers) is implemented. Maximum parsimony is based on an extension of the "minimizing deep coalescences" criterion to phylogenetic networks, whereas maximum likelihood and Bayesian inference are based on the multispecies network coalescent. All methods allow for multiple individuals per species. As computing the likelihood of a phylogenetic network is computationally hard, PhyloNet allows for evaluation and inference of networks using a pseudolikelihood measure. PhyloNet summarizes the results of the various analyzes and generates phylogenetic networks in the extended Newick format that is readily viewable by existing visualization software.
Collapse
Affiliation(s)
| | | | | | - Luay Nakhleh
- Computer Science.,BioSciences, Rice University, 6100 Main Street, Houston, TX 77005, USA
| |
Collapse
|
98
|
Degnan JH. Modeling Hybridization Under the Network Multispecies Coalescent. Syst Biol 2018; 67:786-799. [PMID: 29846734 PMCID: PMC6101600 DOI: 10.1093/sysbio/syy040] [Citation(s) in RCA: 59] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2017] [Revised: 05/13/2018] [Accepted: 05/16/2018] [Indexed: 11/13/2022] Open
Abstract
Simultaneously modeling hybridization and the multispecies coalescent is becoming increasingly common, and inference of species networks in this context is now implemented in several software packages. This article addresses some of the conceptual issues and decisions to be made in this modeling, including whether or not to use branch lengths and issues with model identifiability. This article is based on a talk given at a Spotlight Session at Evolution 2017 meeting in Portland, Oregon. This session included several talks about modeling hybridization and gene flow in the presence of incomplete lineage sorting. Other talks given at this meeting are also included in this special issue of Systematic Biology.
Collapse
Affiliation(s)
- James H Degnan
- Department of Mathematics and Statistics, University of New Mexico, Albuquerque, NM 87131, USA
| |
Collapse
|
99
|
Finding a most parsimonious or likely tree in a network with respect to an alignment. J Math Biol 2018; 78:527-547. [PMID: 30121824 PMCID: PMC6437133 DOI: 10.1007/s00285-018-1282-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2017] [Revised: 05/16/2018] [Indexed: 11/22/2022]
Abstract
Phylogenetic networks are often constructed by merging multiple conflicting phylogenetic signals into a directed acyclic graph. It is interesting to explore whether a network constructed in this way induces biologically-relevant phylogenetic signals that were not present in the input. Here we show that, given a multiple alignment A for a set of taxa X and a rooted phylogenetic network N whose leaves are labelled by X, it is NP-hard to locate a most parsimonious phylogenetic tree displayed by N (with respect to A) even when the level of N—the maximum number of reticulation nodes within a biconnected component—is 1 and A contains only 2 distinct states. (If, additionally, gaps are allowed the problem becomes APX-hard.) We also show that under the same conditions, and assuming a simple binary symmetric model of character evolution, finding a most likely tree displayed by the network is NP-hard. These negative results contrast with earlier work on parsimony in which it is shown that if A consists of a single column the problem is fixed parameter tractable in the level. We conclude with a discussion of why, despite the NP-hardness, both the parsimony and likelihood problem can likely be well-solved in practice.
Collapse
|
100
|
Exploring the Tiers of Rooted Phylogenetic Network Space Using Tail Moves. Bull Math Biol 2018; 80:2177-2208. [PMID: 29948885 PMCID: PMC6061524 DOI: 10.1007/s11538-018-0452-0] [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: 08/17/2017] [Accepted: 06/04/2018] [Indexed: 11/25/2022]
Abstract
Popular methods for exploring the space of rooted phylogenetic trees use rearrangement moves such as rooted Nearest Neighbour Interchange (rNNI) and rooted Subtree Prune and Regraft (rSPR). Recently, these moves were generalized to rooted phylogenetic networks, which are a more suitable representation of reticulate evolutionary histories, and it was shown that any two rooted phylogenetic networks of the same complexity are connected by a sequence of either rSPR or rNNI moves. Here, we show that this is possible using only tail moves, which are a restricted version of rSPR moves on networks that are more closely related to rSPR moves on trees. The connectedness still holds even when we restrict to distance-1 tail moves (a localized version of tail moves). Moreover, we give bounds on the number of (distance-1) tail moves necessary to turn one network into another, which in turn yield new bounds for rSPR, rNNI and SPR (i.e. the equivalent of rSPR on unrooted networks). The upper bounds are constructive, meaning that we can actually find a sequence with at most this length for any pair of networks. Finally, we show that finding a shortest sequence of tail or rSPR moves is NP-hard.
Collapse
|