1
|
Veltman MA, Anthoons B, Schrøder-Nielsen A, Gravendeel B, de Boer HJ. Orchidinae-205: A new genome-wide custom bait set for studying the evolution, systematics, and trade of terrestrial orchids. Mol Ecol Resour 2024; 24:e13986. [PMID: 38899721 DOI: 10.1111/1755-0998.13986] [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: 06/15/2023] [Revised: 05/16/2024] [Accepted: 05/30/2024] [Indexed: 06/21/2024]
Abstract
Terrestrial orchids are a group of genetically understudied, yet culturally and economically important plants. The Orchidinae tribe contains many species that produce edible tubers that are used for the production of traditional delicacies collectively called 'salep'. Overexploitation of wild orchids in the Eastern Mediterranean and Western Asia threatens to drive many of these species to extinction, but cost-effective tools for monitoring their trade are currently lacking. Here we present a custom bait kit for target enrichment and sequencing of 205 novel genetic markers that are tailored to phylogenomic applications in Orchidinae s.l. A subset of 31 markers capture genes putatively involved in the production of glucomannan, a water-soluble polysaccharide that gives salep its distinctive properties. We tested the kit on 73 taxa native to the area, demonstrating universally high locus recovery irrespective of species identity, that exceeds the total sequence length obtained with alternative kits currently available. Phylogenetic inference with concatenation and coalescent approaches was robust and showed high levels of support for most clades, including some which were previously unresolved. Resolution for hybridizing and recently radiated lineages remains difficult, but could be further improved by analysing multiple haplotypes and the non-exonic sequences captured by our kit, with the promise to shed new light on the evolution of enigmatic taxa with a complex speciation history. Offering a step-up from traditional barcoding and universal markers, the genome-wide custom loci targeted by Orchidinae-205 are a valuable new resource to study the evolution, systematics and trade of terrestrial orchids.
Collapse
Affiliation(s)
- Margaretha A Veltman
- Natural History Museum, Oslo, Norway
- Naturalis Biodiversity Center, Leiden, Netherlands
| | | | | | - Barbara Gravendeel
- Naturalis Biodiversity Center, Leiden, Netherlands
- Radboud Institute for Biological and Environmental Sciences, Radboud University, Nijmegen, Netherlands
| | | |
Collapse
|
2
|
Peñafiel Loaiza N, Chafe AH, Moraes R M, Oleas NH, Roncal J. Genotyping-by-sequencing informs conservation of Andean palms sources of non-timber forest products. Evol Appl 2024; 17:e13765. [PMID: 39091352 PMCID: PMC11291087 DOI: 10.1111/eva.13765] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2023] [Revised: 06/30/2024] [Accepted: 07/18/2024] [Indexed: 08/04/2024] Open
Abstract
Conservation and sustainable management of lineages providing non-timber forest products are imperative under the current global biodiversity loss. Most non-timber forest species, however, lack genomic studies that characterize their intraspecific variation and evolutionary history, which inform species' conservation practices. Contrary to many lineages in the Andean biodiversity hotspot that exhibit high diversification, the genus Parajubaea (Arecaceae) has only three species despite the genus' origin 22 million years ago. Two of the three palm species, P. torallyi and P. sunkha, are non-timber forest species endemic to the Andes of Bolivia and are listed as IUCN endangered. The third species, P. cocoides, is a vulnerable species with unknown wild populations. We investigated the evolutionary relationships of Parajubaea species and the genetic diversity and structure of wild Bolivian populations. Sequencing of five low-copy nuclear genes (3753 bp) challenged the hypothesis that P. cocoides is a cultigen that originated from the wild Bolivian species. We further obtained up to 15,134 de novo single-nucleotide polymorphism markers by genotyping-by-sequencing of 194 wild Parajubaea individuals. Our total DNA sequencing effort rejected the taxonomic separation of the two Bolivian species. As expected for narrow endemic species, we observed low genetic diversity, but no inbreeding signal. We found three genetic clusters shaped by geographic distance, which we use to propose three management units. Different percentages of missing genotypic data did not impact the genetic structure of populations. We use the management units to recommend in situ conservation by creating new protected areas, and ex situ conservation through seed collection.
Collapse
Affiliation(s)
- Nicolás Peñafiel Loaiza
- Department of BiologyMemorial University of NewfoundlandSt. John'sNewfoundland and LabradorCanada
- Present address:
Chone y BabahoyoLojaEcuador
| | - Abigail H. Chafe
- Department of BiologyMemorial University of NewfoundlandSt. John'sNewfoundland and LabradorCanada
| | - Mónica Moraes R
- Herbario Nacional de Bolivia, Instituto de EcologíaUniversidad Mayor de San AndrésLa PazBolivia
| | - Nora H. Oleas
- Centro de Investigación de la Biodiversidad y Cambio Climático – BioCamb e Ingeniería en Biodiversidad y Recursos Genéticos, Facultad de Ciencias de Medio AmbienteUniversidad IndoaméricaQuitoEcuador
| | - Julissa Roncal
- Department of BiologyMemorial University of NewfoundlandSt. John'sNewfoundland and LabradorCanada
| |
Collapse
|
3
|
Williamson JL, Gyllenhaal EF, Bauernfeind SM, Bautista E, Baumann MJ, Gadek CR, Marra PP, Ricote N, Valqui T, Bozinovic F, Singh ND, Witt CC. Extreme elevational migration spurred cryptic speciation in giant hummingbirds. Proc Natl Acad Sci U S A 2024; 121:e2313599121. [PMID: 38739790 PMCID: PMC11126955 DOI: 10.1073/pnas.2313599121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 03/19/2024] [Indexed: 05/16/2024] Open
Abstract
The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how "outlier" species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world's largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.
Collapse
Affiliation(s)
- Jessie L. Williamson
- Museum of Southwestern Biology, University of New Mexico, Albuquerque, NM87131
- Department of Biology, University of New Mexico, Albuquerque, NM87131
- Cornell Lab of Ornithology, Cornell University, Ithaca, NY14850
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY14850
| | - Ethan F. Gyllenhaal
- Museum of Southwestern Biology, University of New Mexico, Albuquerque, NM87131
- Department of Biology, University of New Mexico, Albuquerque, NM87131
| | | | - Emil Bautista
- Centro de Ornitología y Biodiversidad, Lima15064, Peru
| | - Matthew J. Baumann
- Museum of Southwestern Biology, University of New Mexico, Albuquerque, NM87131
| | - Chauncey R. Gadek
- Museum of Southwestern Biology, University of New Mexico, Albuquerque, NM87131
- Department of Biology, University of New Mexico, Albuquerque, NM87131
- Environmental Stewardship, Los Alamos National Laboratory, Los Alamos, NM87545
| | - Peter P. Marra
- The Earth Commons Institute, Department of Biology, McCourt School of Public Policy, Georgetown University, Washington, DC20057
| | - Natalia Ricote
- Facultad de Artes Liberales, Departamento de Ciencias, Universidad Adolfo Ibáñez, Santiago7941169, Chile
| | - Thomas Valqui
- Centro de Ornitología y Biodiversidad, Lima15064, Peru
- Facultad de Ciencias Forestales, Universidad Nacional Agraria La Molina, Lima15024, Peru
| | - Francisco Bozinovic
- Departamento de Ecología, Center of Applied Ecology and Sustainability, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Santiago3542000, Chile
| | - Nadia D. Singh
- Department of Biology, Institute of Ecology and Evolution, University of Oregon, Eugene, OR97403
| | - Christopher C. Witt
- Museum of Southwestern Biology, University of New Mexico, Albuquerque, NM87131
- Department of Biology, University of New Mexico, Albuquerque, NM87131
| |
Collapse
|
4
|
Bolívar‐Leguizamón SD, Bocalini F, Silveira LF, Bravo GA. The role of biogeographical barriers on the historical dynamics of passerine birds with a circum-Amazonian distribution. Ecol Evol 2024; 14:e10860. [PMID: 38450322 PMCID: PMC10915597 DOI: 10.1002/ece3.10860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 01/03/2024] [Accepted: 01/08/2024] [Indexed: 03/08/2024] Open
Abstract
Common distributional patterns have provided the foundations of our knowledge of Neotropical biogeography. A distinctive pattern is the "circum-Amazonian distribution", which surrounds Amazonia across the forested lowlands south and east of the basin, the Andean foothills, the Venezuelan Coastal Range, and the Tepuis. The underlying evolutionary and biogeographical mechanisms responsible for this widespread pattern of avian distribution have yet to be elucidated. Here, we test the effects of biogeographical barriers in four species in the passerine family Thamnophilidae by performing comparative demographic analyses of genome-scale data. Specifically, we used flanking regions of ultraconserved regions to estimate population historical parameters and genealogical trees and tested demographic models reflecting contrasting biogeographical scenarios explaining the circum-Amazonian distribution. We found that taxa with circum-Amazonian distribution have at least two main phylogeographical clusters: (1) Andes, often extending into Central America and the Tepuis; and (2) the remaining of their distribution. These clusters are connected through corridors along the Chaco-Cerrado and southeastern Amazonia, allowing gene flow between Andean and eastern South American populations. Demographic histories are consistent with Pleistocene climatic fluctuations having a strong influence on the diversification history of circum-Amazonian taxa, Refugia played a crucial role, enabling both phenotypic and genetic differentiation, yet maintaining substantial interconnectedness to keep considerable levels of gene flow during different dry/cool and warm/humid periods. Additionally, steep environmental gradients appear to play a critical role in maintaining both genetic and phenotypic structure.
Collapse
Affiliation(s)
- Sergio D. Bolívar‐Leguizamón
- Seção de AvesMuseu de Zoologia da Universidade de São PauloSão PauloBrazil
- Laboratório de Zoologia de Vertebrados, Departamento de Ciências Biológicas, Escola Superior de Agricultura “Luiz de Queiroz” –ESALQ–Universidade de São PauloPiracicabaBrazil
| | - Fernanda Bocalini
- Seção de AvesMuseu de Zoologia da Universidade de São PauloSão PauloBrazil
| | - Luís F. Silveira
- Seção de AvesMuseu de Zoologia da Universidade de São PauloSão PauloBrazil
| | - Gustavo A. Bravo
- Seção de AvesMuseu de Zoologia da Universidade de São PauloSão PauloBrazil
- Sección de Ornitología, Colecciones Biológicas, Instituto de Investigación de Recursos Biológicos Alexander von HumboldtClaustro de San AgustínVilla de Leyva, BoyacáColombia
- Museum of Comparative Zoology and Department of Organismic and Evolutionary BiologyHarvard UniversityCambridgeMassachusettsUSA
| |
Collapse
|
5
|
Garrick RC. Genetic signatures of lineage fusion closely resemble population decline. Ecol Evol 2023; 13:e10725. [PMID: 37964788 PMCID: PMC10641302 DOI: 10.1002/ece3.10725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 10/24/2023] [Accepted: 10/27/2023] [Indexed: 11/16/2023] Open
Abstract
Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F S and Z nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.
Collapse
Affiliation(s)
- Ryan C. Garrick
- Department of BiologyUniversity of MississippiOxfordMississippiUSA
| |
Collapse
|
6
|
Yu PL, Fulton JC, Hudson OH, Huguet-Tapia JC, Brawner JT. Next-generation fungal identification using target enrichment and Nanopore sequencing. BMC Genomics 2023; 24:581. [PMID: 37784013 PMCID: PMC10544392 DOI: 10.1186/s12864-023-09691-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 09/21/2023] [Indexed: 10/04/2023] Open
Abstract
BACKGROUND Rapid and accurate pathogen identification is required for disease management. Compared to sequencing entire genomes, targeted sequencing may be used to direct sequencing resources to genes of interest for microbe identification and mitigate the low resolution that single-locus molecular identification provides. This work describes a broad-spectrum fungal identification tool developed to focus high-throughput Nanopore sequencing on genes commonly employed for disease diagnostics and phylogenetic inference. RESULTS Orthologs of targeted genes were extracted from 386 reference genomes of fungal species spanning six phyla to identify homologous regions that were used to design the baits used for enrichment. To reduce the cost of producing probes without diminishing the phylogenetic power, DNA sequences were first clustered, and then consensus sequences within each cluster were identified to produce 26,000 probes that targeted 114 genes. To test the efficacy of our probes, we applied the technique to three species representing Ascomycota and Basidiomycota fungi. The efficiency of enrichment, quantified as mean target coverage over the mean genome-wide coverage, ranged from 200 to 300. Furthermore, enrichment of long reads increased the depth of coverage across the targeted genes and into non-coding flanking sequence. The assemblies generated from enriched samples provided well-resolved phylogenetic trees for taxonomic assignment and molecular identification. CONCLUSIONS Our work provides data to support the utility of targeted Nanopore sequencing for fungal identification and provides a platform that may be extended for use with other phytopathogens.
Collapse
Affiliation(s)
- Pei-Ling Yu
- Department of Plant Pathology, University of Florida, Gainesville, FL, 32611, USA
| | - James C Fulton
- Department of Plant Pathology, University of Florida, Gainesville, FL, 32611, USA
- Florida Department of Agriculture and Consumer Services, Division of Plant Industry, Gainesville, FL, 32608, USA
| | - Owen H Hudson
- Department of Plant Pathology, University of Florida, Gainesville, FL, 32611, USA
| | - Jose C Huguet-Tapia
- Department of Plant Pathology, University of Florida, Gainesville, FL, 32611, USA
| | - Jeremy T Brawner
- Department of Plant Pathology, University of Florida, Gainesville, FL, 32611, USA.
| |
Collapse
|
7
|
Martchenko D, Shafer ABA. Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study. Heredity (Edinb) 2023; 131:273-281. [PMID: 37532838 PMCID: PMC10539292 DOI: 10.1038/s41437-023-00643-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2022] [Revised: 07/22/2023] [Accepted: 07/22/2023] [Indexed: 08/04/2023] Open
Abstract
Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.
Collapse
Affiliation(s)
- Daria Martchenko
- Environmental and Life Sciences Graduate Program, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada.
| | - Aaron B A Shafer
- Environmental and Life Sciences Graduate Program, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada
- Department of Forensics & Environmental and Life Sciences Graduate Program, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada
| |
Collapse
|
8
|
Gijsbers JC, Englebert N, Prata KE, Pichon M, Dinesen Z, Brunner R, Eyal G, González-Zapata FL, Kahng SE, Latijnhouwers KRW, Muir P, Radice VZ, Sánchez JA, Vermeij MJA, Hoegh-Guldberg O, Jacobs SJ, Bongaerts P. Global phylogenomic assessment of Leptoseris and Agaricia reveals substantial undescribed diversity at mesophotic depths. BMC Biol 2023; 21:147. [PMID: 37365558 DOI: 10.1186/s12915-023-01630-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 05/23/2023] [Indexed: 06/28/2023] Open
Abstract
BACKGROUND Mesophotic coral communities are increasingly gaining attention for the unique biological diversity they host, exemplified by the numerous mesophotic fish species that continue to be discovered. In contrast, many of the photosynthetic scleractinian corals observed at mesophotic depths are assumed to be depth-generalists, with very few species characterised as mesophotic-specialists. This presumed lack of a specialised community remains largely untested, as phylogenetic studies on corals have rarely included mesophotic samples and have long suffered from resolution issues associated with traditional sequence markers. RESULTS Here, we used reduced-representation genome sequencing to conduct a phylogenomic assessment of the two dominant mesophotic genera of plating corals in the Indo-Pacific and Western Atlantic, respectively, Leptoseris and Agaricia. While these genome-wide phylogenies broadly corroborated the morphological taxonomy, they also exposed deep divergences within the two genera and undescribed diversity across the current taxonomic species. Five of the eight focal species consisted of at least two sympatric and genetically distinct lineages, which were consistently detected across different methods. CONCLUSIONS The repeated observation of genetically divergent lineages associated with mesophotic depths highlights that there may be many more mesophotic-specialist coral species than currently acknowledged and that an urgent assessment of this largely unstudied biological diversity is warranted.
Collapse
Affiliation(s)
- J C Gijsbers
- California Academy of Sciences, San Francisco, CA, 94118, USA.
| | - N Englebert
- Global Change Institute, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - K E Prata
- California Academy of Sciences, San Francisco, CA, 94118, USA
- School of Biological Sciences, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - M Pichon
- Biodiversity Section, Queensland Museum, Townsville, 4810, Australia
| | - Z Dinesen
- Centre for Biodiversity and Conservation Science, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - R Brunner
- Global Change Institute, The University of Queensland, St Lucia, QLD, 4072, Australia
- ARC Centre of Excellence for Coral Reef Studies, James Cook University, Townsville, QLD, 4811, Australia
| | - G Eyal
- School of Biological Sciences, The University of Queensland, St Lucia, QLD, 4072, Australia
- ARC Centre of Excellence for Coral Reef Studies, The University of Queensland, St Lucia, QLD, 4072, Australia
- The Mina & Everard Goodman Faculty of Life Sciences, Bar-Ilan University, 5290002, Ramat Gan, Israel
| | - F L González-Zapata
- Laboratorio de Biología Molecular Marina (BIOMMAR), Departamento de Ciencias Biológicas, Facultad de Ciencias, Universidad de Los Andes, 111711, Bogotá, Colombia
| | - S E Kahng
- Department of Oceanography, University of Hawaii at Manoa, 1000 Pope Road, Honolulu, HI, 96822, USA
| | - K R W Latijnhouwers
- CARMABI Foundation, Piscaderabaai Z/N, PO Box 2090, Willemstad, Curaçao
- Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Science Park 700, 1098 XH, Amsterdam, The Netherlands
| | - P Muir
- Global Change Institute, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - V Z Radice
- Global Change Institute, The University of Queensland, St Lucia, QLD, 4072, Australia
- Department of Biological Sciences, Old Dominion University, Norfolk, VA, 23529, USA
| | - J A Sánchez
- Laboratorio de Biología Molecular Marina (BIOMMAR), Departamento de Ciencias Biológicas, Facultad de Ciencias, Universidad de Los Andes, 111711, Bogotá, Colombia
| | - M J A Vermeij
- CARMABI Foundation, Piscaderabaai Z/N, PO Box 2090, Willemstad, Curaçao
- Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Science Park 700, 1098 XH, Amsterdam, The Netherlands
| | - O Hoegh-Guldberg
- Global Change Institute, The University of Queensland, St Lucia, QLD, 4072, Australia
- School of Biological Sciences, The University of Queensland, St Lucia, QLD, 4072, Australia
- ARC Centre of Excellence for Coral Reef Studies, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - S J Jacobs
- California Academy of Sciences, San Francisco, CA, 94118, USA
| | - P Bongaerts
- California Academy of Sciences, San Francisco, CA, 94118, USA.
- Global Change Institute, The University of Queensland, St Lucia, QLD, 4072, Australia.
- CARMABI Foundation, Piscaderabaai Z/N, PO Box 2090, Willemstad, Curaçao.
| |
Collapse
|
9
|
Stiller J, Wilson NG, Rouse GW. Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture. BMC Biol 2023; 21:129. [PMID: 37248474 DOI: 10.1186/s12915-023-01628-9] [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: 01/26/2023] [Accepted: 05/16/2023] [Indexed: 05/31/2023] Open
Abstract
BACKGROUND Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture. RESULTS Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia. CONCLUSIONS We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.
Collapse
Affiliation(s)
- Josefin Stiller
- Scripps Institution of Oceanography, University of California San Diego, La Jolla, 92093 , USA.
- Centre for Biodiversity Genomics, University of Copenhagen, 2100, Copenhagen, Denmark.
| | - Nerida G Wilson
- Scripps Institution of Oceanography, University of California San Diego, La Jolla, 92093 , USA
- Research & Collections, Western Australian Museum, Perth, Western Australia, 6106, Australia
- School of Biological Sciences, University of Western Australia, Perth, Western Australia, 6009, Australia
| | - Greg W Rouse
- Scripps Institution of Oceanography, University of California San Diego, La Jolla, 92093 , USA.
| |
Collapse
|
10
|
Spaulding F, McLaughlin JF, Cheek RG, McCracken KG, Glenn TC, Winker K. Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex. Mol Phylogenet Evol 2023; 182:107733. [PMID: 36801373 PMCID: PMC10092703 DOI: 10.1016/j.ympev.2023.107733] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 01/31/2023] [Accepted: 02/09/2023] [Indexed: 02/18/2023]
Abstract
The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ∼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.
Collapse
Affiliation(s)
- Fern Spaulding
- University of Alaska Museum, University of Alaska Fairbanks, Fairbanks, AK, USA; Department of Biology and Wildlife, University of Alaska Fairbanks, Fairbanks, AK, USA.
| | - Jessica F McLaughlin
- Department of Environmental Science, Policy, and Management, University of California Berkeley, Berkeley, CA, USA
| | - Rebecca G Cheek
- Graduate Degree Program in Ecology, Department of Biology, Colorado State University, Fort Collins, CO, USA
| | - Kevin G McCracken
- University of Alaska Museum, University of Alaska Fairbanks, Fairbanks, AK, USA; Department of Biology, University of Miami, Coral Gables, FL, USA
| | - Travis C Glenn
- Department of Environmental Health Science, University of Georgia, Athens, GA, USA
| | - Kevin Winker
- University of Alaska Museum, University of Alaska Fairbanks, Fairbanks, AK, USA; Department of Biology and Wildlife, University of Alaska Fairbanks, Fairbanks, AK, USA
| |
Collapse
|
11
|
Bartoš O, Bohlen J, Šlechtová VB, Kočí J, Röslein J, Janko K. Sequence capture: Obsolete or irreplaceable? A thorough validation across phylogenetic distances and its applicability to hybrids and allopolyploids. Mol Ecol Resour 2023. [PMID: 37122140 DOI: 10.1111/1755-0998.13806] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 04/05/2023] [Accepted: 04/12/2023] [Indexed: 05/02/2023]
Abstract
As whole-genome sequencing has become pervasive, some have suggested that reduced genomic representation approaches, for example, sequence capture, are becoming obsolete. In the present study, we argue that these techniques still provide excellent tools in terms of price and quality of data as well as in their ability to provide markers with specific features, as required, for example, in phylogenomics. A potential drawback of the wide-scale application of reduced representation approaches could be their drop in efficiency with increasing phylogenetic distance from the reference species. While some studies have focused on the degree and performance of reduced representation techniques in such situations, to our knowledge, none of them evaluated their applicability to inter-specific hybrids and polyploids. This highlights a significant gap in current knowledge since there is increasing evidence for the frequent occurrence of natural hybrids and polyploids, as well as for the major importance of both phenomena in evolution. The main aim of the present study was to carry out a thorough validation of SEQcap applicability to (1) a set of non-model taxa with a wide range of phylogenetic relatedness and (2) inter-specific hybrids of various ploidies and genomic compositions. Considering the latter point, we especially focused on mechanisms causing allelic bias and consequent allelic dropout, as these could have confounding effects with respect to the evolutionary genomic dynamics of hybrids, especially in asexuals, which virtually reproduce as a frozen F1 generation.
Collapse
Affiliation(s)
- Oldřich Bartoš
- Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Libechov, Czech Republic
- Department of Zoology, Faculty of Science, Charles University, Prague, Czech Republic
| | - Jörg Bohlen
- Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Libechov, Czech Republic
| | - Vendula Bohlen Šlechtová
- Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Libechov, Czech Republic
| | - Jan Kočí
- Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Libechov, Czech Republic
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| | - Jan Röslein
- Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Libechov, Czech Republic
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| | - Karel Janko
- Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, The Czech Academy of Sciences, Libechov, Czech Republic
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| |
Collapse
|
12
|
Oury N, Noël C, Mona S, Aurelle D, Magalon H. From genomics to integrative species delimitation? The case study of the Indo-Pacific Pocillopora corals. Mol Phylogenet Evol 2023; 184:107803. [PMID: 37120114 DOI: 10.1016/j.ympev.2023.107803] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 04/06/2023] [Accepted: 04/24/2023] [Indexed: 05/01/2023]
Abstract
With the advent of genomics, sequencing thousands of loci from hundreds of individuals now appears feasible at reasonable costs, allowing complex phylogenies to be resolved. This is particularly relevant for cnidarians, for which insufficient data is available due to the small number of currently available markers and obscures species boundaries. Difficulties in inferring gene trees and morphological incongruences further blur the study and conservation of these organisms. Yet, can genomics alone be used to delimit species? Here, focusing on the coral genus Pocillopora, whose colonies play key roles in Indo-Pacific reef ecosystems but have challenged taxonomists for decades, we explored and discussed the usefulness of multiple criteria (genetics, morphology, biogeography and symbiosis ecology) to delimit species of this genus. Phylogenetic inferences, clustering approaches and species delimitation methods based on genome-wide single-nucleotide polymorphisms (SNP) were first used to resolve Pocillopora phylogeny and propose genomic species hypotheses from 356 colonies sampled across the Indo-Pacific (western Indian Ocean, tropical southwestern Pacific and south-east Polynesia). These species hypotheses were then compared to other lines of evidence based on genetic, morphology, biogeography and symbiont associations. Out of 21 species hypotheses delimited by genomics, 13 were strongly supported by all approaches, while six could represent either undescribed species or nominal species that have been synonymised incorrectly. Altogether, our results support (1) the obsolescence of macromorphology (i.e., overall colony and branches shape) but the relevance of micromorphology (i.e., corallite structures) to refine Pocillopora species boundaries, (2) the relevance of the mtORF (coupled with other markers in some cases) as a diagnostic marker of most species, (3) the requirement of molecular identification when species identity of colonies is absolutely necessary to interpret results, as morphology can blur species identification in the field, and (4) the need for a taxonomic revision of the genus Pocillopora. These results give new insights into the usefulness of multiple criteria for resolving Pocillopora, and more widely, scleractinian species boundaries, and will ultimately contribute to the taxonomic revision of this genus and the conservation of its species.
Collapse
Affiliation(s)
- Nicolas Oury
- UMR ENTROPIE (Université de La Réunion, IRD, IFREMER, Université de Nouvelle-Calédonie, CNRS), Université de La Réunion, St Denis, La Réunion, France; Laboratoire Cogitamus, Paris, France.
| | - Cyril Noël
- IFREMER - IRSI - Service de Bioinformatique (SeBiMER), Plouzané, France
| | - Stefano Mona
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, EPHE-PSL, Université PSL, CNRS, SU, UA, Paris, France; EPHE, PSL Research University, Paris, France; Laboratoire d'Excellence CORAIL, Perpignan, France
| | - Didier Aurelle
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, EPHE-PSL, Université PSL, CNRS, SU, UA, Paris, France; Aix Marseille Université, Université de Toulon, CNRS, IRD, MIO, Marseille, France
| | - Hélène Magalon
- UMR ENTROPIE (Université de La Réunion, IRD, IFREMER, Université de Nouvelle-Calédonie, CNRS), Université de La Réunion, St Denis, La Réunion, France; Laboratoire Cogitamus, Paris, France; Laboratoire d'Excellence CORAIL, Perpignan, France
| |
Collapse
|
13
|
Ortiz-Sepulveda CM, Genete M, Blassiau C, Godé C, Albrecht C, Vekemans X, Van Bocxlaer B. Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms. Mol Ecol Resour 2023; 23:659-679. [PMID: 36349833 DOI: 10.1111/1755-0998.13735] [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: 11/14/2021] [Revised: 10/09/2022] [Accepted: 10/19/2022] [Indexed: 11/10/2022]
Abstract
Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.
Collapse
Affiliation(s)
| | - Mathieu Genete
- CNRS, Univ. Lille, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | | | - Cécile Godé
- CNRS, Univ. Lille, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Christian Albrecht
- Department of Animal Ecology and Systematics, Justus Liebig University, D-35392 Giessen, Germany.,Department of Biology, Mbarara University of Science and Technology, Mbarara, Uganda
| | - Xavier Vekemans
- CNRS, Univ. Lille, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | | |
Collapse
|
14
|
Chambers EA, Tarvin RD, Santos JC, Ron SR, Betancourth-Cundar M, Hillis DM, Matz MV, Cannatella DC. 2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics. Ecol Evol 2023; 13:e9842. [PMID: 36911313 PMCID: PMC9994478 DOI: 10.1002/ece3.9842] [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: 01/24/2023] [Revised: 02/02/2023] [Accepted: 02/03/2023] [Indexed: 03/10/2023] Open
Abstract
Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.
Collapse
Affiliation(s)
- E Anne Chambers
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA.,Department of Environmental Science, Policy, and Management and Museum of Vertebrate Zoology University of California Berkeley Berkeley California USA
| | - Rebecca D Tarvin
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA.,Department of Integrative Biology and Museum of Vertebrate Zoology University of California Berkeley Berkeley California USA
| | - Juan C Santos
- Department of Biological Sciences St John's University New York New York USA
| | - Santiago R Ron
- Museo de Zoología, Escuela de Ciencias Biológicas Pontificia Universidad Católica del Ecuador Quito Ecuador
| | | | - David M Hillis
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA
| | - Mikhail V Matz
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA
| | - David C Cannatella
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA
| |
Collapse
|
15
|
Stull GW. Evolutionary origins of the eastern North American-Mesoamerican floristic disjunction: Current status and future prospects. AMERICAN JOURNAL OF BOTANY 2023; 110:1-11. [PMID: 36794648 DOI: 10.1002/ajb2.16142] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 01/25/2023] [Accepted: 01/25/2023] [Indexed: 05/11/2023]
Abstract
Biogeographic disjunction patterns, where multiple taxa are shared between isolated geographic areas, represent excellent systems for investigating the historical assembly of modern biotas and fundamental biological processes such as speciation, diversification, niche evolution, and evolutionary responses to climate change. Studies on plant genera disjunct across the northern hemisphere, particularly between eastern North America (ENA) and eastern Asia (EAS), have yielded tremendous insight on the geologic history and assembly of rich temperate floras. However, one of the most prevalent disjunction patterns involving ENA forests has been largely overlooked: that of taxa disjunct between ENA and cloud forests of Mesoamerica (MAM), with examples including Acer saccharum, Liquidambar styraciflua, Cercis canadensis, Fagus grandifolia, and Epifagus virginiana. Despite the remarkable nature of this disjunction pattern, which has been recognized for over 75 years, there have been few recent efforts to empirically examine its evolutionary and ecological origins. Here I synthesize previous systematic, paleobotanical, phylogenetic, and phylogeographic studies to establish what is known about this disjunction pattern to provide a roadmap for future research. I argue that this disjunction pattern, and the evolution and fossil record of the Mexican flora more broadly, represents a key missing piece in the broader puzzle of northern hemisphere biogeography. I also suggest that the ENA-MAM disjunction represents an excellent system for examining fundamental questions about how traits and life history strategies mediate plant evolutionary responses to climate change and for predicting how broadleaf temperate forests will respond to the ongoing climatic pressures of the Anthropocene.
Collapse
Affiliation(s)
- Gregory W Stull
- Department of Botany, National Museum of Natural History, Smithsonian Institution, Washington, D.C., 20013, USA
| |
Collapse
|
16
|
Johnson O, Ribas CC, Aleixo A, Naka LN, Harvey MG, Brumfield RT. Amazonian birds in more dynamic habitats have less population genetic structure and higher gene flow. Mol Ecol 2023; 32:2186-2205. [PMID: 36798996 DOI: 10.1111/mec.16886] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 02/07/2023] [Accepted: 02/09/2023] [Indexed: 02/18/2023]
Abstract
Understanding the factors that govern variation in genetic structure across species is key to the study of speciation and population genetics. Genetic structure has been linked to several aspects of life history, such as foraging strategy, habitat association, migration distance, and dispersal ability, all of which might influence dispersal and gene flow. Comparative studies of population genetic data from species with differing life histories provide opportunities to tease apart the role of dispersal in shaping gene flow and population genetic structure. Here, we examine population genetic data from sets of bird species specialized on a series of Amazonian habitat types hypothesized to filter for species with dramatically different dispersal abilities: stable upland forest, dynamic floodplain forest, and highly dynamic riverine islands. Using genome-wide markers, we show that habitat type has a significant effect on population genetic structure, with species in upland forest, floodplain forest, and riverine islands exhibiting progressively lower levels of structure. Although morphological traits used as proxies for individual-level dispersal ability did not explain this pattern, population genetic measures of gene flow are elevated in species from more dynamic riverine habitats. Our results suggest that the habitat in which a species occurs drives the degree of population genetic structuring via its impact on long-term fluctuations in levels of gene flow, with species in highly dynamic habitats having particularly elevated gene flow. These differences in genetic variation across taxa specialized in distinct habitats may lead to disparate responses to environmental change or habitat-specific diversification dynamics over evolutionary time scales.
Collapse
Affiliation(s)
- Oscar Johnson
- Department of Biological Sciences and Museum of Natural Science, Louisiana State University, Baton Rouge, Louisiana, USA
| | - Camila C Ribas
- Instituto Nacional de Pesquisas da Amazônia (INPA), Manaus, Amazonas, Brazil
| | - Alexandre Aleixo
- Museu Paraense Emílio Goeldi (MPEG), Belém, Pará, Brazil.,Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland.,Instituto Tecnológico Vale, Belém, Brazil
| | - Luciano N Naka
- Laboratório de Ecologia & Evolução de Aves, Departamento de Zoologia, Universidade Federal de Pernambuco, Recife, Brazil
| | - Michael G Harvey
- Department of Biological Sciences, University of Texas at El Paso, El Paso, Texas, USA
| | - Robb T Brumfield
- Department of Biological Sciences and Museum of Natural Science, Louisiana State University, Baton Rouge, Louisiana, USA
| |
Collapse
|
17
|
Luna LW, Naka LN, Thom G, Knowles LL, Sawakuchi AO, Aleixo A, Ribas CC. Late Pleistocene landscape changes and habitat specialization as promoters of population genomic divergence in Amazonian floodplain birds. Mol Ecol 2023; 32:214-228. [PMID: 36261866 DOI: 10.1111/mec.16741] [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: 06/16/2022] [Revised: 09/28/2022] [Accepted: 10/06/2022] [Indexed: 12/29/2022]
Abstract
Although vicariant processes are expected to leave similar genomic signatures among codistributed taxa, ecological traits such as habitat and stratum can influence genetic divergence within species. Here, we combined landscape history and habitat specialization to understand the historical and ecological factors responsible for current levels of genetic divergence in three species of birds specialized in seasonally flooded habitats in muddy rivers and which are widespread in the Amazon basin but have isolated populations in the Rio Branco. Populations of the white-bellied spinetail (Mazaria propinqua), lesser wagtail-tyrant (Stigmatura napensis) and bicolored conebill (Conirostrum bicolor) are currently isolated in the Rio Branco by the black-waters of the lower Rio Negro, offering a unique opportunity to test the effect of river colour as a barrier to gene flow. We used ultraconserved elements (UCEs) to test alternative hypotheses of population history in a comparative phylogeographical approach by modelling genetic structure, demographic history and testing for shared divergence time among codistributed taxa. Our analyses revealed that (i) all three populations from the Rio Branco floodplains are genetically distinct from other populations along the Amazon River floodplains; (ii) these divergences are the result of at least two distinct events, consistent with species habitat specialization; and (iii) the most likely model of population evolution includes lower population connectivity during the Late Pleistocene transition (~250,000 years ago), with gene flow being completely disrupted after the Last Glacial Maximum (~21,000 years ago). Our findings highlight how landscape evolution modulates population connectivity in habitat specialist species and how organisms can have different responses to the same historical processes of environmental change, depending on their habitat affinity.
Collapse
Affiliation(s)
- Leilton Willians Luna
- Department of Ecosystem Science and Management, Pennsylvania State University, University Park, Pennsylvania, USA.,Department of Zoology, Universidade Federal do Pará/Emílio Goeldi Museum, Belém, Brazil
| | - Luciano Nicolas Naka
- Laboratory of Avian Ecology and Evolution, Department of Zoology, Universidade Federal do Pernambuco, Recife, Brazil
| | - Gregory Thom
- Museum of Natural Science, Louisiana State University, Baton Rouge, Louisiana, USA
| | - Laura Lacey Knowles
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, Michigan, USA
| | | | - Alexandre Aleixo
- Department of Zoology, Universidade Federal do Pará/Emílio Goeldi Museum, Belém, Brazil.,Instituto Tecnológico Vale, Belém, Brazil.,Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland
| | | |
Collapse
|
18
|
Dietz L, Eberle J, Mayer C, Kukowka S, Bohacz C, Baur H, Espeland M, Huber BA, Hutter C, Mengual X, Peters RS, Vences M, Wesener T, Willmott K, Misof B, Niehuis O, Ahrens D. Standardized nuclear markers improve and homogenize species delimitation in Metazoa. Methods Ecol Evol 2022. [DOI: 10.1111/2041-210x.14041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Lars Dietz
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Jonas Eberle
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
- Paris‐Lodron‐University Salzburg Austria
| | - Christoph Mayer
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Sandra Kukowka
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Claudia Bohacz
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Hannes Baur
- Naturhistorisches Museum Bern/Institute of Ecology and Evolution University of Bern Bern Switzerland
| | - Marianne Espeland
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Bernhard A. Huber
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Carl Hutter
- Museum of Natural Sciences and Department of Biological Sciences Louisiana State University Baton Rouge Louisiana USA
| | - Ximo Mengual
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Ralph S. Peters
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Miguel Vences
- Technische Universität Braunschweig Braunschweig Germany
| | - Thomas Wesener
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| | - Keith Willmott
- Florida Museum of Natural History University of Florida Gainesville Florida USA
| | - Bernhard Misof
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
- Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Oliver Niehuis
- Abt. Evolutionsbiologie und Ökologie, Institut für Biologie I, Albert‐Ludwigs‐Universität Freiburg Freiburg im Breisgau Germany
| | - Dirk Ahrens
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB) Bonn Germany
| |
Collapse
|
19
|
Ma J, Deng Y, Zhang M, Yu J. The role of multi-omics in the diagnosis of COVID-19 and the prediction of new therapeutic targets. Virulence 2022; 13:1101-1110. [PMID: 35801633 PMCID: PMC9272836 DOI: 10.1080/21505594.2022.2092941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022] Open
Abstract
The global pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has led to more than 170 million confirmed cases in 223 countries and regions, claiming 3,872,457 lives. Some patients with COVID-19 have mild clinical symptoms despite severe respiratory failure, which greatly increases the difficulty of diagnosis and treatment. It is therefore necessary to identify biological characteristics of SARS-CoV-2, screen novel diagnostic and prognostic biomarkers, as well as to explore potential therapeutic targets for COVID-19. In this comprehensive review, we discuss the current published literature on COVID-19. We find that the comprehensive application of genomics, transcriptomics, proteomics and metabolomics is becoming increasingly important in the treatment of COVID-19. Multi-omics analysis platforms are expected to revolutionize the diagnosis and classification of COVID-19. This review aims to provide a reference for diagnosis, surveillance and clinical decision making related to COVID-19.
Collapse
Affiliation(s)
- Jianli Ma
- Department of Radiation Oncology, Shandong University Cancer Center, Jinan, Shandong Province, People's Republic of China
| | - Yuwei Deng
- Department of Breast Medical Oncology, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang Province, People's Republic of China
| | - Minghui Zhang
- Department of Respiratory Medical Oncology, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang Province, People's Republic of China
| | - Jinming Yu
- Department of Radiation Oncology, Shandong University Cancer Center, Jinan, Shandong Province, People's Republic of China
| |
Collapse
|
20
|
Chaturvedi S, Gompert Z, Feder JL, Osborne OG, Muschick M, Riesch R, Soria-Carrasco V, Nosil P. Climatic similarity and genomic background shape the extent of parallel adaptation in Timema stick insects. Nat Ecol Evol 2022; 6:1952-1964. [PMID: 36280782 PMCID: PMC7613875 DOI: 10.1038/s41559-022-01909-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Accepted: 09/13/2022] [Indexed: 12/15/2022]
Abstract
Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.
Collapse
Affiliation(s)
- Samridhi Chaturvedi
- Department of Integrative Biology, University of California, Berkeley, CA, USA.
- Department of Biology and Ecology Center, Utah State University, Logan, UT, USA.
| | - Zachariah Gompert
- Department of Biology and Ecology Center, Utah State University, Logan, UT, USA.
| | - Jeffrey L Feder
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, USA
| | - Owen G Osborne
- Molecular Ecology and Evolution Bangor, Environment Centre Wales, School of Natural Sciences, Bangor University, Bangor, UK
| | - Moritz Muschick
- Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
- Department of Fish Ecology and Evolution, Eawag, Swiss Federal Institute for Aquatic Science and Technology, Kastanienbaum, Switzerland
| | - Rüdiger Riesch
- Department of Biological Sciences, Royal Holloway University of London, Egham, UK
| | | | - Patrik Nosil
- Department of Biology and Ecology Center, Utah State University, Logan, UT, USA
- CEFE, Univ. Montpellier, CNRS, EPHE, IRD, Univ. Paul Valéry Montpellier 3, Montpellier, France
| |
Collapse
|
21
|
Nunes R, Storer C, Doleck T, Kawahara AY, Pierce NE, Lohman DJ. Predictors of sequence capture in a large-scale anchored phylogenomics project. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.943361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
Abstract
Next-generation sequencing (NGS) technologies have revolutionized phylogenomics by decreasing the cost and time required to generate sequence data from multiple markers or whole genomes. Further, the fragmented DNA of biological specimens collected decades ago can be sequenced with NGS, reducing the need for collecting fresh specimens. Sequence capture, also known as anchored hybrid enrichment, is a method to produce reduced representation libraries for NGS sequencing. The technique uses single-stranded oligonucleotide probes that hybridize with pre-selected regions of the genome that are sequenced via NGS, culminating in a dataset of numerous orthologous loci from multiple taxa. Phylogenetic analyses using these sequences have the potential to resolve deep and shallow phylogenetic relationships. Identifying the factors that affect sequence capture success could save time, money, and valuable specimens that might be destructively sampled despite low likelihood of sequencing success. We investigated the impacts of specimen age, preservation method, and DNA concentration on sequence capture (number of captured sequences and sequence quality) while accounting for taxonomy and extracted tissue type in a large-scale butterfly phylogenomics project. This project used two probe sets to extract 391 loci or a subset of 13 loci from over 6,000 butterfly specimens. We found that sequence capture is a resilient method capable of amplifying loci in samples of varying age (0–111 years), preservation method (alcohol, papered, pinned), and DNA concentration (0.020 ng/μl - 316 ng/ul). Regression analyses demonstrate that sequence capture is positively correlated with DNA concentration. However, sequence capture and DNA concentration are negatively correlated with sample age and preservation method. Our findings suggest that sequence capture projects should prioritize the use of alcohol-preserved samples younger than 20 years old when available. In the absence of such specimens, dried samples of any age can yield sequence data, albeit with returns that diminish with increasing age.
Collapse
|
22
|
Spies I, Tarpey C, Kristiansen T, Fisher M, Rohan S, Hauser L. Genomic differentiation in Pacific cod using
P
ool‐
S
eq. Evol Appl 2022; 15:1907-1924. [PMID: 36426128 PMCID: PMC9679252 DOI: 10.1111/eva.13488] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Revised: 09/05/2022] [Accepted: 09/21/2022] [Indexed: 11/28/2022] Open
Abstract
Patterns of genetic differentiation across the genome can provide insight into selective forces driving adaptation. We used pooled whole genome sequencing, gene annotation, and environmental covariates to evaluate patterns of genomic differentiation and to investigate mechanisms responsible for divergence among proximate Pacific cod (Gadus macrocephalus) populations from the Bering Sea and Aleutian Islands and more distant Washington Coast cod. Samples were taken from eight spawning locations, three of which were replicated to estimate consistency in allele frequency estimation. A kernel smoothing moving weighted average of relative divergence (FST) identified 11 genomic islands of differentiation between the Aleutian Islands and Bering Sea samples. In some islands of differentiation, there was also elevated absolute divergence (dXY) and evidence for selection, despite proximity and potential for gene flow. Similar levels of absolute divergence (dXY) but roughly double the relative divergence (FST) were observed between the distant Bering Sea and Washington Coast samples. Islands of differentiation were much smaller than the four large inversions among Atlantic cod ecotypes. Islands of differentiation between the Bering Sea and Aleutian Island were associated with SNPs from five vision system genes, which can be associated with feeding, predator avoidance, orientation, and socialization. We hypothesize that islands of differentiation between Pacific cod from the Bering Sea and Aleutian Islands provide evidence for adaptive differentiation despite gene flow in this commercially important marine species.
Collapse
Affiliation(s)
- Ingrid Spies
- Resource Ecology and Fisheries Management Division Alaska Fisheries Science Center Seattle Washington USA
| | - Carolyn Tarpey
- School of Aquatic and Fishery Sciences University of Washington Seattle Washington USA
| | | | - Mary Fisher
- School of Aquatic and Fishery Sciences University of Washington Seattle Washington USA
| | - Sean Rohan
- Resource Assessment and Conservation Engineering Division Alaska Fisheries Science Center Seattle Washington USA
| | - Lorenz Hauser
- Resource Ecology and Fisheries Management Division Alaska Fisheries Science Center Seattle Washington USA
| |
Collapse
|
23
|
Pacheco C, Lobo D, Silva P, Álvares F, García EJ, Castro D, Layna JF, López-Bao JV, Godinho R. Assessing the performance of historical skins and bones for museomics using wolf specimens as a case study. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.970249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Advances in the field of museomics have promoted a high sampling demand for natural history collections (NHCs), eventually resulting in damage to invaluable resources to understand historical biodiversity. It is thus essential to achieve a consensus about which historical tissues present the best sources of DNA. In this study, we evaluated the performance of different historical tissues from Iberian wolf NHCs in genome-wide assessments. We targeted three tissues—bone (jaw and femur), maxilloturbinal bone, and skin—that have been favored by traditional taxidermy practices for mammalian carnivores. Specifically, we performed shotgun sequencing and target capture enrichment for 100,000 single nucleotide polymorphisms (SNPs) selected from the commercial Canine HD BeadChip across 103 specimens from 1912 to 2005. The performance of the different tissues was assessed using metrics based on endogenous DNA content, uniquely high-quality mapped reads after capture, and enrichment proportions. All samples succeeded as DNA sources, regardless of their collection year or sample type. Skin samples yielded significantly higher amounts of endogenous DNA compared to both bone types, which yielded equivalent amounts. There was no evidence for a direct effect of tissue type on capture efficiency; however, the number of genotyped SNPs was strictly associated with the starting amount of endogenous DNA. Evaluation of genotyping accuracy for distinct minimum read depths across tissue types showed a consistent overall low genotyping error rate (<7%), even at low (3x) coverage. We recommend the use of skins as reliable and minimally destructive sources of endogenous DNA for whole-genome and target enrichment approaches in mammalian carnivores. In addition, we provide a new 100,000 SNP capture array validated for historical DNA (hDNA) compatible to the Canine HD BeadChip for high-quality DNA. The increasing demand for NHCs as DNA sources should encourage the generation of genomic datasets comparable among studies.
Collapse
|
24
|
Penaud B, Laurent B, Milhes M, Noüs C, Ehrenmann F, Dutech C. SNP4OrphanSpecies: A bioinformatics pipeline to isolate molecular markers for studying genetic diversity of orphan species. Biodivers Data J 2022; 10:e85587. [PMID: 36761595 PMCID: PMC9848450 DOI: 10.3897/bdj.10.e85587] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 06/23/2022] [Indexed: 11/12/2022] Open
Abstract
Background For several decades, an increase in disease or pest emergences due to anthropogenic introduction or environmental changes has been recorded. This increase leads to serious threats to the genetic and species diversity of numerous ecosystems. Many of these events involve species with poor or no genomic resources (called here "orphan species"). This lack of resources is a serious limitation to our understanding of the origin of emergent populations, their ability to adapt to new environments and to predict future consequences to biodiversity. Analyses of genetic diversity are an efficient method to obtain this information rapidly, but require available polymorphic genetic markers. New information We developed a generic bioinformatics pipeline to rapidly isolate such markers with the goal for the pipeline to be applied in studies of invasive taxa from different taxonomic groups, with a special focus on forest fungal pathogens and insect pests. This pipeline is based on: 1) an automated de novo genome assembly obtained from shotgun whole genome sequencing using paired-end Illumina technology; 2) the isolation of single-copy genes conserved in species related to the studied emergent organisms; 3) primer development for multiplexed short sequences obtained from these conserved genes. Previous studies have shown that intronic regions of these conserved genes generally contain several single nucleotide polymorphisms within species. The pipeline's functionality was evaluated with sequenced genomes of five invasive or expanding pathogen and pest species in Europe (Armillariaostoyae (Romagn.) Herink 1973, Bursaphelenchusxylophilus Steiner & Buhrer 1934, Sphaeropsissapinea (fr.) Dicko & B. Sutton 1980, Erysiphealphitoides (Griffon & Maubl.) U. Braun & S. Takam. 2000, Thaumetopoeapityocampa Denis & Schiffermüller, 1775). We successfully isolated several pools of one hundred short gene regions for each assembled genome, which can be amplified in multiplex. The bioinformatics pipeline is user-friendly and requires little computational resources. This easy-to-set-up and run method for genetic marker identification will be useful for numerous laboratories studying biological invasions, but with limited resources and expertise in bioinformatics.
Collapse
Affiliation(s)
- Benjamin Penaud
- BIOGECO, INRAE, Univ. Bordeaux, 33610 Cestas, FranceBIOGECO, INRAE, Univ. Bordeaux33610 CestasFrance
| | - Benoit Laurent
- BIOGECO, INRAE, Univ. Bordeaux, 33610 Cestas, FranceBIOGECO, INRAE, Univ. Bordeaux33610 CestasFrance
| | - Marine Milhes
- INRAE, US 1426, GeT-PlaGe, Genotoul, Castanet-Tolosan, FranceINRAE, US 1426, GeT-PlaGe, GenotoulCastanet-TolosanFrance
| | - Camille Noüs
- Laboratoire Cogitamus, Bordeaux, FranceLaboratoire CogitamusBordeauxFrance
| | - François Ehrenmann
- BIOGECO, INRAE, Univ. Bordeaux, 33610 Cestas, FranceBIOGECO, INRAE, Univ. Bordeaux33610 CestasFrance
| | - Cyril Dutech
- BIOGECO, INRAE, Univ. Bordeaux, 33610 Cestas, FranceBIOGECO, INRAE, Univ. Bordeaux33610 CestasFrance
| |
Collapse
|
25
|
Smith BT, Merwin J, Provost KL, Thom G, Brumfield RT, Ferreira M, Mauck Iii WM, Moyle RG, Wright T, Joseph L. Phylogenomic analysis of the parrots of the world distinguishes artifactual from biological sources of gene tree discordance. Syst Biol 2022; 72:228-241. [PMID: 35916751 DOI: 10.1093/sysbio/syac055] [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: 04/26/2021] [Revised: 02/22/2022] [Accepted: 07/22/2022] [Indexed: 11/14/2022] Open
Abstract
Gene tree discordance is expected in phylogenomic trees and biological processes are often invoked to explain it. However, heterogeneous levels of phylogenetic signal among individuals within datasets may cause artifactual sources of topological discordance. We examined how the information content in tips and subclades impacts topological discordance in the parrots (Order: Psittaciformes), a diverse and highly threatened clade of nearly 400 species. Using ultraconserved elements from 96% of the clade's species-level diversity, we estimated concatenated and species trees for 382 ingroup taxa. We found that discordance among tree topologies was most common at nodes dating between the late Miocene and Pliocene, and often at the taxonomic level of genus. Accordingly, we used two metrics to characterize information content in tips and assess the degree to which conflict between trees was being driven by lower quality samples. Most instances of topological conflict and non-monophyletic genera in the species tree could be objectively identified using these metrics. For subclades still discordant after tip-based filtering, we used a machine learning approach to determine whether phylogenetic signal or noise was the more important predictor of metrics supporting the alternative topologies. We found that when signal favored one of the topologies, noise was the most important variable in poorly performing models that favored the alternative topology. In sum, we show that artifactual sources of gene tree discordance, which are likely a common phenomenon in many datasets, can be distinguished from biological sources by quantifying the information content in each tip and modeling which factors support each topology.
Collapse
Affiliation(s)
- Brian Tilston Smith
- Department of Ornithology, American Museum of Natural History, Central Park West at 79th Street, New York, NY 10024, USA
| | - Jon Merwin
- Department of Ornithology, Academy of Natural Sciences of Drexel University, 1900 Benjamin Franklin Parkway, Philadelphia, PA 19103, USA.,Department of Biodiversity, Earth, and Environmental Science, Drexel University, Philadelphia, PA 19103, USA
| | - Kaiya L Provost
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, 318 W. 12th Avenue, Columbus, OH 43210, USA
| | - Gregory Thom
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Robb T Brumfield
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Mateus Ferreira
- Centro de Estudos da Biodiversidade, Universidade Federal de Roraima, Av. Cap. Ene Garcez, 2413, Boa Vista, RR, Brazil
| | - William M Mauck Iii
- Department of Ornithology, American Museum of Natural History, Central Park West at 79th Street, New York, NY 10024, USA
| | - Robert G Moyle
- Department of Ecology and Evolutionary Biology and Biodiversity Institute, University of Kansas, 1345 Jayhawk Blvd., Lawrence, KS 66045, USA
| | - Timothy Wright
- Department of Biology, New Mexico State University, Las Cruces, NM, 88003, USA
| | - Leo Joseph
- Australian National Wildlife Collection, National Research Collections Australia, CSIRO, GPO Box 1700, Canberra, ACT, 2601, Australia
| |
Collapse
|
26
|
Zhou W, Jenny Xiang QY. Phylogenomics and Biogeography of Castanea (Chestnut) and Hamamelis (Witch-hazel) - Choosing between RAD-seq and Hyb-Seq Approaches. Mol Phylogenet Evol 2022; 176:107592. [DOI: 10.1016/j.ympev.2022.107592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2021] [Revised: 06/18/2022] [Accepted: 07/20/2022] [Indexed: 10/31/2022]
|
27
|
Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations. ORG DIVERS EVOL 2022. [DOI: 10.1007/s13127-022-00568-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Abstract
Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.
Collapse
|
28
|
Salter JF, Hosner PA, Tsai WLE, McCormack JE, Braun EL, Kimball RT, Brumfield RT, Faircloth BC. Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae). Mol Phylogenet Evol 2022; 175:107559. [PMID: 35803448 DOI: 10.1016/j.ympev.2022.107559] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 05/19/2022] [Accepted: 05/26/2022] [Indexed: 01/22/2023]
Abstract
As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.
Collapse
Affiliation(s)
- Jessie F Salter
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA.
| | - Peter A Hosner
- Natural History Museum of Denmark, Center for Global Mountain Biodiversity, and Center for Macroecology, Evolution, and Climate, University of Copenhagen, Copenhagen, Denmark; Department of Biology, University of Florida, Gainesville, FL, USA
| | - Whitney L E Tsai
- Moore Laboratory of Biology, Occidental College, Los Angeles, CA, USA
| | - John E McCormack
- Moore Laboratory of Biology, Occidental College, Los Angeles, CA, USA; Biology Department, Occidental College, Los Angeles, CA, USA
| | - Edward L Braun
- Department of Biology, University of Florida, Gainesville, FL, USA
| | | | - Robb T Brumfield
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA
| | - Brant C Faircloth
- Museum of Natural Science and Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA
| |
Collapse
|
29
|
Ciaccio E, Debray A, Hedin M. Phylogenomics of paleoendemic lampshade spiders (Araneae, Hypochilidae, Hypochilus), with the description of a new species from montane California. Zookeys 2022; 1086:163-204. [PMID: 35221748 PMCID: PMC8873193 DOI: 10.3897/zookeys.1086.77190] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Accepted: 01/18/2022] [Indexed: 12/31/2022] Open
Abstract
Hypochilus is a relictual lineage of Nearctic spiders distributed disjunctly across the United States in three montane regions (California, southern Rocky Mountains, southern Appalachia). Phylogenetic resolution of species relationships in Hypochilus has been challenging, and conserved morphology coupled with extreme genetic divergence has led to uncertain species limits in some complexes. Here, Hypochilus interspecies relationships have been reconstructed and cryptic speciation more critically evaluated using a combination of ultraconserved elements, mitochondrial CO1 by-catch, and morphology. Phylogenomic data strongly support the monophyly of regional clades and support a ((California, Appalachia), southern Rocky Mountains) topology. In Appalachia, five species are resolved as four lineages (H.thorelli Marx, 1888 and H.coylei Platnick, 1987 are clearly sister taxa), but the interrelationships of these four lineages remain unresolved. The Appalachian species H.pococki Platnick, 1987 is recovered as monophyletic but is highly genetically structured at the nuclear level. While algorithmic analyses of nuclear data indicate many species (e.g., all H.pococki populations as species), male morphology instead reveals striking stasis. Within the California clade, nuclear and mitochondrial lineages of H.petrunkevitchi Gertsch, 1958 correspond directly to drainage basins of the southern Sierra Nevada, with H.bernardino Catley, 1994 nested within H.petrunkevitchi and sister to the southernmost basin populations. Combining nuclear, mitochondrial, geographical, and morphological evidence a new species from the Tule River and Cedar Creek drainages is described, Hypochilusxomotesp. nov. We also emphasize the conservation issues that face several microendemic, habitat-specialized species in this remarkable genus.
Collapse
Affiliation(s)
- Erik Ciaccio
- Department of Biology, San Diego State University, San Diego, California, USA San Diego State University San Diego United States of America.,Department of Entomology, Plant Pathology and Nematology, University of Idaho, Idaho, USA University of Idaho Idaho United States of America
| | - Andrew Debray
- Department of Biology, San Diego State University, San Diego, California, USA San Diego State University San Diego United States of America.,Nano PharmaSolutions Inc., San Diego, California, USA Nano PharmaSolutions Inc. San Diego United States of America
| | - Marshal Hedin
- Department of Biology, San Diego State University, San Diego, California, USA San Diego State University San Diego United States of America
| |
Collapse
|
30
|
Incompatibility Phylogenetic Signals between Double-Digest Restriction Site-Associated DNA Sequencing and Plastid Genomes in Chinese Curcuma (Zingiberaceae)—A Recent Qinghai–Tibetan Plateau Diversification Genera. FORESTS 2022. [DOI: 10.3390/f13020280] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Curcuma is of high economic value, credited to its medicinal, edible, and ornamental properties, which possess all signatures of adaptability, and rapid radiation, especially species of Curcuma (Chinese Curcuma, a recent Qinghai–Tibetan Plateau diversification genera) scattered in China. However, little is known about the incongruent phylogenetic signals within this genera from different inheritance patterns that will militate against the further development of this genera. In this research, we applied complete chloroplast genome data together with double-digest restriction site-associated DNA sequencing data (ddRAD-seq) strategy to investigate phylogenetic signals of Chinese Curcuma species, clustering using two RAD analysis pipelines (STACKS and pyRAD). Phylogenetic trees were obtained from each locus based on the maximum likelihood (ML) and multispecies coalescent (BEAST) methods. For visual comparison, multi-method and different datasets were used to infer the phylogeny. We discovered inconsistent relationships for the Chinese Curcuma with varying degrees of support using different methods and datasets.
Collapse
|
31
|
McLean BS, Bell KC, Cook JA. SNP-based Phylogenomic Inference in Holarctic Ground Squirrels (Urocitellus). Mol Phylogenet Evol 2022; 169:107396. [PMID: 35031463 DOI: 10.1016/j.ympev.2022.107396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Revised: 12/02/2021] [Accepted: 12/08/2021] [Indexed: 11/24/2022]
Abstract
Resolution of rapid evolutionary radiations requires harvesting maximal signal from phylogenomic datasets. However, studies of non-model clades often target conserved loci that are characterized by reduced information content, which can negatively affect gene tree precision and species tree accuracy. Single nucleotide polymorphism (SNP)-based methods are an underutilized but potentially valuable tool for estimating phylogeny and divergence times because they do not rely on resolved gene trees, allowing information from many or all variant loci to be leveraged in species tree reconstruction. We evaluated the utility of SNP-based methods in resolving phylogeny of Holarctic ground squirrels (Urocitellus), a radiation that has been difficult to disentangle, even in prior phylogenomic studies. We inferred phylogeny from a dataset of >3,000 ultraconserved element loci (UCEs) using two methods (SNAPP, SVDquartets) and compared our results with a new mitogenome phylogeny. We also systematically evaluated how phasing of UCEs improves per-locus information content, and inference of topology and other parameters within each of these SNP-based methods. Phasing improved topological resolution and branch length estimation at shallow levels (within species complexes), but less so at deeper levels, likely reflecting true uncertainty due to ancestral polymorphisms segregating in these rapidly diverging lineages. We resolved several key clades in Urocitellus and present targeted opportunities for future phylogenomic inquiry. Our results extend the roadmap for use of SNPs to address vertebrate radiations and support comparative analyses at multiple temporal scales.
Collapse
Affiliation(s)
- Bryan S McLean
- University of North Carolina Greensboro, Department of Biology, Greensboro, NC 27402 USA.
| | - Kayce C Bell
- Natural History Museum of Los Angeles County, Department of Mammalogy, Los Angeles, CA 90007 USA.
| | - Joseph A Cook
- University of New Mexico, Department of Biology and Museum of Southwestern Biology, Albuquerque, NM 87131 USA.
| |
Collapse
|
32
|
Ribeiro CVR, Oliveira LP, Batista R, De Sousa M. UCEasy: A software package for automating and simplifying the analysis of ultraconserved elements (UCEs). Biodivers Data J 2021; 9:e78132. [PMID: 34934383 PMCID: PMC8683391 DOI: 10.3897/bdj.9.e78132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Accepted: 12/09/2021] [Indexed: 11/25/2022] Open
Abstract
Background The use of Ultraconserved Elements (UCEs) as genetic markers in phylogenomics has become popular and has provided promising results. Although UCE data can be easily obtained from targeted enriched sequencing, the protocol for in silico analysis of UCEs consist of the execution of heterogeneous and complex tools, a challenge for scientists without training in bioinformatics. Developing tools with the adoption of best practices in research software can lessen this problem by improving the execution of computational experiments, thus promoting better reproducibility. New information We present UCEasy, an easy-to-install and easy-to-use software package with a simple command line interface that facilitates the computational analysis of UCEs from sequencing samples, following the best practices of research software. UCEasy is a wrapper that standardises, automates and simplifies the quality control of raw reads, assembly and extraction and alignment of UCEs, generating at the end a data matrix with different levels of completeness that can be used to infer phylogenetic trees. We demonstrate the functionalities of UCEasy by reproducing the published results of phylogenomic studies of the bird genus Turdus (Aves) and of Adephaga families (Coleoptera) containing genomic datasets to efficiently extract UCEs.
Collapse
Affiliation(s)
- Caio V R Ribeiro
- Coordenação de Ciência da Computação, Centro Universitário do Estado do Pará (CESUPA), Belém, Brazil Coordenação de Ciência da Computação, Centro Universitário do Estado do Pará (CESUPA) Belém Brazil
| | - Lucas P Oliveira
- Instituto de Computação, Universidade Estadual de Campinas (UNICAMP), Campinas, Brazil Instituto de Computação, Universidade Estadual de Campinas (UNICAMP) Campinas Brazil
| | - Romina Batista
- Instituto Nacional de Pesquisas da Amazônia (INPA), Manaus, Brazil Instituto Nacional de Pesquisas da Amazônia (INPA) Manaus Brazil.,Gothenburg Global Biodiversity Centre, Gothenburg, Sweden Gothenburg Global Biodiversity Centre Gothenburg Sweden
| | - Marcos De Sousa
- Museu Paraense Emílio Goeldi (MPEG), Belém, Brazil Museu Paraense Emílio Goeldi (MPEG) Belém Brazil.,Coordenação de Ciência da Computação, Centro Universitário do Estado do Pará (CESUPA), Belém, Brazil Coordenação de Ciência da Computação, Centro Universitário do Estado do Pará (CESUPA) Belém Brazil
| |
Collapse
|
33
|
Dong X, Yi W, Zheng C, Zhu X, Wang S, Xue H, Ye Z, Bu W. Species delimitation of rice seed bugs complex: Insights from mitochondrial genomes and ddRAD‐seq data. ZOOL SCR 2021. [DOI: 10.1111/zsc.12523] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Xue Dong
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| | - Wenbo Yi
- Department of Biology Xinzhou Teachers University Xinzhou China
| | - Chenguang Zheng
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| | - Xiuxiu Zhu
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| | - Shujing Wang
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| | - Huaijun Xue
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| | - Zhen Ye
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| | - Wenjun Bu
- Institute of Entomology College of Life Sciences Nankai University Tianjin China
| |
Collapse
|
34
|
Loera-Sánchez M, Studer B, Kölliker R. A multispecies amplicon sequencing approach for genetic diversity assessments in grassland plant species. Mol Ecol Resour 2021; 22:1725-1745. [PMID: 34918474 PMCID: PMC9305562 DOI: 10.1111/1755-0998.13577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 12/08/2021] [Accepted: 12/10/2021] [Indexed: 11/30/2022]
Abstract
Grasslands are widespread and economically relevant ecosystems at the basis of sustainable roughage production. Plant genetic diversity (PGD; i.e., within‐species diversity) is related to many beneficial effects on the ecosystem functioning of grasslands. The monitoring of PGD in temperate grasslands is complicated by the multiplicity of species present and by a shortage of methods for large‐scale assessments. However, the continuous advancement of high‐throughput DNA sequencing approaches has improved the prospects of broad, multispecies PGD monitoring. Among them, amplicon sequencing stands out as a robust and cost‐effective method. Here, we report a set of 12 multispecies primer pairs that can be used for high‐throughput PGD assessments in multiple grassland plant species. The target loci were selected and tested in two phases: a “discovery phase” based on a sequence capture assay (611 nuclear loci assessed in 16 grassland plant species), which resulted in the selection of 11 loci; and a “validation phase”, in which the selected loci were targeted and sequenced using multispecies primers in test populations of Dactylis glomerata L., Lolium perenne L., Festuca pratensis Huds., Trifolium pratense L. and T. repens L. The multispecies amplicons had nucleotide diversities per species from 5.19 × 10−3 to 1.29 × 10−2, which is in the range of flowering‐related genes but slightly lower than pathogen resistance genes. We conclude that the methodology, the DNA sequence resources, and the primer pairs reported in this study provide the basis for large‐scale, multispecies PGD monitoring in grassland plants.
Collapse
Affiliation(s)
- Miguel Loera-Sánchez
- Molecular Plant Breeding, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
| | - Bruno Studer
- Molecular Plant Breeding, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
| | - Roland Kölliker
- Molecular Plant Breeding, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
| |
Collapse
|
35
|
Nürnberger B, Baird SJE, Čížková D, Bryjová A, Mudd AB, Blaxter ML, Szymura JM. A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata). G3 (BETHESDA, MD.) 2021; 11:6353606. [PMID: 34849761 PMCID: PMC8664441 DOI: 10.1093/g3journal/jkab286] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/16/2021] [Accepted: 08/16/2021] [Indexed: 12/20/2022]
Abstract
Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.
Collapse
Affiliation(s)
- Beate Nürnberger
- Research Facility Studenec, Institute of Vertebrate Biology, Czech Academy of Sciences, 603 65 Brno, Czech Republic
| | - Stuart J E Baird
- Research Facility Studenec, Institute of Vertebrate Biology, Czech Academy of Sciences, 603 65 Brno, Czech Republic
| | - Dagmar Čížková
- Research Facility Studenec, Institute of Vertebrate Biology, Czech Academy of Sciences, 603 65 Brno, Czech Republic
| | - Anna Bryjová
- Research Facility Studenec, Institute of Vertebrate Biology, Czech Academy of Sciences, 603 65 Brno, Czech Republic
| | - Austin B Mudd
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, 94720 CA, USA
| | - Mark L Blaxter
- Tree of Life Programme, Wellcome Sanger Institute, Hinxton, Cambridge CB10 1SA, UK
| | - Jacek M Szymura
- Department of Comparative Anatomy, Jagiellonian University, 30-387 Kraków, Poland
| |
Collapse
|
36
|
Nikolakis ZL, Orton RW, Crother BI. Fine‐scale population structure within an Eastern Nearctic snake complex (
Pituophis melanoleucus
). ZOOL SCR 2021. [DOI: 10.1111/zsc.12522] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Zachary L. Nikolakis
- Department of Biology Southeastern Louisiana University Hammond Louisina USA
- Department of Biology University of Texas at Arlington Arlington Texas USA
| | - Richard W. Orton
- Department of Biology University of Texas at Arlington Arlington Texas USA
| | - Brian I. Crother
- Department of Biology Southeastern Louisiana University Hammond Louisina USA
| |
Collapse
|
37
|
How challenging RADseq data turned out to favor coalescent-based species tree inference. A case study in Aichryson (Crassulaceae). Mol Phylogenet Evol 2021; 167:107342. [PMID: 34785384 DOI: 10.1016/j.ympev.2021.107342] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Revised: 07/05/2021] [Accepted: 10/29/2021] [Indexed: 12/24/2022]
Abstract
Analysing multiple genomic regions while incorporating detection and qualification of discordance among regions has become standard for understanding phylogenetic relationships. In plants, which usually have comparatively large genomes, this is feasible by the combination of reduced-representation library (RRL) methods and high-throughput sequencing enabling the cost effective acquisition of genomic data for thousands of loci from hundreds of samples. One popular RRL method is RADseq. A major disadvantage of established RADseq approaches is the rather short fragment and sequencing range, leading to loci of little individual phylogenetic information. This issue hampers the application of coalescent-based species tree inference. The modified RADseq protocol presented here targets ca. 5,000 loci of 300-600nt length, sequenced with the latest short-read-sequencing (SRS) technology, has the potential to overcome this drawback. To illustrate the advantages of this approach we use the study group Aichryson Webb & Berthelott (Crassulaceae), a plant genus that diversified on the Canary Islands. The data analysis approach used here aims at a careful quality control of the long loci dataset. It involves an informed selection of thresholds for accurate clustering, a thorough exploration of locus properties, such as locus length, coverage and variability, to identify potential biased data and a comparative phylogenetic inference of filtered datasets, accompanied by an evaluation of resulting BS support, gene and site concordance factor values, to improve overall resolution of the resulting phylogenetic trees. The final dataset contains variable loci with an average length of 373nt and facilitates species tree estimation using a coalescent-based summary approach. Additional improvements brought by the approach are critically discussed.
Collapse
|
38
|
Thom G, Gehara M, Smith BT, Miyaki CY, do Amaral FR. Microevolutionary dynamics show tropical valleys are deeper for montane birds of the Atlantic Forest. Nat Commun 2021; 12:6269. [PMID: 34725329 PMCID: PMC8560783 DOI: 10.1038/s41467-021-26537-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Accepted: 10/08/2021] [Indexed: 11/18/2022] Open
Abstract
Tropical mountains hold more biodiversity than their temperate counterparts, and this disparity is often associated with the latitudinal climatic gradient. However, distinguishing the impact of latitude versus the background effects of species history and traits is challenging due to the evolutionary distance between tropical and temperate assemblages. Here, we test whether microevolutionary processes are linked to environmental variation across a sharp latitudinal transition in 21 montane birds of the southern Atlantic Forest in Brazil. We find that effective dispersal within populations in the tropical mountains is lower and genomic differentiation is better predicted by the current environmental complexity of the region than within the subtropical populations. The concordant response of multiple co-occurring populations is consistent with spatial climatic variability as a major process driving population differentiation. Our results provide evidence for how a narrow latitudinal gradient can shape microevolutionary processes and contribute to broader scale biodiversity patterns. There are many hypotheses for why the tropics are more biodiverse than higher latitudes. Phylogenomic analyses of 21 montane birds finds that tropical birds disperse less and have more genetically structured populations than their counterparts at higher latitudes, possibly due to a larger elevational climate gradient in the tropics
Collapse
Affiliation(s)
- Gregory Thom
- Department of Ornithology, American Museum of Natural History, Central Park West at 79th Street, New York, NY, 10024, USA. .,Departamento de Genética e Biologia Evolutiva, Universidade de São Paulo, Rua do Matão, 277, Cidade Universitária, São Paulo, SP, 05508-090, Brazil.
| | - Marcelo Gehara
- Sackler Institute for Comparative Genomics, American Museum of Natural History, New York, NY, 10024, USA.,Department of Earth and Environmental Sciences, Rutgers University Newark, 195 University Ave, Newark, NJ, 07102, USA
| | - Brian Tilston Smith
- Department of Ornithology, American Museum of Natural History, Central Park West at 79th Street, New York, NY, 10024, USA
| | - Cristina Y Miyaki
- Departamento de Genética e Biologia Evolutiva, Universidade de São Paulo, Rua do Matão, 277, Cidade Universitária, São Paulo, SP, 05508-090, Brazil
| | - Fábio Raposo do Amaral
- Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Rua Prof. Artur Riedel, 275, Jardim Eldorado, Diadema, SP, CEP 09972-270, Brazil
| |
Collapse
|
39
|
Ortiz D, Pekár S, Dianat M. Phylogenomics and loci dropout patterns of deeply diverged Zodarion ant-eating spiders suggest a high potential of RAD-seq for genus-level spider phylogenetics. Cladistics 2021; 38:320-334. [PMID: 34699083 DOI: 10.1111/cla.12493] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/02/2021] [Indexed: 11/28/2022] Open
Abstract
RAD sequencing yields large amounts of genome-wide data at a relatively low cost and without requiring previous taxon-specific information, making it ideal for evolutionary studies of highly diversified and neglected organisms. However, concerns about information decay with phylogenetic distance have discouraged its use for assessing supraspecific relationships. Here, using Double Digest Restriction Associated DNA (ddRAD) data, we perform the first deep-level approach to the phylogeny of Zodarion, a highly diversified spider genus. We explore the impact of loci and taxon filtering across concatenated and multispecies coalescent reconstruction methods and investigate the patterns of information dropout in reference to both the time of divergence and the mitochondrial divergence between taxa. We found that relaxed loci-filtering and nested taxon-filtering strategies maximized the amount of molecular information and improved phylogenetic inference. As expected, there was a clear pattern of allele dropout towards deeper time and mitochondrial divergences, but the phylogenetic signal remained strong throughout the phylogeny. Therefore, we inferred topologies that were almost fully resolved, highly supported, and noticeably congruent between setups and inference methods, which highlights overall inconsistency in the taxonomy of Zodarion. Because Zodarion appears to be among the oldest and most mitochondrially diversified spider genera, our results suggest that ddRAD data show high potential for inferring intra-generic relationships across spiders and probably also in other taxonomic groups.
Collapse
Affiliation(s)
- David Ortiz
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czechia
| | - Stano Pekár
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czechia
| | - Malahat Dianat
- Department of Botany and Zoology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czechia
| |
Collapse
|
40
|
Perez MF, Bonatelli IAS, Romeiro-Brito M, Franco FF, Taylor NP, Zappi DC, Moraes EM. Coalescent-based species delimitation meets deep learning: Insights from a highly fragmented cactus system. Mol Ecol Resour 2021; 22:1016-1028. [PMID: 34669256 DOI: 10.1111/1755-0998.13534] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2021] [Revised: 09/16/2021] [Accepted: 10/12/2021] [Indexed: 11/26/2022]
Abstract
Delimiting species boundaries is a major goal in evolutionary biology. An increasing volume of literature has focused on the challenges of investigating cryptic diversity within complex evolutionary scenarios of speciation, including gene flow and demographic fluctuations. New methods based on model selection, such as approximate Bayesian computation, approximate likelihoods, and machine learning are promising tools arising in this field. Here, we introduce a framework for species delimitation using the multispecies coalescent model coupled with a deep learning algorithm based on convolutional neural networks (CNNs). We compared this strategy with a similar ABC approach. We applied both methods to test species boundary hypotheses based on current and previous taxonomic delimitations as well as genetic data (sequences from 41 loci) in Pilosocereus aurisetus, a cactus species complex with a sky-island distribution and taxonomic uncertainty. To validate our method, we also applied the same strategy on data from widely accepted species from the genus Drosophila. The results show that our CNN approach has a high capacity to distinguish among the simulated species delimitation scenarios, with higher accuracy than ABC. For the cactus data set, a splitter hypothesis without gene flow showed the highest probability in both CNN and ABC approaches, a result agreeing with previous taxonomic classifications and in line with the sky-island distribution and low dispersal of P. aurisetus. Our results highlight the cryptic diversity within the P. aurisetus complex and show that CNNs are a promising approach for distinguishing complex evolutionary histories, even outperforming the accuracy of other model-based approaches such as ABC.
Collapse
Affiliation(s)
- Manolo F Perez
- Departamento de Biologia, Universidade Federal de São Carlos, Sorocaba, Brazil.,Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, Brazil
| | - Isabel A S Bonatelli
- Departamento de Biologia, Universidade Federal de São Carlos, Sorocaba, Brazil.,Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Diadema, Brazil
| | | | - Fernando F Franco
- Departamento de Biologia, Universidade Federal de São Carlos, Sorocaba, Brazil
| | | | - Daniela C Zappi
- Programa de Pós Graduação em Botânica, Instituto de Ciências Biológicas, Universidade de Brasília, Brasília, Brazil
| | - Evandro M Moraes
- Departamento de Biologia, Universidade Federal de São Carlos, Sorocaba, Brazil
| |
Collapse
|
41
|
Silva SM, Ribas CC, Aleixo A. Recent population differentiation in the habitat specialist Glossy Antshrike (Aves: Thamnophilidae) across Amazonian seasonally flooded forests. Ecol Evol 2021; 11:11826-11838. [PMID: 34522344 PMCID: PMC8427616 DOI: 10.1002/ece3.7951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Revised: 07/11/2021] [Accepted: 07/13/2021] [Indexed: 11/30/2022] Open
Abstract
We assessed population structure and the spatio-temporal pattern of diversification in the Glossy Antshrike Sakesphorus luctuosus (Aves, Thamnophilidae) to understand the processes shaping the evolutionary history of Amazonian floodplains and address unresolved taxonomic controversies surrounding its species limits. By targeting ultraconserved elements (UCEs) from 32 specimens of S. luctuosus, we identified independent lineages and estimated their differentiation, divergence times, and migration rates. We also estimated current and past demographic histories for each recovered lineage. We found evidence confirming that S. luctuosus consists of a single species, comprising at least four populations, with some highly admixed individuals and overall similar levels of migration between populations. We confirmed the differentiation of the Araguaia River basin population (S. l. araguayae) and gathered circumstantial evidence indicating that the taxon S. hagmanni may represent a highly introgressed population between three distinct phylogroups of S. luctuosus. Divergences between populations occurred during the last 1.2 mya. Signs of population expansions were detected for populations attributed to subspecies S. l. luctuosus, but not for the S. l. araguayae population. Our results support that S. luctuosus has had a complex population history, resulting from a high dependence on southeastern "clear water" seasonally flooded habitats and their availability through time. Spatial and demographic expansions toward the western "white water" flooded forests might be related to recent changes in connectivity and availability of these habitats. Our study reinforces the view that isolation due to absence of suitable habitat has been an important driver of population differentiation within Amazonian flooded forests, but also that differences between várzeas ("white water" floodplains, mostly in southwestern Amazonia) and igapós ("clear water" floodplains, especially located in the east) should be further explored as drivers of micro-evolution for terrestrial species.
Collapse
Affiliation(s)
- Sofia Marques Silva
- Research Centre in Biodiversity and Genetic ResourcesCIBIO/InBIOVairãoPortugal
- Department of ZoologyMuseu Paraense Emílio GoeldiBelémBrazil
| | | | - Alexandre Aleixo
- Department of ZoologyMuseu Paraense Emílio GoeldiBelémBrazil
- Finnish Museum of Natural HistoryUniversity of HelsinkiHelsinkiFinland
| |
Collapse
|
42
|
Ferrer Obiol J, James HF, Chesser RT, Bretagnolle V, González-Solís J, Rozas J, Riutort M, Welch AJ. Integrating Sequence Capture and Restriction Site-Associated DNA Sequencing to Resolve Recent Radiations of Pelagic Seabirds. Syst Biol 2021; 70:976-996. [PMID: 33512506 PMCID: PMC8357341 DOI: 10.1093/sysbio/syaa101] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2020] [Revised: 11/13/2020] [Accepted: 12/15/2020] [Indexed: 01/01/2023] Open
Abstract
The diversification of modern birds has been shaped by a number of radiations. Rapid diversification events make reconstructing the evolutionary relationships among taxa challenging due to the convoluted effects of incomplete lineage sorting (ILS) and introgression. Phylogenomic data sets have the potential to detect patterns of phylogenetic incongruence, and to address their causes. However, the footprints of ILS and introgression on sequence data can vary between different phylogenomic markers at different phylogenetic scales depending on factors such as their evolutionary rates or their selection pressures. We show that combining phylogenomic markers that evolve at different rates, such as paired-end double-digest restriction site-associated DNA (PE-ddRAD) and ultraconserved elements (UCEs), allows a comprehensive exploration of the causes of phylogenetic discordance associated with short internodes at different timescales. We used thousands of UCE and PE-ddRAD markers to produce the first well-resolved phylogeny of shearwaters, a group of medium-sized pelagic seabirds that are among the most phylogenetically controversial and endangered bird groups. We found that phylogenomic conflict was mainly derived from high levels of ILS due to rapid speciation events. We also documented a case of introgression, despite the high philopatry of shearwaters to their breeding sites, which typically limits gene flow. We integrated state-of-the-art concatenated and coalescent-based approaches to expand on previous comparisons of UCE and RAD-Seq data sets for phylogenetics, divergence time estimation, and inference of introgression, and we propose a strategy to optimize RAD-Seq data for phylogenetic analyses. Our results highlight the usefulness of combining phylogenomic markers evolving at different rates to understand the causes of phylogenetic discordance at different timescales. [Aves; incomplete lineage sorting; introgression; PE-ddRAD-Seq; phylogenomics; radiations; shearwaters; UCEs.].
Collapse
Affiliation(s)
- Joan Ferrer Obiol
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, Universitat de Barcelona, Barcelona, Catalonia, Spain
- Institut de Recerca de la Biodiversitat (IRBio), Barcelona, Catalonia, Spain
| | - Helen F James
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, USA
| | - R Terry Chesser
- Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, USA
- U.S. Geological Survey, Patuxent Wildlife Research Center, Laurel, MD, USA
| | - Vincent Bretagnolle
- Centre d’Études Biologiques de Chizé, CNRS & La Rochelle Université, 79360, Villiers en Bois, France
| | - Jacob González-Solís
- Institut de Recerca de la Biodiversitat (IRBio), Barcelona, Catalonia, Spain
- Departament de Biologia Evolutiva, Ecologia i Ciències Ambientals, Facultat de Biologia, Universitat de Barcelona, Barcelona, Catalonia, Spain
| | - Julio Rozas
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, Universitat de Barcelona, Barcelona, Catalonia, Spain
- Institut de Recerca de la Biodiversitat (IRBio), Barcelona, Catalonia, Spain
| | - Marta Riutort
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, Universitat de Barcelona, Barcelona, Catalonia, Spain
- Institut de Recerca de la Biodiversitat (IRBio), Barcelona, Catalonia, Spain
| | | |
Collapse
|
43
|
Raposo do Amaral F, Thom G, Lima-Ribeiro MS, Alvarado-Serrano DF, Montesanti JAC, Pellegrino KCM, Miyaki CY, Hickerson MJ, Maldonado-Coelho M. Rugged relief and climate promote isolation and divergence between two neotropical cold-associated birds. Evolution 2021; 75:2371-2387. [PMID: 34375460 DOI: 10.1111/evo.14318] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Revised: 06/21/2021] [Accepted: 07/19/2021] [Indexed: 01/15/2023]
Abstract
The role of historical factors in establishing patterns of diversity in tropical mountains is of interest to understand the buildup of megadiverse biotas. In these regions, the historical processes of range fragmentation and contraction followed by dispersal are thought to be mediated by the interplay between rugged relief (complex topography) and climate fluctuations and likely explain most of the dynamics of diversification in plants and animals. Although empirical studies addressing the interaction between climate and topography have provided invaluable insights into population divergence and speciation patterns in tropical montane organisms, a more detailed and robust test of such processes in an explicit spatio-temporal framework is still lacking. Consequently, our ability to gain insights into historical range shifts over time and the genomic footprint left by them is limited. Here, we used niche modeling and subgenomic population-level datasets to explore the evolution of two species of warbling finches (genus Microspingus) disjunctly distributed across the Montane Atlantic Forest, a Neotropical region with complex geological and environmental histories. Population structure inferences suggest a scenario of three genetically differentiated populations, which are congruent with both geography and phenotypic variation. Demographic simulations support asynchronous isolation of these populations as recently as ∼40,000 years ago, relatively stable population sizes over recent time, and past gene flow subsequent to divergence. Throughout the last 800,000 years, niche models predicted extensive expansion into lowland areas with increasing overlap of species distributions during glacial periods, with prominent retractions and isolation into higher altitudes during interglacials, which are in line with signs of introgression of currently isolated populations. These results support a dual role of cyclical climatic changes: population divergence and persistence in mountain tops during warm periods followed by periods of expansion and admixture in lower elevations during cold periods. Our results underscore the role of the interplay between landscape and climate as an important mechanism in the evolution of the Neotropical montane biota.
Collapse
Affiliation(s)
- Fábio Raposo do Amaral
- Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Rua Professor Artur Riedel, 275, Diadema, SP, 09972-270, Brazil
| | - Gregory Thom
- Department of Ornithology, Division of Vertebrate Zoology, American Museum of Natural History, Central Park West at 79th Street, New York, NY, 10024-5192, USA
| | - Matheus S Lima-Ribeiro
- Departamento de Ciências Biológicas, Universidade Federal de Jataí, CP 03, Jataí, GO, 75804-020, Brazil
| | - Diego F Alvarado-Serrano
- Department of Biological Sciences, Ohio Center for Ecology and Evolutionary Studies, Ohio University, Life Sciences Building R219, Athens, OH, 45701, USA
| | - Julia A C Montesanti
- Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Rua Professor Artur Riedel, 275, Diadema, SP, 09972-270, Brazil
| | - Katia C M Pellegrino
- Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Rua Professor Artur Riedel, 275, Diadema, SP, 09972-270, Brazil
| | - Cristina Y Miyaki
- Departamento de Genética e Biologia Evolutiva, Universidade de São Paulo, Rua do Matão, 277, Cidade Universitária, São Paulo, SP, 05508-090, Brazil
| | - Michael J Hickerson
- Department of Biology, City College of New York, 160, Convent Avenue, New York, NY, 10031, USA
| | - Marcos Maldonado-Coelho
- Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Rua Professor Artur Riedel, 275, Diadema, SP, 09972-270, Brazil.,Department of Biology, Lund University, Lund, SE-223 62, Sweden
| |
Collapse
|
44
|
Esquerré D, Keogh JS, Demangel D, Morando M, Avila LJ, Sites JW, Ferri-Yáñez F, Leaché AD. Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards. Syst Biol 2021; 71:286-300. [PMID: 34259868 DOI: 10.1093/sysbio/syab058] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Revised: 06/25/2021] [Accepted: 06/30/2021] [Indexed: 01/09/2023] Open
Abstract
Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.
Collapse
Affiliation(s)
- Damien Esquerré
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, ACT, Australia
| | - J Scott Keogh
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, ACT, Australia
| | | | - Mariana Morando
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales (IPEEC- CONICET), Puerto Madryn, Chubut, Argentina
| | - Luciano J Avila
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales (IPEEC- CONICET), Puerto Madryn, Chubut, Argentina
| | - Jack W Sites
- Department of Biology and M.L. Bean Life Science Museum, Brigham Young University, Provo, Utah, USA
| | - Francisco Ferri-Yáñez
- Departamento de Biogeografía y Cambio Global, Museo Nacional de Ciencias Naturales, CSIC & Laboratorio Internacional en Cambio Global CSIC-PUC (LINCGlobal), Calle José Gutiérrez Abascal, 2, 28006, Madrid, Spain
| | - Adam D Leaché
- Department of Biology & Burke Museum of Natural History and Culture, University of Washington, Seattle, Washington, USA
| |
Collapse
|
45
|
Smith BT, Mauck WM, Benz BW, Andersen MJ. Uneven Missing Data Skew Phylogenomic Relationships within the Lories and Lorikeets. Genome Biol Evol 2021; 12:1131-1147. [PMID: 32470111 PMCID: PMC7486955 DOI: 10.1093/gbe/evaa113] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/26/2020] [Indexed: 01/21/2023] Open
Abstract
The resolution of the Tree of Life has accelerated with advances in DNA sequencing technology. To achieve dense taxon sampling, it is often necessary to obtain DNA from historical museum specimens to supplement modern genetic samples. However, DNA from historical material is generally degraded, which presents various challenges. In this study, we evaluated how the coverage at variant sites and missing data among historical and modern samples impacts phylogenomic inference. We explored these patterns in the brush-tongued parrots (lories and lorikeets) of Australasia by sampling ultraconserved elements in 105 taxa. Trees estimated with low coverage characters had several clades where relationships appeared to be influenced by whether the sample came from historical or modern specimens, which were not observed when more stringent filtering was applied. To assess if the topologies were affected by missing data, we performed an outlier analysis of sites and loci, and a data reduction approach where we excluded sites based on data completeness. Depending on the outlier test, 0.15% of total sites or 38% of loci were driving the topological differences among trees, and at these sites, historical samples had 10.9× more missing data than modern ones. In contrast, 70% data completeness was necessary to avoid spurious relationships. Predictive modeling found that outlier analysis scores were correlated with parsimony informative sites in the clades whose topologies changed the most by filtering. After accounting for biased loci and understanding the stability of relationships, we inferred a more robust phylogenetic hypothesis for lories and lorikeets.
Collapse
Affiliation(s)
- Brian Tilston Smith
- Department of Ornithology, American Museum of Natural History, New York, New York
| | - William M Mauck
- Department of Ornithology, American Museum of Natural History, New York, New York.,New York Genome Center, New York, New York
| | - Brett W Benz
- Museum of Zoology and Department of Ecology and Evolutionary Biology, University of Michigan
| | - Michael J Andersen
- Department of Biology and Museum of Southwestern Biology, University of New Mexico
| |
Collapse
|
46
|
Zhou W, Soghigian J, Xiang QYJ. A New Pipeline for Removing Paralogs in Target Enrichment Data. Syst Biol 2021; 71:410-425. [PMID: 34146111 PMCID: PMC8974407 DOI: 10.1093/sysbio/syab044] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 06/04/2021] [Accepted: 06/12/2021] [Indexed: 12/30/2022] Open
Abstract
Target enrichment (such as Hyb-Seq) is a well-established high throughput sequencing
method that has been increasingly used for phylogenomic studies. Unfortunately, current
widely used pipelines for analysis of target enrichment data do not have a vigorous
procedure to remove paralogs in target enrichment data. In this study, we develop a
pipeline we call Putative Paralogs Detection (PPD) to better address putative paralogs
from enrichment data. The new pipeline is an add-on to the existing HybPiper pipeline, and
the entire pipeline applies criteria in both sequence similarity and heterozygous sites at
each locus in the identification of paralogs. Users may adjust the thresholds of sequence
identity and heterozygous sites to identify and remove paralogs according to the level of
phylogenetic divergence of their group of interest. The new pipeline also removes highly
polymorphic sites attributed to errors in sequence assembly and gappy regions in the
alignment. We demonstrated the value of the new pipeline using empirical data generated
from Hyb-Seq and the Angiosperms353 kit for two woody genera Castanea
(Fagaceae, Fagales) and Hamamelis (Hamamelidaceae, Saxifragales).
Comparisons of data sets showed that the PPD identified many more putative paralogs than
the popular method HybPiper. Comparisons of tree topologies and divergence times showed
evident differences between data from HybPiper and data from our new PPD pipeline. We
further evaluated the accuracy and error rates of PPD by BLAST mapping of putative
paralogous and orthologous sequences to a reference genome sequence of Castanea
mollissima. Compared to HybPiper alone, PPD identified substantially more
paralogous gene sequences that mapped to multiple regions of the reference genome (31
genes for PPD compared with 4 genes for HybPiper alone). In conjunction with HybPiper,
paralogous genes identified by both pipelines can be removed resulting in the construction
of more robust orthologous gene data sets for phylogenomic and divergence time analyses.
Our study demonstrates the value of Hyb-Seq with data derived from the Angiosperms353
probe set for elucidating species relationships within a genus, and argues for the
importance of additional steps to filter paralogous genes and poorly aligned regions
(e.g., as occur through assembly errors), such as our new PPD pipeline described in this
study. [Angiosperms353; Castanea; divergence time;
Hamamelis; Hyb-Seq, paralogs, phylogenomics.]
Collapse
Affiliation(s)
- Wenbin Zhou
- Department of Plant and Microbial Biology, North Carolina State University, Raleigh, NC 27965, USA
| | - John Soghigian
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC 27965, USA
| | - Qiu-Yun Jenny Xiang
- Department of Plant and Microbial Biology, North Carolina State University, Raleigh, NC 27965, USA
| |
Collapse
|
47
|
Comparison of sequence-capture and ddRAD approaches in resolving species and populations in hexacorallian anthozoans. Mol Phylogenet Evol 2021; 163:107233. [PMID: 34139346 DOI: 10.1016/j.ympev.2021.107233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2020] [Revised: 05/25/2021] [Accepted: 06/11/2021] [Indexed: 10/21/2022]
Abstract
Genome-level sequencing is the next step in understanding species-level relationships within Anthozoa (soft corals, anemones, stony corals, and their kin) as morphological and PCR-directed (single-locus) sequencing methods often fall short of differentiating species. The sea anemone genus Metridium is a common northern temperate sea anemone whose species are difficult to differentiate using morphology alone. Here we use Metridium as a case study to confirm the low level of information available in six loci for species differentiation commonly sequenced for Actiniaria and explore and compare the efficacy of ddRAD and sequence-capture methods in species-level systematics and biogeographic studies. We produce phylogenetic trees from concatenated datasets and perform DAPC and STRUCTURE analyses using SNP data. The six conventional loci are not able to consistently differentiate species within Metridium. The sequence-capture dataset resulted in high support and resolution for both current species and relationships between geographic areas. The ddRAD datasets displayed ambiguity among species, and support between major geographic groupings was not as high as the sequence-capture datasets. The level of resolution and support resulting from the sequence-capture data, combined with the ability to add additional individuals and expand beyond the genus Metridium over time, emphasizes the utility of sequence-capture methods for both systematics and future biogeographic studies within anthozoans. We discuss the strengths and weaknesses of the genomic approaches in light of our findings and suggest potential implications for the biogeography of Metridium based on our sampling.
Collapse
|
48
|
Esquerré D, Donnellan SC, Pavón-Vázquez CJ, Fenker J, Keogh JS. Phylogeography, historical demography and systematics of the world's smallest pythons (Pythonidae, Antaresia). Mol Phylogenet Evol 2021; 161:107181. [PMID: 33892100 DOI: 10.1016/j.ympev.2021.107181] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Revised: 04/06/2021] [Accepted: 04/15/2021] [Indexed: 11/18/2022]
Abstract
Advances from empirical studies in phylogeography, systematics and species delimitation highlight the importance of integrative approaches for quantifying taxonomic diversity. Genomic data have greatly improved our ability to discern both systematic diversity and evolutionary history. Here we combine analyses of mitochondrial DNA sequences, thousands of genome-wide SNPs and linear and geometric morphometrics on Antaresia, a clade of four currently recognised dwarf pythons from Australia and New Guinea (Antaresia childreni, A. stimsoni, A. maculosa and A. perthensis). Our integrative analyses of phylogenetics, population structure, species delimitation, historical demography and morphometrics revealed that the true evolutionary diversity is not well reflected in the current appraisal of the diversity of the group. We find that Antaresia childreni and A. stimsoni comprise a widespread network of populations connected by gene flow and without evidence of species-level divergence among them. However, A. maculosa shows considerable genetic structuring which leads us to recognise two subspecies in northeastern Australia and a new species in Torres Strait and New Guinea. These two contrasting cases of over and under estimation of diversity, respectively, illustrate the power of thorough integrative approaches into understanding evolution of biodiversity. Furthermore, our analyses of historical demographic patterns highlight the importance of the Kimberley, Pilbara and Cape York as origins of biodiversity in Australia.
Collapse
Affiliation(s)
- Damien Esquerré
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, ACT 0200, Australia.
| | | | - Carlos J Pavón-Vázquez
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, ACT 0200, Australia
| | - Jéssica Fenker
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, ACT 0200, Australia
| | - J Scott Keogh
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Canberra, ACT 0200, Australia
| |
Collapse
|
49
|
Terraneo TI, Benzoni F, Arrigoni R, Baird AH, Mariappan KG, Forsman ZH, Wooster MK, Bouwmeester J, Marshell A, Berumen ML. Phylogenomics of Porites from the Arabian Peninsula. Mol Phylogenet Evol 2021; 161:107173. [PMID: 33813021 DOI: 10.1016/j.ympev.2021.107173] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 03/25/2021] [Accepted: 03/29/2021] [Indexed: 11/16/2022]
Abstract
The advent of high throughput sequencing technologies provides an opportunity to resolve phylogenetic relationships among closely related species. By incorporating hundreds to thousands of unlinked loci and single nucleotide polymorphisms (SNPs), phylogenomic analyses have a far greater potential to resolve species boundaries than approaches that rely on only a few markers. Scleractinian taxa have proved challenging to identify using traditional morphological approaches and many groups lack an adequate set of molecular markers to investigate their phylogenies. Here, we examine the potential of Restriction-site Associated DNA sequencing (RADseq) to investigate phylogenetic relationships and species limits within the scleractinian coral genus Porites. A total of 126 colonies were collected from 16 localities in the seas surrounding the Arabian Peninsula and ascribed to 12 nominal and two unknown species based on their morphology. Reference mapping was used to retrieve and compare nearly complete mitochondrial genomes, ribosomal DNA, and histone loci. De novo assembly and reference mapping to the P. lobata coral transcriptome were compared and used to obtain thousands of genome-wide loci and SNPs. A suite of species discovery methods (phylogenetic, ordination, and clustering analyses) and species delimitation approaches (coalescent-based, species tree, and Bayesian Factor delimitation) suggested the presence of eight molecular lineages, one of which included six morphospecies. Our phylogenomic approach provided a fully supported phylogeny of Porites from the Arabian Peninsula, suggesting the power of RADseq data to solve the species delineation problem in this speciose coral genus.
Collapse
Affiliation(s)
- Tullia I Terraneo
- Red Sea Research Centre, Division of Biological and Environmental Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia; ARC Centre of Excellence for Coral Reef Studies, James Cook University, Townsville 4811, QLD, Australia.
| | - Francesca Benzoni
- Red Sea Research Centre, Division of Biological and Environmental Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Roberto Arrigoni
- Red Sea Research Centre, Division of Biological and Environmental Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia; European Commission, Joint Research Centre (JRC), Ispra, Italy; Department of Biology and Evolution of Marine Organisms (BEOM), Stazione Zoologica Anton Dohrn Napoli, Villa Comunale, 80121 Napoli, Italy
| | - Andrew H Baird
- ARC Centre of Excellence for Coral Reef Studies, James Cook University, Townsville 4811, QLD, Australia
| | - Kiruthiga G Mariappan
- Red Sea Research Centre, Division of Biological and Environmental Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Zac H Forsman
- Hawaii Institute of Marine Biology, Kaneohe 96744, HI, USA
| | - Michael K Wooster
- Red Sea Research Centre, Division of Biological and Environmental Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | | | - Alyssa Marshell
- Department of Marine Science and Fisheries, College of Agricultural and Marine Sciences, Sultan Qaboos University, Muscat, Oman
| | - Michael L Berumen
- Red Sea Research Centre, Division of Biological and Environmental Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
50
|
Andersen MJ, McCullough JM, Gyllenhaal EF, Mapel XM, Haryoko T, Jønsson KA, Joseph L. Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds. Mol Ecol 2021; 30:2087-2103. [PMID: 33615597 PMCID: PMC8252742 DOI: 10.1111/mec.15856] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Revised: 02/02/2021] [Accepted: 02/12/2021] [Indexed: 01/03/2023]
Abstract
Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.
Collapse
Affiliation(s)
- Michael J Andersen
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Jenna M McCullough
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Ethan F Gyllenhaal
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Xena M Mapel
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA.,Animal Genomics, ETH Zürich, Lindau, Switzerland
| | - Tri Haryoko
- Museum Zoologicum Bogoriense, Research Centre for Biology, Indonesian Institute of Sciences (LIPI), Cibinong, Indonesia
| | - Knud A Jønsson
- Natural History Museum of Denmark, University of Copenhagen, Copenhagen Ø, Denmark
| | - Leo Joseph
- Australian National Wildlife Collection, CSIRO National Research Collections, Canberra, Australian Capital Territory, Australia
| |
Collapse
|