1
|
Bartoš O, Chmel M, Swierczková I. The overlooked evolutionary dynamics of 16S rRNA revises its role as the "gold standard" for bacterial species identification. Sci Rep 2024; 14:9067. [PMID: 38643216 PMCID: PMC11032355 DOI: 10.1038/s41598-024-59667-3] [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: 11/13/2023] [Accepted: 04/12/2024] [Indexed: 04/22/2024] Open
Abstract
The role of 16S rRNA has been and largely remains crucial for the identification of microbial organisms. Although 16S rRNA could certainly be described as one of the most studied sequences ever, the current view of it remains somewhat ambiguous. While some consider 16S rRNA to be a variable marker with resolution power down to the strain level, others consider them to be living fossils that carry information about the origin of domains of cellular life. We show that 16S rRNA is clearly an evolutionarily very rigid sequence, making it a largely unique and irreplaceable marker, but its applicability beyond the genus level is highly limited. Interestingly, it seems that the evolutionary rigidity is not driven by functional constraints of the sequence (RNA-protein interactions), but rather results from the characteristics of the host organism. Our results suggest that, at least in some lineages, Horizontal Gene Transfer (HGT) within genera plays an important role for the evolutionary non-dynamics (stasis) of 16S rRNA. Such genera exhibit an apparent lack of diversification at the 16S rRNA level in comparison to the rest of a genome. However, why it is limited specifically and solely to 16S rRNA remains enigmatic.
Collapse
Affiliation(s)
- Oldřich Bartoš
- Military Health Institute, Military Medical Agency, 16200, Prague, Czech Republic.
| | - Martin Chmel
- Military Health Institute, Military Medical Agency, 16200, Prague, Czech Republic
- Department of Infectious Diseases, First Faculty of Medicine, Charles University and Military University Hospital Prague, 12108, Prague, Czech Republic
| | - Iva Swierczková
- Military Health Institute, Military Medical Agency, 16200, Prague, Czech Republic
| |
Collapse
|
2
|
Xiao B, Rey-lglesia A, Yuan J, Hu J, Song S, Hou Y, Chen X, Germonpré M, Bao L, Wang S, Taogetongqimuge, Valentinovna LL, Lister AM, Lai X, Sheng G. Relationships of Late Pleistocene giant deer as revealed by Sinomegaceros mitogenomes from East Asia. iScience 2023; 26:108406. [PMID: 38047074 PMCID: PMC10690636 DOI: 10.1016/j.isci.2023.108406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 09/26/2023] [Accepted: 11/03/2023] [Indexed: 12/05/2023] Open
Abstract
The giant deer, widespread in northern Eurasia during the Late Pleistocene, have been classified as western Megaloceros and eastern Sinomegaceros through morphological studies. While Megaloceros's evolutionary history has been unveiled through mitogenomes, Sinomegaceros remains molecularly unexplored. Herein, we generated mitogenomes of giant deer from East Asia. We find that, in contrast to the morphological differences between Megaloceros and Sinomegaceros, they are mixed in the mitochondrial phylogeny, and Siberian specimens suggest a range contact or overlap between these two groups. Meanwhile, one deep divergent clade and another surviving until 20.1 thousand years ago (ka) were detected in northeastern China, the latter implying this area as a potential refugium during the Last Glacial Maximum (LGM). Moreover, stable isotope analyses indicate correlations between climate-introduced vegetation changes and giant deer extinction. Our study demonstrates the genetic relationship between eastern and western giant deer and explores the promoters of their extirpation in northern East Asia.
Collapse
Affiliation(s)
- Bo Xiao
- State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
- School of Earth Sciences, China University of Geosciences, Wuhan 430074, China
| | - Alba Rey-lglesia
- Globe Institute, University of Copenhagen, Copenhagen, 1350 Copenhagen K, Denmark
| | - Junxia Yuan
- State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
- Faculty of Materials Science and Chemistry, China University of Geosciences, Wuhan 430078, China
| | - Jiaming Hu
- State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
- School of Earth Sciences, China University of Geosciences, Wuhan 430074, China
| | - Shiwen Song
- State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
- School of Environmental Studies, China University of Geosciences, Wuhan 430078, China
| | - Yamei Hou
- Key Laboratory of Vertebrate Evolution and Human Origins, Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, Beijing 100044, China
| | - Xi Chen
- Department of Cultural Heritage and Museology, Nanjing Normal University, Nanjing 210046, China
| | - Mietje Germonpré
- Royal Belgian Institute of Natural Sciences, 1000 Brussels, Belgium
| | - Lei Bao
- Ordos Institute of Cultural Relics and Archaeology, Ordos 017010, China
| | | | | | - Lbova Liudmila Valentinovna
- Graduate School of International Relations, Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Grazhdansky Av., 28, Russia
| | | | - Xulong Lai
- State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
- School of Earth Sciences, China University of Geosciences, Wuhan 430074, China
| | - Guilian Sheng
- State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
- School of Environmental Studies, China University of Geosciences, Wuhan 430078, China
| |
Collapse
|
3
|
Comte A, Tricou T, Tannier E, Joseph J, Siberchicot A, Penel S, Allio R, Delsuc F, Dray S, de Vienne DM. PhylteR: Efficient Identification of Outlier Sequences in Phylogenomic Datasets. Mol Biol Evol 2023; 40:msad234. [PMID: 37879113 PMCID: PMC10655845 DOI: 10.1093/molbev/msad234] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Revised: 09/29/2023] [Accepted: 10/18/2023] [Indexed: 10/27/2023] Open
Abstract
In phylogenomics, incongruences between gene trees, resulting from both artifactual and biological reasons, can decrease the signal-to-noise ratio and complicate species tree inference. The amount of data handled today in classical phylogenomic analyses precludes manual error detection and removal. However, a simple and efficient way to automate the identification of outliers from a collection of gene trees is still missing. Here, we present PhylteR, a method that allows rapid and accurate detection of outlier sequences in phylogenomic datasets, i.e. species from individual gene trees that do not follow the general trend. PhylteR relies on DISTATIS, an extension of multidimensional scaling to 3 dimensions to compare multiple distance matrices at once. In PhylteR, these distance matrices extracted from individual gene phylogenies represent evolutionary distances between species according to each gene. On simulated datasets, we show that PhylteR identifies outliers with more sensitivity and precision than a comparable existing method. We also show that PhylteR is not sensitive to ILS-induced incongruences, which is a desirable feature. On a biological dataset of 14,463 genes for 53 species previously assembled for Carnivora phylogenomics, we show (i) that PhylteR identifies as outliers sequences that can be considered as such by other means, and (ii) that the removal of these sequences improves the concordance between the gene trees and the species tree. Thanks to the generation of numerous graphical outputs, PhylteR also allows for the rapid and easy visual characterization of the dataset at hand, thus aiding in the precise identification of errors. PhylteR is distributed as an R package on CRAN and as containerized versions (docker and singularity).
Collapse
Affiliation(s)
- Aurore Comte
- French Institute of Bioinformatics (IFB)—South Green Bioinformatics Platform, Bioversity, CIRAD, INRAE, IRD, Montpellier, France
- IRD, CIRAD, INRAE, Institut Agro, PHIM Plant Health Institute, Montpellier University, Montpellier, France
| | - Théo Tricou
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
| | - Eric Tannier
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
- Centre de Recherches Inria de Lyon, Villeurbanne, France
| | - Julien Joseph
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
| | - Aurélie Siberchicot
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
| | - Simon Penel
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
| | - Rémi Allio
- CBGP, INRAE, CIRAD, IRD, Montpellier SupAgro, Univ. Montpellier, Montpellier, France
| | | | - Stéphane Dray
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
| | - Damien M de Vienne
- Université de Lyon, Université Lyon 1, UMR CNRS 5558 Laboratoire de Biométrie et Biologie Évolutive, Villeurbanne, France
| |
Collapse
|
4
|
DeRaad DA, McCullough JM, DeCicco LH, Hime PM, Joseph L, Andersen MJ, Moyle RG. Mitonuclear discordance results from incomplete lineage sorting, with no detectable evidence for gene flow, in a rapid radiation of Todiramphus kingfishers. Mol Ecol 2023; 32:4844-4862. [PMID: 37515525 DOI: 10.1111/mec.17080] [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/28/2022] [Revised: 03/15/2023] [Accepted: 06/19/2023] [Indexed: 07/31/2023]
Abstract
Many organisms possess multiple discrete genomes (i.e. nuclear and organellar), which are inherited separately and may have unique and even conflicting evolutionary histories. Phylogenetic reconstructions from these discrete genomes can yield different patterns of relatedness, a phenomenon known as cytonuclear discordance. In many animals, mitonuclear discordance (i.e. discordant evolutionary histories between the nuclear and mitochondrial genomes) has been widely documented, but its causes are often considered idiosyncratic and inscrutable. We show that a case of mitonuclear discordance in Todiramphus kingfishers can be explained by extensive genome-wide incomplete lineage sorting (ILS), likely a result of the explosive diversification history of this genus. For these kingfishers, quartet frequencies reveal that the nuclear genome is dominated by discordant topologies, with none of the internal branches in our consensus nuclear tree recovered in >50% of genome-wide gene trees. Meanwhile, a lack of inter-species shared ancestry, non-significant pairwise tests for gene flow, and little evidence for meaningful migration edges between species, leads to the conclusion that gene flow cannot explain the mitonuclear discordance we observe. This lack of evidence for gene flow combined with evidence for extensive genome-wide gene tree discordance, a hallmark of ILS, leads us to conclude that the mitonuclear discordance we observe likely results from ILS, specifically deep coalescence of the mitochondrial genome. Based on this case study, we hypothesize that similar demographic histories in other 'great speciator' taxa across the Indo-Pacific likely predispose these groups to high levels of ILS and high likelihoods of mitonuclear discordance.
Collapse
Affiliation(s)
- Devon A DeRaad
- Biodiversity Institute and Natural History Museum, University of Kansas, Lawrence, Kansas, USA
| | - Jenna M McCullough
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Lucas H DeCicco
- Biodiversity Institute and Natural History Museum, University of Kansas, Lawrence, Kansas, USA
| | - Paul M Hime
- Biodiversity Institute and Natural History Museum, University of Kansas, Lawrence, Kansas, USA
| | - Leo Joseph
- Australian National Wildlife Collection, CSIRO National Research Collections Australia, Canberra, Australian Capital Territory, Australia
| | - Michael J Andersen
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Robert G Moyle
- Biodiversity Institute and Natural History Museum, University of Kansas, Lawrence, Kansas, USA
| |
Collapse
|
5
|
Contradictory Phylogenetic Signals in the Laurasiatheria Anomaly Zone. Genes (Basel) 2022; 13:genes13050766. [PMID: 35627151 PMCID: PMC9141728 DOI: 10.3390/genes13050766] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 04/12/2022] [Accepted: 04/21/2022] [Indexed: 02/04/2023] Open
Abstract
Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.
Collapse
|
6
|
Sharaf A, Elateek S. Orthology Prediction and Phylogenetic Analysis Methods in Plants. Methods Mol Biol 2022; 2512:3-22. [PMID: 35817996 DOI: 10.1007/978-1-0716-2429-6_1] [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] [Indexed: 06/15/2023]
Abstract
In this chapter, we outline a pipeline for ortholog prediction and phylogenetic analysis in plants. This computational pipeline uses algorithms from different software to enable bioinformatic-beginner biologists to predict orthologs that can be shared with many distinct plant nonmodel and model species and identify gene loss events. Prediction of orthologs allows (1) investigation of the evolutionary relationships of plant genomes, (2) discovery of their origin, function, and (3) the impact of their adaptability to the environment.We developed a pipeline to fit, not only eukaryote but also prokaryote organisms, with small or large genomes. All results acquired from the orthologs predication will enable phylogenetic tree construction, using gene and species (phylogenomic) phylogeny approaches.
Collapse
Affiliation(s)
- Abdoallah Sharaf
- Institute of Plant Molecular Biology, Biology Centre, Czech Academy of Sciences, České Budějovice, Czech Republic.
- Genetic Department, Faculty of Agriculture, Ain Shams University, Cairo, Egypt.
| | - Sawsan Elateek
- Genetic Department, Faculty of Agriculture, Ain Shams University, Cairo, Egypt
| |
Collapse
|
7
|
Literman R, Schwartz R. Genome-Scale Profiling Reveals Noncoding Loci Carry Higher Proportions of Concordant Data. Mol Biol Evol 2021; 38:2306-2318. [PMID: 33528497 PMCID: PMC8136493 DOI: 10.1093/molbev/msab026] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Many evolutionary relationships remain controversial despite whole-genome sequencing data. These controversies arise, in part, due to challenges associated with accurately modeling the complex phylogenetic signal coming from genomic regions experiencing distinct evolutionary forces. Here, we examine how different regions of the genome support or contradict well-established relationships among three mammal groups using millions of orthologous parsimony-informative biallelic sites (PIBS) distributed across primate, rodent, and Pecora genomes. We compared PIBS concordance percentages among locus types (e.g. coding sequences (CDS), introns, intergenic regions), and contrasted PIBS utility over evolutionary timescales. Sites derived from noncoding sequences provided more data and proportionally more concordant sites compared with those from CDS in all clades. CDS PIBS were also predominant drivers of tree incongruence in two cases of topological conflict. PIBS derived from most locus types provided surprisingly consistent support for splitting events spread across the timescales we examined, although we find evidence that CDS and intronic PIBS may, respectively and to a limited degree, inform disproportionately about older and younger splits. In this era of accessible wholegenome sequence data, these results:1) suggest benefits to more intentionally focusing on noncoding loci as robust data for tree inference and 2) reinforce the importance of accurate modeling, especially when using CDS data.
Collapse
Affiliation(s)
- Robert Literman
- Department of Biological Sciences, University of Rhode Island, South Kingstown, RI, USA.,Center for Food Safety and Applied Nutrition, Office of Regulatory Science, U.S. Food and Drug Administration, College Park, MD, USA
| | - Rachel Schwartz
- Department of Biological Sciences, University of Rhode Island, South Kingstown, RI, USA
| |
Collapse
|
8
|
Cerca J, Rivera-Colón AG, Ferreira MS, Ravinet M, Nowak MD, Catchen JM, Struck TH. Incomplete lineage sorting and ancient admixture, and speciation without morphological change in ghost-worm cryptic species. PeerJ 2021; 9:e10896. [PMID: 33614296 PMCID: PMC7879940 DOI: 10.7717/peerj.10896] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 01/13/2021] [Indexed: 12/14/2022] Open
Abstract
Morphologically similar species, that is cryptic species, may be similar or quasi-similar owing to the deceleration of morphological evolution and stasis. While the factors underlying the deceleration of morphological evolution or stasis in cryptic species remain unknown, decades of research in the field of paleontology on punctuated equilibrium have originated clear hypotheses. Species are expected to remain morphologically identical in scenarios of shared genetic variation, such as hybridization and incomplete lineage sorting, or in scenarios where bottlenecks reduce genetic variation and constrain the evolution of morphology. Here, focusing on three morphologically similar Stygocapitella species, we employ a whole-genome amplification method (WGA) coupled with double-digestion restriction-site associated DNA sequencing (ddRAD) to reconstruct the evolutionary history of the species complex. We explore population structure, use population-level statistics to determine the degree of connectivity between populations and species, and determine the most likely demographic scenarios which generally reject for recent hybridization. We find that the combination of WGA and ddRAD allowed us to obtain genomic-level data from microscopic eukaryotes (∼1 millimetre) opening up opportunities for those working with population genomics and phylogenomics in such taxa. The three species share genetic variance, likely from incomplete lineage sorting and ancient admixture. We speculate that the degree of shared variation might underlie morphological similarity in the Atlantic species complex.
Collapse
Affiliation(s)
- José Cerca
- Department of Environmental Science, Policy, and Management, University of California, University of California, Berkeley, Berkeley, CA, United States of America
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
- Natural History Museum, University of Oslo, Oslo, Norway
| | - Angel G. Rivera-Colón
- Department of Evolution, Ecology, and Behavior, University of Illinois at Urbana-Champaign, Urbana Champaign, IL, United States of America
| | - Mafalda S. Ferreira
- Division of Biological Sciences, University of Montana, Missoula, MT, United States of America
- Departamento de Biologia, Universidade do Porto, Porto, Porto, Portugal
- CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Porto, Porto, Portugal
| | - Mark Ravinet
- School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Oslo, Norway
| | | | - Julian M. Catchen
- Department of Evolution, Ecology, and Behavior, University of Illinois at Urbana-Champaign, Urbana Champaign, IL, United States of America
| | | |
Collapse
|
9
|
Jiang X, Edwards SV, Liu L. The Multispecies Coalescent Model Outperforms Concatenation Across Diverse Phylogenomic Data Sets. Syst Biol 2021; 69:795-812. [PMID: 32011711 PMCID: PMC7302055 DOI: 10.1093/sysbio/syaa008] [Citation(s) in RCA: 35] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Revised: 12/24/2019] [Accepted: 01/02/2020] [Indexed: 11/30/2022] Open
Abstract
A statistical framework of model comparison and model validation is essential to resolving the debates over concatenation and coalescent models in phylogenomic data analysis. A set of statistical tests are here applied and developed to evaluate and compare the adequacy of substitution, concatenation, and multispecies coalescent (MSC) models across 47 phylogenomic data sets collected across tree of life. Tests for substitution models and the concatenation assumption of topologically congruent gene trees suggest that a poor fit of substitution models, rejected by 44% of loci, and concatenation models, rejected by 38% of loci, is widespread. Logistic regression shows that the proportions of GC content and informative sites are both negatively correlated with the fit of substitution models across loci. Moreover, a substantial violation of the concatenation assumption of congruent gene trees is consistently observed across six major groups (birds, mammals, fish, insects, reptiles, and others, including other invertebrates). In contrast, among those loci adequately described by a given substitution model, the proportion of loci rejecting the MSC model is 11%, significantly lower than those rejecting the substitution and concatenation models. Although conducted on reduced data sets due to computational constraints, Bayesian model validation and comparison both strongly favor the MSC over concatenation across all data sets; the concatenation assumption of congruent gene trees rarely holds for phylogenomic data sets with more than 10 loci. Thus, for large phylogenomic data sets, model comparisons are expected to consistently and more strongly favor the coalescent model over the concatenation model. We also found that loci rejecting the MSC have little effect on species tree estimation. Our study reveals the value of model validation and comparison in phylogenomic data analysis, as well as the need for further improvements of multilocus models and computational tools for phylogenetic inference. [Bayes factor; Bayesian model validation; coalescent prior; congruent gene trees; independent prior; Metazoa; posterior predictive simulation.]
Collapse
Affiliation(s)
- Xiaodong Jiang
- Department of Statistics, University of Georgia, 310 Herty Drive, Athens, GA 30602, USA
| | - Scott V Edwards
- Department of Organismic and Evolutionary Biology and Museum of Comparative Zoology, Harvard, 26 Oxford Street, Cambridge, MA 02138, USA
| | - Liang Liu
- Department of Statistics, University of Georgia, 310 Herty Drive, Athens, GA 30602, USA.,Institute of Bioinformatics, University of Georgia, 120 Green Street, Athens, GA 30602, USA
| |
Collapse
|
10
|
Van Dam MH, Henderson JB, Esposito L, Trautwein M. Genomic Characterization and Curation of UCEs Improves Species Tree Reconstruction. Syst Biol 2020; 70:307-321. [PMID: 32750133 PMCID: PMC7875437 DOI: 10.1093/sysbio/syaa063] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2019] [Revised: 07/26/2020] [Accepted: 07/29/2020] [Indexed: 12/12/2022] Open
Abstract
Ultraconserved genomic elements (UCEs) are generally treated as independent loci in phylogenetic analyses. The identification pipeline for UCE probes does not require prior knowledge of genetic identity, only selecting loci that are highly conserved, single copy, without repeats, and of a particular length. Here, we characterized UCEs from 11 phylogenomic studies across the animal tree of life, from birds to marine invertebrates. We found that within vertebrate lineages, UCEs are mostly intronic and intergenic, while in invertebrates, the majority are in exons. We then curated four different sets of UCE markers by genomic category from five different studies including: birds, mammals, fish, Hymenoptera (ants, wasps, and bees), and Coleoptera (beetles). Of genes captured by UCEs, we find that many are represented by two or more UCEs, corresponding to nonoverlapping segments of a single gene. We considered these UCEs to be nonindependent, merged all UCEs that belonged to a particular gene, constructed gene and species trees, and then evaluated the subsequent effect of merging cogenic UCEs on gene and species tree reconstruction. Average bootstrap support for merged UCE gene trees was significantly improved across all data sets apparently driven by the increase in loci length. Additionally, we conducted simulations and found that gene trees generated from merged UCEs were more accurate than those generated by unmerged UCEs. As loci length improves gene tree accuracy, this modest degree of UCE characterization and curation impacts downstream analyses and demonstrates the advantages of incorporating basic genomic characterizations into phylogenomic analyses. [Anchored hybrid enrichment; ants; ASTRAL; bait capture; carangimorph; Coleoptera; conserved nonexonic elements; exon capture; gene tree; Hymenoptera; mammal; phylogenomic markers; songbird; species tree; ultraconserved elements; weevils.]
Collapse
Affiliation(s)
- Matthew H Van Dam
- Entomology Department, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA.,Center for Comparative Genomics, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA
| | - James B Henderson
- Center for Comparative Genomics, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA
| | - Lauren Esposito
- Entomology Department, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA.,Center for Comparative Genomics, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA
| | - Michelle Trautwein
- Entomology Department, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA.,Center for Comparative Genomics, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, CA 94118, USA
| |
Collapse
|
11
|
Abstract
Knowing phylogenetic relationships among species is fundamental for many studies in biology. An accurate phylogenetic tree underpins our understanding of the major transitions in evolution, such as the emergence of new body plans or metabolism, and is key to inferring the origin of new genes, detecting molecular adaptation, understanding morphological character evolution and reconstructing demographic changes in recently diverged species. Although data are ever more plentiful and powerful analysis methods are available, there remain many challenges to reliable tree building. Here, we discuss the major steps of phylogenetic analysis, including identification of orthologous genes or proteins, multiple sequence alignment, and choice of substitution models and inference methodologies. Understanding the different sources of errors and the strategies to mitigate them is essential for assembling an accurate tree of life.
Collapse
|
12
|
Swanson MT, Oliveros CH, Esselstyn JA. A phylogenomic rodent tree reveals the repeated evolution of masseter architectures. Proc Biol Sci 2020; 286:20190672. [PMID: 31064307 DOI: 10.1098/rspb.2019.0672] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
Understanding the number of times a trait has evolved is a necessary foundation for comprehending its potential relationships with selective regimes, developmental constraints and evolutionary diversification. Rodents make up over 40% of extant mammalian species, and their ecological and evolutionary success has been partially attributed to the increase in biting efficiency that resulted from a forward shift of one or two portions of the masseter muscle from the zygomatic arch onto the rostrum. This forward shift has occurred in three discrete ways, but the number of times it has occurred has never been explicitly quantified. We estimated an ultrametric phylogeny, the first to include all rodent families, using thousands of ultraconserved elements. We examined support for evolutionary relationships among the five rodent suborders and then incorporated relevant fossils, fitted models of character evolution, and used stochastic character mapping to determine that a portion of the masseter muscle has moved forward onto the rostrum at least seven times (with one reversal) during the approximately 70 Myr history of rodents. Combined, the repeated evolution of this key innovation, its increasing prevalence through time, and the species diversity of clades with this character underscores the adaptive value of improved biting efficiency and the relative ease with which some advantageous traits arise.
Collapse
Affiliation(s)
- Mark T Swanson
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University , Baton Rouge, LA , USA
| | - Carl H Oliveros
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University , Baton Rouge, LA , USA
| | - Jacob A Esselstyn
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University , Baton Rouge, LA , USA
| |
Collapse
|
13
|
Islam M, Sarker K, Das T, Reaz R, Bayzid MS. STELAR: a statistically consistent coalescent-based species tree estimation method by maximizing triplet consistency. BMC Genomics 2020; 21:136. [PMID: 32039704 PMCID: PMC7011378 DOI: 10.1186/s12864-020-6519-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2019] [Accepted: 01/20/2020] [Indexed: 12/14/2022] Open
Abstract
Background Species tree estimation is frequently based on phylogenomic approaches that use multiple genes from throughout the genome. However, estimating a species tree from a collection of gene trees can be complicated due to the presence of gene tree incongruence resulting from incomplete lineage sorting (ILS), which is modelled by the multi-species coalescent process. Maximum likelihood and Bayesian MCMC methods can potentially result in accurate trees, but they do not scale well to large datasets. Results We present STELAR (Species Tree Estimation by maximizing tripLet AgReement), a new fast and highly accurate statistically consistent coalescent-based method for estimating species trees from a collection of gene trees. We formalized the constrained triplet consensus (CTC) problem and showed that the solution to the CTC problem is a statistically consistent estimate of the species tree under the multi-species coalescent (MSC) model. STELAR is an efficient dynamic programming based solution to the CTC problem which is highly accurate and scalable. We evaluated the accuracy of STELAR in comparison with SuperTriplets, which is an alternate fast and highly accurate triplet-based supertree method, and with MP-EST and ASTRAL – two of the most popular and accurate coalescent-based methods. Experimental results suggest that STELAR matches the accuracy of ASTRAL and improves on MP-EST and SuperTriplets. Conclusions Theoretical and empirical results (on both simulated and real biological datasets) suggest that STELAR is a valuable technique for species tree estimation from gene tree distributions.
Collapse
Affiliation(s)
- Mazharul Islam
- Department of Computer Science and Engineering, Bangladesh University of Engineering and Technology, Dhaka, 1205, Bangladesh
| | - Kowshika Sarker
- Department of Computer Science and Engineering, Bangladesh University of Engineering and Technology, Dhaka, 1205, Bangladesh
| | - Trisha Das
- Department of Computer Science and Engineering, Bangladesh University of Engineering and Technology, Dhaka, 1205, Bangladesh
| | - Rezwana Reaz
- Department of Computer Science, The University of Texas at Austin, Texas, 78712, USA
| | - Md Shamsuzzoha Bayzid
- Department of Computer Science and Engineering, Bangladesh University of Engineering and Technology, Dhaka, 1205, Bangladesh.
| |
Collapse
|
14
|
Li G, Figueiró HV, Eizirik E, Murphy WJ. Recombination-Aware Phylogenomics Reveals the Structured Genomic Landscape of Hybridizing Cat Species. Mol Biol Evol 2020; 36:2111-2126. [PMID: 31198971 PMCID: PMC6759079 DOI: 10.1093/molbev/msz139] [Citation(s) in RCA: 69] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Current phylogenomic approaches implicitly assume that the predominant phylogenetic signal within a genome reflects the true evolutionary history of organisms, without assessing the confounding effects of postspeciation gene flow that can produce a mosaic of phylogenetic signals that interact with recombinational variation. Here, we tested the validity of this assumption with a phylogenomic analysis of 27 species of the cat family, assessing local effects of recombination rate on species tree inference and divergence time estimation across their genomes. We found that the prevailing phylogenetic signal within the autosomes is not always representative of the most probable speciation history, due to ancient hybridization throughout felid evolution. Instead, phylogenetic signal was concentrated within regions of low recombination, and notably enriched within large X chromosome recombination cold spots that exhibited recurrent patterns of strong genetic differentiation and selective sweeps across mammalian orders. By contrast, regions of high recombination were enriched for signatures of ancient gene flow, and these sequences inflated crown-lineage divergence times by ∼40%. We conclude that existing phylogenomic approaches to infer the Tree of Life may be highly misleading without considering the genomic architecture of phylogenetic signal relative to recombination rate and its interplay with historical hybridization.
Collapse
Affiliation(s)
- Gang Li
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX
| | - Henrique V Figueiró
- PUCRS, Escola de Ciências, Laboratory of Genomics and Molecular Biology, Porto Alegre, Brazil.,INCT-EECBio, Brazil
| | - Eduardo Eizirik
- PUCRS, Escola de Ciências, Laboratory of Genomics and Molecular Biology, Porto Alegre, Brazil.,INCT-EECBio, Brazil
| | - William J Murphy
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX
| |
Collapse
|
15
|
Upham NS, Esselstyn JA, Jetz W. Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation. PLoS Biol 2019; 17:e3000494. [PMID: 31800571 PMCID: PMC6892540 DOI: 10.1371/journal.pbio.3000494] [Citation(s) in RCA: 471] [Impact Index Per Article: 94.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2019] [Accepted: 10/24/2019] [Indexed: 12/18/2022] Open
Abstract
Big, time-scaled phylogenies are fundamental to connecting evolutionary processes to modern biodiversity patterns. Yet inferring reliable phylogenetic trees for thousands of species involves numerous trade-offs that have limited their utility to comparative biologists. To establish a robust evolutionary timescale for all approximately 6,000 living species of mammals, we developed credible sets of trees that capture root-to-tip uncertainty in topology and divergence times. Our "backbone-and-patch" approach to tree building applies a newly assembled 31-gene supermatrix to two levels of Bayesian inference: (1) backbone relationships and ages among major lineages, using fossil node or tip dating, and (2) species-level "patch" phylogenies with nonoverlapping in-groups that each correspond to one representative lineage in the backbone. Species unsampled for DNA are either excluded ("DNA-only" trees) or imputed within taxonomic constraints using branch lengths drawn from local birth-death models ("completed" trees). Joining time-scaled patches to backbones results in species-level trees of extant Mammalia with all branches estimated under the same modeling framework, thereby facilitating rate comparisons among lineages as disparate as marsupials and placentals. We compare our phylogenetic trees to previous estimates of mammal-wide phylogeny and divergence times, finding that (1) node ages are broadly concordant among studies, and (2) recent (tip-level) rates of speciation are estimated more accurately in our study than in previous "supertree" approaches, in which unresolved nodes led to branch-length artifacts. Credible sets of mammalian phylogenetic history are now available for download at http://vertlife.org/phylosubsets, enabling investigations of long-standing questions in comparative biology.
Collapse
Affiliation(s)
- Nathan S. Upham
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Center for Biodiversity & Global Change, Yale University, New Haven, Connecticut, United States of America
| | - Jacob A. Esselstyn
- Department of Biological Sciences and Museum of Natural Science, Louisiana State University, Baton Rouge, Louisiana, United States of America
| | - Walter Jetz
- Department of Ecology & Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
- Center for Biodiversity & Global Change, Yale University, New Haven, Connecticut, United States of America
| |
Collapse
|
16
|
Zhang ML, Li ML, Ayoola AO, Murphy RW, Wu DD, Shao Y. Conserved sequences identify the closest living relatives of primates. Zool Res 2019; 40:532-540. [PMID: 31393097 PMCID: PMC6822925 DOI: 10.24272/j.issn.2095-8137.2019.057] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Elucidating the closest living relatives of extant primates is essential for fully understanding important biological processes related to the genomic and phenotypic evolution of primates, especially of humans. However, the phylogenetic placement of these primate relatives remains controversial, with three primary hypotheses currently espoused based on morphological and molecular evidence. In the present study, we used two algorithms to analyze differently partitioned genomic datasets consisting of 45.4 Mb of conserved non-coding elements and 393 kb of concatenated coding sequences to test these hypotheses. We assessed different genomic histories and compared with other molecular studies found solid support for colugos being the closest living relatives of primates. Our phylogeny showed Cercopithecinae to have low levels of nucleotide divergence, especially for Papionini, and gibbons to have a high rate of divergence. The MCMCtree comprehensively updated divergence dates of early evolution of Primatomorpha and Primates.
Collapse
Affiliation(s)
- Mei-Ling Zhang
- Department of Acute Infectious Diseases Control and Prevention, Yunnan Provincial Centre for Disease Control and Prevention, Kunming Yunnan 650022, China
| | - Ming-Li Li
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China.,Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming Yunnan 650223, China
| | - Adeola Oluwakemi Ayoola
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China.,Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming Yunnan 650223, China
| | - Robert W Murphy
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China.,Centre for Biodiversity and Conservation Biology, Royal Ontario Museum, Toronto M5S 2C6, Canada
| | - Dong-Dong Wu
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China, E-mail:
| | - Yong Shao
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China, E-mail:
| |
Collapse
|
17
|
Dapper AL, Payseur BA. Molecular evolution of the meiotic recombination pathway in mammals. Evolution 2019; 73:2368-2389. [PMID: 31579931 DOI: 10.1111/evo.13850] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Accepted: 09/07/2019] [Indexed: 02/06/2023]
Abstract
Meiotic recombination shapes evolution and helps to ensure proper chromosome segregation in most species that reproduce sexually. Recombination itself evolves, with species showing considerable divergence in the rate of crossing-over. However, the genetic basis of this divergence is poorly understood. Recombination events are produced via a complicated, but increasingly well-described, cellular pathway. We apply a phylogenetic comparative approach to a carefully selected panel of genes involved in the processes leading to crossovers-spanning double-strand break formation, strand invasion, the crossover/non-crossover decision, and resolution-to reconstruct the evolution of the recombination pathway in eutherian mammals and identify components of the pathway likely to contribute to divergence between species. Eleven recombination genes, predominantly involved in the stabilization of homologous pairing and the crossover/non-crossover decision, show evidence of rapid evolution and positive selection across mammals. We highlight TEX11 and associated genes involved in the synaptonemal complex and the early stages of the crossover/non-crossover decision as candidates for the evolution of recombination rate. Evolutionary comparisons to MLH1 count, a surrogate for the number of crossovers, reveal a positive correlation between genome-wide recombination rate and the rate of evolution at TEX11 across the mammalian phylogeny. Our results illustrate the power of viewing the evolution of recombination from a pathway perspective.
Collapse
Affiliation(s)
- Amy L Dapper
- Laboratory of Genetics, University of Wisconsin, Madison, Wisconsin, 53706.,Department of Biological Sciences, Mississippi State University, Mississippi, 39762
| | - Bret A Payseur
- Laboratory of Genetics, University of Wisconsin, Madison, Wisconsin, 53706
| |
Collapse
|
18
|
Hamilton CA, St Laurent RA, Dexter K, Kitching IJ, Breinholt JW, Zwick A, Timmermans MJTN, Barber JR, Kawahara AY. Phylogenomics resolves major relationships and reveals significant diversification rate shifts in the evolution of silk moths and relatives. BMC Evol Biol 2019; 19:182. [PMID: 31533606 PMCID: PMC6751749 DOI: 10.1186/s12862-019-1505-1] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2019] [Accepted: 08/29/2019] [Indexed: 03/13/2023] Open
Abstract
BACKGROUND Silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying spectacular forms and diverse ecological traits, relatively little attention has been given to understanding their evolution and drivers of their diversity. To begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, with a newly developed DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections. RESULTS The well-supported tree was overall consistent with prior morphological and molecular studies, although some taxa were misplaced. The bombycid Arotros Schaus was formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths). CONCLUSIONS Our study establishes a backbone for future evolutionary, comparative, and taxonomic studies of Bombycoidea. We postulate that the rate shifts identified are due to the well-documented bat-moth "arms race". Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.
Collapse
Affiliation(s)
- C A Hamilton
- Florida Museum of Natural History, University of Florida, Gainesville, FL, 32611, USA.
- Department of Entomology, Plant Pathology & Nematology, University of Idaho, Moscow, ID, 83844, USA.
| | - R A St Laurent
- Florida Museum of Natural History, University of Florida, Gainesville, FL, 32611, USA
| | - K Dexter
- Florida Museum of Natural History, University of Florida, Gainesville, FL, 32611, USA
| | - I J Kitching
- Department of Life Sciences, Natural History Museum, Cromwell Road, London, SW7 5BD, UK
| | - J W Breinholt
- Florida Museum of Natural History, University of Florida, Gainesville, FL, 32611, USA
- RAPiD Genomics, 747 SW 2nd Avenue #314, Gainesville, FL, 32601, USA
| | - A Zwick
- Australian National Insect Collection, CSIRO, Clunies Ross St, Acton, ACT, Canberra, 2601, Australia
| | - M J T N Timmermans
- Department of Natural Sciences, Middlesex University, The Burroughs, London, NW4 4BT, UK
| | - J R Barber
- Department of Biological Sciences, Boise State University, Boise, ID, 83725, USA
| | - A Y Kawahara
- Florida Museum of Natural History, University of Florida, Gainesville, FL, 32611, USA.
| |
Collapse
|
19
|
Mendes FK, Livera AP, Hahn MW. The perils of intralocus recombination for inferences of molecular convergence. Philos Trans R Soc Lond B Biol Sci 2019; 374:20180244. [PMID: 31154973 DOI: 10.1098/rstb.2018.0244] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Accurate inferences of convergence require that the appropriate tree topology be used. If there is a mismatch between the tree a trait has evolved along and the tree used for analysis, then false inferences of convergence ('hemiplasy') can occur. To avoid problems of hemiplasy when there are high levels of gene tree discordance with the species tree, researchers have begun to construct tree topologies from individual loci. However, due to intralocus recombination, even locus-specific trees may contain multiple topologies within them. This implies that the use of individual tree topologies discordant with the species tree can still lead to incorrect inferences about molecular convergence. Here, we examine the frequency with which single exons and single protein-coding genes contain multiple underlying tree topologies, in primates and Drosophila, and quantify the effects of hemiplasy when using trees inferred from individual loci. In both clades, we find that there are most often multiple diagnosable topologies within single exons and whole genes, with 91% of Drosophila protein-coding genes containing multiple topologies. Because of this underlying topological heterogeneity, even using trees inferred from individual protein-coding genes results in 25% and 38% of substitutions falsely labelled as convergent in primates and Drosophila, respectively. While constructing local trees can reduce the problem of hemiplasy, our results suggest that it will be difficult to completely avoid false inferences of convergence. We conclude by suggesting several ways forward in the analysis of convergent evolution, for both molecular and morphological characters. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.
Collapse
Affiliation(s)
- Fábio K Mendes
- 1 Department of Computer Science, The University of Auckland , Auckland 1010 , New Zealand.,2 Department of Biology, Indiana University , Bloomington, IN 47405 , USA
| | - Andrew P Livera
- 2 Department of Biology, Indiana University , Bloomington, IN 47405 , USA
| | - Matthew W Hahn
- 2 Department of Biology, Indiana University , Bloomington, IN 47405 , USA.,3 Department of Computer Science, Indiana University , Bloomington, IN 47405 , USA
| |
Collapse
|
20
|
Hirano T, Saito T, Tsunamoto Y, Koseki J, Ye B, Do VT, Miura O, Suyama Y, Chiba S. Enigmatic incongruence between mtDNA and nDNA revealed by multi-locus phylogenomic analyses in freshwater snails. Sci Rep 2019; 9:6223. [PMID: 30996240 PMCID: PMC6470147 DOI: 10.1038/s41598-019-42682-0] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2018] [Accepted: 04/05/2019] [Indexed: 02/08/2023] Open
Abstract
Phylogenetic incongruence has frequently been encountered among different molecular markers. Recent progress in molecular phylogenomics has provided detailed and important information for evolutionary biology and taxonomy. Here we focused on the freshwater viviparid snails (Cipangopaludina chinensis chinensis and C. c. laeta) of East Asia. We conducted phylogenetic analyses and divergence time estimation using two mitochondrial markers. We also performed population genetic analyses using genome-wide SNPs. We investigated how and which phylogenetic patterns reflect shell morphology. The results showed these two species could be separated into four major mitochondrial clades, whereas the nuclear clusters supported two groups. The phylogenetic patterns of both mtDNA and nDNA largely reflected the geographical distribution. Shell morphology reflected the phylogenetic clusters based on nDNA. The findings also showed these two species diversified in the Pliocene to early Pleistocene era, and occurred introgressive hybridisation. The results also raise the taxonomic issue of the two species.
Collapse
Affiliation(s)
- Takahiro Hirano
- Department of Biological Sciences, University of Idaho, Moscow, Idaho, USA.
| | - Takumi Saito
- Graduate school of Life Sciences, Tohoku University, Miyagi, Japan
| | - Yoshihiro Tsunamoto
- Kawatabi Field Science Center, Graduate School of Agricultural Science, Tohoku University, Miyagi, Japan
| | - Joichiro Koseki
- Graduate school of Life Sciences, Tohoku University, Miyagi, Japan
| | - Bin Ye
- Graduate school of Life Sciences, Tohoku University, Miyagi, Japan
- Agricultural Experiment Station, Zhejiang University, Hangzhou, China
| | - Van Tu Do
- Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Osamu Miura
- Faculty of Agriculture and Marine Science, Kochi University, Kochi, Japan
| | - Yoshihisa Suyama
- Kawatabi Field Science Center, Graduate School of Agricultural Science, Tohoku University, Miyagi, Japan
| | - Satoshi Chiba
- Graduate school of Life Sciences, Tohoku University, Miyagi, Japan
- Center for Northeast Asian Studies, Tohoku University, Miyagi, Japan
| |
Collapse
|
21
|
Advances in Computational Methods for Phylogenetic Networks in the Presence of Hybridization. BIOINFORMATICS AND PHYLOGENETICS 2019. [DOI: 10.1007/978-3-030-10837-3_13] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
22
|
Mongiardino Koch N, Coppard SE, Lessios HA, Briggs DEG, Mooi R, Rouse GW. A phylogenomic resolution of the sea urchin tree of life. BMC Evol Biol 2018; 18:189. [PMID: 30545284 PMCID: PMC6293586 DOI: 10.1186/s12862-018-1300-4] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Accepted: 11/19/2018] [Indexed: 12/21/2022] Open
Abstract
BACKGROUND Echinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved. RESULTS We performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We show that phylogenetic signal for novel resolutions of these lineages is strong and distributed throughout the genome, and fail to recover systematic biases as drivers of our results. CONCLUSIONS Our investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.
Collapse
Affiliation(s)
| | - Simon E. Coppard
- Department of Biology, Hamilton College, Clinton, NY USA
- Smithsonian Tropical Research Institute, Balboa, Panama
| | | | - Derek E. G. Briggs
- Department of Geology and Geophysics, Yale University, New Haven, CT USA
- Peabody Museum of Natural History, Yale University, New Haven, CT USA
| | - Rich Mooi
- Department of Invertebrate Zoology and Geology, California Academy of Sciences, San Francisco, CA USA
| | - Greg W. Rouse
- Scripps Institution of Oceanography, UC San Diego, La Jolla, CA USA
| |
Collapse
|
23
|
Gossmann TI, Bockwoldt M, Diringer L, Schwarz F, Schumann VF. Evidence for Strong Fixation Bias at 4-fold Degenerate Sites Across Genes in the Great Tit Genome. Front Ecol Evol 2018. [DOI: 10.3389/fevo.2018.00203] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
24
|
Abstract
This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408–427 Ma and 399–414 Ma, respectively, with the interval of 9–13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous–Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.
Collapse
Affiliation(s)
- Naoko Takezaki
- Life Science Research Center, Kagawa University, Kitagun, Kagawa, Japan
| |
Collapse
|
25
|
Stange M, Sánchez-Villagra MR, Salzburger W, Matschiner M. Bayesian Divergence-Time Estimation with Genome-Wide Single-Nucleotide Polymorphism Data of Sea Catfishes (Ariidae) Supports Miocene Closure of the Panamanian Isthmus. Syst Biol 2018; 67:681-699. [PMID: 29385552 PMCID: PMC6005153 DOI: 10.1093/sysbio/syy006] [Citation(s) in RCA: 98] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Accepted: 01/23/2018] [Indexed: 02/03/2023] Open
Abstract
The closure of the Isthmus of Panama has long been considered to be one of the best defined biogeographic calibration points for molecular divergence-time estimation. However, geological and biological evidence has recently cast doubt on the presumed timing of the initial isthmus closure around 3 Ma but has instead suggested the existence of temporary land bridges as early as the Middle or Late Miocene. The biological evidence supporting these earlier land bridges was based either on only few molecular markers or on concatenation of genome-wide sequence data, an approach that is known to result in potentially misleading branch lengths and divergence times, which could compromise the reliability of this evidence. To allow divergence-time estimation with genomic data using the more appropriate multispecies coalescent (MSC) model, we here develop a new method combining the single-nucleotide polymorphism-based Bayesian species-tree inference of the software SNAPP with a molecular clock model that can be calibrated with fossil or biogeographic constraints. We validate our approach with simulations and use our method to reanalyze genomic data of Neotropical army ants (Dorylinae) that previously supported divergence times of Central and South American populations before the isthmus closure around 3 Ma. Our reanalysis with the MSC model shifts all of these divergence times to ages younger than 3 Ma, suggesting that the older estimates supporting the earlier existence of temporary land bridges were artifacts resulting at least partially from the use of concatenation. We then apply our method to a new restriction-site associated DNA-sequencing data set of Neotropical sea catfishes (Ariidae) and calibrate their species tree with extensive information from the fossil record. We identify a series of divergences between groups of Caribbean and Pacific sea catfishes around 10 Ma, indicating that processes related to the emergence of the isthmus led to vicariant speciation already in the Late Miocene, millions of years before the final isthmus closure.
Collapse
Affiliation(s)
- Madlen Stange
- Department of Palaeontology and Museum, University of Zurich, Karl-Schmid-Strasse 4, 8006 Zurich, Switzerland.,Zoological Institute, University of Basel, Vesalgasse 1, 4051 Basel, Switzerland
| | - Marcelo R Sánchez-Villagra
- Department of Palaeontology and Museum, University of Zurich, Karl-Schmid-Strasse 4, 8006 Zurich, Switzerland
| | - Walter Salzburger
- Zoological Institute, University of Basel, Vesalgasse 1, 4051 Basel, Switzerland.,Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, 0316 Oslo, Norway
| | - Michael Matschiner
- Zoological Institute, University of Basel, Vesalgasse 1, 4051 Basel, Switzerland.,Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, 0316 Oslo, Norway
| |
Collapse
|
26
|
Mongiardino Koch N, Gauthier JA. Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata. PLoS One 2018; 13:e0202729. [PMID: 30133514 PMCID: PMC6105018 DOI: 10.1371/journal.pone.0202729] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2018] [Accepted: 08/08/2018] [Indexed: 12/23/2022] Open
Abstract
Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.
Collapse
Affiliation(s)
- Nicolás Mongiardino Koch
- Department of Geology and Geophysics, Yale University, New Haven, Connecticut, United States of America
| | - Jacques A. Gauthier
- Department of Geology and Geophysics, Yale University, New Haven, Connecticut, United States of America
- Yale Peabody Museum of Natural History, New Haven, Connecticut, United States of America
| |
Collapse
|
27
|
Spalink D, Stoffel K, Walden GK, Hulse-Kemp AM, Hill TA, Van Deynze A, Bohs L. Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae). Mol Phylogenet Evol 2018; 126:293-302. [PMID: 29702214 DOI: 10.1016/j.ympev.2018.04.030] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2017] [Revised: 04/20/2018] [Accepted: 04/20/2018] [Indexed: 11/17/2022]
Abstract
The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.
Collapse
Affiliation(s)
- Daniel Spalink
- Department of Biology, University of Utah, Salt Lake City, UT, USA.
| | - Kevin Stoffel
- Department of Plant Sciences, University of California, Davis, CA, USA
| | - Genevieve K Walden
- Department of Biology, University of Utah, Salt Lake City, UT, USA; Plant Pest Diagnostics Center, California Department of Food and Agriculture, 3294 Meadowview Road, Sacramento, CA 95832-1448 USA
| | - Amanda M Hulse-Kemp
- Department of Plant Sciences, University of California, Davis, CA, USA; USDA-ARS, Genomics and Bioinformatics Research Unit, Raleigh, NC, USA
| | - Theresa A Hill
- Department of Plant Sciences, University of California, Davis, CA, USA
| | - Allen Van Deynze
- Department of Plant Sciences, University of California, Davis, CA, USA
| | - Lynn Bohs
- Department of Biology, University of Utah, Salt Lake City, UT, USA
| |
Collapse
|
28
|
Springer MS, Gatesy J. Delimiting Coalescence Genes (C-Genes) in Phylogenomic Data Sets. Genes (Basel) 2018; 9:genes9030123. [PMID: 29495400 PMCID: PMC5867844 DOI: 10.3390/genes9030123] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Revised: 02/02/2018] [Accepted: 02/19/2018] [Indexed: 02/07/2023] Open
Abstract
coalescence methods have emerged as a popular alternative for inferring species trees with large genomic datasets, because these methods explicitly account for incomplete lineage sorting. However, statistical consistency of summary coalescence methods is not guaranteed unless several model assumptions are true, including the critical assumption that recombination occurs freely among but not within coalescence genes (c-genes), which are the fundamental units of analysis for these methods. Each c-gene has a single branching history, and large sets of these independent gene histories should be the input for genome-scale coalescence estimates of phylogeny. By contrast, numerous studies have reported the results of coalescence analyses in which complete protein-coding sequences are treated as c-genes even though exons for these loci can span more than a megabase of DNA. Empirical estimates of recombination breakpoints suggest that c-genes may be much shorter, especially when large clades with many species are the focus of analysis. Although this idea has been challenged recently in the literature, the inverse relationship between c-gene size and increased taxon sampling in a dataset-the 'recombination ratchet'-is a fundamental property of c-genes. For taxonomic groups characterized by genes with long intron sequences, complete protein-coding sequences are likely not valid c-genes and are inappropriate units of analysis for summary coalescence methods unless they occur in recombination deserts that are devoid of incomplete lineage sorting (ILS). Finally, it has been argued that coalescence methods are robust when the no-recombination within loci assumption is violated, but recombination must matter at some scale because ILS, a by-product of recombination, is the raison d'etre for coalescence methods. That is, extensive recombination is required to yield the large number of independently segregating c-genes used to infer a species tree. If coalescent methods are powerful enough to infer the correct species tree for difficult phylogenetic problems in the anomaly zone, where concatenation is expected to fail because of ILS, then there should be a decreasing probability of inferring the correct species tree using longer loci with many intralocus recombination breakpoints (i.e., increased levels of concatenation).
Collapse
Affiliation(s)
- Mark S Springer
- Department of Evolution, Ecology, and Organismal Biology, University of California, Riverside, CA 92521, USA.
| | - John Gatesy
- Division of Vertebrate Zoology and Sackler Institute for Comparative Genomics, American Museum of Natural History, New York, NY 10024, USA.
| |
Collapse
|
29
|
|
30
|
Esselstyn JA, Oliveros CH, Swanson MT, Faircloth BC. Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements. Genome Biol Evol 2017; 9:2308-2321. [PMID: 28934378 PMCID: PMC5604124 DOI: 10.1093/gbe/evx168] [Citation(s) in RCA: 79] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/25/2017] [Indexed: 12/21/2022] Open
Abstract
The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.
Collapse
Affiliation(s)
- Jacob A. Esselstyn
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Carl H. Oliveros
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Mark T. Swanson
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Brant C. Faircloth
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| |
Collapse
|
31
|
Esselstyn JA, Oliveros CH, Swanson MT, Faircloth BC. Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements. Genome Biol Evol 2017. [PMID: 28934378 DOI: 10.1093/gbe/evx168)] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.
Collapse
Affiliation(s)
- Jacob A Esselstyn
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Carl H Oliveros
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Mark T Swanson
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Brant C Faircloth
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge
| |
Collapse
|