1
|
Ortiz D, Pekár S, Dianat M. A consequential one-night stand: Episodic historical hybridization leads to mitochondrial takeover in sympatric desert ant-eating spiders. Mol Phylogenet Evol 2024; 199:108167. [PMID: 39103025 DOI: 10.1016/j.ympev.2024.108167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Revised: 07/03/2024] [Accepted: 08/01/2024] [Indexed: 08/07/2024]
Abstract
Disentangling the genomic intricacies underlying speciation and the causes of discordance between sources of evidence can offer remarkable insights into evolutionary dynamics. The ant-eating spider Zodarion nitidum, found across the Middle East and Egypt, displays yellowish and blackish morphs that co-occur sympatrically. These morphs additionally differ in behavioral and physiological features and show complete pre-mating reproductive isolation. In contrast, they possess similar sexual features and lack distinct differences in their mitochondrial DNA. We analyzed both Z. nitidum morphs and outgroups using genome-wide and additional mitochondrial DNA data. The genomic evidence indicated that Yellow and Black are reciprocally independent lineages without signs of recent admixture. Interestingly, the sister group of Yellow is not Black but Z. luctuosum, a morphologically distinct species. Genomic gene flow analyses pinpointed an asymmetric nuclear introgression event, with Yellow contributing nearly 5 % of its genome to Black roughly 320,000 years ago, intriguingly aligning with the independently estimated origin of the mitochondrial DNA of Black. We conclude that the blackish and yellowish morphs of Z. nitidum are long-diverged distinct species, and that the ancient and modest genomic introgression event registered resulted in a complete mitochondrial takeover of Black by Yellow. This investigation underscores the profound long-term effects that even modest hybridization events can have on the genome of organisms. It also exemplifies the utility of phylogenetic networks for estimating historical events and how integrating independent lines of evidence can increase the reliability of such estimations.
Collapse
Affiliation(s)
- David Ortiz
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czechia.
| | - Stano Pekár
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czechia
| | - Malahat Dianat
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czechia; Institute of Vertebrate Biology, Czech Academy of Sciences, Brno, Czechia
| |
Collapse
|
2
|
Smith ML, Hahn MW. Selection leads to false inferences of introgression using popular methods. Genetics 2024; 227:iyae089. [PMID: 38805070 DOI: 10.1093/genetics/iyae089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2023] [Revised: 10/28/2023] [Accepted: 05/21/2024] [Indexed: 05/29/2024] Open
Abstract
Detecting introgression between closely related populations or species is a fundamental objective in evolutionary biology. Existing methods for detecting migration and inferring migration rates from population genetic data often assume a neutral model of evolution. Growing evidence of the pervasive impact of selection on large portions of the genome across diverse taxa suggests that this assumption is unrealistic in most empirical systems. Further, ignoring selection has previously been shown to negatively impact demographic inferences (e.g. of population size histories). However, the impacts of biologically realistic selection on inferences of migration remain poorly explored. Here, we simulate data under models of background selection, selective sweeps, balancing selection, and adaptive introgression. We show that ignoring selection sometimes leads to false inferences of migration in popularly used methods that rely on the site frequency spectrum. Specifically, balancing selection and some models of background selection result in the rejection of isolation-only models in favor of isolation-with-migration models and lead to elevated estimates of migration rates. BPP, a method that analyzes sequence data directly, showed false positives for all conditions at recent divergence times, but balancing selection also led to false positives at medium-divergence times. Our results suggest that such methods may be unreliable in some empirical systems, such that new methods that are robust to selection need to be developed.
Collapse
Affiliation(s)
- Megan L Smith
- Department of Biological Sciences, Mississippi State University, Starkville, MS 39762, USA
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
- Department of Computer Science, Indiana University, Bloomington, IN 47405, USA
| |
Collapse
|
3
|
Pavón-Vázquez CJ, Rana Q, Farleigh K, Crispo E, Zeng M, Liliah J, Mulcahy D, Ascanio A, Jezkova T, Leaché AD, Flouri T, Yang Z, Blair C. Gene Flow and Isolation in the Arid Nearctic Revealed by Genomic Analyses of Desert Spiny Lizards. Syst Biol 2024; 73:323-342. [PMID: 38190300 DOI: 10.1093/sysbio/syae001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Revised: 12/18/2023] [Accepted: 01/05/2024] [Indexed: 01/10/2024] Open
Abstract
The opposing forces of gene flow and isolation are two major processes shaping genetic diversity. Understanding how these vary across space and time is necessary to identify the environmental features that promote diversification. The detection of considerable geographic structure in taxa from the arid Nearctic has prompted research into the drivers of isolation in the region. Several geographic features have been proposed as barriers to gene flow, including the Colorado River, Western Continental Divide (WCD), and a hypothetical Mid-Peninsular Seaway in Baja California. However, recent studies suggest that the role of barriers in genetic differentiation may have been overestimated when compared to other mechanisms of divergence. In this study, we infer historical and spatial patterns of connectivity and isolation in Desert Spiny Lizards (Sceloporus magister) and Baja Spiny Lizards (Sceloporus zosteromus), which together form a species complex composed of parapatric lineages with wide distributions in arid western North America. Our analyses incorporate mitochondrial sequences, genomic-scale data, and past and present climatic data to evaluate the nature and strength of barriers to gene flow in the region. Our approach relies on estimates of migration under the multispecies coalescent to understand the history of lineage divergence in the face of gene flow. Results show that the S. magister complex is geographically structured, but we also detect instances of gene flow. The WCD is a strong barrier to gene flow, while the Colorado River is more permeable. Analyses yield conflicting results for the catalyst of differentiation of peninsular lineages in S. zosteromus. Our study shows how large-scale genomic data for thoroughly sampled species can shed new light on biogeography. Furthermore, our approach highlights the need for the combined analysis of multiple sources of evidence to adequately characterize the drivers of divergence.
Collapse
Affiliation(s)
- Carlos J Pavón-Vázquez
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY 11201, USA
- Facultad de Estudios Superiores Iztacala, Universidad Nacional Autónoma de México, Colonia Los Reyes Ixtacala, Tlalnepantla, Estado de México, C.P. 54090, México
| | - Qaantah Rana
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY 11201, USA
| | - Keaka Farleigh
- Department of Biology, Miami University, Oxford, OH 45056, USA
| | - Erika Crispo
- Department of Biology, Pace University, One Pace Plaza, New York, NY 10038, USA
| | - Mimi Zeng
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY 11201, USA
| | - Jeevanie Liliah
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY 11201, USA
| | - Daniel Mulcahy
- Collection Future, Museum für Naturkunde, Leibniz-Institute for Evolution and Biodiversity Science, Berlin 10115, Germany
| | - Alfredo Ascanio
- Department of Biology, Miami University, Oxford, OH 45056, USA
| | - Tereza Jezkova
- Department of Biology, Miami University, Oxford, OH 45056, USA
| | - Adam D Leaché
- Department of Biology & Burke Museum of Natural History and Culture, University of Washington, Seattle, WA 98195, USA
| | - Tomas Flouri
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - 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
| |
Collapse
|
4
|
Braichenko S, Borges R, Kosiol C. Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion. Mol Biol Evol 2024; 41:msae138. [PMID: 38980178 PMCID: PMC11272101 DOI: 10.1093/molbev/msae138] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 04/19/2024] [Accepted: 07/06/2024] [Indexed: 07/10/2024] Open
Abstract
The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.
Collapse
Affiliation(s)
- Svitlana Braichenko
- Centre for Biological Diversity, School of Biology, University of St Andrews, Fife KY16 9TH, UK
- Institute of Genetics and Cancer, University of Edinburgh, Edinburgh EH4 2XU, UK
| | - Rui Borges
- Institut für Populationsgenetik, Vetmeduni Vienna, Wien 1210, Austria
| | - Carolin Kosiol
- Centre for Biological Diversity, School of Biology, University of St Andrews, Fife KY16 9TH, UK
| |
Collapse
|
5
|
Pang XX, Zhang DY. Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information. Syst Biol 2024; 73:207-222. [PMID: 38224495 PMCID: PMC11129598 DOI: 10.1093/sysbio/syad077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Revised: 12/17/2023] [Accepted: 12/27/2023] [Indexed: 01/17/2024] Open
Abstract
In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.
Collapse
Affiliation(s)
- Xiao-Xu Pang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
| | - Da-Yong Zhang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
6
|
Rosser N, Seixas F, Queste LM, Cama B, Mori-Pezo R, Kryvokhyzha D, Nelson M, Waite-Hudson R, Goringe M, Costa M, Elias M, Mendes Eleres de Figueiredo C, Freitas AVL, Joron M, Kozak K, Lamas G, Martins ARP, McMillan WO, Ready J, Rueda-Muñoz N, Salazar C, Salazar P, Schulz S, Shirai LT, Silva-Brandão KL, Mallet J, Dasmahapatra KK. Hybrid speciation driven by multilocus introgression of ecological traits. Nature 2024; 628:811-817. [PMID: 38632397 PMCID: PMC11041799 DOI: 10.1038/s41586-024-07263-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 03/01/2024] [Indexed: 04/19/2024]
Abstract
Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.
Collapse
Affiliation(s)
- Neil Rosser
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA.
- Department of Biology, University of York, York, UK.
| | - Fernando Seixas
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
| | | | - Bruna Cama
- Department of Biology, University of York, York, UK
| | - Ronald Mori-Pezo
- URKU Estudios Amazónicos, Tarapoto, Perú
- Universidad Nacional Autónoma de Alto Amazona, Yurimaguas, Perú
| | - Dmytro Kryvokhyzha
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
- Department of Clinical Sciences, Lund University Diabetes Centre, Malmö, Sweden
| | | | | | - Matt Goringe
- Department of Biology, University of York, York, UK
| | | | - Marianne Elias
- Institut Systématique, Evolution, Biodiversité, UMR 7205 MNHN-CNRS-EPHE-UPMC Sorbonne Universités, Muséum National d'Histoire Naturelle, Paris, France
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Clarisse Mendes Eleres de Figueiredo
- Institute for Biological Sciences, Federal University of Pará (UFPA), Belém, Brazil
- Centre for Advanced Studies of Biodiversity (CEABIO), Belém, Brazil
| | - André Victor Lucci Freitas
- Departamento de Biologia Animal and Museu de Diversidade Biológica, Instituto de Biologia, Universidade Estadual de Campinas, São Paulo, Brazil
| | - Mathieu Joron
- Centre d'Ecologie Fonctionnelle et Evolutive, UMR 5175 CNRS, Université de Montpellier-Université Paul Valéry Montpellier-EPHE, Montpellier, France
| | - Krzysztof Kozak
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Gerardo Lamas
- Museo de Historia Natural, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | | | - W Owen McMillan
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Jonathan Ready
- Institute for Biological Sciences, Federal University of Pará (UFPA), Belém, Brazil
- Centre for Advanced Studies of Biodiversity (CEABIO), Belém, Brazil
| | - Nicol Rueda-Muñoz
- Biology Program, Faculty of Natural Sciences, Universidad del Rosario, Bogotá, Colombia
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences, Universidad del Rosario, Bogotá, Colombia
| | - Patricio Salazar
- Ecology and Evolutionary Biology, School of Biosciences, University of Sheffield, Sheffield, UK
| | - Stefan Schulz
- Institut für Organische Chemie, Technische Universität Braunschweig, Braunschweig, Germany
| | - Leila T Shirai
- Departamento de Biologia Animal and Museu de Diversidade Biológica, Instituto de Biologia, Universidade Estadual de Campinas, São Paulo, Brazil
| | - Karina L Silva-Brandão
- Leibniz Institute for the Analysis of Biodiversity Change, Museum de Natur Hamburg Zoology, Hamburg, Germany
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA.
| | - Kanchon K Dasmahapatra
- Department of Biology, University of York, York, UK
- Leverhulme Centre for Anthropocene Biodiversity, Department of Biology, University of York, York, UK
| |
Collapse
|
7
|
Sgarlata GM, Rasolondraibe E, Salmona J, Le Pors B, Ralantoharijaona T, Rakotonanahary A, Jan F, Manzi S, Iribar A, Zaonarivelo JR, Volasoa Andriaholinirina N, Rasoloharijaona S, Chikhi L. The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species. Mol Phylogenet Evol 2024; 193:107997. [PMID: 38128795 DOI: 10.1016/j.ympev.2023.107997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2023] [Revised: 12/06/2023] [Accepted: 12/18/2023] [Indexed: 12/23/2023]
Abstract
Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.
Collapse
Affiliation(s)
| | - Emmanuel Rasolondraibe
- Département de Biologie Animale et Ecologie, Faculté des Sciences, Université de Mahajanga, Mahajanga, Madagascar.
| | - Jordi Salmona
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande, 6, 2780-156 Oeiras, Portugal; Centre de Recherche sur la Biodiversité et l'Environnement (CRBE),Université de Toulouse, CNRS, IRD, Toulouse INP, Université Toulouse 3 -Paul Sabatier (UT3), Toulouse, France.
| | - Barbara Le Pors
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande, 6, 2780-156 Oeiras, Portugal
| | - Tantely Ralantoharijaona
- Département de Biologie Animale et Ecologie, Faculté des Sciences, Université de Mahajanga, Mahajanga, Madagascar
| | - Ando Rakotonanahary
- Département de Biologie Animale et Ecologie, Faculté des Sciences, Université de Mahajanga, Mahajanga, Madagascar.
| | - Fabien Jan
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande, 6, 2780-156 Oeiras, Portugal
| | - Sophie Manzi
- Centre de Recherche sur la Biodiversité et l'Environnement (CRBE),Université de Toulouse, CNRS, IRD, Toulouse INP, Université Toulouse 3 -Paul Sabatier (UT3), Toulouse, France.
| | - Amaia Iribar
- Centre de Recherche sur la Biodiversité et l'Environnement (CRBE),Université de Toulouse, CNRS, IRD, Toulouse INP, Université Toulouse 3 -Paul Sabatier (UT3), Toulouse, France.
| | - John Rigobert Zaonarivelo
- Département des Sciences de la Nature et de l'Environnement, Université d'Antsiranana, 201 Antsiranana, Madagascar.
| | | | - Solofonirina Rasoloharijaona
- Département de Biologie Animale et Ecologie, Faculté des Sciences, Université de Mahajanga, Mahajanga, Madagascar
| | - Lounès Chikhi
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande, 6, 2780-156 Oeiras, Portugal; Centre de Recherche sur la Biodiversité et l'Environnement (CRBE),Université de Toulouse, CNRS, IRD, Toulouse INP, Université Toulouse 3 -Paul Sabatier (UT3), Toulouse, France.
| |
Collapse
|
8
|
Sánchez KI, Recknagel H, Elmer KR, Avila LJ, Morando M. Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group). Evolution 2024; 78:716-733. [PMID: 38262697 DOI: 10.1093/evolut/qpae009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Revised: 12/20/2023] [Accepted: 01/17/2024] [Indexed: 01/25/2024]
Abstract
Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.
Collapse
Affiliation(s)
- Kevin I Sánchez
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales, Consejo Nacional de Investigaciones Científicas y Técnicas (IPEEC-CONICET), Puerto Madryn, Chubut, Argentina
| | - Hans Recknagel
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Kathryn R Elmer
- School of Biodiversity, One Health and Veterinary Medicine, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, United Kingdom
| | - Luciano J Avila
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales, Consejo Nacional de Investigaciones Científicas y Técnicas (IPEEC-CONICET), Puerto Madryn, Chubut, Argentina
| | - Mariana Morando
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales, Consejo Nacional de Investigaciones Científicas y Técnicas (IPEEC-CONICET), Puerto Madryn, Chubut, Argentina
- Departamento de Biología y Ambiente, Universidad Nacional de la Patagonia San Juan Bosco, Sede Puerto Madryn, Puerto Madryn, Chubut, Argentina
| |
Collapse
|
9
|
Wang X, Ye X, Li R, Zai X, Hu M, Wang S, Ren H, Jin Y, Xu J, Yue J. Spatio-temporal spread and evolution of Lassa virus in West Africa. BMC Infect Dis 2024; 24:314. [PMID: 38486143 PMCID: PMC10941413 DOI: 10.1186/s12879-024-09200-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 03/06/2024] [Indexed: 03/17/2024] Open
Abstract
BACKGROUND Lassa fever is a hemorrhagic disease caused by Lassa virus (LASV), which has been classified by the World Health Organization as one of the top infectious diseases requiring prioritized research. Previous studies have provided insights into the classification and geographic characteristics of LASV lineages. However, the factor of the distribution and evolution characteristics and phylodynamics of the virus was still limited. METHODS To enhance comprehensive understanding of LASV, we employed phylogenetic analysis, reassortment and recombination detection, and variation evaluation utilizing publicly available viral genome sequences. RESULTS The results showed the estimated the root of time of the most recent common ancestor (TMRCA) for large (L) segment was approximately 634 (95% HPD: [385879]), whereas the TMRCA for small (S) segment was around 1224 (95% HPD: [10301401]). LASV primarily spread from east to west in West Africa through two routes, and in route 2, the virus independently spread to surrounding countries through Liberia, resulting in a wider spread of LASV. From 1969 to 2018, the effective population size experienced two significant increased, indicating the enhanced genetic diversity of LASV. We also found the evolution rate of L segment was faster than S segment, further results showed zinc-binding protein had the fastest evolution rate. Reassortment events were detected in multiple lineages including sub-lineage IIg, while recombination events were observed within lineage V. Significant amino acid changes in the glycoprotein precursor of LASV were identified, demonstrating sequence diversity among lineages in LASV. CONCLUSION This study comprehensively elucidated the transmission and evolution of LASV in West Africa, providing detailed insights into reassortment events, recombination events, and amino acid variations.
Collapse
Affiliation(s)
- Xia Wang
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China
- Medical College of Guizhou University, Guiyang, 550025, China
| | - Xianwei Ye
- Medical College of Guizhou University, Guiyang, 550025, China
- Guizhou Provincial People's Hospital, Guiyang, 550002, China
| | - Ruihua Li
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China
| | - Xiaodong Zai
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China
| | - Mingda Hu
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China
| | - Shaoyan Wang
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China
| | - Hongguang Ren
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China
| | - Yuan Jin
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China.
| | - Junjie Xu
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China.
| | - Junjie Yue
- Laboratory of Advanced Biotechnology & State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, 100071, China.
| |
Collapse
|
10
|
Morita K, Saito T, Uechi T, Sawada N, Miura O. Out of the ancient lake: Multiple riverine colonizations and diversification of the freshwater snails in the genus Semisulcospira around Lake Biwa. Mol Phylogenet Evol 2024; 191:107987. [PMID: 38081401 DOI: 10.1016/j.ympev.2023.107987] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2023] [Revised: 11/20/2023] [Accepted: 12/07/2023] [Indexed: 12/18/2023]
Abstract
Ancient lakes are a hotspot of biodiversity. Freshwater species often experience spectacular species radiation after colonizing lakes from riverine habitats. Therefore, the relationship between the fauna of the ancient lakes and the surrounding riverine system has a special significance in understanding their origin and evolutionary history. The study of ancient lake species often focused on the lake colonization of riverine species. In contrast, far less attention has been placed on the reverse direction: the riverine colonization of the lake species, despite its importance in disentangling their complex evolutionary history. The freshwater snails in the genus Semisulcospira involve endemic groups that radiated in the ancient Lake Biwa. Using genetics and fossil records, we inferred that the ancestors of these lake-endemic Semisulcospira snails historically colonized the riverine habitats at least three times during the Middle Pleistocene. Each colonization resulted in the formation of a new lineage that was genetically and morphologically distinct from other lineages. Further, one of these colonizations was followed by hybridization with a cosmopolitan riverine species, which potentially facilitated the population persistence of the colonizers in the new environment. Despite their complex histories, all these colonizers were currently grouped within a single species, Semisulcospira kurodai, suggesting cryptic diversity in this species. This study highlights the significance of the riverine colonizations of the lake species to fully understand the diversification history of freshwater fauna in and around the ancient lakes.
Collapse
Affiliation(s)
- Kohei Morita
- Faculty of Agriculture and Marine Science, Kochi University, 200 Monobe, Nankoku, Kochi 783-8502 Japan
| | - Takumi Saito
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
| | - Takeru Uechi
- Major in Environmental Management, Graduate School of Agriculture, Kindai University, 3327-204 Nakamachi, Nara 631-8505, Japan
| | - Naoto Sawada
- Department of Zoology, Graduate School of Science, Kyoto University, Kitashirakawa-oiwakecho, Sakyo, Kyoto, Kyoto 606-8502 Japan
| | - Osamu Miura
- Faculty of Agriculture and Marine Science, Kochi University, 200 Monobe, Nankoku, Kochi 783-8502 Japan.
| |
Collapse
|
11
|
Thawornwattana Y, Seixas F, Yang Z, Mallet J. Major patterns in the introgression history of Heliconius butterflies. eLife 2023; 12:RP90656. [PMID: 38108819 PMCID: PMC10727504 DOI: 10.7554/elife.90656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2023] Open
Abstract
Gene flow between species, although usually deleterious, is an important evolutionary process that can facilitate adaptation and lead to species diversification. It also makes estimation of species relationships difficult. Here, we use the full-likelihood multispecies coalescent (MSC) approach to estimate species phylogeny and major introgression events in Heliconius butterflies from whole-genome sequence data. We obtain a robust estimate of species branching order among major clades in the genus, including the 'melpomene-silvaniform' group, which shows extensive historical and ongoing gene flow. We obtain chromosome-level estimates of key parameters in the species phylogeny, including species divergence times, present-day and ancestral population sizes, as well as the direction, timing, and intensity of gene flow. Our analysis leads to a phylogeny with introgression events that differ from those obtained in previous studies. We find that Heliconius aoede most likely represents the earliest-branching lineage of the genus and that 'silvaniform' species are paraphyletic within the melpomene-silvaniform group. Our phylogeny provides new, parsimonious histories for the origins of key traits in Heliconius, including pollen feeding and an inversion involved in wing pattern mimicry. Our results demonstrate the power and feasibility of the full-likelihood MSC approach for estimating species phylogeny and key population parameters despite extensive gene flow. The methods used here should be useful for analysis of other difficult species groups with high rates of introgression.
Collapse
Affiliation(s)
| | - Fernando Seixas
- Department of Organismic and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College LondonLondonUnited Kingdom
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| |
Collapse
|
12
|
San Jose M, Doorenweerd C, Geib S, Barr N, Dupuis JR, Leblanc L, Kauwe A, Morris KY, Rubinoff D. Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex. Mol Phylogenet Evol 2023; 188:107892. [PMID: 37524217 DOI: 10.1016/j.ympev.2023.107892] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 07/07/2023] [Accepted: 07/28/2023] [Indexed: 08/02/2023]
Abstract
As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.
Collapse
Affiliation(s)
- Michael San Jose
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA.
| | - Camiel Doorenweerd
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA
| | - Scott Geib
- Tropical Crop and Commodity Protection Research Unit, Daniel K Inouye U.S. Pacific Basin Agricultural Center, USDA Agricultural Research Services, Hilo, HI, USA
| | - Norman Barr
- United States Department of Agriculture, Animal and Plant Health Inspection Service, Plant Protection and Quarantine, Science & Technology, Insect Management and Molecular Diagnostics Laboratory, 22675 N. Moorefield Road, Edinburg, TX 78541, USA
| | - Julian R Dupuis
- University of Kentucky, Department of Entomology, S225 Ag Science Center North, 1100 South Limestone, Lexington, KY, 40546-0091, USA
| | - Luc Leblanc
- University of Idaho, Department of Entomology, Plant Pathology and Nematology, 875 Perimeter Drive, MS2329, Moscow, ID, 83844-2329, USA
| | - Angela Kauwe
- Tropical Crop and Commodity Protection Research Unit, Daniel K Inouye U.S. Pacific Basin Agricultural Center, USDA Agricultural Research Services, Hilo, HI, USA
| | - Kimberley Y Morris
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA; Tropical Crop and Commodity Protection Research Unit, Daniel K Inouye U.S. Pacific Basin Agricultural Center, USDA Agricultural Research Services, Hilo, HI, USA
| | - Daniel Rubinoff
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA
| |
Collapse
|
13
|
Flouri T, Jiao X, Huang J, Rannala B, Yang Z. Efficient Bayesian inference under the multispecies coalescent with migration. Proc Natl Acad Sci U S A 2023; 120:e2310708120. [PMID: 37871206 PMCID: PMC10622872 DOI: 10.1073/pnas.2310708120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2023] [Accepted: 08/15/2023] [Indexed: 10/25/2023] Open
Abstract
Analyses of genome sequence data have revealed pervasive interspecific gene flow and enriched our understanding of the role of gene flow in speciation and adaptation. Inference of gene flow using genomic data requires powerful statistical methods. Yet current likelihood-based methods involve heavy computation and are feasible for small datasets only. Here, we implement the multispecies-coalescent-with-migration model in the Bayesian program bpp, which can be used to test for gene flow and estimate migration rates, as well as species divergence times and population sizes. We develop Markov chain Monte Carlo algorithms for efficient sampling from the posterior, enabling the analysis of genome-scale datasets with thousands of loci. Implementation of both introgression and migration models in the same program allows us to test whether gene flow occurred continuously over time or in pulses. Analyses of genomic data from Anopheles mosquitoes demonstrate rich information in typical genomic datasets about the mode and rate of gene flow.
Collapse
Affiliation(s)
- Tomáš Flouri
- Department of Genetics, Evolution, and Environment, University College London, LondonWC1E 6BT, United Kingdom
| | - Xiyun Jiao
- Department of Statistics and Data Science, China Southern University of Science and Technology, Shenzhen518055, China
| | - Jun Huang
- Department of Intelligent Medical Engineering, School of Biomedical Engineering, Capital Medical University, Beijing100069, China
| | - Bruce Rannala
- Department of Evolution and Ecology, University of California, Davis, CA95616
| | - Ziheng Yang
- Department of Genetics, Evolution, and Environment, University College London, LondonWC1E 6BT, United Kingdom
| |
Collapse
|
14
|
Sánchez KI, Diaz Huesa EG, Breitman MF, Avila LJ, Sites JW, Morando M. Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade). Syst Biol 2023; 72:739-752. [PMID: 37097104 DOI: 10.1093/sysbio/syad019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 03/28/2023] [Accepted: 04/11/2023] [Indexed: 04/26/2023] Open
Abstract
In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].
Collapse
Affiliation(s)
- Kevin I Sánchez
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales, Consejo Nacional de Investigaciones Científicas y Técnicas (IPEEC-CONICET), Puerto Madryn, U9120ACD, Argentina
| | - Emilce G Diaz Huesa
- Instituto de Diversidad y Evolución Austral, Consejo Nacional de Investigaciones Científicas y Técnicas (IDEAus-CONICET), Puerto Madryn, U9120ACD, Argentina
| | - María F Breitman
- Department of Biology and Environmental Science, Auburn University at Montgomery, Montgomery, 36117, USA
| | - Luciano J Avila
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales, Consejo Nacional de Investigaciones Científicas y Técnicas (IPEEC-CONICET), Puerto Madryn, U9120ACD, Argentina
| | - Jack W Sites
- Department of Biology, Austin Peay State University, Clarksville, 37044, USA
| | - Mariana Morando
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales, Consejo Nacional de Investigaciones Científicas y Técnicas (IPEEC-CONICET), Puerto Madryn, U9120ACD, Argentina
- Universidad Nacional de la Patagonia San Juan Bosco (UNPSJB), Puerto Madryn, U9120ACD, Argentina
| |
Collapse
|
15
|
Mcguire JA, Huang X, Reilly SB, Iskandar DT, Wang-Claypool CY, Werning S, Chong RA, Lawalata SZS, Stubbs AL, Frederick JH, Brown RM, Evans BJ, Arifin U, Riyanto A, Hamidy A, Arida E, Koo MS, Supriatna J, Andayani N, Hall R. Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation. Syst Biol 2023; 72:885-911. [PMID: 37074804 PMCID: PMC10405571 DOI: 10.1093/sysbio/syad020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2022] [Revised: 03/14/2023] [Accepted: 04/13/2023] [Indexed: 04/20/2023] Open
Abstract
The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].
Collapse
Affiliation(s)
- Jimmy A Mcguire
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
| | - Xiaoting Huang
- College of Marine Life Sciences, Ocean University of China, No. 5 Yushan Road, Qindao, Shandong, 266003, PR China
| | - Sean B Reilly
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, CA 95060, USA
| | - Djoko T Iskandar
- School of Life Sciences and Technology, Institut Teknologi Bandung, Bandung, Indonesia
| | - Cynthia Y Wang-Claypool
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
| | - Sarah Werning
- Department of Anatomy, Des Moines University, 3200 Grand Avenue, Des Moines, IA 50312-4198, USA
| | - Rebecca A Chong
- Department of Biology, University of Hawaii at Manoa, Honolulu, HI 96822, USA
| | - Shobi Z S Lawalata
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
- United in Diversity Foundation, Jalan Hayam Wuruk, Jakarta, Indonesia
| | - Alexander L Stubbs
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
| | - Jeffrey H Frederick
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
- Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
| | - Rafe M Brown
- Biodiversity Institute and Department of Ecology and Evolutionary Biology, 1345 Jayhawk Blvd., University of Kansas, Lawrence, KS 66045, USA
| | - Ben J Evans
- Biology Department, McMaster University, Hamilton, Ontario, Canada
| | - Umilaela Arifin
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
- School of Life Sciences and Technology, Institut Teknologi Bandung, Bandung, Indonesia
- Center for Taxonomy and Morphology, Zoologisches Museum Hamburg, Leibniz Institute for the Analysis of Biodiversity Change, Martin-Luther-King-Platz 3, R230 20146 Hamburg, Germany
| | - Awal Riyanto
- Laboratory of Herpetology, Museum Zoologicum Bogoriense, Research Center for Biosystematics and Evolution, National Research and Innovation Agency of Indonesia (BRIN), Cibinong 16911, Indonesia
| | - Amir Hamidy
- Laboratory of Herpetology, Museum Zoologicum Bogoriense, Research Center for Biosystematics and Evolution, National Research and Innovation Agency of Indonesia (BRIN), Cibinong 16911, Indonesia
| | - Evy Arida
- Research Center for Applied Zoology, National Research and Innovation Agency of Indonesia (BRIN), Cibinong 16911, Indonesia
| | - Michelle S Koo
- Museum of Vertebrate Zoology, University of California, Berkeley, CA 94720, USA
| | - Jatna Supriatna
- Department of Biology, Institute for Sustainable Earth and Resources (I-SER), Gedung Laboratorium Multidisiplin, and Research Center for Climate Change (RCCC-UI), Gedung Laboratorium Multidisiplin, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok 16424, Indonesia
| | - Noviar Andayani
- Department of Biology, Institute for Sustainable Earth and Resources (I-SER), Gedung Laboratorium Multidisiplin, and Research Center for Climate Change (RCCC-UI), Gedung Laboratorium Multidisiplin, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok 16424, Indonesia
| | - Robert Hall
- SE Asia Research Group (SEARG), Department of Earth Sciences, Royal Holloway University of London, Egham, Surrey TW20 0EX, UK
| |
Collapse
|
16
|
Tiley GP, Flouri T, Jiao X, Poelstra JW, Xu B, Zhu T, Rannala B, Yoder AD, Yang Z. Estimation of species divergence times in presence of cross-species gene flow. Syst Biol 2023; 72:820-836. [PMID: 36961245 PMCID: PMC10405360 DOI: 10.1093/sysbio/syad015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 03/22/2023] [Indexed: 03/25/2023] Open
Abstract
Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent].
Collapse
Affiliation(s)
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London, UK
| | - Xiyun Jiao
- Department of Genetics, Evolution and Environment, University College London, London, UK
- Department of Statistics and Data Science, China Southern University of Science and Technology, Shenzhen, Guangdong, China
| | | | - Bo Xu
- Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China
| | - Tianqi Zhu
- National Center for Mathematics and Interdisciplinary Sciences, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, China
- Key Laboratory of Random Complex Structures and Data Science, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, China
| | - Bruce Rannala
- Department of Evolution and Ecology, University of California, Davis, Davis, CA, USA
| | - Anne D Yoder
- Department of Biology, Duke University, Durham, NC, USA
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, UK
| |
Collapse
|
17
|
Thawornwattana Y, Huang J, Flouri T, Mallet J, Yang Z. Inferring the Direction of Introgression Using Genomic Sequence Data. Mol Biol Evol 2023; 40:msad178. [PMID: 37552932 PMCID: PMC10439365 DOI: 10.1093/molbev/msad178] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 08/01/2023] [Accepted: 08/02/2023] [Indexed: 08/10/2023] Open
Abstract
Genomic data are informative about the history of species divergence and interspecific gene flow, including the direction, timing, and strength of gene flow. However, gene flow in opposite directions generates similar patterns in multilocus sequence data, such as reduced sequence divergence between the hybridizing species. As a result, inference of the direction of gene flow is challenging. Here, we investigate the information about the direction of gene flow present in genomic sequence data using likelihood-based methods under the multispecies-coalescent-with-introgression model. We analyze the case of two species, and use simulation to examine cases with three or four species. We find that it is easier to infer gene flow from a small population to a large one than in the opposite direction, and easier to infer inflow (gene flow from outgroup species to an ingroup species) than outflow (gene flow from an ingroup species to an outgroup species). It is also easier to infer gene flow if there is a longer time of separate evolution between the initial divergence and subsequent introgression. When introgression is assumed to occur in the wrong direction, the time of introgression tends to be correctly estimated and the Bayesian test of gene flow is often significant, while estimates of introgression probability can be even greater than the true probability. We analyze genomic sequences from Heliconius butterflies to demonstrate that typical genomic datasets are informative about the direction of interspecific gene flow, as well as its timing and strength.
Collapse
Affiliation(s)
| | - Jun Huang
- School of Biomedical Engineering, Capital Medical University, Beijing 100069, P.R. China
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| |
Collapse
|
18
|
Ji J, Jackson DJ, Leaché AD, Yang Z. Power of Bayesian and Heuristic Tests to Detect Cross-Species Introgression with Reference to Gene Flow in the Tamias quadrivittatus Group of North American Chipmunks. Syst Biol 2023; 72:446-465. [PMID: 36504374 PMCID: PMC10275556 DOI: 10.1093/sysbio/syac077] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 11/15/2022] [Accepted: 12/01/2022] [Indexed: 10/25/2023] Open
Abstract
In the past two decades, genomic data have been widely used to detect historical gene flow between species in a variety of plants and animals. The Tamias quadrivittatus group of North America chipmunks, which originated through a series of rapid speciation events, are known to undergo massive amounts of mitochondrial introgression. Yet in a recent analysis of targeted nuclear loci from the group, no evidence for cross-species introgression was detected, indicating widespread cytonuclear discordance. The study used the heuristic method HYDE to detect gene flow, which may suffer from low power. Here we use the Bayesian method implemented in the program BPP to re-analyze these data. We develop a Bayesian test of introgression, calculating the Bayes factor via the Savage-Dickey density ratio using the Markov chain Monte Carlo (MCMC) sample under the model of introgression. We take a stepwise approach to constructing an introgression model by adding introgression events onto a well-supported binary species tree. The analysis detected robust evidence for multiple ancient introgression events affecting the nuclear genome, with introgression probabilities reaching 63%. We estimate population parameters and highlight the fact that species divergence times may be seriously underestimated if ancient cross-species gene flow is ignored in the analysis. We examine the assumptions and performance of HYDE and demonstrate that it lacks power if gene flow occurs between sister lineages or if the mode of gene flow does not match the assumed hybrid-speciation model with symmetrical population sizes. Our analyses highlight the power of likelihood-based inference of cross-species gene flow using genomic sequence data. [Bayesian test; BPP; chipmunks; introgression; MSci; multispecies coalescent; Savage-Dickey density ratio.].
Collapse
Affiliation(s)
- Jiayi Ji
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Donavan J Jackson
- Department of Biology and Burke Museum of Natural History and Culture, University of Washington, Box 351800, Seattle, WA 98195-1800, USA
| | - Adam D Leaché
- Department of Biology and Burke Museum of Natural History and Culture, University of Washington, Box 351800, Seattle, WA 98195-1800, USA
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| |
Collapse
|
19
|
Lopes F, Oliveira LR, Beux Y, Kessler A, Cárdenas-Alayza S, Majluf P, Páez-Rosas D, Chaves J, Crespo E, Brownell RL, Baylis AMM, Sepúlveda M, Franco-Trecu V, Loch C, Robertson BC, Peart CR, Wolf JBW, Bonatto SL. Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator. SCIENCE ADVANCES 2023; 9:eadf6601. [PMID: 37134171 PMCID: PMC10156116 DOI: 10.1126/sciadv.adf6601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.
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, PUCRS, Porto Alegre, Brazil
- Laboratório de Ecologia de Mamíferos, Universidade do Vale do Rio dos Sinos, São Leopoldo, Brazil
- Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland
| | - Larissa R Oliveira
- Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland
- Grupo de Estudos de Mamíferos Aquáticos do Rio Grande do Sul (GEMARS), Torres, Brazil
| | - Yago Beux
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, PUCRS, Porto Alegre, Brazil
| | - Amanda Kessler
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, PUCRS, Porto Alegre, Brazil
| | - Susana Cárdenas-Alayza
- Centro para la Sostenibilidad Ambiental, Universidad Peruana Cayetano Heredia, Lima, Peru
- Departamento de Ciencias Biológicas y Fisiológicas, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Patricia Majluf
- Centro para la Sostenibilidad Ambiental, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Diego Páez-Rosas
- Colegio de Ciencias Biológicas y Ambientales, COCIBA, Universidad San Francisco de Quito, Quito, Ecuador
- Dirección del Parque Nacional Galápagos, Oficina Técnica San Cristobal, Islas Galápagos, Ecuador
| | - Jaime Chaves
- Colegio de Ciencias Biológicas y Ambientales, COCIBA, Universidad San Francisco de Quito, Quito, Ecuador
- Galapagos Science Center, Puerto Baquerizo Moreno, Ecuador
- Department of Biology, San Francisco State University, 1800 Holloway Ave, San Francisco, CA, USA
| | - Enrique Crespo
- Laboratório de Mamíferos Marinos, CESIMAR - CCT CENPAT, CONICET, Puerto Madryn, Argentina
| | - Robert L Brownell
- Southwest Fisheries Science Center, NOAA Fisheries, La Jolla, CA, USA
| | | | - Maritza Sepúlveda
- Centro de Investigación y Gestión de Recursos Naturales (CIGREN), Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
| | - Valentina Franco-Trecu
- Departamento de Ecología y Evolución, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Carolina Loch
- Sir John Walsh Research Institute, Faculty of Dentistry, University of Otago, Dunedin, New Zealand
| | | | - Claire R Peart
- Division of Evolutionary Biology, LMU Munich, München, Germany
| | - Jochen B W Wolf
- Division of Evolutionary Biology, LMU Munich, München, Germany
| | - Sandro L Bonatto
- Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul, PUCRS, Porto Alegre, Brazil
| |
Collapse
|
20
|
Yi H, Dong S, Yang L, Wang J, Kidner C, Kang M. Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns. Mol Phylogenet Evol 2023; 184:107801. [PMID: 37088242 DOI: 10.1016/j.ympev.2023.107801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Revised: 04/07/2023] [Accepted: 04/18/2023] [Indexed: 04/25/2023]
Abstract
Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.
Collapse
Affiliation(s)
- Huiqin Yi
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China; South China National Botanical Garden, Guangzhou 510650, China
| | - Shiying Dong
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China; South China National Botanical Garden, Guangzhou 510650, China
| | - Lihua Yang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China; South China National Botanical Garden, Guangzhou 510650, China
| | - Jing Wang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China; South China National Botanical Garden, Guangzhou 510650, China
| | - Catherine Kidner
- Institute of Molecular Plant Sciences, University of Edinburgh, Daniel Rutherford Building Max Born Crescent, The King's Buildings, Edinburgh EH9 3BF, UK; Royal Botanic Garden Edinburgh, 20a Inverleith Row, Edinburgh EH3 5LR, UK
| | - Ming Kang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China; South China National Botanical Garden, Guangzhou 510650, China.
| |
Collapse
|
21
|
Ortiz D, Pekár S, Bryjová A. Gene flow assessment helps to distinguish strong genomic structure from speciation in an Iberian ant-eating spider. Mol Phylogenet Evol 2023; 180:107682. [PMID: 36574825 DOI: 10.1016/j.ympev.2022.107682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 12/01/2022] [Accepted: 12/21/2022] [Indexed: 12/26/2022]
Abstract
Although genomic data is boosting our understanding of evolution, we still lack a solid framework to perform reliable genome-based species delineation. This problem is especially critical in the case of phylogeographically structured organisms, with allopatric populations showing similar divergence patterns as species. Here, we assess the species limits and phylogeography of Zodarion alacre, an ant-eating spider widely distributed across the Iberian Peninsula. We first performed species delimitation based on genome-wide data and then validated these results using additional evidence. A commonly employed species delimitation strategy detected four distinct lineages with almost no admixture, which present allopatric distributions. These lineages showed ecological differentiation but no clear morphological differentiation, and evidence of introgression in a mitochondrial barcode. Phylogenomic networks found evidence of substantial gene flow between lineages. Finally, phylogeographic methods highlighted remarkable isolation by distance and detected evidence of range expansion from south-central Portugal to central-north Spain. We conclude that despite their deep genomic differentiation, the lineages of Z. alacre do not show evidence of complete speciation. Our results likely shed light on why Zodarion is among the most diversified spider genera despite its limited distribution and support the use of gene flow evidence to inform species boundaries.
Collapse
Affiliation(s)
- David Ortiz
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czechia.
| | - Stano Pekár
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Brno, Czechia
| | - Anna Bryjová
- Institute of Vertebrate Biology of the Czech Academy of Sciences, Brno, Czechia
| |
Collapse
|
22
|
Edwards SV, Tonini JFR, Mcinerney N, Welch C, Beerli P. Multilocus phylogeography, population genetics and niche evolution of Australian brown and black-tailed treecreepers (Aves: Climacteris). Biol J Linn Soc Lond 2023. [DOI: 10.1093/biolinnean/blac144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Abstract
The Carpentarian barrier across north-eastern Australia is a major biogeographic barrier and a generator of biodiversity within the Australian Monsoonal Tropics. Here we present a continent-wide analysis of mitochondrial (control region) and autosomal (14 anonymous loci) sequence and indel variation and niche modelling of brown and black-tailed treecreepers (Climacteris picumnus and Climacteris melanurus), a clade with a classic distribution on either side of the Carpentarian barrier. mtDNA control region sequences exhibited reciprocal monophyly and strong differentiation (Fst = 0.91), and revealed a signature of a recent selective sweep in C. picumnus. A variety of tests support an isolation-with-migration model of divergence, albeit with low levels of gene flow across the Carpentarian barrier and a divergence time between species of ~1.7–2.8 Mya. Palaeoecological niche models show that both range size as measured by available habitat and estimated historical population sizes of both species declined in the past ~600 kyr and that the area of interspecific range overlap was never historically large, perhaps decreasing opportunities for extensive gene flow. The relatively long divergence time and low opportunity for gene flow may have facilitated speciation more so than in other co-distributed bird taxa across the Australian Monsoonal Tropics.
Collapse
Affiliation(s)
- Scott V Edwards
- Museum of Comparative Zoology, Harvard University , Cambridge, MA 02138 , USA
- Department of Organismic and Evolutionary Biology, Harvard University , Cambridge, MA 02138 , USA
| | - João F R Tonini
- Museum of Comparative Zoology, Harvard University , Cambridge, MA 02138 , USA
- Department of Organismic and Evolutionary Biology, Harvard University , Cambridge, MA 02138 , USA
- Department of Biology, University of Richmond , Richmond, VA 23217 , USA
| | - Nancy Mcinerney
- Smithsonian's National Zoo and Conservation Biology Institute , NW, Washington, DC 20008 , USA
| | - Corey Welch
- Department of Biology and Burke Museum, University of Washington , Seattle, WA 98195 , USA
- STEM Scholars Program, Student Innovation Center, Iowa State University , Ames, IA 50011 , USA
| | - Peter Beerli
- Department of Scientific Computing, Florida State University, Florida State University , Tallahassee, FL 32306 , USA
| |
Collapse
|
23
|
Zozaya SM, Teasdale LC, Tedeschi LG, Higgie M, Hoskin CJ, Moritz C. Initiation of speciation across multiple dimensions in a rock-restricted, tropical lizard. Mol Ecol 2023; 32:680-695. [PMID: 36394360 PMCID: PMC10099344 DOI: 10.1111/mec.16787] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 11/02/2022] [Accepted: 11/08/2022] [Indexed: 11/18/2022]
Abstract
Population isolation and concomitant genetic divergence, resulting in strong phylogeographical structure, is a core aspect of speciation initiation. If and how speciation then proceeds and ultimately completes depends on multiple factors that mediate reproductive isolation, including divergence in genomes, ecology and mating traits. Here we explored these multiple dimensions in two young (Plio-Pleistocene) species complexes of gekkonid lizards (Heteronotia) from the Kimberley-Victoria River regions of tropical Australia. Using mitochondrial DNA screening and exon capture phylogenomics, we show that the rock-restricted Heteronotia planiceps exhibits exceptional fine-scale phylogeographical structure compared to the codistributed habitat generalist Heteronotia binoei. This indicates pervasive population isolation and persistence in the rock-specialist, and thus a high rate of speciation initiation across this geographically complex region, with levels of genomic divergence spanning the "grey zone" of speciation. Proximal lineages of H. planiceps were often separated by different rock substrates, suggesting a potential role for ecological isolation; however, phylogenetic incongruence and historical introgression were inferred between one such pair. Ecomorphological divergence among lineages within both H. planiceps and H. binoei was limited, except that limestone-restricted lineages of H. planiceps tended to be larger than rock-generalists. By contrast, among-lineage divergence in the chemical composition of epidermal pore secretions (putative mating trait) exceeded ecomorphology in both complexes, but with less trait overlap among lineages in H. planiceps. This system-particularly the rock-specialist H. planiceps-highlights the role of multidimensional divergence during incipient speciation, with divergence in genomes, ecomorphology and chemical signals all at play at very fine spatial scales.
Collapse
Affiliation(s)
- Stephen M Zozaya
- Research School of Biology, Australian National University, Australian Capital Territory, Canberra, Australia
| | - Luisa C Teasdale
- Research School of Biology, Australian National University, Australian Capital Territory, Canberra, Australia.,Department of Molecular Biology, Max Planck Institute for Developmental Biology, Tübingen, Germany
| | - Leonardo G Tedeschi
- Research School of Biology, Australian National University, Australian Capital Territory, Canberra, Australia
| | - Megan Higgie
- College of Science and Engineering, James Cook University, Queensland, Townsville, Australia
| | - Conrad J Hoskin
- College of Science and Engineering, James Cook University, Queensland, Townsville, Australia
| | - Craig Moritz
- Research School of Biology, Australian National University, Australian Capital Territory, Canberra, Australia
| |
Collapse
|
24
|
Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles. Animals (Basel) 2023; 13:ani13030471. [PMID: 36766360 PMCID: PMC9913427 DOI: 10.3390/ani13030471] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2022] [Revised: 01/13/2023] [Accepted: 01/15/2023] [Indexed: 02/01/2023] Open
Abstract
Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.
Collapse
|
25
|
Allman ES, Baños H, Mitchell JD, Rhodes JA. The tree of blobs of a species network: identifiability under the coalescent. J Math Biol 2022; 86:10. [PMID: 36472708 PMCID: PMC10062380 DOI: 10.1007/s00285-022-01838-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Revised: 08/31/2022] [Accepted: 11/17/2022] [Indexed: 12/12/2022]
Abstract
Inference of species networks from genomic data under the Network Multispecies Coalescent Model is currently severely limited by heavy computational demands. It also remains unclear how complicated networks can be for consistent inference to be possible. As a step toward inferring a general species network, this work considers its tree of blobs, in which non-cut edges are contracted to nodes, so only tree-like relationships between the taxa are shown. An identifiability theorem, that most features of the unrooted tree of blobs can be determined from the distribution of gene quartet topologies, is established. This depends upon an analysis of gene quartet concordance factors under the model, together with a new combinatorial inference rule. The arguments for this theoretical result suggest a practical algorithm for tree of blobs inference, to be fully developed in a subsequent work.
Collapse
Affiliation(s)
- Elizabeth S Allman
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK, 99775, USA
| | - Hector Baños
- Department of Biochemistry and Molecular Biology, Faculty of Medicine, Dalhousie University, Halifax, NS, Canada
- Department of Mathematics and Statistics, Faculty of Science, Dalhousie University, Halifax, NS, Canada
| | - Jonathan D Mitchell
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK, 99775, USA
- School of Natural Sciences (Mathematics), University of Tasmania, Hobart, TAS, 7001, Australia
- ARC Centre of Excellence for Plant Success in Nature and Agriculture, University of Tasmania, Hobart, TAS, 7001, Australia
| | - John A Rhodes
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK, 99775, USA.
| |
Collapse
|
26
|
Huang J, Thawornwattana Y, Flouri T, Mallet J, Yang Z. Inference of Gene Flow between Species under Misspecified Models. Mol Biol Evol 2022; 39:6783212. [PMID: 36317198 PMCID: PMC9729068 DOI: 10.1093/molbev/msac237] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
Genomic sequence data provide a rich source of information about the history of species divergence and interspecific hybridization or introgression. Despite recent advances in genomics and statistical methods, it remains challenging to infer gene flow, and as a result, one may have to estimate introgression rates and times under misspecified models. Here we use mathematical analysis and computer simulation to examine estimation bias and issues of interpretation when the model of gene flow is misspecified in analysis of genomic datasets, for example, if introgression is assigned to the wrong lineages. In the case of two species, we establish a correspondence between the migration rate in the continuous migration model and the introgression probability in the introgression model. When gene flow occurs continuously through time but in the analysis is assumed to occur at a fixed time point, common evolutionary parameters such as species divergence times are surprisingly well estimated. However, the time of introgression tends to be estimated towards the recent end of the period of continuous gene flow. When introgression events are assigned incorrectly to the parental or daughter lineages, introgression times tend to collapse onto species divergence times, with introgression probabilities underestimated. Overall, our analyses suggest that the simple introgression model is useful for extracting information concerning between-specific gene flow and divergence even when the model may be misspecified. However, for reliable inference of gene flow it is important to include multiple samples per species, in particular, from hybridizing species.
Collapse
Affiliation(s)
| | | | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, United Kingdom
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138
| | | |
Collapse
|
27
|
Borges R, Boussau B, Höhna S, Pereira RJ, Kosiol C. Polymorphism‐aware estimation of species trees and evolutionary forces from genomic sequences with
RevBayes. Methods Ecol Evol 2022. [DOI: 10.1111/2041-210x.13980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Rui Borges
- Institut für Populationsgenetik, Vetmeduni Vienna Wien Austria
| | - Bastien Boussau
- Université de Lyon, Université Claude Bernard Lyon 1 Villeurbanne France
| | - Sebastian Höhna
- GeoBio‐Center, Ludwig‐Maximilians‐Universität München Munich Germany
- Department of Earth and Environmental Sciences, Paleontology & Geobiology Ludwig‐Maximilians‐Universität München Munich Germany
| | - Ricardo J. Pereira
- Division of Evolutionary Biology, Department of Biology II Ludwig‐Maximilians‐Universität München Martinsried Germany
| | - Carolin Kosiol
- Institut für Populationsgenetik, Vetmeduni Vienna Wien Austria
- Centre for Biological Diversity University of St Andrews St Andrews UK
| |
Collapse
|
28
|
Legried B, Terhorst J. Rates of convergence in the two-island and isolation-with-migration models. Theor Popul Biol 2022; 147:16-27. [PMID: 36007782 DOI: 10.1016/j.tpb.2022.08.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2021] [Revised: 08/10/2022] [Accepted: 08/11/2022] [Indexed: 11/25/2022]
Abstract
A number of powerful demographic inference methods have been developed in recent years, with the goal of fitting rich evolutionary models to genetic data obtained from many populations. In this paper we investigate the statistical performance of these methods in the specific case where there is continuous migration between populations. Compared with earlier work, migration significantly complicates the theoretical analysis and requires new techniques. We employ the theories of phase-type distributions and concentration of measure in order to study the two-island and isolation-with-migration models, resulting in both upper and lower bounds on rates of convergence for parametric estimators in migration models. For the upper bounds, we consider inferring rates of coalescent and migration on the basis of directly observing pairwise coalescent times, and, more realistically, when (conditionally) Poisson-distributed mutations dropped on latent trees are observed. We complement these upper bounds with information-theoretic lower bounds which establish a limit, in terms of sample size, below which inference is effectively impossible.
Collapse
Affiliation(s)
- Brandon Legried
- Department of Statistics, University of Michigan, United States of America
| | - Jonathan Terhorst
- Department of Statistics, University of Michigan, United States of America.
| |
Collapse
|
29
|
Barley AJ, Nieto-Montes de Oca A, Manríquez-Morán NL, Thomson RC. The evolutionary network of whiptail lizards reveals predictable outcomes of hybridization. Science 2022; 377:773-777. [PMID: 35951680 DOI: 10.1126/science.abn1593] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Hybridization between diverging lineages is associated with the generation and loss of species diversity, introgression, adaptation, and changes in reproductive mode, but it is unknown when and why it results in these divergent outcomes. We estimate a comprehensive evolutionary network for the largest group of unisexual vertebrates and use it to understand the evolutionary outcomes of hybridization. Our results show that rates of introgression between species decrease with time since divergence and suggest that species must attain a threshold of evolutionary divergence before hybridization results in transitions to unisexuality. Rates of hybridization also predict genome-wide patterns of genetic diversity in whiptail lizards. These results distinguish among models for hybridization that have not previously been tested and suggest that the evolutionary outcomes can be predictable.
Collapse
Affiliation(s)
- Anthony J Barley
- Department of Evolution and Ecology, University of California, Davis, Davis, CA 95616, USA.,School of Life Sciences, University of Hawai'i, Honolulu, HI 96822, USA
| | - Adrián Nieto-Montes de Oca
- Laboratorio de Herpetología and Museo de Zoología Alfonso L. Herrera, Departamento de Biología Evolutiva, Facultad de Ciencias, Universidad Nacional Autónoma de México, Ciudad Universitaria, Alcadía Coyoacán, Ciudad de México, México
| | - Norma L Manríquez-Morán
- Laboratorio de Sistemática Molecular, Centro de Investigaciones Biológicas, Universidad Autónoma del Estado de Hidalgo, Colonia Carboneras, Mineral de la Reforma, Hidalgo, México
| | - Robert C Thomson
- School of Life Sciences, University of Hawai'i, Honolulu, HI 96822, USA
| |
Collapse
|
30
|
Poelstra JW, Montero BK, Lüdemann J, Yang Z, Rakotondranary SJ, Hohenlohe P, Stetter N, Ganzhorn JU, Yoder AD. RADseq data reveal a lack of admixture in a mouse lemur contact zone contrary to previous microsatellite results. Proc Biol Sci 2022; 289:20220596. [PMID: 35946151 PMCID: PMC9364002 DOI: 10.1098/rspb.2022.0596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Microsatellites have been a workhorse of evolutionary genetic studies for decades and are still commonly in use for estimating signatures of genetic diversity at the population and species level across a multitude of taxa. Yet, the very high mutation rate of these loci is a double-edged sword, conferring great sensitivity at shallow levels of analysis (e.g. paternity analysis) but yielding considerable uncertainty for deeper evolutionary comparisons. For the present study, we used reduced representation genome-wide data (restriction site-associated DNA sequencing (RADseq)) to test for patterns of interspecific hybridization previously characterized using microsatellite data in a contact zone between two closely related mouse lemur species in Madagascar (Microcebus murinus and Microcebus griseorufus). We revisit this system by examining populations in, near, and far from the contact zone, including many of the same individuals that had previously been identified as hybrids with microsatellite data. Surprisingly, we find no evidence for admixed nuclear ancestry. Instead, re-analyses of microsatellite data and simulations suggest that previously inferred hybrids were false positives and that the program NewHybrids can be particularly sensitive to erroneously inferring hybrid ancestry. Combined with results from coalescent-based analyses and evidence for local syntopic co-occurrence, we conclude that the two mouse lemur species are in fact completely reproductively isolated, thus providing a new understanding of the evolutionary rate whereby reproductive isolation can be achieved in a primate.
Collapse
Affiliation(s)
- Jelmer W. Poelstra
- Department of Biology, Duke University, Durham, NC 27708, USA,Molecular and Cellular Imaging Center, Ohio State University, Wooster, OH 44691, USA
| | - B. Karina Montero
- Institute of Zoology, Department of Animal Ecology and Conservation, Universität Hamburg, Hamburg, 20146, Germany
| | - Jan Lüdemann
- Institute of Zoology, Department of Animal Ecology and Conservation, Universität Hamburg, Hamburg, 20146, Germany
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, UK
| | - S. Jacques Rakotondranary
- Institute of Zoology, Department of Animal Ecology and Conservation, Universität Hamburg, Hamburg, 20146, Germany,Anthropobiologie et Développement Durable, Faculté des Sciences, Université d'Antananarivo, PO Box 906, Antananarivo 101, Madagascar
| | - Paul Hohenlohe
- Institute for Bioinformatics and Evolutionary Studies, Department of Biological Sciences, University of Idaho, Moscow, ID 83844, USA
| | - Nadine Stetter
- Institute of Zoology, Department of Animal Ecology and Conservation, Universität Hamburg, Hamburg, 20146, Germany,Bernhard Nocht Institute for Tropical Medicine, 20359 Hamburg, Germany
| | - Jörg U. Ganzhorn
- Institute of Zoology, Department of Animal Ecology and Conservation, Universität Hamburg, Hamburg, 20146, Germany
| | - Anne D. Yoder
- Department of Biology, Duke University, Durham, NC 27708, USA
| |
Collapse
|
31
|
Flouri T, Huang J, Jiao X, Kapli P, Rannala B, Yang Z. Bayesian phylogenetic inference using relaxed-clocks and the multispecies coalescent. Mol Biol Evol 2022; 39:6652437. [PMID: 35907248 PMCID: PMC9366188 DOI: 10.1093/molbev/msac161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes–Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed.
Collapse
Affiliation(s)
- Tomáš Flouri
- Department of Genetics, Evolution, and Environment, University College London, Gower Street, London WC1E 6BT, UK
| | - Jun Huang
- Department of Genetics, Evolution, and Environment, University College London, Gower Street, London WC1E 6BT, UK.,School of Biomedical Engineering, Capital Medical University, Beijing, 100069, China
| | - Xiyun Jiao
- Department of Genetics, Evolution, and Environment, University College London, Gower Street, London WC1E 6BT, UK.,Department of Statistics and Data Science, China Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Paschalia Kapli
- Department of Genetics, Evolution, and Environment, 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, Evolution, and Environment, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
32
|
Pang XX, Zhang DY. Impact of Ghost Introgression on Coalescent-based Species Tree Inference and Estimation of Divergence Time. Syst Biol 2022; 72:35-49. [PMID: 35799362 DOI: 10.1093/sysbio/syac047] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 06/25/2022] [Accepted: 07/05/2022] [Indexed: 11/15/2022] Open
Abstract
The species studied in any evolutionary investigation generally constitute a small proportion of all the species currently existing or that have gone extinct. It is therefore likely that introgression, which is widespread across the tree of life, involves "ghosts," i.e., unsampled, unknown, or extinct lineages. However, the impact of ghost introgression on estimations of species trees has rarely been studied and is poorly understood. Here, we use mathematical analysis and simulations to examine the robustness of species tree methods based on the multispecies coalescent model to introgression from a ghost or extant lineage. We found that many results originally obtained for introgression between extant species can easily be extended to ghost introgression, such as the strongly interactive effects of incomplete lineage sorting (ILS) and introgression on the occurrence of anomalous gene trees (AGTs). The relative performance of the summary species tree method (ASTRAL) and the full-likelihood method (*BEAST) varies under different introgression scenarios, with the former being more robust to gene flow between non-sister species whereas the latter performing better under certain conditions of ghost introgression. When an outgroup ghost (defined as a lineage that diverged before the most basal species under investigation) acts as the donor of the introgressed genes, the time of root divergence among the investigated species generally was overestimated, whereas ingroup introgression, as commonly perceived, can only lead to underestimation. In many cases of ingroup introgression that may or may not involve ghost lineages, the stronger the ILS, the higher the accuracy achieved in estimating the time of root divergence, although the topology of the species tree is more prone to be biased by the effect of introgression.
Collapse
Affiliation(s)
- Xiao-Xu Pang
- State Key Laboratory of Earth Surface Processes and Resource Ecology and Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
| | - Da-Yong Zhang
- State Key Laboratory of Earth Surface Processes and Resource Ecology and Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
33
|
Cuesta-Porta V, Barrera-Ruiz UM, Cibrián-Tovar D, Pujade-Villar J. Discovery of the male of Loxaulus hyalinus, and implications for the sympatric species Loxaulus laeta (Hymenoptera: Cynipidae: Cynipini): an integrative taxonomical case of species delimitation. J NAT HIST 2022. [DOI: 10.1080/00222933.2022.2077150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Víctor Cuesta-Porta
- Facultat de Biologia, Departament de Biologia Evolutiva, Universitat de Barcelona, Barcelona, Catalonia, Spain
| | | | - David Cibrián-Tovar
- División de Ciencias Forestales, Universidad Autónoma Chapingo, Chapingo, Mexico
| | - Juli Pujade-Villar
- Facultat de Biologia, Departament de Biologia Evolutiva, Universitat de Barcelona, Barcelona, Catalonia, Spain
| |
Collapse
|
34
|
Kong S, Pons JC, Kubatko L, Wicke K. Classes of explicit phylogenetic networks and their biological and mathematical significance. J Math Biol 2022; 84:47. [PMID: 35503141 DOI: 10.1007/s00285-022-01746-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Revised: 01/18/2022] [Accepted: 03/31/2022] [Indexed: 11/24/2022]
Abstract
The evolutionary relationships among organisms have traditionally been represented using rooted phylogenetic trees. However, due to reticulate processes such as hybridization or lateral gene transfer, evolution cannot always be adequately represented by a phylogenetic tree, and rooted phylogenetic networks that describe such complex processes have been introduced as a generalization of rooted phylogenetic trees. In fact, estimating rooted phylogenetic networks from genomic sequence data and analyzing their structural properties is one of the most important tasks in contemporary phylogenetics. Over the last two decades, several subclasses of rooted phylogenetic networks (characterized by certain structural constraints) have been introduced in the literature, either to model specific biological phenomena or to enable tractable mathematical and computational analyses. In the present manuscript, we provide a thorough review of these network classes, as well as provide a biological interpretation of the structural constraints underlying these networks where possible. In addition, we discuss how imposing structural constraints on the network topology can be used to address the scalability and identifiability challenges faced in the estimation of phylogenetic networks from empirical data.
Collapse
Affiliation(s)
- Sungsik Kong
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, USA
| | - Joan Carles Pons
- Department of Mathematics and Computer Science, University of the Balearic Islands, Palma, 07122, Spain
| | - Laura Kubatko
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, USA.,Department of Statistics, The Ohio State University, Columbus, OH, USA
| | - Kristina Wicke
- Department of Mathematics, The Ohio State University, Columbus, OH, USA.
| |
Collapse
|
35
|
Yang Z, Flouri T. Estimation of Cross-Species Introgression Rates Using Genomic Data Despite Model Unidentifiability. Mol Biol Evol 2022; 39:msac083. [PMID: 35417543 PMCID: PMC9087891 DOI: 10.1093/molbev/msac083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Full-likelihood implementations of the multispecies coalescent with introgression (MSci) model treat genealogical fluctuations across the genome as a major source of information to infer the history of species divergence and gene flow using multilocus sequence data. However, MSci models are known to have unidentifiability issues, whereby different models or parameters make the same predictions about the data and cannot be distinguished by the data. Previous studies of unidentifiability have focused on heuristic methods based on gene trees and do not make an efficient use of the information in the data. Here we study the unidentifiability of MSci models under the full-likelihood methods. We characterize the unidentifiability of the bidirectional introgression (BDI) model, which assumes that gene flow occurs in both directions. We derive simple rules for arbitrary BDI models, which create unidentifiability of the label-switching type. In general, an MSci model with k BDI events has 2k unidentifiable modes or towers in the posterior, with each BDI event between sister species creating within-model parameter unidentifiability and each BDI event between nonsister species creating between-model unidentifiability. We develop novel algorithms for processing Markov chain Monte Carlo samples to remove label-switching problems and implement them in the bpp program. We analyze real and synthetic data to illustrate the utility of the BDI models and the new algorithms. We discuss the unidentifiability of heuristic methods and provide guidelines for the use of MSci models to infer gene flow using genomic data.
Collapse
Affiliation(s)
- Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E 6BT, UK
| | | |
Collapse
|
36
|
Hurt C, Hildreth P, Williams C. A genomic perspective on the conservation status of the endangered Nashville crayfish (Faxonius shoupi). CONSERV GENET 2022. [DOI: 10.1007/s10592-022-01438-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
37
|
Zhu T, Flouri T, Yang Z. A simulation study to examine the impact of recombination on phylogenomic inferences under the multispecies coalescent model. Mol Ecol 2022; 31:2814-2829. [PMID: 35313033 PMCID: PMC9321900 DOI: 10.1111/mec.16433] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 01/25/2022] [Accepted: 02/28/2022] [Indexed: 11/28/2022]
Affiliation(s)
- Tianqi Zhu
- Institute of Applied Mathematics Academy of Mathematics and Systems Science Chinese Academy of Sciences Beijing 100190 China
- Key Laboratory of Random Complex Structures and Data Science, Academy of Mathematics and Systems Science, Chinese Academy of Sciences Beijing 100190 China
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment University College London London WC1E 6BT UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment University College London London WC1E 6BT UK
| |
Collapse
|
38
|
Sanderson MJ, Búrquez A, Copetti D, McMahon MM, Zeng Y, Wojciechowski MF. Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis. Syst Biol 2022; 71:1178-1194. [PMID: 35244183 DOI: 10.1093/sysbio/syac017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2020] [Revised: 02/18/2022] [Accepted: 02/25/2022] [Indexed: 11/14/2022] Open
Abstract
Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.
Collapse
Affiliation(s)
- Michael J Sanderson
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | - Alberto Búrquez
- Instituto de Ecología, Unidad Hermosillo, Universidad Nacional Autónoma de México, Hermosillo, Sonora, Mexico
| | - Dario Copetti
- Arizona Genomics Institute, School of Plant Sciences, University of Arizona, Tucson, AZ, 85721 USA
| | | | - Yichao Zeng
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | | |
Collapse
|
39
|
Panzeri D, Guidi Nissim W, Labra M, Grassi F. Revisiting the Domestication Process of African Vigna Species (Fabaceae): Background, Perspectives and Challenges. PLANTS (BASEL, SWITZERLAND) 2022; 11:plants11040532. [PMID: 35214865 PMCID: PMC8879845 DOI: 10.3390/plants11040532] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/25/2022] [Revised: 02/12/2022] [Accepted: 02/13/2022] [Indexed: 05/14/2023]
Abstract
Legumes are one of the most economically important and biodiverse families in plants recognised as the basis to develop functional foods. Among these, the Vigna genus stands out as a good representative because of its relatively recent African origin as well as its outstanding potential. Africa is a great biodiversity centre in which a great number of species are spread, but only three of them, Vigna unguiculata, Vigna subterranea and Vigna vexillata, were successfully domesticated. This review aims at analysing and valorising these species by considering the perspective of human activity and what effects it exerts. For each species, we revised the origin history and gave a focus on where, when and how many times domestication occurred. We provided a brief summary of bioactive compounds naturally occurring in these species that are fundamental for human wellbeing. The great number of wild lineages is a key point to improve landraces since the domestication process caused a loss of gene diversity. Their genomes hide a precious gene pool yet mostly unexplored, and genes lost during human activity can be recovered from the wild lineages and reintroduced in cultivated forms through modern technologies. Finally, we describe how all this information is game-changing to the design of future crops by domesticating de novo.
Collapse
Affiliation(s)
- Davide Panzeri
- Department of Biotechnology and Bioscience, University of Milan-Bicocca, Piazza della Scienza 2, 20126 Milano, Italy; (W.G.N.); (M.L.)
- Correspondence: (D.P.); (F.G.)
| | - Werther Guidi Nissim
- Department of Biotechnology and Bioscience, University of Milan-Bicocca, Piazza della Scienza 2, 20126 Milano, Italy; (W.G.N.); (M.L.)
- Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, Viale delle Idee 30, 50019 Sesto Fiorentino, Italy
| | - Massimo Labra
- Department of Biotechnology and Bioscience, University of Milan-Bicocca, Piazza della Scienza 2, 20126 Milano, Italy; (W.G.N.); (M.L.)
| | - Fabrizio Grassi
- Department of Biotechnology and Bioscience, University of Milan-Bicocca, Piazza della Scienza 2, 20126 Milano, Italy; (W.G.N.); (M.L.)
- Correspondence: (D.P.); (F.G.)
| |
Collapse
|
40
|
Hibbins MS, Hahn MW. Phylogenomic approaches to detecting and characterizing introgression. Genetics 2022; 220:iyab173. [PMID: 34788444 PMCID: PMC9208645 DOI: 10.1093/genetics/iyab173] [Citation(s) in RCA: 51] [Impact Index Per Article: 25.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Accepted: 10/02/2021] [Indexed: 12/26/2022] Open
Abstract
Phylogenomics has revealed the remarkable frequency with which introgression occurs across the tree of life. These discoveries have been enabled by the rapid growth of methods designed to detect and characterize introgression from whole-genome sequencing data. A large class of phylogenomic methods makes use of data across species to infer and characterize introgression based on expectations from the multispecies coalescent. These methods range from simple tests, such as the D-statistic, to model-based approaches for inferring phylogenetic networks. Here, we provide a detailed overview of the various signals that different modes of introgression are expected leave in the genome, and how current methods are designed to detect them. We discuss the strengths and pitfalls of these approaches and identify areas for future development, highlighting the different signals of introgression, and the power of each method to detect them. We conclude with a discussion of current challenges in inferring introgression and how they could potentially be addressed.
Collapse
Affiliation(s)
- Mark S Hibbins
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
- Department of Computer Science, Indiana University, Bloomington, IN 47405, USA
| |
Collapse
|
41
|
Suvorov A, Kim BY, Wang J, Armstrong EE, Peede D, D'Agostino ERR, Price DK, Waddell P, Lang M, Courtier-Orgogozo V, David JR, Petrov D, Matute DR, Schrider DR, Comeault AA. Widespread introgression across a phylogeny of 155 Drosophila genomes. Curr Biol 2022; 32:111-123.e5. [PMID: 34788634 PMCID: PMC8752469 DOI: 10.1016/j.cub.2021.10.052] [Citation(s) in RCA: 91] [Impact Index Per Article: 45.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Revised: 09/29/2021] [Accepted: 10/22/2021] [Indexed: 01/12/2023]
Abstract
Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.
Collapse
Affiliation(s)
- Anton Suvorov
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA.
| | - Bernard Y Kim
- Department of Biology, Stanford University, Stanford, CA, USA
| | - Jeremy Wang
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | | | - David Peede
- Department of Biology, University of North Carolina, Chapel Hill, NC 27599, USA
| | | | - Donald K Price
- School of Life Sciences, University of Nevada, Las Vegas, NV 89119, USA
| | - Peter Waddell
- School of Fundamental Sciences, Massey University, Palmerston North 4442, New Zealand
| | - Michael Lang
- CNRS, Institut Jacques Monod, Université de Paris, Paris 75013, France
| | | | - Jean R David
- Laboratoire Evolution, Génomes, Comportement, Ecologie (EGCE) CNRS, IRD, Univ. Paris-sud, Université Paris-Saclay, Gif sur Yvette 91190, France; Institut de Systématique, Evolution, Biodiversité, CNRS, MNHN, UPMC, EPHE, Muséum National d'Histoire Naturelle, Sorbonne Universités, Paris 75005, France
| | - Dmitri Petrov
- Department of Biology, Stanford University, Stanford, CA, USA
| | - Daniel R Matute
- Department of Biology, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Daniel R Schrider
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Aaron A Comeault
- Molecular Ecology & Evolution Group, School of Natural Sciences, Bangor University, Bangor, Gwynedd LL57 2DGA, UK.
| |
Collapse
|
42
|
Jiao X, Flouri T, Yang Z. Multispecies coalescent and its applications to infer species phylogenies and cross-species gene flow. Natl Sci Rev 2022; 8:nwab127. [PMID: 34987842 PMCID: PMC8692950 DOI: 10.1093/nsr/nwab127] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Revised: 07/10/2021] [Accepted: 07/11/2021] [Indexed: 02/06/2023] Open
Abstract
Multispecies coalescent (MSC) is the extension of the single-population coalescent model to multiple species. It integrates the phylogenetic process of species divergences and the population genetic process of coalescent, and provides a powerful framework for a number of inference problems using genomic sequence data from multiple species, including estimation of species divergence times and population sizes, estimation of species trees accommodating discordant gene trees, inference of cross-species gene flow and species delimitation. In this review, we introduce the major features of the MSC model, discuss full-likelihood and heuristic methods of species tree estimation and summarize recent methodological advances in inference of cross-species gene flow. We discuss the statistical and computational challenges in the field and research directions where breakthroughs may be likely in the next few years.
Collapse
Affiliation(s)
- Xiyun Jiao
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| |
Collapse
|
43
|
Borges R, Boussau B, Szöllősi GJ, Kosiol C. Nucleotide Usage Biases Distort Inferences of the Species Tree. Genome Biol Evol 2022; 14:6496956. [PMID: 34983052 PMCID: PMC8829901 DOI: 10.1093/gbe/evab290] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/27/2021] [Indexed: 12/15/2022] Open
Abstract
Despite the importance of natural selection in species’ evolutionary history, phylogenetic methods that take into account population-level processes typically ignore selection. The assumption of neutrality is often based on the idea that selection occurs at a minority of loci in the genome and is unlikely to compromise phylogenetic inferences significantly. However, genome-wide processes like GC-bias and some variation segregating at the coding regions are known to evolve in the nearly neutral range. As we are now using genome-wide data to estimate species trees, it is natural to ask whether weak but pervasive selection is likely to blur species tree inferences. We developed a polymorphism-aware phylogenetic model tailored for measuring signatures of nucleotide usage biases to test the impact of selection in the species tree. Our analyses indicate that although the inferred relationships among species are not significantly compromised, the genetic distances are systematically underestimated in a node-height-dependent manner: that is, the deeper nodes tend to be more underestimated than the shallow ones. Such biases have implications for molecular dating. We dated the evolutionary history of 30 worldwide fruit fly populations, and we found signatures of GC-bias considerably affecting the estimated divergence times (up to 23%) in the neutral model. Our findings call for the need to account for selection when quantifying divergence or dating species evolution.
Collapse
Affiliation(s)
- Rui Borges
- Institut für Populationsgenetik, Vetmeduni Vienna, Wien, Austria
| | - Bastien Boussau
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS UMR 5558, LBBE, Villeurbanne, France
| | - Gergely J Szöllősi
- Department of Biological Physics, Eötvös University, Budapest , Hungary.,MTA-ELTE "Lendület" Evolutionary Genomics Research Group, Budapest, Hungary.,Evolutionary Systems Research Group, Centre for Ecological Research, Hungarian Academy of Sciences, Tihany, Hungary
| | - Carolin Kosiol
- Institut für Populationsgenetik, Vetmeduni Vienna, Wien, Austria.,Centre for Biological Diversity, University of St Andrews, St Andrews, United Kingdom
| |
Collapse
|
44
|
Abstract
Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]
Collapse
|
45
|
Thawornwattana Y, Seixas FA, Yang Z, Mallet J. OUP accepted manuscript. Syst Biol 2022; 71:1159-1177. [PMID: 35169847 PMCID: PMC9366460 DOI: 10.1093/sysbio/syac009] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Revised: 02/01/2022] [Accepted: 02/08/2022] [Indexed: 11/21/2022] Open
Abstract
Introgressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species make estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ the multispecies coalescent (MSC) model with and without gene flow to infer the species phylogeny and cross-species introgression events using genomic data from six members of the erato-sara clade of Heliconius butterflies. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of Heliconius telesiphe, a likely hybrid speciation origin for Heliconius hecalesia, and gene flow between the sister species Heliconius erato and Heliconius himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the MSC makes it possible to extract rich historical information of species divergence and gene flow from genomic data. [3s; bpp; gene flow; Heliconius; hybrid speciation; introgression; inversion; multispecies coalescent]
Collapse
Affiliation(s)
- Yuttapong Thawornwattana
- Correspondence to be sent to: Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA; E-mail: ; (Y.T. and J.M.); Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK; E-mail: (Z.Y.)
| | - Fernando A Seixas
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | - Ziheng Yang
- Correspondence to be sent to: Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA; E-mail: ; (Y.T. and J.M.); Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK; E-mail: (Z.Y.)
| | - James Mallet
- Correspondence to be sent to: Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA; E-mail: ; (Y.T. and J.M.); Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK; E-mail: (Z.Y.)
| |
Collapse
|
46
|
Finger N, Farleigh K, Bracken JT, Leaché AD, François O, Yang Z, Flouri T, Charran T, Jezkova T, Williams DA, Blair C. Genome-scale data reveal deep lineage divergence and a complex demographic history in the Texas horned lizard (Phrynosoma cornutum) throughout the southwestern and central US. Genome Biol Evol 2021; 14:6443127. [PMID: 34849831 PMCID: PMC8735750 DOI: 10.1093/gbe/evab260] [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] [Accepted: 11/12/2021] [Indexed: 12/03/2022] Open
Abstract
The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find that divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.
Collapse
Affiliation(s)
- Nicholas Finger
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY, 11201, USA
| | - Keaka Farleigh
- Department of Biology, Miami University, 501 E High St, Oxford, OH, 45056, USA
| | - Jason T Bracken
- Department of Biology, Miami University, 501 E High St, Oxford, OH, 45056, USA
| | - Adam D Leaché
- Department of Biology & Burke Museum of Natural History and Culture, University of Washington, Seattle, WA, 98195, USA
| | - Olivier François
- Faculty of Medicine, University Grenoble-Alpes, TIMC-IMAG UMR 5525, Grenoble, La Tronche, F38706, France 38000
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, Darwin Building, Gower Street, London, WC1E 6BT, UK
| | - Tomas Flouri
- Department of Genetics, Evolution and Environment, University College London, Darwin Building, Gower Street, London, WC1E 6BT, UK
| | - Tristan Charran
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY, 11201, USA
| | - Tereza Jezkova
- Department of Biology, Miami University, 501 E High St, Oxford, OH, 45056, USA
| | - Dean A Williams
- Department of Biology, Texas Christian University, 2800 S University Dr, Fort Worth, TX, 76129, USA
| | - 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
| |
Collapse
|
47
|
Abstract
Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.
Collapse
Affiliation(s)
- Nathaniel B Edelman
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138, USA; .,Current affiliation: Yale Institute for Biospheric Studies and Yale School of the Environment, Yale University, New Haven, Connecticut 06511, USA;
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138, USA;
| |
Collapse
|
48
|
Calderón-Acevedo CA, Bagley JC, Muchhala N. Genome-wide ultraconserved elements resolve phylogenetic relationships and biogeographic history among Neotropical leaf-nosed bats in the genus Anoura (Phyllostomidae). Mol Phylogenet Evol 2021; 167:107356. [PMID: 34774763 DOI: 10.1016/j.ympev.2021.107356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Revised: 10/26/2021] [Accepted: 11/08/2021] [Indexed: 10/19/2022]
Abstract
AnouraGray, 1838 are Neotropical nectarivorous bats and the most speciose genus within the phyllostomid subfamily Glossophaginae. However, Anoura species limits remain debated, and phylogenetic relationships remain poorly known, because previous studies used limited Anoura taxon sampling or focused primarily on higher-level relationships. Here, we conduct the first phylogenomic study of Anoura by analyzing 2039 genome-wide ultraconserved elements (UCEs) sequenced for 42 individuals from 8 Anoura species/lineages plus two outgroups. Overall, our results based on UCEs resolved relationships in the genus and supported (1) the monophyly of small-bodied Anoura species (previously genus Lonchoglossa); (2) monotypic status of A. caudifer; and (3) nested positions of "A. carishina", A. caudifer aequatoris, and A. geoffroyi peruana specimens within A. latidens, A. caudifer and A. geoffroyi, respectively (suggesting that these taxa are not distinct species). Additionally, (4) phylogenetic networks allowing reticulate edges did not explain gene tree discordance better than the species tree (without introgression), indicating that a coalescent model accounting for discordance solely through incomplete lineage sorting fit our data well. Sensitivity analyses indicated that our species tree results were not adversely affected by varying taxon sampling across loci. Tree calibration and Bayesian coalescent analyses dated the onset of diversification within Anoura to around ∼ 6-9 million years ago in the Miocene, with extant species diverging mainly within the past ∼ 4 million years. We inferred a historical biogeographical scenario for Anoura of parapatric speciation fragmenting the range of a wide-ranging ancestral lineage centered in the Central to Northern Andes, along with Pliocene-Pleistocene dispersal or founder event speciation in Amazonia and the Brazilian Atlantic forest during the last ∼ 2.5 million years.
Collapse
Affiliation(s)
- Camilo A Calderón-Acevedo
- Department of Biology, University of Missouri-St. Louis, One University Blvd., 223 Research Bldg., St. Louis, MO 63121, USA; Department of Earth and Environmental Science, Rutgers University, 195 University Ave., Boyden Hall 433, Newark, NJ, 07102 USA.
| | - Justin C Bagley
- Department of Biology, University of Missouri-St. Louis, One University Blvd., 223 Research Bldg., St. Louis, MO 63121, USA; Department of Biology, Jacksonville State University, 242 Martin Hall, 700 Pelham Rd North, Jacksonville, AL 36265, USA; Department of Biology, Virginia Commonwealth University, 1000 W Cary St., Suite 126, Richmond, VA 23284, USA.
| | - Nathan Muchhala
- Department of Biology, University of Missouri-St. Louis, One University Blvd., 223 Research Bldg., St. Louis, MO 63121, USA.
| |
Collapse
|
49
|
Mirarab S, Nakhleh L, Warnow T. Multispecies Coalescent: Theory and Applications in Phylogenetics. ANNUAL REVIEW OF ECOLOGY, EVOLUTION, AND SYSTEMATICS 2021. [DOI: 10.1146/annurev-ecolsys-012121-095340] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Species tree estimation is a basic part of many biological research projects, ranging from answering basic evolutionary questions (e.g., how did a group of species adapt to their environments?) to addressing questions in functional biology. Yet, species tree estimation is very challenging, due to processes such as incomplete lineage sorting, gene duplication and loss, horizontal gene transfer, and hybridization, which can make gene trees differ from each other and from the overall evolutionary history of the species. Over the last 10–20 years, there has been tremendous growth in methods and mathematical theory for estimating species trees and phylogenetic networks, and some of these methods are now in wide use. In this survey, we provide an overview of the current state of the art, identify the limitations of existing methods and theory, and propose additional research problems and directions.
Collapse
Affiliation(s)
- Siavash Mirarab
- Electrical and Computer Engineering Department, University of California, San Diego, La Jolla, California 92093, USA
| | - Luay Nakhleh
- Department of Computer Science, Rice University, Houston, Texas 77005, USA
| | - Tandy Warnow
- Department of Computer Science, University of Illinois Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
50
|
Tihelka E, Cai C, Giacomelli M, Lozano-Fernandez J, Rota-Stabelli O, Huang D, Engel MS, Donoghue PCJ, Pisani D. The evolution of insect biodiversity. Curr Biol 2021; 31:R1299-R1311. [PMID: 34637741 DOI: 10.1016/j.cub.2021.08.057] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Insects comprise over half of all described animal species. Together with the Protura (coneheads), Collembola (springtails) and Diplura (two-pronged bristletails), insects form the Hexapoda, a terrestrial arthropod lineage characterised by possessing six legs. Exponential growth of genome-scale data for the hexapods has substantially altered our understanding of the origin and evolution of insect biodiversity. Phylogenomics has provided a new framework for reconstructing insect evolutionary history, resolving their position among the arthropods and some long-standing internal controversies such as the placement of the termites, twisted-winged insects, lice and fleas. However, despite the greatly increased size of phylogenomic datasets, contentious relationships among key insect clades remain unresolved. Further advances in insect phylogeny cannot rely on increased depth and breadth of genome and taxon sequencing. Improved modelling of the substitution process is fundamental to countering tree-reconstruction artefacts, while gene content, modelling of duplications and deletions, and comparative morphology all provide complementary lines of evidence to test hypotheses emerging from the analysis of sequence data. Finally, the integration of molecular and morphological data is key to the incorporation of fossil species within insect phylogeny. The emerging integrated framework of insect evolution will help explain the origins of insect megadiversity in terms of the evolution of their body plan, species diversity and ecology. Future studies of insect phylogeny should build upon an experimental, hypothesis-driven approach where the robustness of hypotheses generated is tested against increasingly realistic evolutionary models as well as complementary sources of phylogenetic evidence.
Collapse
Affiliation(s)
- Erik Tihelka
- School of Earth Sciences, University of Bristol, Bristol, UK; State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology, and Centre for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China.
| | - Chenyang Cai
- School of Earth Sciences, University of Bristol, Bristol, UK; State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology, and Centre for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China.
| | | | - Jesus Lozano-Fernandez
- School of Biological Sciences, University of Bristol, Bristol, UK; Institute of Evolutionary Biology (CSIC-UPF), Barcelona, Spain
| | - Omar Rota-Stabelli
- Research and Innovation Centre, Fondazione Edmund Mach, 38010 San Michele all Adige, Italy; Center Agriculture Food Environment, University of Trento, 38010 San Michele all Adige, Italy
| | - Diying Huang
- State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology, and Centre for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China
| | - Michael S Engel
- Division of Entomology, Natural History Museum, University of Kansas, Lawrence, KS, USA; Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, USA
| | | | - Davide Pisani
- School of Earth Sciences, University of Bristol, Bristol, UK; School of Biological Sciences, University of Bristol, Bristol, UK.
| |
Collapse
|