1
|
Ren H, Wong TKF, Minh BQ, Lanfear R. MixtureFinder: Estimating DNA Mixture Models for Phylogenetic Analyses. Mol Biol Evol 2025; 42:msae264. [PMID: 39715360 PMCID: PMC11704958 DOI: 10.1093/molbev/msae264] [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: 03/19/2024] [Revised: 11/26/2024] [Accepted: 12/19/2024] [Indexed: 12/25/2024] Open
Abstract
In phylogenetic studies, both partitioned models and mixture models are used to account for heterogeneity in molecular evolution among the sites of DNA sequence alignments. Partitioned models require the user to specify the grouping of sites into subsets, and then assume that each subset of sites can be modeled by a single common process. Mixture models do not require users to prespecify subsets of sites, and instead calculate the likelihood of every site under every model, while co-estimating the model weights and parameters. While much research has gone into the optimization of partitioned models by merging user-specified subsets, there has been less attention paid to the optimization of mixture models for DNA sequence alignments. In this study, we first ask whether a key assumption of partitioned models-that each user-specified subset can be modeled by a single common process-is supported by the data. Having shown that this is not the case, we then design, implement, test, and apply an algorithm, MixtureFinder, to select the optimum number of classes for a mixture model of Q-matrices for the standard models of DNA sequence evolution. We show this algorithm performs well on simulated and empirical datasets and suggest that it may be useful for future empirical studies. MixtureFinder is available in IQ-TREE2, and a tutorial for using MixtureFinder can be found here: http://www.iqtree.org/doc/Complex-Models#mixture-models.
Collapse
Affiliation(s)
- Huaiyan Ren
- School of Computing, College of Engineering, Computing and Cybernetics, Australian National University, Canberra, ACT 2600, Australia
- Ecology and Evolution, Research School of Biology, College of Science, Australian National University, Canberra, ACT 2600, Australia
| | - Thomas K F Wong
- School of Computing, College of Engineering, Computing and Cybernetics, Australian National University, Canberra, ACT 2600, Australia
- Ecology and Evolution, Research School of Biology, College of Science, Australian National University, Canberra, ACT 2600, Australia
| | - Bui Quang Minh
- School of Computing, College of Engineering, Computing and Cybernetics, Australian National University, Canberra, ACT 2600, Australia
| | - Robert Lanfear
- Ecology and Evolution, Research School of Biology, College of Science, Australian National University, Canberra, ACT 2600, Australia
| |
Collapse
|
2
|
Ferranti DA, Delwiche CF. Investigating the evolution of green algae with a large transcriptomic data set. JOURNAL OF PHYCOLOGY 2024; 60:1406-1419. [PMID: 39404089 DOI: 10.1111/jpy.13509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 08/29/2024] [Accepted: 09/04/2024] [Indexed: 12/28/2024]
Abstract
The colonization of land by plants approximately 450-500 million years ago (Mya) is one of the most important events in the history of life on Earth. Land plants, hereafter referred to as "embryophytes," comprise the foundation of every terrestrial biome, making them an essential lineage for the origin and maintenance of biodiversity. The embryophytes form a monophyletic clade within one of the two major phyla of the green algae (Viridiplantae), the Streptophyta. Estimates from fossil data and molecular clock analyses suggest the Streptophyte algae (Charophytes) diverged from the other main phylum of green algae, the Chlorophyta, as much as 1500 Mya. Here we present a phylogenetic analysis using transcriptomic and genomic data of 62 green algae and embryophyte operational taxonomic units, 31 of which were assembled de novo for this project. We have focused on identifying the charophyte lineage that is sister to embryophytes, and show that the Zygnematophyceae have the strongest support, followed by the Charophyceae. Furthermore, we have examined amino acid and codon usage across the tree and determined these data broadly follow the phylogenetic tree. We concluded by searching the data set for protein domains and gene families known to be important in embryophytes. Many of these domains and genes have homologous sequences in the charophyte lineages, giving insight into the processes that underlay the colonization of the land by plants. This provides new insights into green algal diversification, identifies previously unknown attributes of genome evolution within the group, and shows how functional mechanisms have evolved over time.
Collapse
Affiliation(s)
- David A Ferranti
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, Maryland, USA
| | - Charles F Delwiche
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, Maryland, USA
| |
Collapse
|
3
|
Schrago CG, Mello B. Challenges in Assembling the Dated Tree of Life. Genome Biol Evol 2024; 16:evae229. [PMID: 39475308 PMCID: PMC11523137 DOI: 10.1093/gbe/evae229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/15/2024] [Indexed: 11/02/2024] Open
Abstract
The assembly of a comprehensive and dated Tree of Life (ToL) remains one of the most formidable challenges in evolutionary biology. The complexity of life's history, involving both vertical and horizontal transmission of genetic information, defies its representation by a simple bifurcating phylogeny. With the advent of genome and metagenome sequencing, vast amounts of data have become available. However, employing this information for phylogeny and divergence time inference has introduced significant theoretical and computational hurdles. This perspective addresses some key methodological challenges in assembling the dated ToL, namely, the identification and classification of homologous genes, accounting for gene tree-species tree mismatch due to population-level processes along with duplication, loss, and horizontal gene transfer, and the accurate dating of evolutionary events. Ultimately, the success of this endeavor requires new approaches that integrate knowledge databases with optimized phylogenetic algorithms capable of managing complex evolutionary models.
Collapse
Affiliation(s)
- Carlos G Schrago
- Department of Genetics, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil
| | - Beatriz Mello
- Department of Genetics, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil
| |
Collapse
|
4
|
Mello B, Schrago CG. Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock. Genome Biol Evol 2024; 16:evae199. [PMID: 39332907 PMCID: PMC11430275 DOI: 10.1093/gbe/evae199] [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] [Accepted: 09/08/2024] [Indexed: 09/29/2024] Open
Abstract
Relaxing the molecular clock using models of how substitution rates change across lineages has become essential for addressing evolutionary problems. The diversity of rate evolution models and their implementations are substantial, and studies have demonstrated their impact on divergence time estimates can be as significant as that of calibration information. In this review, we trace the development of rate evolution models from the proposal of the molecular clock concept to the development of sophisticated Bayesian and non-Bayesian methods that handle rate variation in phylogenies. We discuss the various approaches to modeling rate evolution, provide a comprehensive list of available software, and examine the challenges and advancements of the prevalent Bayesian framework, contrasting them to faster non-Bayesian methods. Lastly, we offer insights into potential advancements in the field in the era of big data.
Collapse
Affiliation(s)
- Beatriz Mello
- Department of Genetics, Federal University of Rio de Janeiro, Rio de Janeiro, RJ 21941-617, Brazil
| | - Carlos G Schrago
- Department of Genetics, Federal University of Rio de Janeiro, Rio de Janeiro, RJ 21941-617, Brazil
| |
Collapse
|
5
|
Khurana MP, Scheidwasser-Clow N, Penn MJ, Bhatt S, Duchêne DA. The Limits of the Constant-rate Birth-Death Prior for Phylogenetic Tree Topology Inference. Syst Biol 2024; 73:235-246. [PMID: 38153910 PMCID: PMC11129600 DOI: 10.1093/sysbio/syad075] [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: 02/06/2023] [Revised: 12/20/2023] [Accepted: 12/27/2023] [Indexed: 12/30/2023] Open
Abstract
Birth-death models are stochastic processes describing speciation and extinction through time and across taxa and are widely used in biology for inference of evolutionary timescales. Previous research has highlighted how the expected trees under the constant-rate birth-death (crBD) model tend to differ from empirical trees, for example, with respect to the amount of phylogenetic imbalance. However, our understanding of how trees differ between the crBD model and the signal in empirical data remains incomplete. In this Point of View, we aim to expose the degree to which the crBD model differs from empirically inferred phylogenies and test the limits of the model in practice. Using a wide range of topology indices to compare crBD expectations against a comprehensive dataset of 1189 empirically estimated trees, we confirm that crBD model trees frequently differ topologically compared with empirical trees. To place this in the context of standard practice in the field, we conducted a meta-analysis for a subset of the empirical studies. When comparing studies that used Bayesian methods and crBD priors with those that used other non-crBD priors and non-Bayesian methods (i.e., maximum likelihood methods), we do not find any significant differences in tree topology inferences. To scrutinize this finding for the case of highly imbalanced trees, we selected the 100 trees with the greatest imbalance from our dataset, simulated sequence data for these tree topologies under various evolutionary rates, and re-inferred the trees under maximum likelihood and using the crBD model in a Bayesian setting. We find that when the substitution rate is low, the crBD prior results in overly balanced trees, but the tendency is negligible when substitution rates are sufficiently high. Overall, our findings demonstrate the general robustness of crBD priors across a broad range of phylogenetic inference scenarios but also highlight that empirically observed phylogenetic imbalance is highly improbable under the crBD model, leading to systematic bias in data sets with limited information content.
Collapse
Affiliation(s)
- Mark P Khurana
- Section of Epidemiology, Department of Public Health, University of Copenhagen, 1352 Copenhagen, Denmark
| | - Neil Scheidwasser-Clow
- Section of Epidemiology, Department of Public Health, University of Copenhagen, 1352 Copenhagen, Denmark
| | - Matthew J Penn
- Department of Statistics, University of Oxford, OX1 3LB, Oxford, UK
| | - Samir Bhatt
- Section of Epidemiology, Department of Public Health, University of Copenhagen, 1352 Copenhagen, Denmark
- MRC Centre for Global Infectious Disease Analysis, School of Public Health, Imperial College London, SW7 2AZ, London, UK
| | - David A Duchêne
- Centre for Evolutionary Hologenomics, University of Copenhagen, 1352 Copenhagen, Denmark
| |
Collapse
|
6
|
Chen LQ, Li X, Yao X, Li DZ, Barrett C, dePamphilis CW, Yu WB. Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae. PLANT MOLECULAR BIOLOGY 2024; 114:40. [PMID: 38622367 DOI: 10.1007/s11103-024-01440-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 03/09/2024] [Indexed: 04/17/2024]
Abstract
Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.
Collapse
Affiliation(s)
- Li-Qiong Chen
- Center for Integrative Conservation & Yunnan Key Laboratory for the Conservation of Tropical Rainforests and Asian Elephants, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla, Yunnan, 666303, China
| | - Xin Li
- Center for Integrative Conservation & Yunnan Key Laboratory for the Conservation of Tropical Rainforests and Asian Elephants, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla, Yunnan, 666303, China
- Division of BiologicalScience, Graduate School of Science and Technology, Nara Institute of Science and Technology, Ikoma, Nara, 630-0192, Japan
| | - Xin Yao
- Center for Integrative Conservation & Yunnan Key Laboratory for the Conservation of Tropical Rainforests and Asian Elephants, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla, Yunnan, 666303, China
| | - De-Zhu Li
- Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
| | - Craig Barrett
- Department of Biology, West Virginia University, Morgantown, West Virginia, 26506, USA
| | - Claude W dePamphilis
- Department of Biology, The Pennsylvania State University, University Park, State College, Pennsylvania, 16802, USA
| | - Wen-Bin Yu
- Center for Integrative Conservation & Yunnan Key Laboratory for the Conservation of Tropical Rainforests and Asian Elephants, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla, Yunnan, 666303, China.
- Southeast Asia Biodiversity Research Institute, Chinese Academy of Sciences, Mengla, Yunnan, 666303, China.
| |
Collapse
|
7
|
Kuo LY, Su HJ, Koubínová D, Xie PJ, Whitehouse C, Ebihara A, Grant JR. Organellar phylogenomics of Ophioglossaceae fern genera. FRONTIERS IN PLANT SCIENCE 2024; 14:1294716. [PMID: 38288414 PMCID: PMC10823028 DOI: 10.3389/fpls.2023.1294716] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Accepted: 12/27/2023] [Indexed: 01/31/2024]
Abstract
Previous phylogenies showed conflicting relationships among the subfamilies and genera within the fern family Ophioglossaceae. However, their classification remains unsettled where contrasting classifications recognize four to 15 genera. Since these treatments are mostly based on phylogenetic evidence using limited, plastid-only loci, a phylogenomic understanding is actually necessary to provide conclusive insight into the systematics of the genera. In this study, we have therefore compiled datasets with the broadest sampling of Ophioglossaceae genera to date, including all fifteen currently recognized genera, especially for the first time the South African endemic genus Rhizoglossum. Notably, our comprehensive phylogenomic matrix is based on both plastome and mitogenome genes. Inferred from the coding sequences of 83 plastid and 37 mitochondrial genes, a strongly supported topology for these subfamilies is presented, and is established by analyses using different partitioning approaches and substitution models. At the generic level, most relationships are well resolved except for few within the subfamily Ophioglossoideae. With this new phylogenomic scheme, key morphological and genomic changes were further identified along this backbone. In addition, we confirmed numerous horizontally transferred (HGT) genes in the genera Botrypus, Helminthostachys, Mankyua, Sahashia, and Sceptridium. These HGT genes are most likely located in mitogenomes and are predominately donated from angiosperm Santalales or non-Ophioglossaceae ferns. By our in-depth searches of the organellar genomes, we also provided phylogenetic overviews for the plastid and mitochondrial MORFFO genes found in these Ophioglossaceae ferns.
Collapse
Affiliation(s)
- Li-Yaung Kuo
- Institute of Molecular & Cellular Biology, National Tsing Hua University, Hsinchu, Taiwan
| | - Huei-Jiun Su
- Department of Earth and Life Sciences, University of Taipei, Taipei, Taiwan
| | - Darina Koubínová
- University of Neuchâtel, Laboratory of Evolutionary Genetics, Neuchâtel, Switzerland
| | - Pei-Jun Xie
- Institute of Molecular & Cellular Biology, National Tsing Hua University, Hsinchu, Taiwan
| | | | - Atsushi Ebihara
- Department of Botany, National Museum of Nature and Science, Tsukuba, Japan
| | - Jason R. Grant
- University of Neuchâtel, Laboratory of Evolutionary Genetics, Neuchâtel, Switzerland
| |
Collapse
|
8
|
Bowman J, Enard D, Lynch VJ. Phylogenomics reveals an almost perfect polytomy among the almost ungulates ( Paenungulata). BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.07.570590. [PMID: 38106080 PMCID: PMC10723481 DOI: 10.1101/2023.12.07.570590] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.
Collapse
Affiliation(s)
- Jacob Bowman
- Department of Biological Sciences, University at Buffalo, SUNY, 551 Cooke Hall, Buffalo, NY, USA
| | - David Enard
- Department of Ecology and Evolutionary Biology. University of Arizona, Tucson, AZ, USA
| | - Vincent J. Lynch
- Department of Biological Sciences, University at Buffalo, SUNY, 551 Cooke Hall, Buffalo, NY, USA
| |
Collapse
|
9
|
Jacques F, Bolivar P, Pietras K, Hammarlund EU. Roadmap to the study of gene and protein phylogeny and evolution-A practical guide. PLoS One 2023; 18:e0279597. [PMID: 36827278 PMCID: PMC9955684 DOI: 10.1371/journal.pone.0279597] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2022] [Accepted: 12/12/2022] [Indexed: 02/25/2023] Open
Abstract
Developments in sequencing technologies and the sequencing of an ever-increasing number of genomes have revolutionised studies of biodiversity and organismal evolution. This accumulation of data has been paralleled by the creation of numerous public biological databases through which the scientific community can mine the sequences and annotations of genomes, transcriptomes, and proteomes of multiple species. However, to find the appropriate databases and bioinformatic tools for respective inquiries and aims can be challenging. Here, we present a compilation of DNA and protein databases, as well as bioinformatic tools for phylogenetic reconstruction and a wide range of studies on molecular evolution. We provide a protocol for information extraction from biological databases and simple phylogenetic reconstruction using probabilistic and distance methods, facilitating the study of biodiversity and evolution at the molecular level for the broad scientific community.
Collapse
Affiliation(s)
- Florian Jacques
- Lund University Cancer Centre, Department of Laboratory Medicine, Lund University, Lund, Sweden
- Lund Stem Cell Center, Department of Laboratory Medicine, Lund University, Lund, Sweden
| | - Paulina Bolivar
- Lund University Cancer Centre, Department of Laboratory Medicine, Lund University, Lund, Sweden
| | - Kristian Pietras
- Lund University Cancer Centre, Department of Laboratory Medicine, Lund University, Lund, Sweden
| | - Emma U. Hammarlund
- Lund University Cancer Centre, Department of Laboratory Medicine, Lund University, Lund, Sweden
- Lund Stem Cell Center, Department of Laboratory Medicine, Lund University, Lund, Sweden
| |
Collapse
|
10
|
Giacomelli M, Rossi ME, Lozano-Fernandez J, Feuda R, Pisani D. Resolving tricky nodes in the tree of life through amino acid recoding. iScience 2022; 25:105594. [PMID: 36458253 PMCID: PMC9706708 DOI: 10.1016/j.isci.2022.105594] [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/17/2021] [Revised: 09/05/2022] [Accepted: 11/11/2022] [Indexed: 11/17/2022] Open
Abstract
Genomic data allowed a detailed resolution of the Tree of Life, but "tricky nodes" such as the root of the animals remain unresolved. Genome-scale datasets are heterogeneous as genes and species are exposed to different pressures, and this can negatively impacts phylogenetic accuracy. We use simulated genomic-scale datasets and show that recoding amino acid data improves accuracy when the model does not account for the compositional heterogeneity of the amino acid alignment. We apply our findings to three datasets addressing the root of the animal tree, where the debate centers on whether sponges (Porifera) or comb jellies (Ctenophora) represent the sister of all other animals. We show that results from empirical data follow predictions from simulations and suggest that, at the least in phylogenies inferred from amino acid sequences, a placement of the ctenophores as sister to all the other animals is best explained as a tree reconstruction artifact.
Collapse
Affiliation(s)
- Mattia Giacomelli
- Bristol Palaeobiology Group, School of Biological Sciences, University of Bristol, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ, UK
| | - Maria Eleonora Rossi
- Bristol Palaeobiology Group, School of Earth Sciences, University of Bristol, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
| | - Jesus Lozano-Fernandez
- Department of Genetics, Microbiology and Statistics, & Biodiversity Research Institute (IRBio), Faculty of Biology, University of Barcelona, Barcelona, Spain
| | - Roberto Feuda
- Department of Genetics and Genome Biology, University of Leicester, Leicester, UK
| | - Davide Pisani
- Bristol Palaeobiology Group, School of Biological Sciences, University of Bristol, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ, UK
- Bristol Palaeobiology Group, School of Earth Sciences, University of Bristol, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
| |
Collapse
|
11
|
Barrett CF, Santee MV, Fama NM, Freudenstein JV, Simon SJ, Sinn BT. Lineage and role in integrative taxonomy of a heterotrophic orchid complex. Mol Ecol 2022; 31:4762-4781. [PMID: 35837745 PMCID: PMC9452484 DOI: 10.1111/mec.16617] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2022] [Revised: 07/06/2022] [Accepted: 07/11/2022] [Indexed: 11/30/2022]
Abstract
Lineage-based species definitions applying coalescent approaches to species delimitation have become increasingly popular. Yet, the application of these methods and the recognition of lineage-only definitions have recently been questioned. Species delimitation criteria that explicitly consider both lineages and evidence for ecological role shifts provide an opportunity to incorporate ecologically meaningful data from multiple sources in studies of species boundaries. Here, such criteria were applied to a problematic group of mycoheterotrophic orchids, the Corallorhiza striata complex, analysing genomic, morphological, phenological, reproductive-mode, niche, and fungal host data. A recently developed method for generating genomic polymorphism data-ISSRseq-demonstrates evidence for four distinct lineages, including a previously unidentified lineage in the Coast Ranges and Cascades of California and Oregon, USA. There is divergence in morphology, phenology, reproductive mode, and fungal associates among the four lineages. Integrative analyses, conducted in population assignment and redundancy analysis frameworks, provide evidence of distinct genomic lineages and a similar pattern of divergence in the extended data, albeit with weaker signal. However, none of the extended data sets fully satisfy the condition of a significant role shift, which requires evidence of fixed differences. The four lineages identified in the current study are recognized at the level of variety, short of comprising different species. This study represents the most comprehensive application of lineage + role to date and illustrates the advantages of such an approach.
Collapse
Affiliation(s)
- Craig F. Barrett
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, West Virginia, USA 26506
| | - Mathilda V. Santee
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, West Virginia, USA 26506
| | - Nicole M. Fama
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, West Virginia, USA 26506
| | - John V. Freudenstein
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, 1315 Kinnear Rd., Columbus, Ohio, USA 43212
| | - Sandra J. Simon
- Department of Biology, West Virginia University Institute of Technology, Beckley, WV, USA
| | - Brandon T. Sinn
- Department of Biology and Earth Science, Otterbein University, Westerville, OH, USA
- Department of Botany and Ecology, University of Latvia, Jelgavas iela 1, Riga, LV-1004, Latvia
| |
Collapse
|
12
|
Lemarcis T, Fedosov AE, Kantor YI, Abdelkrim J, Zaharias P, Puillandre N. Neogastropod (Mollusca, Gastropoda) phylogeny: A step forward with mitogenomes. ZOOL SCR 2022; 51:550-561. [PMID: 36245672 PMCID: PMC9544082 DOI: 10.1111/zsc.12552] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Revised: 05/24/2022] [Accepted: 06/13/2022] [Indexed: 11/28/2022]
Abstract
The Neogastropoda (Mollusca, Gastropoda) encompass more than 15,000 described species of marine predators, including several model organisms in toxinology, embryology and physiology. However, their phylogenetic relationships remain mostly unresolved and their classification unstable. We took advantage of the many mitogenomes published in GenBank to produce a new molecular phylogeny of the neogastropods. We completed the taxon sampling by using an in-house bioinformatic pipeline to retrieve mitochondrial genes from 13 transcriptomes, corresponding to five families not represented in GenBank, for a final dataset of 113 taxa. Because mitogenomic data are prone to reconstruction artefacts, eight different evolutionary models were applied to reconstruct phylogenetic trees with IQTREE, RAxML and MrBayes. If the over-parametrization of some models produced trees with aberrant internal long branches, the global topology of the trees remained stable over models and softwares, and several relationships were revealed or found supported here for the first time. However, even if our dataset encompasses 60% of the valid families of neogastropods, some key taxa are missing and should be added in the future before proposing a revision of the classification of the neogastropods. Our study also demonstrates that even complex models struggle to satisfactorily handle the evolutionary history of mitogenomes, still leading to long-branch attractions in phylogenetic trees. Other approaches, such as reduced-genome strategies, must be envisaged to fully resolve the neogastropod phylogeny.
Collapse
Affiliation(s)
- Thomas Lemarcis
- Institut Systématique Evolution Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, EPHESorbonne Université, Université des AntillesParisFrance
| | - Alexander E. Fedosov
- Institut Systématique Evolution Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, EPHESorbonne Université, Université des AntillesParisFrance
- A. N. Severtsov Institute of Ecology and EvolutionRussian Academy of SciencesMoscowRussia
| | - Yuri I. Kantor
- Institut Systématique Evolution Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, EPHESorbonne Université, Université des AntillesParisFrance
- A. N. Severtsov Institute of Ecology and EvolutionRussian Academy of SciencesMoscowRussia
| | - Jawad Abdelkrim
- UAR 2700 'Acquisition et Analyse de Données pour l'Histoire naturelle', Service d'Analyse de Données, CNRSMuséum National d'Histoire Naturelle, Sorbonne UniversitésParis CedexFrance
| | - Paul Zaharias
- Institut Systématique Evolution Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, EPHESorbonne Université, Université des AntillesParisFrance
| | - Nicolas Puillandre
- Institut Systématique Evolution Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, EPHESorbonne Université, Université des AntillesParisFrance
| |
Collapse
|
13
|
Ge X, Zang H, Ye X, Peng L, Wang B, Lian G, Sun C. Comparative Mitogenomic Analyses of Hydropsychidae Revealing the Novel Rearrangement of Protein-Coding Gene and tRNA (Trichoptera: Annulipalpia). INSECTS 2022; 13:759. [PMID: 36135460 PMCID: PMC9501032 DOI: 10.3390/insects13090759] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/02/2022] [Revised: 08/17/2022] [Accepted: 08/20/2022] [Indexed: 06/16/2023]
Abstract
Gene rearrangement of the mitochondrial genome of insects, especially the rearrangement of protein-coding genes, has long been a hot topic for entomologists. Although mitochondrial gene rearrangement is common within Annulipalpia, protein-coding gene rearrangement is relatively rare. As the largest family in Annulipalpia, the available mitogenomes from Hydropsychidae Curtis, 1835 are scarce, and thus restrict our interpretation of the mitogenome characteristic. In this study, we obtained 19 novel mitogenomes of Hydropsychidae, of which the mitogenomes of the genus Arctopsyche are published for the first time. Coupled with published hydropsychid mitogenome, we analyzed the nucleotide composition evolutionary rates and gene rearrangements of the mitogenomes among subfamilies. As a result, we found two novel gene rearrangement patterns within Hydropsychidae, including rearrangement of protein-coding genes. Meanwhile, our results consider that the protein-coding gene arrangement of Potamyia can be interpreted by the tandem duplication/random loss (TDRL) model. In addition, the phylogenetic relationships within Hydropsychidae constructed by two strategies (Bayesian inference and maximum likelihood) strongly support the monophyly of Arctopscychinae, Diplectroninae, Hydropsychinae, and Macronematinae. Our study provides new insights into the mechanisms and patterns of mitogenome rearrangements in Hydropsychidae.
Collapse
Affiliation(s)
- Xinyu Ge
- Lab of Taxonomy & Aquatic Insects, Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
| | - Haoming Zang
- Lab of Taxonomy & Aquatic Insects, Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
| | - Xiaoyun Ye
- Environmental Monitoring Station of Qingtian County, Lishui 323999, China
| | - Lang Peng
- Lab of Taxonomy & Aquatic Insects, Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
| | - Beixin Wang
- Lab of Taxonomy & Aquatic Insects, Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
| | - Gang Lian
- Zhejiang Key Laboratory of Ecological and Environmental Monitoring, Forewarning and Quality Control, Zhejiang Province Ecological Environment Monitoring Centre, Hangzhou 310012, China
| | - Changhai Sun
- Lab of Taxonomy & Aquatic Insects, Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
| |
Collapse
|
14
|
Pelosi JA, Kim EH, Barbazuk WB, Sessa EB. Phylotranscriptomics Illuminates the Placement of Whole Genome Duplications and Gene Retention in Ferns. FRONTIERS IN PLANT SCIENCE 2022; 13:882441. [PMID: 35909764 PMCID: PMC9330400 DOI: 10.3389/fpls.2022.882441] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2022] [Accepted: 06/16/2022] [Indexed: 05/31/2023]
Abstract
Ferns are the second largest clade of vascular plants with over 10,000 species, yet the generation of genomic resources for the group has lagged behind other major clades of plants. Transcriptomic data have proven to be a powerful tool to assess phylogenetic relationships, using thousands of markers that are largely conserved across the genome, and without the need to sequence entire genomes. We assembled the largest nuclear phylogenetic dataset for ferns to date, including 2884 single-copy nuclear loci from 247 transcriptomes (242 ferns, five outgroups), and investigated phylogenetic relationships across the fern tree, the placement of whole genome duplications (WGDs), and gene retention patterns following WGDs. We generated a well-supported phylogeny of ferns and identified several regions of the fern phylogeny that demonstrate high levels of gene tree-species tree conflict, which largely correspond to areas of the phylogeny that have been difficult to resolve. Using a combination of approaches, we identified 27 WGDs across the phylogeny, including 18 large-scale events (involving more than one sampled taxon) and nine small-scale events (involving only one sampled taxon). Most inferred WGDs occur within single lineages (e.g., orders, families) rather than on the backbone of the phylogeny, although two inferred events are shared by leptosporangiate ferns (excluding Osmundales) and Polypodiales (excluding Lindsaeineae and Saccolomatineae), clades which correspond to the majority of fern diversity. We further examined how retained duplicates following WGDs compared across independent events and found that functions of retained genes were largely convergent, with processes involved in binding, responses to stimuli, and certain organelles over-represented in paralogs while processes involved in transport, organelles derived from endosymbiotic events, and signaling were under-represented. To date, our study is the most comprehensive investigation of the nuclear fern phylogeny, though several avenues for future research remain unexplored.
Collapse
Affiliation(s)
- Jessie A. Pelosi
- Department of Biology, University of Florida, Gainesville, FL, United States
| | - Emily H. Kim
- Department of Microbiology and Cell Science, University of Florida, Gainesville, FL, United States
| | - W. Brad Barbazuk
- Department of Biology, University of Florida, Gainesville, FL, United States
- Genetics Institute, University of Florida, Gainesville, FL, United States
| | - Emily B. Sessa
- Department of Biology, University of Florida, Gainesville, FL, United States
| |
Collapse
|
15
|
Du XY, Kuo LY, Zuo ZY, Li DZ, Lu JM. Structural Variation of Plastomes Provides Key Insight Into the Deep Phylogeny of Ferns. FRONTIERS IN PLANT SCIENCE 2022; 13:862772. [PMID: 35645990 PMCID: PMC9134734 DOI: 10.3389/fpls.2022.862772] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 04/04/2022] [Indexed: 06/02/2023]
Abstract
Structural variation of plastid genomes (plastomes), particularly large inversions and gene losses, can provide key evidence for the deep phylogeny of plants. In this study, we investigated the structural variation of fern plastomes in a phylogenetic context. A total of 127 plastomes representing all 50 recognized families and 11 orders of ferns were sampled, making it the most comprehensive plastomic analysis of fern lineages to date. The samples included 42 novel plastomes of 15 families with a focus on Hymenophyllales and Gleicheniales. We reconstructed a well-supported phylogeny of all extant fern families, detected significant structural synapomorphies, including 9 large inversions, 7 invert repeat region (IR) boundary shifts, 10 protein-coding gene losses, 7 tRNA gene losses or anticodon changes, and 19 codon indels (insertions or deletions) across the deep phylogeny of ferns, particularly on the backbone nodes. The newly identified inversion V5, together with the newly inferred expansion of the IR boundary R5, can be identified as a synapomorphy of a clade composed of Dipteridaceae, Matoniaceae, Schizaeales, and the core leptosporangiates, while a unique inversion V4, together with an expansion of the IR boundary R4, was verified as a synapomorphy of Gleicheniaceae. This structural evidence is in support of our phylogenetic inference, thus providing key insight into the paraphyly of Gleicheniales. The inversions of V5 and V7 together filled the crucial gap regarding how the "reversed" gene orientation in the IR region characterized by most extant ferns (Schizaeales and the core leptosporangiates) evolved from the inferred ancestral type as retained in Equisetales and Osmundales. The tRNA genes trnR-ACG and trnM-CAU were assumed to be relicts of the early-divergent fern lineages but intact in most Polypodiales, particularly in eupolypods; and the loss of the tRNA genes trnR-CCG, trnV-UAC, and trnR-UCU in fern plastomes was much more prevalent than previously thought. We also identified several codon indels in protein-coding genes within the core leptosporangiates, which may be identified as synapomorphies of specific families or higher ranks. This study provides an empirical case of integrating structural and sequence information of plastomes to resolve deep phylogeny of plants.
Collapse
Affiliation(s)
- Xin-Yu Du
- Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Li-Yaung Kuo
- Institute of Molecular and Cellular Biology, National Tsing Hua University, Hsinchu, Taiwan
| | - Zheng-Yu Zuo
- Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - De-Zhu Li
- Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Jin-Mei Lu
- Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| |
Collapse
|
16
|
Dong W, Li E, Liu Y, Xu C, Wang Y, Liu K, Cui X, Sun J, Suo Z, Zhang Z, Wen J, Zhou S. Phylogenomic approaches untangle early divergences and complex diversifications of the olive plant family. BMC Biol 2022; 20:92. [PMID: 35468824 PMCID: PMC9040247 DOI: 10.1186/s12915-022-01297-0] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Accepted: 04/13/2022] [Indexed: 12/27/2022] Open
Abstract
BACKGROUND Deep-branching phylogenetic relationships are often difficult to resolve because phylogenetic signals are obscured by the long history and complexity of evolutionary processes, such as ancient introgression/hybridization, polyploidization, and incomplete lineage sorting (ILS). Phylogenomics has been effective in providing information for resolving both deep- and shallow-scale relationships across all branches of the tree of life. The olive family (Oleaceae) is composed of 25 genera classified into five tribes with tribe Oleeae consisting of four subtribes. Previous phylogenetic analyses showed that ILS and/or hybridization led to phylogenetic incongruence in the family. It was essential to distinguish phylogenetic signal conflicts, and explore mechanisms for the uncertainties concerning relationships of the olive family, especially at the deep-branching nodes. RESULTS We used the whole plastid genome and nuclear single nucleotide polymorphism (SNP) data to infer the phylogenetic relationships and to assess the variation and rates among the main clades of the olive family. We also used 2608 and 1865 orthologous nuclear genes to infer the deep-branching relationships among tribes of Oleaceae and subtribes of tribe Oleeae, respectively. Concatenated and coalescence trees based on the plastid genome, nuclear SNPs and multiple nuclear genes suggest events of ILS and/or ancient introgression during the diversification of Oleaceae. Additionally, there was extreme heterogeneity in the substitution rates across the tribes. Furthermore, our results supported that introgression/hybridization, rather than ILS, is the main factor for phylogenetic discordance among the five tribes of Oleaceae. The tribe Oleeae is supported to have originated via ancient hybridization and polyploidy, and its most likely parentages are the ancestral lineage of Jasmineae or its sister group, which is a "ghost lineage," and Forsythieae. However, ILS and ancient introgression are mainly responsible for the phylogenetic discordance among the four subtribes of tribe Oleeae. CONCLUSIONS This study showcases that using multiple sequence datasets (plastid genomes, nuclear SNPs and thousands of nuclear genes) and diverse phylogenomic methods such as data partition, heterogeneous models, quantifying introgression via branch lengths (QuIBL) analysis, and species network analysis can facilitate untangling long and complex evolutionary processes of ancient introgression, paleopolyploidization, and ILS.
Collapse
Affiliation(s)
- Wenpan Dong
- Laboratory of Systematic Evolution and Biogeography of Woody Plants, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, 100083, China.
| | - Enze Li
- Laboratory of Systematic Evolution and Biogeography of Woody Plants, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, 100083, China
| | - Yanlei Liu
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China
| | - Chao Xu
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China
| | - Yushuang Wang
- Laboratory of Systematic Evolution and Biogeography of Woody Plants, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, 100083, China
| | - Kangjia Liu
- Laboratory of Systematic Evolution and Biogeography of Woody Plants, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, 100083, China
| | - Xingyong Cui
- Laboratory of Systematic Evolution and Biogeography of Woody Plants, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, 100083, China
| | - Jiahui Sun
- State Key Laboratory Breeding Base of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, 100700, China.
| | - Zhili Suo
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China
| | - Zhixiang Zhang
- Laboratory of Systematic Evolution and Biogeography of Woody Plants, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, 100083, China
| | - Jun Wen
- Department of Botany, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20013-7012, USA.
| | - Shiliang Zhou
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China
| |
Collapse
|
17
|
Krug PJ, Caplins SA, Algoso K, Thomas K, Valdés ÁA, Wade R, Wong NLWS, Eernisse DJ, Kocot KM. Phylogenomic resolution of the root of Panpulmonata, a hyperdiverse radiation of gastropods: new insight into the evolution of air breathing. Proc Biol Sci 2022; 289:20211855. [PMID: 35382597 PMCID: PMC8984808 DOI: 10.1098/rspb.2021.1855] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2021] [Accepted: 02/21/2022] [Indexed: 11/12/2022] Open
Abstract
Transitions to terrestriality have been associated with major animal radiations including land snails and slugs in Stylommatophora (>20 000 described species), the most successful lineage of 'pulmonates' (a non-monophyletic assemblage of air-breathing gastropods). However, phylogenomic studies have failed to robustly resolve relationships among traditional pulmonates and affiliated marine lineages that comprise clade Panpulmonata (Mollusca, Gastropoda), especially two key taxa: Sacoglossa, a group including photosynthetic sea slugs, and Siphonarioidea, intertidal limpet-like snails with a non-contractile pneumostome (narrow opening to a vascularized pallial cavity). To clarify the evolutionary history of the panpulmonate radiation, we performed phylogenomic analyses on datasets of up to 1160 nuclear protein-coding genes for 110 gastropods, including 40 new transcriptomes for Sacoglossa and Siphonarioidea. All 18 analyses recovered Sacoglossa as the sister group to a clade we named Pneumopulmonata, within which Siphonarioidea was sister to the remaining lineages in most analyses. Comparative modelling indicated shifts to marginal habitat (estuarine, mangrove and intertidal zones) preceded and accelerated the evolution of a pneumostome, present in the pneumopulmonate ancestor along with a one-sided plicate gill. These findings highlight key intermediate stages in the evolution of air-breathing snails, supporting the hypothesis that adaptation to marginal zones played an important role in major sea-to-land transitions.
Collapse
Affiliation(s)
- Patrick J. Krug
- Department of Biological Sciences, California State University, Los Angeles, CA 90032-8201, USA
| | | | - Krisha Algoso
- Department of Biological Sciences, California State University, Los Angeles, CA 90032-8201, USA
| | - Kanique Thomas
- Department of Biological Sciences, California State University, Los Angeles, CA 90032-8201, USA
| | - Ángel A. Valdés
- Department of Biological Sciences, California State Polytechnic University, Pomona, CA 91768, USA
| | - Rachael Wade
- Department of Botany, University of British Columbia, Vancouver, BC, Canada V6T 1Z4
| | - Nur Leena W. S. Wong
- International Institute of Aquaculture and Aquatic Sciences, Universiti Putra Malaysia, Selangor, Malaysia
| | - Douglas J. Eernisse
- Department of Biological Science, California State University, Fullerton, CA 92834, USA
| | - Kevin M. Kocot
- Department of Biological Sciences and Alabama Museum of Natural History, The University of Alabama, Tuscaloosa, AL 35487, USA
| |
Collapse
|
18
|
Hou Z, Ma X, Shi X, Li X, Yang L, Xiao S, De Clerck O, Leliaert F, Zhong B. Phylotranscriptomic insights into a Mesoproterozoic-Neoproterozoic origin and early radiation of green seaweeds (Ulvophyceae). Nat Commun 2022; 13:1610. [PMID: 35318329 PMCID: PMC8941102 DOI: 10.1038/s41467-022-29282-9] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Accepted: 03/09/2022] [Indexed: 01/09/2023] Open
Abstract
The Ulvophyceae, a major group of green algae, is of particular evolutionary interest because of its remarkable morphological and ecological diversity. Its phylogenetic relationships and diversification timeline, however, are still not fully resolved. In this study, using an extensive nuclear gene dataset, we apply coalescent- and concatenation-based approaches to reconstruct the phylogeny of the Ulvophyceae and to explore the sources of conflict in previous phylogenomic studies. The Ulvophyceae is recovered as a paraphyletic group, with the Bryopsidales being a sister group to the Chlorophyceae, and the remaining taxa forming a clade (Ulvophyceae sensu stricto). Molecular clock analyses with different calibration strategies emphasize the large impact of fossil calibrations, and indicate a Meso-Neoproterozoic origin of the Ulvophyceae (sensu stricto), earlier than previous estimates. The results imply that ulvophyceans may have had a profound influence on oceanic redox structures and global biogeochemical cycles at the Mesoproterozoic-Neoproterozoic transition. “Ulvophyceae is a remarkably morphologically and ecologically diverse clade of green algae. Here, the authors reconstruct the Ulvophyceae phylogeny, showing that these algae originated earlier than expected and may have influenced biogeochemical cycles at the Mesoproterozoic-Neoproterozoic transition.”
Collapse
Affiliation(s)
- Zheng Hou
- College of Life Sciences, Nanjing Normal University, Nanjing, China
| | - Xiaoya Ma
- College of Life Sciences, Nanjing Normal University, Nanjing, China
| | - Xuan Shi
- College of Life Sciences, Nanjing Normal University, Nanjing, China
| | - Xi Li
- College of Life Sciences, Nanjing Normal University, Nanjing, China
| | - Lingxiao Yang
- College of Life Sciences, Nanjing Normal University, Nanjing, China
| | - Shuhai Xiao
- Department of Geosciences and Global Change Center, Virginia Tech, Blacksburg, VA, USA
| | - Olivier De Clerck
- Phycology Research Group and Center for Molecular Phylogenetics and Evolution, Ghent University, Ghent, Belgium
| | - Frederik Leliaert
- Phycology Research Group and Center for Molecular Phylogenetics and Evolution, Ghent University, Ghent, Belgium.,Meise Botanic Garden, Meise, Belgium
| | - Bojian Zhong
- College of Life Sciences, Nanjing Normal University, Nanjing, China.
| |
Collapse
|
19
|
Samigullin T, Logacheva M, Terentieva E, Degtjareva G, Pimenov M, Valiejo-Roman C. Plastid Phylogenomic Analysis of Tordylieae Tribe (Apiaceae, Apioideae). PLANTS (BASEL, SWITZERLAND) 2022; 11:709. [PMID: 35270181 PMCID: PMC8912408 DOI: 10.3390/plants11050709] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 02/28/2022] [Accepted: 03/01/2022] [Indexed: 11/16/2022]
Abstract
Based on the nrDNA ITS sequence data, the Tordylieae tribe is recognized as monophyletic with three major lineages: the subtribe Tordyliinae, the Cymbocarpum clade, and the Lefebvrea clade. Recent phylogenomic investigations showed incongruence between the nuclear and plastid genome evolution in the tribe. To assess phylogenetic relations and structure evolution of plastomes in Tordylieae, we generated eleven complete plastome sequences using the genome skimming approach and compared them with the available data from this tribe and close relatives. Newly assembled plastomes had lengths ranging from 141,148 to 150,103 base pairs and contained 122-127 genes, including 79-82 protein-coding genes, 35-37 tRNAs, and 8 rRNAs. We observed substantial differences in the inverted repeat length and gene content, accompanied by a complex picture of multiple JLA and JLB shifts. In concatenated phylogenetic analyses, Tordylieae plastomes formed at least three not closely related lineages with plastomes of the Lefebvrea clade as a sister group to plastomes from the Selineae tribe. The newly obtained data have increased our knowledge on the range of plastome variability in Apiaceae.
Collapse
Affiliation(s)
- Tahir Samigullin
- Department of Evolutionary Biochemistry, A. N. Belozersky Research Institute of Physicochemical Biology, Lomonosov Moscow State University, Leninskie Gory 1–40, 119992 Moscow, Russia; (M.L.); (C.V.-R.)
| | - Maria Logacheva
- Department of Evolutionary Biochemistry, A. N. Belozersky Research Institute of Physicochemical Biology, Lomonosov Moscow State University, Leninskie Gory 1–40, 119992 Moscow, Russia; (M.L.); (C.V.-R.)
- Center of Life Sciences, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, Bld. 1, 121205 Moscow, Russia
| | - Elena Terentieva
- Botanical Garden, Faculty of Biology, Lomonosov Moscow State University, Leninskie Gory 1/12, 119992 Moscow, Russia; (E.T.); (G.D.); (M.P.)
| | - Galina Degtjareva
- Botanical Garden, Faculty of Biology, Lomonosov Moscow State University, Leninskie Gory 1/12, 119992 Moscow, Russia; (E.T.); (G.D.); (M.P.)
| | - Michael Pimenov
- Botanical Garden, Faculty of Biology, Lomonosov Moscow State University, Leninskie Gory 1/12, 119992 Moscow, Russia; (E.T.); (G.D.); (M.P.)
| | - Carmen Valiejo-Roman
- Department of Evolutionary Biochemistry, A. N. Belozersky Research Institute of Physicochemical Biology, Lomonosov Moscow State University, Leninskie Gory 1–40, 119992 Moscow, Russia; (M.L.); (C.V.-R.)
| |
Collapse
|
20
|
Zheng B, Han Y, Yuan R, Liu J, van Achterberg C, Tang P, Chen X. Comparative Mitochondrial Genomics of 104 Darwin Wasps (Hymenoptera: Ichneumonidae) and Its Implication for Phylogeny. INSECTS 2022; 13:insects13020124. [PMID: 35206698 PMCID: PMC8874624 DOI: 10.3390/insects13020124] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/26/2021] [Revised: 01/17/2022] [Accepted: 01/21/2022] [Indexed: 02/04/2023]
Abstract
Simple Summary Nearly a hundred mitochondrial genomes of ichneumonid wasps are newly reported. Comparative mitogenomics of 104 mitochondrial genomes representing 33 subfamilies of Ichneumonidae, as well as its implications for phylogeny, were studied. We found that the mitochondrial genomes of ichneumonid wasps were highly conserved in their base composition and had low evolutionary rates, but were diverse in gene order. There are 38 types of gene rearrangement events in 104 ichneumonid mitochondrial genomes, of which 30 novel rearrangement types (R3-6, R8-R10, R12-R15, R17-R18, R20-R35 and R38) and a hot spot rearrangement around R1, with a shuffled tRNA cluster trnW-trnY-trnC and trnI-trnQ-trnM, were detected. The relationships among these subfamilies are firstly discussed based on mitochondrial genomes at a large scale. We suggest five subfamily groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively. Abstract Ichneumonidae is one of the largest families of insects with a mega-diversity of specialized morphological and biological characteristics. We newly sequenced 92 mitochondrial genomes of ichneumonid wasps and found that they have a conserved base composition and a lower evolutionary rate than that of other families of parasitic Hymenoptera. There are 38 types of gene order in the ichneumonid mitochondrial genome, with 30 novel types identified in 104 ichneumonids. We also found that the rearrangement events occur more frequently in Ophioniformes than in Ichneumoniformes and Pimpliformes. Furthermore, the higher Ophioniformes and their relative lineages shared the transposition of trnL2 to trnI-trnQ-trnM tRNA cluster. We confirmed five higher-level groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively. The results will improve our understanding of the diversity and evolution of Ichneumonidae.
Collapse
Affiliation(s)
- Boying Zheng
- State Key Laboratory of Rice Biology, Zhejiang University, Hangzhou 310058, China; (B.Z.); (Y.H.); (R.Y.); (C.v.A.); (X.C.)
- Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
| | - Yuanyuan Han
- State Key Laboratory of Rice Biology, Zhejiang University, Hangzhou 310058, China; (B.Z.); (Y.H.); (R.Y.); (C.v.A.); (X.C.)
- Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
| | - Ruizhong Yuan
- State Key Laboratory of Rice Biology, Zhejiang University, Hangzhou 310058, China; (B.Z.); (Y.H.); (R.Y.); (C.v.A.); (X.C.)
- Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
| | - Jingxian Liu
- Department of Entomology, College of Plant Protection, South China Agricultural University, Guangzhou 510642, China;
| | - Cornelis van Achterberg
- State Key Laboratory of Rice Biology, Zhejiang University, Hangzhou 310058, China; (B.Z.); (Y.H.); (R.Y.); (C.v.A.); (X.C.)
- Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
| | - Pu Tang
- State Key Laboratory of Rice Biology, Zhejiang University, Hangzhou 310058, China; (B.Z.); (Y.H.); (R.Y.); (C.v.A.); (X.C.)
- Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
- Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou 310058, China
- Zhejiang Provincial Key Laboratory of Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou 310058, China
- Correspondence:
| | - Xuexin Chen
- State Key Laboratory of Rice Biology, Zhejiang University, Hangzhou 310058, China; (B.Z.); (Y.H.); (R.Y.); (C.v.A.); (X.C.)
- Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
- Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou 310058, China
- Zhejiang Provincial Key Laboratory of Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
21
|
Baker CM, Buckman-Young RS, Costa CS, Giribet G. Phylogenomic Analysis of Velvet Worms (Onychophora) Uncovers an Evolutionary Radiation in the Neotropics. Mol Biol Evol 2021; 38:5391-5404. [PMID: 34427671 PMCID: PMC8662635 DOI: 10.1093/molbev/msab251] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
Onychophora ("velvet worms") are charismatic soil invertebrates known for their status as a "living fossil," their phylogenetic affiliation to arthropods, and their distinctive biogeographic patterns. However, several aspects of their internal phylogenetic relationships remain unresolved, limiting our understanding of the group's evolutionary history, particularly with regard to changes in reproductive mode and dispersal ability. To address these gaps, we used RNA sequencing and phylogenomic analysis of transcriptomes to reconstruct the evolutionary relationships and infer divergence times within the phylum. We recovered a fully resolved and well-supported phylogeny for the circum-Antarctic family Peripatopsidae, which retains signals of Gondwanan vicariance and showcases the evolutionary lability of reproductive mode in the family. Within the Neotropical clade of Peripatidae, though, we found that amino acid-translated sequence data masked nearly all phylogenetic signal, resulting in highly unstable and poorly supported relationships. Analyses using nucleotide sequence data were able to resolve many more relationships, though we still saw discordant phylogenetic signal between genes, probably indicative of a rapid, mid-Cretaceous radiation in the group. Finally, we hypothesize that the unique reproductive mode of placentotrophic viviparity found in all Neotropical peripatids may have facilitated the multiple inferred instances of over-water dispersal and establishment on oceanic islands.
Collapse
Affiliation(s)
- Caitlin M Baker
- Museum of Comparative Zoology, Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
| | - Rebecca S Buckman-Young
- Museum of Comparative Zoology, Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
| | - Cristiano S Costa
- Laboratório de Sistemática e Taxonomia de Artrópodes Terrestres, Departamento de Biologia e Zoologia, Instituto de Biociências, Universidade Federal de Mato Grosso, Cuiabá, Brazil
| | - Gonzalo Giribet
- Museum of Comparative Zoology, Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
| |
Collapse
|
22
|
Benavides LR, Jiang C, Giribet G. Mimopidae is the sister group to all other scolopendromorph centipedes (Chilopoda, Scolopendromorpha): a phylotranscriptomic approach. ORG DIVERS EVOL 2021. [DOI: 10.1007/s13127-021-00502-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
23
|
Abdullah, Mehmood F, Rahim A, Heidari P, Ahmed I, Poczai P. Comparative plastome analysis of Blumea, with implications for genome evolution and phylogeny of Asteroideae. Ecol Evol 2021; 11:7810-7826. [PMID: 34188853 PMCID: PMC8216946 DOI: 10.1002/ece3.7614] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Accepted: 04/09/2021] [Indexed: 12/25/2022] Open
Abstract
The genus Blumea (Asteroideae, Asteraceae) comprises about 100 species, including herbs, shrubs, and small trees. Previous studies have been unable to resolve taxonomic issues and the phylogeny of the genus Blumea due to the low polymorphism of molecular markers. Therefore, suitable polymorphic regions need to be identified. Here, we de novo assembled plastomes of the three Blumea species B. oxyodonta, B. tenella, and B. balsamifera and compared them with 26 other species of Asteroideae after correction of annotations. These species have quadripartite plastomes with similar gene content, genome organization, and inverted repeat contraction and expansion comprising 113 genes, including 80 protein-coding, 29 transfer RNA, and 4 ribosomal RNA genes. The comparative analysis of codon usage, amino acid frequency, microsatellite repeats, oligonucleotide repeats, and transition and transversion substitutions has revealed high resemblance among the newly assembled species of Blumea. We identified 10 highly polymorphic regions with nucleotide diversity above 0.02, including rps16-trnQ, ycf1, ndhF-rpl32, petN-psbM, and rpl32-trnL, and they may be suitable for the development of robust, authentic, and cost-effective markers for barcoding and inference of the phylogeny of the genus Blumea. Among these highly polymorphic regions, five regions also co-occurred with oligonucleotide repeats and support use of repeats as a proxy for the identification of polymorphic loci. The phylogenetic analysis revealed a close relationship between Blumea and Pluchea within the tribe Inuleae. At tribe level, our phylogeny supports a sister relationship between Astereae and Anthemideae rooted as Gnaphalieae, Calenduleae, and Senecioneae. These results are contradictory to recent studies which reported a sister relationship between "Senecioneae and Anthemideae" and "Astereae and Gnaphalieae" or a sister relationship between Astereae and Gnaphalieae rooted as Calenduleae, Anthemideae, and then Senecioneae using nuclear genome sequences. The conflicting phylogenetic signals observed at the tribal level between plastidt and nuclear genome data require further investigation.
Collapse
Affiliation(s)
- Abdullah
- Department of BiochemistryFaculty of Biological SciencesQuaid‐i‐Azam UniversityIslamabadPakistan
| | - Furrukh Mehmood
- Department of BiochemistryFaculty of Biological SciencesQuaid‐i‐Azam UniversityIslamabadPakistan
| | - Abdur Rahim
- Department of ZoologyPost Graduate College NowsheraAbdul Wali Khan UniversityMardanPakistan
| | - Parviz Heidari
- Faculty of AgricultureShahrood University of TechnologyShahroodIran
| | - Ibrar Ahmed
- Alpha Genomics Private LimitedIslamabadPakistan
| | - Péter Poczai
- Finnish Museum of Natural HistoryUniversity of HelsinkiHelsinkiFinland
- Faculty of Biological and Environmental SciencesUniversity of HelsinkiHelsinkiFinland
| |
Collapse
|
24
|
Li X, Hou Z, Xu C, Shi X, Yang L, Lewis LA, Zhong B. Large Phylogenomic Data sets Reveal Deep Relationships and Trait Evolution in Chlorophyte Green Algae. Genome Biol Evol 2021; 13:6265471. [PMID: 33950183 PMCID: PMC8271138 DOI: 10.1093/gbe/evab101] [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] [Accepted: 05/04/2021] [Indexed: 12/01/2022] Open
Abstract
The chlorophyte green algae (Chlorophyta) are species-rich ancient groups ubiquitous in various habitats with high cytological diversity, ranging from microscopic to macroscopic organisms. However, the deep phylogeny within core Chlorophyta remains unresolved, in part due to the relatively sparse taxon and gene sampling in previous studies. Here we contribute new transcriptomic data and reconstruct phylogenetic relationships of core Chlorophyta based on four large data sets up to 2,698 genes of 70 species, representing 80% of extant orders. The impacts of outgroup choice, missing data, bootstrap-support cutoffs, and model misspecification in phylogenetic inference of core Chlorophyta are examined. The species tree topologies of core Chlorophyta from different analyses are highly congruent, with strong supports at many relationships (e.g., the Bryopsidales and the Scotinosphaerales-Dasycladales clade). The monophyly of Chlorophyceae and of Trebouxiophyceae as well as the uncertain placement of Chlorodendrophyceae and Pedinophyceae corroborate results from previous studies. The reconstruction of ancestral scenarios illustrates the evolution of the freshwater-sea and microscopic–macroscopic transition in the Ulvophyceae, and the transformation of unicellular→colonial→multicellular in the chlorophyte green algae. In addition, we provided new evidence that serine is encoded by both canonical codons and noncanonical TAG code in Scotinosphaerales, and stop-to-sense codon reassignment in the Ulvophyceae has originated independently at least three times. Our robust phylogenetic framework of core Chlorophyta unveils the evolutionary history of phycoplast, cyto-morphology, and noncanonical genetic codes in chlorophyte green algae.
Collapse
Affiliation(s)
- Xi Li
- College of Life Sciences, Nanjing Normal University, China
| | - Zheng Hou
- College of Life Sciences, Nanjing Normal University, China
| | - Chenjie Xu
- College of Life Sciences, Nanjing Normal University, China
| | - Xuan Shi
- College of Life Sciences, Nanjing Normal University, China
| | - Lingxiao Yang
- College of Life Sciences, Nanjing Normal University, China
| | - Louise A Lewis
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, Connecticut, USA
| | - Bojian Zhong
- College of Life Sciences, Nanjing Normal University, China
| |
Collapse
|
25
|
Abstract
Multiple sequence alignment is a core first step in many bioinformatics analyses, and errors in these alignments can have negative consequences for scientific studies. In this article, we review some of the recent literature evaluating multiple sequence alignment methods and identify specific challenges that arise when performing these evaluations. In particular, we discuss the different trends observed in simulation studies and when using biological benchmarks. Overall, we find that multiple sequence alignment, far from being a "solved problem," would benefit from new attention.
Collapse
Affiliation(s)
- Tandy Warnow
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
26
|
Williams TA, Schrempf D, Szöllősi GJ, Cox CJ, Foster PG, Embley TM. Inferring the deep past from molecular data. Genome Biol Evol 2021; 13:6192802. [PMID: 33772552 PMCID: PMC8175050 DOI: 10.1093/gbe/evab067] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/22/2021] [Indexed: 12/17/2022] Open
Abstract
There is an expectation that analyses of molecular sequences might be able to distinguish between alternative hypotheses for ancient relationships, but the phylogenetic methods used and types of data analyzed are of critical importance in any attempt to recover historical signal. Here, we discuss some common issues that can influence the topology of trees obtained when using overly simple models to analyze molecular data that often display complicated patterns of sequence heterogeneity. To illustrate our discussion, we have used three examples of inferred relationships which have changed radically as models and methods of analysis have improved. In two of these examples, the sister-group relationship between thermophilic Thermus and mesophilic Deinococcus, and the position of long-branch Microsporidia among eukaryotes, we show that recovering what is now generally considered to be the correct tree is critically dependent on the fit between model and data. In the third example, the position of eukaryotes in the tree of life, the hypothesis that is currently supported by the best available methods is fundamentally different from the classical view of relationships between major cellular domains. Since heterogeneity appears to be pervasive and varied among all molecular sequence data, and even the best available models can still struggle to deal with some problems, the issues we discuss are generally relevant to phylogenetic analyses. It remains essential to maintain a critical attitude to all trees as hypotheses of relationship that may change with more data and better methods.
Collapse
Affiliation(s)
- Tom A Williams
- School of Biological Sciences, University of Bristol, Bristol BS8 1TQ, United Kingdom
| | - Dominik Schrempf
- Dept. of Biological Physics, Eötvös Loránd University, 1117 Budapest, Hungary
| | - Gergely J Szöllősi
- Dept. of Biological Physics, Eötvös Loránd University, 1117 Budapest, Hungary.,MTA-ELTE "Lendület" Evolutionary Genomics Research Group, 1117 Budapest, Hungary.,Institute of Evolution, Centre for Ecological Research, 1121 Budapest, Hungary
| | - Cymon J Cox
- Centro de Ciências do Mar, Universidade do Algarve, Gambelas, 8005-319 Faro, Portugal
| | - Peter G Foster
- Department of Life Sciences, Natural History Museum, London SW7 5BD, United Kingdom
| | - T Martin Embley
- Biosciences Institute, Centre for Bacterial Cell Biology, Newcastle University, Newcastle upon Tyne NE2 4AX, United Kingdom
| |
Collapse
|
27
|
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol Biol Evol 2021; 37:1530-1534. [PMID: 32011700 PMCID: PMC7182206 DOI: 10.1093/molbev/msaa015] [Citation(s) in RCA: 5912] [Impact Index Per Article: 1478.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
IQ-TREE (http://www.iqtree.org, last accessed February 6, 2020) is a user-friendly and widely used software package for phylogenetic inference using maximum likelihood. Since the release of version 1 in 2014, we have continuously expanded IQ-TREE to integrate a plethora of new models of sequence evolution and efficient computational approaches of phylogenetic inference to deal with genomic data. Here, we describe notable features of IQ-TREE version 2 and highlight the key advantages over other software.
Collapse
Affiliation(s)
- Bui Quang Minh
- Research School of Computer Science, Australian National University, Canberra, ACT, Australia.,Department of Ecology and Evolution, Research School of Biology, Australian National University, Canberra, ACT, Australia
| | - Heiko A Schmidt
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Olga Chernomor
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Dominik Schrempf
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria.,Department of Biological Physics, Eötvös Lórand University, Budapest, Hungary
| | - Michael D Woodhams
- Discipline of Mathematics, University of Tasmania, Hobart, TAS, Australia
| | - Arndt von Haeseler
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria.,Bioinformatics and Computational Biology, Faculty of Computer Science, University of Vienna, Vienna, Austria
| | - Robert Lanfear
- Department of Ecology and Evolution, Research School of Biology, Australian National University, Canberra, ACT, Australia
| |
Collapse
|
28
|
Papetti C, Babbucci M, Dettai A, Basso A, Lucassen M, Harms L, Bonillo C, Heindler FM, Patarnello T, Negrisolo E. Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish. Genome Biol Evol 2021; 13:6133229. [PMID: 33570582 PMCID: PMC7936035 DOI: 10.1093/gbe/evab017] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/22/2021] [Indexed: 12/21/2022] Open
Abstract
The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.
Collapse
Affiliation(s)
- Chiara Papetti
- Department of Biology, University of Padova, Padova 35121,Italy.,Consorzio Nazionale Interuniversitario per le Scienze del Mare (CoNISMa), Roma 00196, Italy
| | - Massimiliano Babbucci
- Department of Comparative Biomedicine and Food Science, University of Padova, Legnaro 35020, Italy
| | - Agnes Dettai
- Institut de Systematique, Evolution, Biodiversité (ISYEB) Muséum national d'Histoire naturelle-CNRS-Sorbonne Université-EPHE, MNHN, Paris 75005, France
| | - Andrea Basso
- Department of Comparative Biomedicine and Food Science, University of Padova, Legnaro 35020, Italy
| | - Magnus Lucassen
- Helmholtz Centre for Polar and Marine Research, Alfred Wegener Institute, Am Handelshafen 12, Bremerhaven 27570, Germany
| | - Lars Harms
- Helmholtz Centre for Polar and Marine Research, Alfred Wegener Institute, Am Handelshafen 12, Bremerhaven 27570, Germany.,Helmholtz Institute for Functional Marine Biodiversity, University of Oldenburg (HIFMB), Ammerlsity of Oldenburg (HIFMOldenburg 26129, Germany
| | - Celine Bonillo
- Service de Systématique Moléculaire, UMS2700 Acquisition et Analyse de Données (2AD), MNHN, Paris 75005, France
| | | | - Tomaso Patarnello
- Department of Comparative Biomedicine and Food Science, University of Padova, Legnaro 35020, Italy
| | - Enrico Negrisolo
- Department of Comparative Biomedicine and Food Science, University of Padova, Legnaro 35020, Italy.,CRIBI Interdepartmental Research Centre for Innovative Biotechnologies, University of Padova, viale G. Colombo 3, Padova 35121, Italy
| |
Collapse
|
29
|
Chang Y, Rochon D, Sekimoto S, Wang Y, Chovatia M, Sandor L, Salamov A, Grigoriev IV, Stajich JE, Spatafora JW. Genome-scale phylogenetic analyses confirm Olpidium as the closest living zoosporic fungus to the non-flagellated, terrestrial fungi. Sci Rep 2021; 11:3217. [PMID: 33547391 PMCID: PMC7865070 DOI: 10.1038/s41598-021-82607-4] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Accepted: 01/19/2021] [Indexed: 12/24/2022] Open
Abstract
The zoosporic obligate endoparasites, Olpidium, hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium, but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium, Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.
Collapse
Affiliation(s)
- Ying Chang
- Department of Botany and Plant Pathology, College of Agricultural Sciences, Oregon State University, Oregon, USA.
| | - D'Ann Rochon
- Agriculture and Agri-Food Canada, Pacific Agri-Food Research Centre, Summerland, BC, Canada
| | - Satoshi Sekimoto
- Department of Botany, University of British Columbia, Vancouver, BC, Canada
- Research and Development Center, Mitsubishi-Chemical Foods Corporation, Yokohama, Japan
| | - Yan Wang
- Department of Microbiology and Plant Pathology, University of California Riverside, Riverside, CA, USA
- Institute for Integrative Genome Biology, University of California, Riverside, CA, USA
- Department of Biological Sciences, University of Toronto Scarborough, Toronto, Canada
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Canada
| | - Mansi Chovatia
- Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Laura Sandor
- Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Asaf Salamov
- Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Igor V Grigoriev
- Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Department of Plant and Microbial Biology, University of California Berkeley, Berkeley, CA, USA
| | - Jason E Stajich
- Department of Microbiology and Plant Pathology, University of California Riverside, Riverside, CA, USA
- Institute for Integrative Genome Biology, University of California, Riverside, CA, USA
| | - Joseph W Spatafora
- Department of Botany and Plant Pathology, College of Agricultural Sciences, Oregon State University, Oregon, USA
| |
Collapse
|
30
|
Varney RM, Brenzinger B, Malaquias MAE, Meyer CP, Schrödl M, Kocot KM. Assessment of mitochondrial genomes for heterobranch gastropod phylogenetics. BMC Ecol Evol 2021; 21:6. [PMID: 33514315 PMCID: PMC7853304 DOI: 10.1186/s12862-020-01728-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Accepted: 11/26/2020] [Indexed: 01/07/2023] Open
Abstract
Background Heterobranchia is a diverse clade of marine, freshwater, and terrestrial gastropod molluscs. It includes such disparate taxa as nudibranchs, sea hares, bubble snails, pulmonate land snails and slugs, and a number of (mostly small-bodied) poorly known snails and slugs collectively referred to as the “lower heterobranchs”. Evolutionary relationships within Heterobranchia have been challenging to resolve and the group has been subject to frequent and significant taxonomic revision. Mitochondrial (mt) genomes can be a useful molecular marker for phylogenetics but, to date, sequences have been available for only a relatively small subset of Heterobranchia. Results To assess the utility of mitochondrial genomes for resolving evolutionary relationships within this clade, eleven new mt genomes were sequenced including representatives of several groups of “lower heterobranchs”. Maximum likelihood analyses of concatenated matrices of the thirteen protein coding genes found weak support for most higher-level relationships even after several taxa with extremely high rates of evolution were excluded. Bayesian inference with the CAT + GTR model resulted in a reconstruction that is much more consistent with the current understanding of heterobranch phylogeny. Notably, this analysis recovered Valvatoidea and Orbitestelloidea in a polytomy with a clade including all other heterobranchs, highlighting these taxa as important to understanding early heterobranch evolution. Also, dramatic gene rearrangements were detected within and between multiple clades. However, a single gene order is conserved across the majority of heterobranch clades. Conclusions Analysis of mitochondrial genomes in a Bayesian framework with the site heterogeneous CAT + GTR model resulted in a topology largely consistent with the current understanding of heterobranch phylogeny. However, mitochondrial genomes appear to be too variable to serve as good phylogenetic markers for robustly resolving a number of deeper splits within this clade.
Collapse
Affiliation(s)
- Rebecca M Varney
- Department of Biological Sciences, The University of Alabama, Campus Box 870344, Tuscaloosa, AL, 35487, USA
| | - Bastian Brenzinger
- SNSB-Bavarian State Collection of Zoology, Münchhausenstr. 21, 81247, München, Germany
| | | | - Christopher P Meyer
- National Museum of Natural History, Smithsonian Institution, 10th St. & Constitution Ave. NW, Washington, D.C., 20560, USA
| | - Michael Schrödl
- SNSB-Bavarian State Collection of Zoology, Münchhausenstr. 21, 81247, München, Germany.,BioGeoCenter LMU (Ludwig Maximillion University Munich), University of Munich, Biozentrum, Großhaderner Str. 2, 82152, Planegg-Martinsried, Germany
| | - Kevin M Kocot
- Department of Biological Sciences, The University of Alabama, Campus Box 870344, Tuscaloosa, AL, 35487, USA. .,Alabama Museum of Natural History, Campus Box 870344, Tuscaloosa, AL, 35487, USA.
| |
Collapse
|
31
|
Warnow T, Mirarab S. Multiple Sequence Alignment for Large Heterogeneous Datasets Using SATé, PASTA, and UPP. Methods Mol Biol 2021; 2231:99-119. [PMID: 33289889 DOI: 10.1007/978-1-0716-1036-7_7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The estimation of very large multiple sequence alignments is a challenging problem that requires special techniques in order to achieve high accuracy. Here we describe two software packages-PASTA and UPP-for constructing alignments on large and ultra-large datasets. Both methods have been able to produce highly accurate alignments on 1,000,000 sequences, and trees computed on these alignments are also highly accurate. PASTA provides the best tree accuracy when the input sequences are all full-length, but UPP provides improved accuracy compared to PASTA and other methods when the input contains a large number of fragmentary sequences. Both methods are available in open source form on GitHub.
Collapse
Affiliation(s)
- Tandy Warnow
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| | - Siavash Mirarab
- Electrical and Computer Engineering, University of California at San Diego, La Jolla, CA, USA
| |
Collapse
|
32
|
Roa-Varón A, Dikow RB, Carnevale G, Tornabene L, Baldwin CC, Li C, Hilton EJ. Confronting Sources of Systematic Error to Resolve Historically Contentious Relationships: A Case Study Using Gadiform Fishes (Teleostei, Paracanthopterygii, Gadiformes). Syst Biol 2020; 70:739-755. [PMID: 33346841 PMCID: PMC8561434 DOI: 10.1093/sysbio/syaa095] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Revised: 11/27/2020] [Accepted: 12/02/2020] [Indexed: 11/14/2022] Open
Abstract
Reliable estimation of phylogeny is central to avoid inaccuracy in downstream macroevolutionary inferences. However, limitations exist in the implementation of concatenated and summary coalescent approaches, and Bayesian and full coalescent inference methods may not yet be feasible for computation of phylogeny using complicated models and large data sets. Here, we explored methodological (e.g., optimality criteria, character sampling, model selection) and biological (e.g., heterotachy, branch length heterogeneity) sources of systematic error that can result in biased or incorrect parameter estimates when reconstructing phylogeny by using the gadiform fishes as a model clade. Gadiformes include some of the most economically important fishes in the world (e.g., Cods, Hakes, and Rattails). Despite many attempts, a robust higher-level phylogenetic framework was lacking due to limited character and taxonomic sampling, particularly from several species-poor families that have been recalcitrant to phylogenetic placement. We compiled the first phylogenomic data set, including 14,208 loci (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$>$\end{document}2.8 M bp) from 58 species representing all recognized gadiform families, to infer a time-calibrated phylogeny for the group. Data were generated with a gene-capture approach targeting coding DNA sequences from single-copy protein-coding genes. Species-tree and concatenated maximum-likelihood (ML) analyses resolved all family-level relationships within Gadiformes. While there were a few differences between topologies produced by the DNA and the amino acid data sets, most of the historically unresolved relationships among gadiform lineages were consistently well resolved with high support in our analyses regardless of the methodological and biological approaches used. However, at deeper levels, we observed inconsistency in branch support estimates between bootstrap and gene and site coefficient factors (gCF, sCF). Despite numerous short internodes, all relationships received unequivocal bootstrap support while gCF and sCF had very little support, reflecting hidden conflict across loci. Most of the gene-tree and species-tree discordance in our study is a result of short divergence times, and consequent lack of informative characters at deep levels, rather than incomplete lineage sorting. We use this phylogeny to establish a new higher-level classification of Gadiformes as a way of clarifying the evolutionary diversification of the order. We recognize 17 families in five suborders: Bregmacerotoidei, Gadoidei, Ranicipitoidei, Merluccioidei, and Macrouroidei (including two subclades). A time-calibrated analysis using 15 fossil taxa suggests that Gadiformes evolved \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$\sim $\end{document}79.5 Ma in the late Cretaceous, but that most extant lineages diverged after the Cretaceous–Paleogene (K-Pg) mass extinction (66 Ma). Our results reiterate the importance of examining phylogenomic analyses for evidence of systematic error that can emerge as a result of unsuitable modeling of biological factors and/or methodological issues, even when data sets are large and yield high support for phylogenetic relationships. [Branch length heterogeneity; Codfishes; commercial fish species; Cretaceous-Paleogene (K-Pg); heterotachy; systematic error; target enrichment.]
Collapse
Affiliation(s)
- Adela Roa-Varón
- National Systematics Laboratory of the National Oceanic Atmospheric Administration Fisheries Service, 10th St. & Constitution Ave. NW, Washington, DC 20560, USA.,Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, 10th St. & Constitution Ave. NW, Washington, DC 20560, USA.,Virginia Institute of Marine Science, William & Mary, Gloucester Point, VA 23062, USA
| | - Rebecca B Dikow
- Data Science Lab, Office of the Chief Information Officer, Smithsonian Institution, 10th St. & Constitution Ave. NW, Washington, DC 20560, USA
| | - Giorgio Carnevale
- Dipartimento di Scienze della Terra, Università degli Studi di Torino, Via Valperga Caluso, 35, 10125 Torino, Italy
| | - Luke Tornabene
- School of Aquatic and Fishery Sciences, University of Washington, 1122 NE Boat Street, Seattle, WA 98105, USA.,Burke Museum of Natural History and Culture, 4300 15th Ave NE, Seattle, WA 98105, USA
| | - Carole C Baldwin
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, 10th St. & Constitution Ave. NW, Washington, DC 20560, USA
| | - Chenhong Li
- Shanghai Universities Key Laboratory of Marine Animal Taxonomy and Evolution, 999 Hucheng Ring Rd, Pudong, Shanghai, China
| | - Eric J Hilton
- Virginia Institute of Marine Science, William & Mary, Gloucester Point, VA 23062, USA
| |
Collapse
|
33
|
Leuchtenberger AF, Crotty SM, Drucks T, Schmidt HA, Burgstaller-Muehlbacher S, von Haeseler A. Distinguishing Felsenstein Zone from Farris Zone Using Neural Networks. Mol Biol Evol 2020; 37:3632-3641. [PMID: 32637998 PMCID: PMC7743852 DOI: 10.1093/molbev/msaa164] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Maximum likelihood and maximum parsimony are two key methods for phylogenetic tree reconstruction. Under certain conditions, each of these two methods can perform more or less efficiently, resulting in unresolved or disputed phylogenies. We show that a neural network can distinguish between four-taxon alignments that were evolved under conditions susceptible to either long-branch attraction or long-branch repulsion. When likelihood and parsimony methods are discordant, the neural network can provide insight as to which tree reconstruction method is best suited to the alignment. When applied to the contentious case of Strepsiptera evolution, our method shows robust support for the current scientific view, that is, it places Strepsiptera with beetles, distant from flies.
Collapse
Affiliation(s)
- Alina F Leuchtenberger
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Stephen M Crotty
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
- School of Mathematical Sciences, University of Adelaide, Adelaide, SA, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, University of Adelaide, Adelaide, SA, Australia
| | - Tamara Drucks
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Heiko A Schmidt
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Sebastian Burgstaller-Muehlbacher
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Arndt von Haeseler
- Center for Integrative Bioinformatics Vienna, Max Perutz Labs, University of Vienna and Medical University of Vienna, Vienna, Austria
- Bioinformatics and Computational Biology, Faculty of Computer Science, University of Vienna, Vienna, Austria
| |
Collapse
|
34
|
A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly. BMC Evol Biol 2020; 20:141. [PMID: 33138771 PMCID: PMC7607883 DOI: 10.1186/s12862-020-01703-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Accepted: 10/19/2020] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics. RESULTS Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS). CONCLUSIONS Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.
Collapse
|
35
|
Le Kim T, Le Sy V. mPartition: A Model-Based Method for Partitioning Alignments. J Mol Evol 2020; 88:641-652. [PMID: 32864711 DOI: 10.1007/s00239-020-09963-z] [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/31/2019] [Accepted: 08/08/2020] [Indexed: 10/23/2022]
Abstract
Maximum likelihood (ML) analysis of nucleotide or amino-acid alignments is widely used to infer evolutionary relationships among species. Computing the likelihood of a phylogenetic tree from such alignments is a complicated task because the evolutionary processes typically vary across sites. A number of studies have shown that partitioning alignments into sub-alignments of sites, where each sub-alignment is analyzed using a different model of evolution (e.g., GTR + I + G), is a sensible strategy. Current partitioning methods group sites into subsets based on the inferred rates of evolution at the sites. However, these do not provide sufficient information to adequately reflect the substitution processes of characters at the sites. Moreover, the site rate-based methods group all invariant sites into one subset, potentially resulting in wrong phylogenetic trees. In this study, we propose a partitioning method, called mPartition, that combines not only the evolutionary rates but also substitution models at sites to partition alignments. Analyses of different partitioning methods on both real and simulated datasets showed that mPartition was better than the other partitioning methods tested. Notably, mPartition overcame the pitfall of grouping all invariant sites into one subset. Using mPartition may lead to increased accuracy of ML-based phylogenetic inference, especially for multiple loci or whole genome datasets.
Collapse
Affiliation(s)
- Thu Le Kim
- University of Engineering and Technology, Vietnam National University Hanoi, 144 Xuan Thuy, Cau Giay, Hanoi, 10000, Vietnam.,Hanoi University of Science and Technology, 1st Dai Co Viet, Hai Ba Trung, Hanoi, 10000, Vietnam
| | - Vinh Le Sy
- University of Engineering and Technology, Vietnam National University Hanoi, 144 Xuan Thuy, Cau Giay, Hanoi, 10000, Vietnam.
| |
Collapse
|
36
|
Smirnov V, Warnow T. Phylogeny Estimation Given Sequence Length Heterogeneity. Syst Biol 2020; 70:268-282. [PMID: 32692823 PMCID: PMC7875441 DOI: 10.1093/sysbio/syaa058] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Revised: 07/14/2020] [Accepted: 07/15/2020] [Indexed: 12/21/2022] Open
Abstract
Phylogeny estimation is a major step in many biological studies, and has many well known challenges. With the dropping cost of sequencing technologies, biologists now have increasingly large datasets available for use in phylogeny estimation. Here we address the challenge of estimating a tree given large datasets with a combination of full-length sequences and fragmentary sequences, which can arise due to a variety of reasons, including sample collection, sequencing technologies, and analytical pipelines. We compare two basic approaches: (1) computing an alignment on the full dataset and then computing a maximum likelihood tree on the alignment, or (2) constructing an alignment and tree on the full length sequences and then using phylogenetic placement to add the remaining sequences (which will generally be fragmentary) into the tree. We explore these two approaches on a range of simulated datasets, each with 1000 sequences and varying in rates of evolution, and two biological datasets. Our study shows some striking performance differences between methods, especially when there is substantial sequence length heterogeneity and high rates of evolution. We find in particular that using UPP to align sequences and RAxML to compute a tree on the alignment provides the best accuracy, substantially outperforming trees computed using phylogenetic placement methods. We also find that FastTree has poor accuracy on alignments containing fragmentary sequences. Overall, our study provides insights into the literature comparing different methods and pipelines for phylogenetic estimation, and suggests directions for future method development. [Phylogeny estimation, sequence length heterogeneity, phylogenetic placement.]
Collapse
Affiliation(s)
- Vladimir Smirnov
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| | - Tandy Warnow
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| |
Collapse
|
37
|
Jermiin LS, Catullo RA, Holland BR. A new phylogenetic protocol: dealing with model misspecification and confirmation bias in molecular phylogenetics. NAR Genom Bioinform 2020; 2:lqaa041. [PMID: 33575594 PMCID: PMC7671319 DOI: 10.1093/nargab/lqaa041] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Revised: 05/18/2020] [Accepted: 06/04/2020] [Indexed: 12/15/2022] Open
Abstract
Molecular phylogenetics plays a key role in comparative genomics and has increasingly significant impacts on science, industry, government, public health and society. In this paper, we posit that the current phylogenetic protocol is missing two critical steps, and that their absence allows model misspecification and confirmation bias to unduly influence phylogenetic estimates. Based on the potential offered by well-established but under-used procedures, such as assessment of phylogenetic assumptions and tests of goodness of fit, we introduce a new phylogenetic protocol that will reduce confirmation bias and increase the accuracy of phylogenetic estimates.
Collapse
Affiliation(s)
- Lars S Jermiin
- CSIRO Land & Water, Canberra, ACT 2601, Australia
- Research School of Biology, Australian National University, Canberra, ACT 2601, Australia
- School of Biology & Environment Science, University College Dublin, Belfield, Dublin 4, Ireland
- Earth Institute, University College Dublin, Belfield, Dublin 4, Ireland
| | - Renee A Catullo
- CSIRO Land & Water, Canberra, ACT 2601, Australia
- Research School of Biology, Australian National University, Canberra, ACT 2601, Australia
- School of Science and Health & Hawkesbury Institute of the Environment, Western Sydney University, Penrith, NSW 2751, Australia
| | - Barbara R Holland
- School of Natural Sciences, University of Tasmania, Hobart, TAS 7001, Australia
| |
Collapse
|
38
|
Zamora J, Ekman S. Phylogeny and character evolution in the Dacrymycetes, and systematics of Unilacrymaceae and Dacryonaemataceae fam. nov. PERSOONIA 2020; 44:161-205. [PMID: 33116340 PMCID: PMC7567964 DOI: 10.3767/persoonia.2020.44.07] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/22/2019] [Accepted: 10/24/2019] [Indexed: 02/07/2023]
Abstract
We present a multilocus phylogeny of the class Dacrymycetes, based on data from the 18S, ITS, 28S, RPB1, RPB2, TEF-1α, 12S, and ATP6 DNA regions, with c. 90 species including the types of most currently accepted genera. A variety of methodological approaches was used to infer phylogenetic relationships among the Dacrymycetes, from a supermatrix strategy using maximum likelihood and Bayesian inference on a concatenated dataset, to coalescence-based calculations, such as quartet-based summary methods of independent single-locus trees, and Bayesian integration of single-locus trees into a species tree under the multispecies coalescent. We evaluate for the first time the taxonomic usefulness of some cytological phenotypic characters, i.e., vacuolar contents (vacuolar bodies and lipid bodies), number of nuclei of recently discharged basidiospores, and pigments, with especial emphasis on carotenoids. These characters, along with several others traditionally used for the taxonomy of this group (basidium shape, presence and morphology of clamp connections, morphology of the terminal cells of cortical/marginal hyphae, presence and degree of ramification of the hyphidia), are mapped on the resulting phylogenies and their evolution through the class Dacrymycetes discussed. Our analyses reveal five lineages that putatively represent five different families, four of which are accepted and named. Three out of these four lineages correspond to previously circumscribed and published families (Cerinomycetaceae, Dacrymycetaceae, and Unilacrymaceae), and one is proposed as the new family Dacryonaemataceae. Provisionally, only a single order, Dacrymycetales, is accepted within the class. Furthermore, the systematics of the two smallest families, Dacryonaemataceae and Unilacrymaceae, are investigated to the species level, using coalescence-based species delimitation on multilocus DNA data, and a detailed morphological study including morphometric analyses of the basidiospores. Three species are accepted in Dacryonaema, the type, Da. rufum, the newly combined Da. macnabbii (basionym Dacrymyces macnabbii), and a new species named Da. macrosporum. Two species are accepted in Unilacryma, the new U. bispora, and the type, U. unispora, the latter treated in a broad sense pending improved sampling across the Holarctic.
Collapse
Affiliation(s)
- J.C. Zamora
- Museum of Evolution, Uppsala University, Norbyvägen 16, SE-75236 Uppsala, Sweden
- Departamento de Biología Vegetal II, Facultad de Farmacia, Universidad Complutense de Madrid, Ciudad Universitaria, plaza de Ramón y Cajal s/n, E-28040, Madrid, Spain
| | - S. Ekman
- Museum of Evolution, Uppsala University, Norbyvägen 16, SE-75236 Uppsala, Sweden
| |
Collapse
|
39
|
Lu C, Huang X, Deng J. The challenge of Coccidae (Hemiptera: Coccoidea) mitochondrial genomes: The case of Saissetia coffeae with novel truncated tRNAs and gene rearrangements. Int J Biol Macromol 2020; 158:854-864. [PMID: 32387610 DOI: 10.1016/j.ijbiomac.2020.04.257] [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: 12/05/2019] [Revised: 04/26/2020] [Accepted: 04/27/2020] [Indexed: 11/23/2022]
Abstract
There have been few reports of complete mitochondrial genomes (mitogenomes) of scale insects, and it has been indicated that complex and novel structures in their mitogenomes may lead to difficulties in sequencing, assembly and annotation. Transfer RNAs (tRNAs) usually possess typical cloverleaf secondary structures, and truncated tRNAs are rarely found in insect mitogenomes. Here, we report a complete Saissetia coffeae mitogenome (15,389 bp) with high A+T content (84.7%) sequenced by next-generation sequencing (NGS) methods. Genes in the mitogenome were annotated, and nine tRNAs were not found using MITOS. Most of the detected tRNAs were significantly truncated without the dihydrouridine (DHU) arm or the TΨC (T) arm. In addition, the 9 "lost" tRNAs containing mismatched base pairs were retrieved based on the tRNA annotation workflow for Coccidae described in our study. The gene arrangement in the Saissetia coffeae mitogenome was significantly different from that in other hemipteran insects. Additionally, Bayesian and maximum likelihood trees based on the mitochondrial genes showed a long branch of the Saissetia lineage, indicating significant nonsynonymous substitutions or high evolutionary rates in the Saissetia lineage. We provide a reference mitogenome for the assembly and annotation of the Coccidae mitogenome and offer insights into the evolution of scale insects.
Collapse
Affiliation(s)
- Congcong Lu
- State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou 350002, China
| | - Xiaolei Huang
- State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou 350002, China
| | - Jun Deng
- State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou 350002, China.
| |
Collapse
|
40
|
Sun X, Yu D, Xie Z, Dong J, Ding Y, Yao H, Greenslade P. Phylomitogenomic analyses on collembolan higher taxa with enhanced taxon sampling and discussion on method selection. PLoS One 2020; 15:e0230827. [PMID: 32282807 PMCID: PMC7153868 DOI: 10.1371/journal.pone.0230827] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 03/09/2020] [Indexed: 12/20/2022] Open
Abstract
Collembola are a basal group of Hexapoda renowned for both unique morphological characters and significant ecological roles. However, a robust and plausible phylogenetic relationship between its deeply divergent lineages has yet to be achieved. We carried out a mitophylogenomic study based on a so far the most comprehensive mitochondrial genome dataset. Our data matrix contained mitogenomes of 31 species from almost all major families of all four orders, with 16 mitogenomes newly sequenced and annotated. We compared the linear arrangements of genes along mitochondria across species. Then we conducted 13 analyses each under a different combination of character coding, partitioning scheme and heterotachy models, and assessed their performance in phylogenetic inference. Several hypothetical tree topologies were also tested. Mitogenomic structure comparison revealed that most species share the same gene order of putative ancestral pancrustacean pattern, while seven species from Onychiuridae, Poduridae and Symphypleona bear different levels of gene rearrangements, indicating phylogenetic signals. Tomoceroidea was robustly recovered for the first time in the presence of all its families and subfamilies. Monophyly of Onychiuroidea was supported using unpartitioned models alleviating LBA. Paronellidae was revealed polyphyletic with two subfamilies inserted independently into Entomobryidae. Although Entomobryomorpha has not been well supported, more than half of the analyses obtained convincing topologies by placing Tomoceroidea within or near remaining Entomobryomorpha. The relationship between elongate-shaped and spherical-shaped collembolans still remained ambiguous, but Neelipleona tend to occupy the basal position in most trees. This study showed that mitochondrial genomes could provide important information for reconstructing the relationships among Collembola when suitable analytical approaches are implemented. Of all the data refining and model selecting schemes used in this study, the combination of nucleotide sequences, partitioning model and exclusion of third codon positions performed better in generating more reliable tree topology and higher node supports than others.
Collapse
Affiliation(s)
- Xin Sun
- J.F. Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Göttingen, Germany
- Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun, China
| | - Daoyuan Yu
- Soil Ecology Lab, College of Resources and Environmental Sciences, Nanjing Agricultural University, Nanjing, China
- * E-mail: ,
| | - Zhijing Xie
- Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Jie Dong
- Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing, China
| | - Yinhuan Ding
- Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing, China
| | - Haifeng Yao
- Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Penelope Greenslade
- Environmental Management, School of Applied and Biomedical Science, Federation University, Ballarat, Victoria, Australia
- Division of Biology, Australian National University, Australian Capital Territory, Australia
| |
Collapse
|
41
|
Sun X, Ding Y, Orr MC, Zhang F. Streamlining universal single-copy orthologue and ultraconserved element design: A case study in Collembola. Mol Ecol Resour 2020; 20. [PMID: 32065730 DOI: 10.1111/1755-0998.13146] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Revised: 01/18/2020] [Accepted: 02/10/2020] [Indexed: 11/27/2022]
Abstract
Genomic data sets are increasingly central to ecological and evolutionary biology, but far fewer resources are available for invertebrates. Powerful new computational tools and the rapidly decreasing cost of Illumina sequencing are beginning to change this, enabling rapid genome assembly and reference marker extraction. We have developed and tested a practical workflow for developing genomic resources in nonmodel groups with real-world data on Collembola (springtails), one of the most dominant soil animals on Earth. We designed universal molecular marker sets, single-copy orthologues (BUSCOs) and ultraconserved elements (UCEs), using three existing and 11 newly generated genomes. Both marker types were tested in silico via marker capture success and phylogenetic performance. The new genomes were assembled with Illumina short reads and 9,585-14,743 protein-coding genes were predicted with ab initio and protein homology evidence. We identified 1,997 benchmarking universal single-copy orthologues (BUSCOs) across 14 genomes and created and assessed a custom BUSCO data set for extracting single-copy genes. We also developed a new UCE probe set containing 46,087 baits targeting 1,885 loci. We successfully captured 1,437-1,865 BUSCOs and 975-1,186 UCEs across 14 genomes. Phylogenomic reconstructions using these markers proved robust, giving new insight on deep-time collembolan relationships. Our study demonstrates the feasibility of generating thousands of universal markers from highly efficient whole-genome sequencing, providing a valuable resource for genome-scale investigations in evolutionary biology and ecology.
Collapse
Affiliation(s)
- Xin Sun
- J. F. Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Göttingen, Germany.,Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing, China
| | - Yinhuan Ding
- Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing, China
| | - Michael C Orr
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Feng Zhang
- Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing, China.,Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
42
|
Williams TA, Cox CJ, Foster PG, Szöllősi GJ, Embley TM. Phylogenomics provides robust support for a two-domains tree of life. Nat Ecol Evol 2020; 4:138-147. [PMID: 31819234 PMCID: PMC6942926 DOI: 10.1038/s41559-019-1040-x] [Citation(s) in RCA: 132] [Impact Index Per Article: 26.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2019] [Accepted: 10/15/2019] [Indexed: 11/09/2022]
Abstract
Hypotheses about the origin of eukaryotic cells are classically framed within the context of a universal 'tree of life' based on conserved core genes. Vigorous ongoing debate about eukaryote origins is based on assertions that the topology of the tree of life depends on the taxa included and the choice and quality of genomic data analysed. Here we have reanalysed the evidence underpinning those claims and apply more data to the question by using supertree and coalescent methods to interrogate >3,000 gene families in archaea and eukaryotes. We find that eukaryotes consistently originate from within the archaea in a two-domains tree when due consideration is given to the fit between model and data. Our analyses support a close relationship between eukaryotes and Asgard archaea and identify the Heimdallarchaeota as the current best candidate for the closest archaeal relatives of the eukaryotic nuclear lineage.
Collapse
Affiliation(s)
- Tom A Williams
- School of Biological Sciences, University of Bristol, Bristol, UK.
| | - Cymon J Cox
- Centro de Ciências do Mar, Universidade do Algarve, Faro, Portugal
| | - Peter G Foster
- Department of Life Sciences, Natural History Museum, London, UK
| | - Gergely J Szöllősi
- MTA-ELTE "Lendület" Evolutionary Genomics Research Group, Budapest, Hungary
- Department of Biological Physics, Eötvös Loránd University, Budapest, Hungary
- Evolutionary Systems Research Group, Centre for Ecological Research, Hungarian Academy of Sciences, Tihany, Hungary
| | - T Martin Embley
- Institute for Cell and Molecular Biosciences, University of Newcastle, Newcastle upon Tyne, UK.
| |
Collapse
|
43
|
Zhang D, Li WX, Zou H, Wu SG, Li M, Jakovlić I, Zhang J, Chen R, Wang G. Homoplasy or plesiomorphy? Reconstruction of the evolutionary history of mitochondrial gene order rearrangements in the subphylum Neodermata. Int J Parasitol 2019; 49:819-829. [PMID: 31401064 DOI: 10.1016/j.ijpara.2019.05.010] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2019] [Revised: 05/15/2019] [Accepted: 05/22/2019] [Indexed: 12/31/2022]
Abstract
Recent mitogenomic studies have exposed a gene order (GO) shared by two classes, four orders and 31 species ('common GO') within the flatworm subphylum Neodermata. There are two possible hypotheses for this phenomenon: convergent evolution (homoplasy) or shared ancestry (plesiomorphy). To test those, we conducted a meta-analysis on all available mitogenomes to infer the evolutionary history of GO in Neodermata. To improve the resolution, we added a newly sequenced mitogenome that exhibited the common GO, Euryhaliotrema johni (Ancyrocephalinae), to the dataset. Phylogenetic analyses conducted on two datasets (nucleotides of all 36 genes and amino acid sequences of 12 protein coding genes) and four algorithms (MrBayes, RAxML, IQ-TREE and PhyloBayes) produced topology instability towards the tips, so ancestral GO reconstructions were conducted using TreeREx and MLGO programs using all eight obtained topologies, plus three unique topologies from previous studies. The results consistently supported the second hypothesis, resolving the common GO as a plesiomorphic ancestral GO for Neodermata, Cestoda, Monopisthocotylea, Cestoda + Trematoda and Cestoda + Trematoda + Monopisthocotylea. This allowed us to trace the evolutionary GO scenarios from each common ancestor to its descendants amongst the Monogenea and Cestoda classes, and propose that the common GO was most likely retained throughout all of the common ancestors, leading to the extant species possessing the common GO. Neodermatan phylogeny inferred from GOs was largely incongruent with all 11 topologies described above, but it did support the mitogenomic dataset in resolving Polyopisthocotylea as the earliest neodermatan branch. Although highly derived GOs might be of some use in resolving isolated taxonomic and phylogenetic uncertainties, we conclude that, due to the discontinuous nature of their evolution, they tend to produce artefactual phylogenetic relationships, which makes them unsuitable for phylogenetic reconstruction in Neodermata. Wider and denser sampling of neodermatan mitogenomic sequences will be needed to infer the evolutionary pathways leading to the observed diversity of GOs with confidence.
Collapse
Affiliation(s)
- Dong Zhang
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, PR China; University of Chinese Academy of Sciences, Beijing, PR China
| | - Wen X Li
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, PR China
| | - Hong Zou
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, PR China
| | - Shan G Wu
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, PR China
| | - Ming Li
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, PR China
| | | | - Jin Zhang
- Bio-Transduction Lab, Wuhan 430075, PR China
| | - Rong Chen
- Bio-Transduction Lab, Wuhan 430075, PR China
| | - Guitang Wang
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, PR China.
| |
Collapse
|
44
|
Zhang D, Zou H, Hua CJ, Li WX, Mahboob S, Al-Ghanim KA, Al-Misned F, Jakovlić I, Wang GT. Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda. Genome Biol Evol 2019; 11:1797-1812. [PMID: 31192351 PMCID: PMC6601869 DOI: 10.1093/gbe/evz121] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/08/2019] [Indexed: 01/04/2023] Open
Abstract
The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.
Collapse
Affiliation(s)
- Dong Zhang
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
- University of Chinese Academy of Sciences, Beijing, P.R. China
| | - Hong Zou
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
| | - Cong-Jie Hua
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
| | - Wen-Xiang Li
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
| | - Shahid Mahboob
- Department of Zoology, College of Science, King Saud University, Riyadh, Saudi Arabia
- Department of Zoology, GC University, Faisalabad, Pakistan
| | | | - Fahad Al-Misned
- Department of Zoology, College of Science, King Saud University, Riyadh, Saudi Arabia
| | | | - Gui-Tang Wang
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, P.R. China
| |
Collapse
|