1
|
Burbrink FT, Myers EA, Pyron RA. Understanding species limits through the formation of phylogeographic lineages. Ecol Evol 2024; 14:e70263. [PMID: 39364037 PMCID: PMC11446989 DOI: 10.1002/ece3.70263] [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: 06/03/2024] [Revised: 08/12/2024] [Accepted: 08/16/2024] [Indexed: 10/05/2024] Open
Abstract
The outcomes of speciation across organismal dimensions (e.g., ecological, genetic, phenotypic) are often assessed using phylogeographic methods. At one extreme, reproductively isolated lineages represent easily delimitable species differing in many or all dimensions, and at the other, geographically distinct genetic segments introgress across broad environmental gradients with limited phenotypic disparity. In the ambiguous gray zone of speciation, where lineages are genetically delimitable but still interacting ecologically, it is expected that these lineages represent species in the context of ontology and the evolutionary species concept when they are maintained over time with geographically well-defined hybrid zones, particularly at the intersection of distinct environments. As a result, genetic structure is correlated with environmental differences and not space alone, and a subset of genes fail to introgress across these zones as underlying genomic differences accumulate. We present a set of tests that synthesize species delimitation with the speciation process. We can thereby assess historical demographics and diversification processes while understanding how lineages are maintained through space and time by exploring spatial and genome clines, genotype-environment interactions, and genome scans for selected loci. Employing these tests in eight lineage-pairs of snakes in North America, we show that six pairs represent 12 "good" species and that two pairs represent local adaptation and regional population structure. The distinct species pairs all have the signature of divergence before or near the mid-Pleistocene, often with low migration, stable hybrid zones of varying size, and a subset of loci showing selection on alleles at the hybrid zone corresponding to transitions between distinct ecoregions. Locally adapted populations are younger, exhibit higher migration, and less ecological differentiation. Our results demonstrate that interacting lineages can be delimited using phylogeographic and population genetic methods that properly integrate spatial, temporal, and environmental data.
Collapse
Affiliation(s)
- Frank T Burbrink
- Department of Herpetology American Museum of Natural History New York New York USA
| | - Edward A Myers
- Department of Herpetology California Academy of Sciences San Francisco California USA
| | - R Alexander Pyron
- Department of Biological Sciences The George Washington University Washington DC USA
| |
Collapse
|
2
|
Liu A, Geraldes A, Taylor EB. Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish. Mol Ecol 2024; 33:e17459. [PMID: 38994921 DOI: 10.1111/mec.17459] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 04/23/2024] [Accepted: 04/30/2024] [Indexed: 07/13/2024]
Abstract
Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.
Collapse
Affiliation(s)
- Amy Liu
- Department of Zoology, University of British Columbia, Vancouver, British Columbia, Canada
- Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Armando Geraldes
- Department of Zoology, University of British Columbia, Vancouver, British Columbia, Canada
- Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Eric B Taylor
- Department of Zoology, University of British Columbia, Vancouver, British Columbia, Canada
- Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
- Beaty Biodiversity Museum, University of British Columbia, Vancouver, British Columbia, Canada
| |
Collapse
|
3
|
Myers EA, Rautsaw RM, Borja M, Jones J, Grünwald CI, Holding ML, Grazziotin F, Parkinson CL. Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus). Syst Biol 2024:syae018. [PMID: 38695290 DOI: 10.1093/sysbio/syae018] [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: 05/10/2023] [Indexed: 08/11/2024] Open
Abstract
Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.
Collapse
Affiliation(s)
- Edward A Myers
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
- Department of Herpetology, California Academy of Sciences, San Francisco, CA 94118, USA
| | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
| | - Miguel Borja
- Facultad de Ciencias Biológicas, Universdad Juárez del Estado de Durango, Av. Universidad s/n. Fracc. Filadelfia, Gómez Palacio, Durango., Mex
| | - Jason Jones
- Herp.mx A.C. C.P. 28989, Villa de Álvarez, Colima, Mexico
| | - Christoph I Grünwald
- Herp.mx A.C. C.P. 28989, Villa de Álvarez, Colima, Mexico
- Biodiversa A.C., Avenida de la Ribera #203, C.P. 45900, Chapala, Jalisco, Mexico
| | - Matthew L Holding
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
- Life Sciences Institute, University of Michigan, Ann Arbor, MI 48109, USA
| | - Felipe Grazziotin
- Laboratório Especial de Coleções Zoológicas, Instituto Butantan, Avenida Vital Brasil, São Paulo, 05503-900, Brazil
| | | |
Collapse
|
4
|
Lyu B, Liu Q, Wu Y, Nguyen TQ, Che J, Nguyen SN, Myers EA, Burbrink FT, Guo P, Wang J. Genomic analysis reveals deep population divergence in the water snake Trimerodytes percarinatus (Serpentes, Natricidae). Ecol Evol 2024; 14:e11278. [PMID: 38628918 PMCID: PMC11019134 DOI: 10.1002/ece3.11278] [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: 09/30/2023] [Revised: 02/26/2024] [Accepted: 04/02/2024] [Indexed: 04/19/2024] Open
Abstract
Although several phylogeographic studies of Asian snakes have been conducted, most have focused on pitvipers, with non-venomous snakes, such as colubrids or natricids, remaining poorly studied. The Chinese keelback water snake (Trimerodytes percarinatus Boulenger) is a widespread, semiaquatic, non-venomous species occurring in China and southeastern Asia. Based on mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, we explored the population genetic structure, genetic diversity, and evolutionary history of this species. MtDNA-based phylogenetic analysis showed that T. percarinatus was composed of five highly supported and geographically structured lineages. SNP-based phylogenetic analysis, principal component analysis, and population structure analysis consistently revealed four distinct, geographically non-overlapping lineages, which was different from the mtDNA-based analysis in topology. Estimation of divergence dates and ancestral area of origin suggest that T. percarinatus originated ~12.68 million years ago (95% highest posterior density: 10.36-15.96 Mya) in a region covering southwestern China and Vietnam. Intraspecific divergence may have been triggered by the Qinghai-Xizang Plateau uplift. Population demographics and ecological niche modeling indicated that the effective population size fluctuated during 0.5 Mya and 0.002 Mya. Based on the data collected here, we also comment on the intraspecific taxonomy of T. percarinatus and question the validity of the subspecies T. p. suriki.
Collapse
Affiliation(s)
- Bing Lyu
- Ministry of Education Key Laboratory for Ecology of Tropical Islands, Key Laboratory of Tropical Animal and Plant Ecology of Hainan Province, College of Life SciencesHainan Normal UniversityHaikouChina
- Faculty of Agriculture, Forestry and Food EngineeringYibin UniversityYibinChina
| | - Qin Liu
- Faculty of Agriculture, Forestry and Food EngineeringYibin UniversityYibinChina
| | - Yayong Wu
- Faculty of Agriculture, Forestry and Food EngineeringYibin UniversityYibinChina
| | - Truong Q. Nguyen
- Institute of Ecology and Biological ResourcesVietnam Academy of Science and TechnologyHanoiVietnam
- Vietnam Academy of Science and TechnologyGraduate University of Science and TechnologyHanoiVietnam
| | - Jing Che
- State Key Laboratory of Genetic Resources and Evolution & Yunnan Key Laboratory of Biodiversity and Ecological Conservation of Gaoligong Mountain, Kunming Institute of ZoologyChinese Academy of SciencesKunmingChina
| | - Sang N. Nguyen
- Institute of Tropical BiologyVietnam Academy of Science and TechnologyHo Chi Minh CityVietnam
| | - Edward A. Myers
- Department of HerpetologyCalifornia Academy of SciencesSan FranciscoCaliforniaUSA
| | - Frank T. Burbrink
- Department of HerpetologyAmerican Museum of Natural HistoryNew YorkNew YorkUSA
| | - Peng Guo
- Faculty of Agriculture, Forestry and Food EngineeringYibin UniversityYibinChina
| | - Jichao Wang
- Ministry of Education Key Laboratory for Ecology of Tropical Islands, Key Laboratory of Tropical Animal and Plant Ecology of Hainan Province, College of Life SciencesHainan Normal UniversityHaikouChina
| |
Collapse
|
5
|
Smith CF, Nikolakis ZL, Perry BW, Schield DR, Meik JM, Saviola AJ, Castoe TA, Parker J, Mackessy SP. The best of both worlds? Rattlesnake hybrid zones generate complex combinations of divergent venom phenotypes that retain high toxicity. Biochimie 2023; 213:176-189. [PMID: 37451532 DOI: 10.1016/j.biochi.2023.07.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Revised: 06/27/2023] [Accepted: 07/10/2023] [Indexed: 07/18/2023]
Abstract
Studying the consequences of hybridization between closely related species with divergent traits can reveal patterns of evolution that shape and maintain extreme trophic adaptations. Snake venoms are an excellent model system for examining the evolutionary and ecological patterns that underlie highly selected polymorphic traits. Here we investigate hybrid venom phenotypes that result from natural introgression between two rattlesnake species that express highly divergent venom phenotypes: Crotalus o. concolor and C. v. viridis. Though not yet documented, interbreeding between these species may lead to novel venom phenotypes with unique activities that break the typical trends of venom composition in rattlesnakes. The characteristics of these unusual phenotypes could unveil the roles of introgression in maintaining patterns of venom composition and variation, including the near ubiquitous dichotomy between neurotoxic or degradative venoms observed across rattlesnakes. We use RADseq data to infer patterns of gene flow and hybrid ancestry between these diverged lineages and link these genetic data with analyses of venom composition, biological activity, and whole animal model toxicity tests to understand the impacts of introgression on venom composition. We find that introgressed populations express admixed venom phenotypes that do not sacrifice biological activity (lethal toxicity) or overall abundance of dominant toxins compared to parental venoms. These hybridized venoms therefore do not represent a trade-off in functionality between the typical phenotypic extremes but instead represent a unique combination of characters whose expression appears limited to the hybrid zone.
Collapse
Affiliation(s)
- Cara F Smith
- School of Biological Sciences, 501 20th Street, University of Northern Colorado, Greeley, CO, 80639, USA
| | - Zachary L Nikolakis
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX, 76019, USA
| | - Blair W Perry
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX, 76019, USA
| | - Drew R Schield
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX, 76019, USA
| | - Jesse M Meik
- Department of Biological Sciences, Tarleton State University, 1333 W. Washington Street, Stephenville, TX, 76402, USA
| | - Anthony J Saviola
- Department of Biochemistry and Molecular Genetics, 12801 East 17th Avenue, University of Colorado Denver, Aurora, CO, 80045, USA
| | - Todd A Castoe
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX, 76019, USA
| | - Joshua Parker
- Fresno City College, 1101 E. University Avenue, Fresno, CA, 93741, USA
| | - Stephen P Mackessy
- School of Biological Sciences, 501 20th Street, University of Northern Colorado, Greeley, CO, 80639, USA.
| |
Collapse
|
6
|
Finger N, Farleigh K, Bracken JT, Leaché AD, François O, Yang Z, Flouri T, Charran T, Jezkova T, Williams DA, Blair C. Genome-scale data reveal deep lineage divergence and a complex demographic history in the Texas horned lizard (Phrynosoma cornutum) throughout the southwestern and central US. Genome Biol Evol 2021; 14:6443127. [PMID: 34849831 PMCID: PMC8735750 DOI: 10.1093/gbe/evab260] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/12/2021] [Indexed: 12/03/2022] Open
Abstract
The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find that divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.
Collapse
Affiliation(s)
- Nicholas Finger
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY, 11201, USA
| | - Keaka Farleigh
- Department of Biology, Miami University, 501 E High St, Oxford, OH, 45056, USA
| | - Jason T Bracken
- Department of Biology, Miami University, 501 E High St, Oxford, OH, 45056, USA
| | - Adam D Leaché
- Department of Biology & Burke Museum of Natural History and Culture, University of Washington, Seattle, WA, 98195, USA
| | - Olivier François
- Faculty of Medicine, University Grenoble-Alpes, TIMC-IMAG UMR 5525, Grenoble, La Tronche, F38706, France 38000
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, Darwin Building, Gower Street, London, WC1E 6BT, UK
| | - Tomas Flouri
- Department of Genetics, Evolution and Environment, University College London, Darwin Building, Gower Street, London, WC1E 6BT, UK
| | - Tristan Charran
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY, 11201, USA
| | - Tereza Jezkova
- Department of Biology, Miami University, 501 E High St, Oxford, OH, 45056, USA
| | - Dean A Williams
- Department of Biology, Texas Christian University, 2800 S University Dr, Fort Worth, TX, 76129, USA
| | - Christopher Blair
- Department of Biological Sciences, New York City College of Technology, The City University of New York, 285 Jay Street, Brooklyn, NY, 11201, USA.,Biology PhD Program, CUNY Graduate Center, 365 5th Ave, New York, NY, 10016, USA
| |
Collapse
|
7
|
Blais BR, Smith BE, Placyk JS, Casper GS, Spellman GM. Phylogeography of the smooth greensnake, Opheodrys vernalis (Squamata: Colubridae): divergent lineages and variable demographics in a widely distributed yet enigmatic species. Biol J Linn Soc Lond 2021. [DOI: 10.1093/biolinnean/blab124] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Abstract
Phylogeographic studies can uncover robust details about the population structure, demographics, and diversity of species. The smooth greensnake, Opheodrys vernalis, is a small, cryptic snake occupying mesic grassland and sparsely wooded habitats. Although O. vernalis has a wide geographical range, many metapopulations are patchy and some are declining. We used mitochondrial DNA and double digest restriction-site associated DNA sequencing to construct the first phylogeographic assessment of O. vernalis. Genomic analysis of 119 individuals (mitochondrial DNA) and a subset of another 45 smooth greensnakes (nuclear DNA; N = 3031 single nucleotide polymorphisms) strongly supports two longitudinally separated lineages, with admixture in the Great Lakes region. Post-Pleistocene secondary contact best explains admixture from populations advancing northwards. Overall, populations expressed low heterozygosity, variable inbreeding rates, and moderate to high differentiation. Disjunct populations in the Rocky Mountains and central Great Plains regions might be contracting relicts, whereas northerly populations in more continuous mesic habitats (e.g., Prairie Pothole region, southern Canada) had signals of population expansion. Broadly, conservation management efforts should be focused on local populations, because habitat connectivity may facilitate gene flow and genetic diversity.
Collapse
Affiliation(s)
- Brian R Blais
- School of Natural Sciences, Black Hills State University, Spearfish, SD, USA
| | - Brian E Smith
- School of Natural Sciences, Black Hills State University, Spearfish, SD, USA
| | - John S Placyk
- Department of Biology, University of Texas at Tyler, 3900 University Boulevard, Tyler, TX, USA
| | - Gary S Casper
- University of Wisconsin-Milwaukee Field Station, Saukville, WI, USA
| | - Garth M Spellman
- Department of Zoology, Denver Museum of Nature & Science, 2001 Colorado Boulevard, Denver, CO, USA
| |
Collapse
|
8
|
Seneci L, Zdenek CN, Bourke LA, Cochran C, Sánchez EE, Neri-Castro E, Bénard-Valle M, Alagón A, Frank N, Fry BG. A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms. Comp Biochem Physiol C Toxicol Pharmacol 2021; 245:109034. [PMID: 33766656 PMCID: PMC8162888 DOI: 10.1016/j.cbpc.2021.109034] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/16/2021] [Revised: 03/02/2021] [Accepted: 03/15/2021] [Indexed: 12/15/2022]
Abstract
What factors influence the evolution of a heavily selected functional trait in a diverse clade? This study adopts rattlesnakes as a model group to investigate the evolutionary history of venom coagulotoxicity in the wider context of phylogenetics, natural history, and biology. Venom-induced clotting of human plasma and fibrinogen was determined and mapped onto the rattlesnake phylogenetic tree to reconstruct the evolution of coagulotoxicity across the group. Our results indicate that venom phenotype is often independent of phylogenetic relationships in rattlesnakes, suggesting the importance of diet and/or other environmental variables in driving venom evolution. Moreover, the striking inter- and intraspecific variability in venom activity on human blood highlights the considerable variability faced by physicians treating envenomation. This study is the most comprehensive effort to date to describe and characterize the evolutionary and biological aspects of coagulotoxins in rattlesnake venom. Further research at finer taxonomic levels is recommended to elucidate patterns of variation within species and lineages.
Collapse
Affiliation(s)
- Lorenzo Seneci
- Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072, Australia; Institute of Biology Leiden (IBL), Leiden University, 2333 BE Leiden, the Netherlands
| | - Christina N Zdenek
- Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072, Australia
| | - Lachlan A Bourke
- Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072, Australia
| | - Chip Cochran
- Department of Earth and Biological Sciences, Loma Linda University, Loma Linda, CA 92350, USA
| | - Elda E Sánchez
- National Natural Toxins Research Center (NNTRC), Department of Chemistry, Texas A&M University-Kingsville, MSC 224, 975 West Avenue B, Kingsville, TX 78363, USA
| | - Edgar Neri-Castro
- Departamento de Medicina Molecular y Bioprocesos, Instituto de Biotecnologa, Universidad Nacional Autónoma de México, Av. Universidad 2001, Cuernavaca, Morelos 62210, Mexico
| | - Melisa Bénard-Valle
- Departamento de Medicina Molecular y Bioprocesos, Instituto de Biotecnologa, Universidad Nacional Autónoma de México, Av. Universidad 2001, Cuernavaca, Morelos 62210, Mexico
| | - Alejandro Alagón
- Departamento de Medicina Molecular y Bioprocesos, Instituto de Biotecnologa, Universidad Nacional Autónoma de México, Av. Universidad 2001, Cuernavaca, Morelos 62210, Mexico
| | | | - Bryan G Fry
- Venom Evolution Lab, School of Biological Sciences, University of Queensland, St Lucia, QLD 4072, Australia.
| |
Collapse
|
9
|
Rautsaw RM, Schramer TD, Acuña R, Arick LN, DiMeo M, Mercier KP, Schrum M, Mason AJ, Margres MJ, Strickland JL, Parkinson CL. Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes. Mol Biol Evol 2021; 38:745-760. [PMID: 33035326 PMCID: PMC7947766 DOI: 10.1093/molbev/msaa266] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.
Collapse
Affiliation(s)
- Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Biology, University of Central Florida, Orlando, FL
| | | | - Rachel Acuña
- Department of Biology, University of Central Florida, Orlando, FL
| | - Lindsay N Arick
- Department of Biology, University of Central Florida, Orlando, FL
| | - Mark DiMeo
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Biology, University of Central Florida, Orlando, FL
| | - Kathryn P Mercier
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Biology, University of Central Florida, Orlando, FL.,Department of Biology, City College of New York, New York, NY.,PhD Program in Biology, The Graduate Center of the City University of New York, New York, NY
| | - Michael Schrum
- Department of Biology, University of Central Florida, Orlando, FL
| | - Andrew J Mason
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Biology, University of Central Florida, Orlando, FL
| | - Mark J Margres
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA.,Department of Integrative Biology, University of South Florida, Tampa, FL
| | - Jason L Strickland
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Biology, University of Central Florida, Orlando, FL.,Department of Biology, University of South Alabama, Mobile, AL
| | - Christopher L Parkinson
- Department of Biological Sciences, Clemson University, Clemson, SC.,Department of Biology, University of Central Florida, Orlando, FL.,Department of Forestry and Environmental Conservation, Clemson University, Clemson, SC
| |
Collapse
|
10
|
Ramírez‐Reyes T, Flores‐Villela O, Piñero D, Lathrop A, Murphy RW. Genomic assessment of the
Phyllodactylus tuberculosus
complex (Reptilia: Phyllodactylidae) in America. ZOOL SCR 2021. [DOI: 10.1111/zsc.12492] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tonatiuh Ramírez‐Reyes
- Posgrado en Ciencias Biológicas Universidad Nacional Autónoma de México México City México
- Museo de Zoología Facultad de Ciencias Universidad Nacional Autónoma de Méxic México City México
- Departamento de Ecología Evolutiva Instituto de EcologíaUniversidad Nacional Autónoma de México México City México
| | - Oscar Flores‐Villela
- Museo de Zoología Facultad de Ciencias Universidad Nacional Autónoma de Méxic México City México
| | - Daniel Piñero
- Departamento de Ecología Evolutiva Instituto de EcologíaUniversidad Nacional Autónoma de México México City México
| | - Amy Lathrop
- Centre for Biodiversity and Conservation Biology Royal Ontario Museum Toronto ON Canada
| | - Robert W. Murphy
- Centre for Biodiversity and Conservation Biology Royal Ontario Museum Toronto ON Canada
| |
Collapse
|
11
|
Humboldt-Paputsachis C, Fernandez Gil P. Análisis morfológico y morfométrico de Tityus (Tityus) sorataensis Kraepelin 1911 (Escorpionida: Buthidae) de dos valles mesotérmicos andinos, Quime y Cheje, La Paz-Bolivia. JOURNAL OF THE SELVA ANDINA RESEARCH SOCIETY 2021. [DOI: 10.36610/j.jsars.2021.120100003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
12
|
Chan KO, Hutter CR, Wood PL, Grismer LL, Das I, Brown RM. Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex. Mol Ecol 2020; 29:3970-3987. [PMID: 32808335 DOI: 10.1111/mec.15603] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2019] [Revised: 07/29/2020] [Accepted: 08/13/2020] [Indexed: 02/06/2023]
Abstract
Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.
Collapse
Affiliation(s)
- Kin O Chan
- Lee Kong Chian National History Museum, Faculty of Science, National University of Singapore, Singapore
| | - Carl R Hutter
- Biodiversity Institute and Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, USA.,Museum of Natural Sciences and Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA
| | - Perry L Wood
- Biodiversity Institute and Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, USA.,Department of Biological Sciences & Museum of Natural History, Auburn University, Auburn, AL, USA
| | - L L Grismer
- Herpetology Laboratory, Department of Biology, La Sierra University, Riverside, CA, USA
| | - Indraneil Das
- Institute of Biodiversity and Environmental Conservation, Universiti Malaysia Sarawak, Kota Samarahan, Sarawak, Malaysia
| | - Rafe M Brown
- Biodiversity Institute and Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, USA
| |
Collapse
|
13
|
Loureiro LO, Engstrom MD, Lim BK. Comparative phylogeography of mainland and insular species of Neotropical molossid bats ( Molossus). Ecol Evol 2020; 10:389-409. [PMID: 31993120 PMCID: PMC6972955 DOI: 10.1002/ece3.5903] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Revised: 11/13/2019] [Accepted: 11/15/2019] [Indexed: 01/01/2023] Open
Abstract
Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype-by-sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.
Collapse
Affiliation(s)
- Livia O. Loureiro
- Department of Natural HistoryRoyal Ontario MuseumTorontoONCanada
- Department of Ecology and Evolutionary BiologyUniversity of TorontoTorontoONCanada
| | - Mark D. Engstrom
- Department of Natural HistoryRoyal Ontario MuseumTorontoONCanada
- Department of Ecology and Evolutionary BiologyUniversity of TorontoTorontoONCanada
| | - Burton K. Lim
- Department of Natural HistoryRoyal Ontario MuseumTorontoONCanada
| |
Collapse
|
14
|
Campbell P, Arévalo L, Martin H, Chen C, Sun S, Rowe AH, Webster MS, Searle JB, Pasch B. Vocal divergence is concordant with genomic evidence for strong reproductive isolation in grasshopper mice ( Onychomys). Ecol Evol 2019; 9:12886-12896. [PMID: 31788222 PMCID: PMC6875671 DOI: 10.1002/ece3.5770] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2019] [Revised: 09/19/2019] [Accepted: 09/24/2019] [Indexed: 12/13/2022] Open
Abstract
Behavioral barriers to gene flow often evolve faster than intrinsic incompatibilities and can eliminate the opportunity for hybridization between interfertile species. While acoustic signal divergence is a common driver of premating isolation in birds and insects, its contribution to speciation in mammals is less studied. Here we characterize the incidence of, and potential barriers to, hybridization among three closely related species of grasshopper mice (genus Onychomys). All three species use long-distance acoustic signals to attract and localize mates; Onychomys arenicola and Onychomys torridus are acoustically similar and morphologically cryptic whereas Onychomys leucogaster is larger and acoustically distinct. We used genotyping-by-sequencing (GBS) to test for evidence of introgression in 227 mice from allopatric and sympatric localities in the western United States and northern Mexico. We conducted laboratory mating trials for all species pairs to assess reproductive compatibility, and recorded vocalizations from O. arenicola and O. torridus in sympatry and allopatry to test for evidence of acoustic character displacement. Hybridization was rare in nature and, contrary to prior evidence for O. torridus/O. arenicola hybrids, only involved O. leucogaster and O. arenicola. In contrast, laboratory crosses between O. torridus and O. arenicola produced litters whereas O. leucogaster and O. arenicola crosses did not. Call fundamental frequency in O. torridus and O. arenicola was indistinguishable in allopatry but significantly differentiated in sympatry, a pattern consistent with reproductive character displacement. These results suggest that assortative mating based on a long-distance signal is an important isolating mechanism between O. torridus and O. arenicola and highlight the importance of behavioral barriers in determining the permeability of species boundaries.
Collapse
Affiliation(s)
- Polly Campbell
- Department of Integrative BiologyOklahoma State UniversityStillwaterOKUSA
- Department of Evolution, Ecology, and Organismal BiologyUniversity of California, RiversideRiversideCAUSA
| | - Lena Arévalo
- Department of Integrative BiologyOklahoma State UniversityStillwaterOKUSA
- Department of Developmental PathologyUniversity of BonnBonnGermany
| | - Heather Martin
- Department of Integrative BiologyOklahoma State UniversityStillwaterOKUSA
| | - Charles Chen
- Department of Biochemistry and Molecular BiologyOklahoma State UniversityStillwaterOKUSA
| | - Shuzhen Sun
- Department of Biochemistry and Molecular BiologyOklahoma State UniversityStillwaterOKUSA
- Department of Forest and Conservation SciencesForest Science CentreThe University of British ColumbiaVancouverBCCanada
| | - Ashlee H. Rowe
- Department of BiologyThe University of OklahomaNormanOKUSA
| | - Michael S. Webster
- Macaulay LibraryCornell Lab of OrnithologyCornell UniversityIthacaNYUSA
- Department of Neurobiology and BehaviorCornell UniversityIthacaNYUSA
| | - Jeremy B. Searle
- Department of Ecology and Evolutionary BiologyCornell UniversityIthacaNYUSA
| | - Bret Pasch
- Department of Biological SciencesNorthern Arizona UniversityFlagstaffAZUSA
| |
Collapse
|
15
|
Methylation content sensitive enzyme ddRAD (MCSeEd): a reference-free, whole genome profiling system to address cytosine/adenine methylation changes. Sci Rep 2019; 9:14864. [PMID: 31619715 PMCID: PMC6795852 DOI: 10.1038/s41598-019-51423-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2019] [Accepted: 09/26/2019] [Indexed: 12/18/2022] Open
Abstract
Methods for investigating DNA methylation nowadays either require a reference genome and high coverage, or investigate only CG methylation. Moreover, no large-scale analysis can be performed for N6-methyladenosine (6 mA) at an affordable price. Here we describe the methylation content sensitive enzyme double-digest restriction-site-associated DNA (ddRAD) technique (MCSeEd), a reduced-representation, reference-free, cost-effective approach for characterizing whole genome methylation patterns across different methylation contexts (e.g., CG, CHG, CHH, 6 mA). MCSeEd can also detect genetic variations among hundreds of samples. MCSeEd is based on parallel restrictions carried out by combinations of methylation insensitive and sensitive endonucleases, followed by next-generation sequencing. Moreover, we present a robust bioinformatic pipeline (available at https://bitbucket.org/capemaster/mcseed/src/master/ ) for differential methylation analysis combined with single nucleotide polymorphism calling without or with a reference genome.
Collapse
|
16
|
Wollenberg Valero KC, Marshall JC, Bastiaans E, Caccone A, Camargo A, Morando M, Niemiller ML, Pabijan M, Russello MA, Sinervo B, Werneck FP, Sites JW, Wiens JJ, Steinfartz S. Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians. Genes (Basel) 2019; 10:genes10090646. [PMID: 31455040 PMCID: PMC6769790 DOI: 10.3390/genes10090646] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Revised: 07/21/2019] [Accepted: 08/05/2019] [Indexed: 12/22/2022] Open
Abstract
In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.
Collapse
Affiliation(s)
| | - Jonathon C Marshall
- Department of Zoology, Weber State University, 1415 Edvalson Street, Dept. 2505, Ogden, UT 84401, USA
| | - Elizabeth Bastiaans
- Department of Biology, State University of New York, College at Oneonta, Oneonta, NY 13820, USA
| | - Adalgisa Caccone
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 06520, USA
| | - Arley Camargo
- Centro Universitario de Rivera, Universidad de la República, Ituzaingó 667, Rivera 40000, Uruguay
| | - Mariana Morando
- Instituto Patagónico para el Estudio de los Ecosistemas Continentales (IPEEC, CENPAT-CONICET) Bv. Brown 2915, Puerto Madryn U9120ACD, Argentina
| | - Matthew L Niemiller
- Department of Biological Sciences, The University of Alabama in Huntsville, Huntsville, AL 35899, USA
| | - Maciej Pabijan
- Department of Comparative Anatomy, Institute of Zoology and Biomedical Research, Jagiellonian University, ul. Gronostajowa 9, 30-387 Kraków, Poland
| | - Michael A Russello
- Department of Biology, University of British Columbia, Okanagan Campus, 3247 University Way, Kelowna, BC V1V 1V7, Canada
| | - Barry Sinervo
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Coastal Biology Building, 130 McAllister Way, Santa Cruz, CA 95060, USA
| | - Fernanda P Werneck
- Programa de Coleções Científicas Biológicas, Coordenação de Biodiversidade, Instituto Nacional de Pesquisas da Amazônia, Manaus 69060-000, Brazil
| | - Jack W Sites
- Department of Biological and Marine Sciences, University of Hull, Cottingham Road, Hull HU6 7RX, UK
| | - John J Wiens
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721, USA
| | - Sebastian Steinfartz
- Molecular Evolution and Systematics of Animals, Institute of Biology, University of Leipzig, Talstrasse 33, 04103 Leipzig, Germany
| |
Collapse
|
17
|
Schield DR, Perry BW, Adams RH, Card DC, Jezkova T, Pasquesi GIM, Nikolakis ZL, Row K, Meik JM, Smith CF, Mackessy SP, Castoe TA. Allopatric divergence and secondary contact with gene flow: a recurring theme in rattlesnake speciation. Biol J Linn Soc Lond 2019. [DOI: 10.1093/biolinnean/blz077] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Abstract
The study of recently diverged lineages whose geographical ranges come into contact can provide insight into the early stages of speciation and the potential roles of reproductive isolation in generating and maintaining species. Such insight can also be important for understanding the strategies and challenges for delimiting species within recently diverged species complexes. Here, we use mitochondrial and nuclear genetic data to study population structure, gene flow and demographic history across a geographically widespread rattlesnake clade, the western rattlesnake species complex (Crotalus cerberus, Crotalus viridis, Crotalus oreganus and relatives), which contains multiple lineages with ranges that overlap geographically or contact one another. We find evidence that the evolutionary history of this group does not conform to a bifurcating tree model and that pervasive gene flow has broadly influenced patterns of present-day genetic diversity. Our results suggest that lineage diversity has been shaped largely by drift and divergent selection in isolation, followed by secondary contact, in which reproductive isolating mechanisms appear weak and insufficient to prevent introgression, even between anciently diverged lineages. The complexity of divergence and secondary contact with gene flow among lineages also provides new context for why delimiting species within this complex has been difficult and contentious historically.
Collapse
Affiliation(s)
- Drew R Schield
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
| | - Blair W Perry
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
| | - Richard H Adams
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
| | - Daren C Card
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
- Museum of Comparative Zoology, Harvard University, Cambridge, MA, USA
| | - Tereza Jezkova
- Department of Biology, Miami University of Ohio, Oxford, OH, USA
| | | | | | - Kristopher Row
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
| | - Jesse M Meik
- Department of Biological Sciences, Tarleton State University, Stephenville, TX, USA
| | - Cara F Smith
- School of Biological Sciences, University of Northern Colorado, Greeley, CO, USA
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, Greeley, CO, USA
| | - Todd A Castoe
- Department of Biology, University of Texas at Arlington, Arlington, TX, USA
| |
Collapse
|
18
|
Adams RH, Schield DR, Castoe TA. Recent Advances in the Inference of Gene Flow from Population Genomic Data. ACTA ACUST UNITED AC 2019. [DOI: 10.1007/s40610-019-00120-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
19
|
Provost KL, Mauck WM, Smith BT. Genomic divergence in allopatric Northern Cardinals of the North American warm deserts is linked to behavioral differentiation. Ecol Evol 2018; 8:12456-12478. [PMID: 30619558 PMCID: PMC6309012 DOI: 10.1002/ece3.4596] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2018] [Revised: 09/06/2018] [Accepted: 09/07/2018] [Indexed: 01/22/2023] Open
Abstract
Biogeographic barriers are considered important in initiating speciation through geographic isolation, but they rarely indiscriminately and completely reduce gene flow across entire communities. Explicitly demonstrating which factors are associated with gene-flow levels across barriers would help elucidate how speciation is initiated and isolation maintained. Here, we investigated the association of behavioral isolation on population differentiation in Northern Cardinals (Cardinalis cardinalis) distributed across the Cochise Filter Barrier, a region of transitional habitat which separates the Sonoran and Chihuahuan deserts of North America. Using genomewide markers, we modeled demographic history by fitting the data to isolation and isolation-with-migration models. The best-fit model indicated that desert populations diverged in the Pleistocene with low, historic, and asymmetric gene flow across the barrier. We then tested behavioral isolation using reciprocal call-broadcast experiments to compare song recognition between deserts, controlling for song dialect changes within deserts. We found that male Northern Cardinals in both deserts were most aggressive to local songs and failed to recognize across-barrier songs. A correlation of genomic differentiation and strong song discrimination is consistent with a model where speciation is initiated across a barrier and maintained by behavioral isolation.
Collapse
Affiliation(s)
- Kaiya L. Provost
- Department of OrnithologyAmerican Museum of Natural HistoryNew YorkNew York
- Department of Ecology, Evolution, and Environmental BiologyColumbia UniversityNew YorkNew York
- Richard Gilder Graduate SchoolAmerican Museum of Natural HistoryNew YorkNew York
| | - William M. Mauck
- Department of OrnithologyAmerican Museum of Natural HistoryNew YorkNew York
- Present address:
New York Genome CenterNew YorkNew York
| | | |
Collapse
|
20
|
Card DC, Perry BW, Adams RH, Schield DR, Young AS, Andrew AL, Jezkova T, Pasquesi GI, Hales NR, Walsh MR, Rochford MR, Mazzotti FJ, Hart KM, Hunter ME, Castoe TA. Novel ecological and climatic conditions drive rapid adaptation in invasive Florida Burmese pythons. Mol Ecol 2018; 27:4744-4757. [DOI: 10.1111/mec.14885] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2018] [Revised: 08/24/2018] [Accepted: 09/14/2018] [Indexed: 12/30/2022]
Affiliation(s)
- Daren C. Card
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Blair W. Perry
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Richard H. Adams
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Drew R. Schield
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Acacia S. Young
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Audra L. Andrew
- Department of Biology The University of Texas at Arlington Arlington Texas
| | | | | | - Nicole R. Hales
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Matthew R. Walsh
- Department of Biology The University of Texas at Arlington Arlington Texas
| | - Michael R. Rochford
- Fort Lauderdale Research & Education Center University of Florida Fort Lauderdale Florida
| | - Frank J. Mazzotti
- Fort Lauderdale Research & Education Center University of Florida Fort Lauderdale Florida
| | - Kristen M. Hart
- U. S. Geological Survey Wetland and Aquatic Research Center Davie Florida
| | - Margaret E. Hunter
- U. S. Geological Survey Wetland and Aquatic Research Center Gainesville Florida
| | - Todd A. Castoe
- Department of Biology The University of Texas at Arlington Arlington Texas
| |
Collapse
|
21
|
Schield DR, Adams RH, Card DC, Corbin AB, Jezkova T, Hales NR, Meik JM, Perry BW, Spencer CL, Smith LL, García GC, Bouzid NM, Strickland JL, Parkinson CL, Borja M, Castañeda-Gaytán G, Bryson RW, Flores-Villela OA, Mackessy SP, Castoe TA. Cryptic genetic diversity, population structure, and gene flow in the Mojave rattlesnake (Crotalus scutulatus). Mol Phylogenet Evol 2018; 127:669-681. [DOI: 10.1016/j.ympev.2018.06.013] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 05/30/2018] [Accepted: 06/06/2018] [Indexed: 10/28/2022]
|
22
|
Rutherford S, Rossetto M, Bragg JG, McPherson H, Benson D, Bonser SP, Wilson PG. Speciation in the presence of gene flow: population genomics of closely related and diverging Eucalyptus species. Heredity (Edinb) 2018; 121:126-141. [PMID: 29632325 DOI: 10.1038/s41437-018-0073-2] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Revised: 03/12/2018] [Accepted: 03/12/2018] [Indexed: 01/22/2023] Open
Abstract
Speciation is a complex process that is fundamental to the origins of biological diversity. While there has been considerable progress in our understanding of speciation, there are still many unanswered questions, especially regarding barriers to gene flow in diverging populations. Eucalyptus is an appropriate system for investigating speciation mechanisms since it comprises species that are rapidly evolving across heterogeneous environments. We examined patterns of genetic variation within and among six closely related Eucalyptus species in subgenus Eucalyptus section Eucalyptus in south-eastern Australia (commonly known as the "green ashes"). We used reduced representation genome sequencing to genotype samples from populations across altitudinal and latitudinal gradients. We found one species, Eucalyptus cunninghamii, to be highly genetically differentiated from the others, and a population of mallees from Mount Banks to be genetically distinct and therefore likely to be a new undescribed species. Only modest levels of differentiation were found between all other species in the study. There was population structure within some species (e.g., E. obstans) corresponding to geographical factors, indicating that vicariance may have played a role in the evolution of the group. Overall, we found that lineages within the green ashes are differentiated to varying extents, from strongly diverged to much earlier stages of the speciation continuum. Furthermore, our results suggest the green ashes represent a group where a range of mechanisms (e.g., reticulate evolution and vicariance) have been operating in concert. These findings not only offer insights into recent speciation mechanisms in Eucalyptus, but also other species complexes.
Collapse
Affiliation(s)
- Susan Rutherford
- Evolution and Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW, Kensington, Sydney, Australia. .,National Herbarium of NSW, Royal Botanic Garden Sydney, Sydney, Australia.
| | - Maurizio Rossetto
- National Herbarium of NSW, Royal Botanic Garden Sydney, Sydney, Australia
| | - Jason G Bragg
- National Herbarium of NSW, Royal Botanic Garden Sydney, Sydney, Australia
| | - Hannah McPherson
- National Herbarium of NSW, Royal Botanic Garden Sydney, Sydney, Australia
| | - Doug Benson
- National Herbarium of NSW, Royal Botanic Garden Sydney, Sydney, Australia
| | - Stephen P Bonser
- Evolution and Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW, Kensington, Sydney, Australia
| | - Peter G Wilson
- National Herbarium of NSW, Royal Botanic Garden Sydney, Sydney, Australia
| |
Collapse
|
23
|
Kearns AM, Restani M, Szabo I, Schrøder-Nielsen A, Kim JA, Richardson HM, Marzluff JM, Fleischer RC, Johnsen A, Omland KE. Genomic evidence of speciation reversal in ravens. Nat Commun 2018; 9:906. [PMID: 29500409 PMCID: PMC5834606 DOI: 10.1038/s41467-018-03294-w] [Citation(s) in RCA: 85] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2016] [Accepted: 02/01/2018] [Indexed: 01/10/2023] Open
Abstract
Many species, including humans, have emerged via complex reticulate processes involving hybridisation. Under certain circumstances, hybridisation can cause distinct lineages to collapse into a single lineage with an admixed mosaic genome. Most known cases of such 'speciation reversal' or 'lineage fusion' involve recently diverged lineages and anthropogenic perturbation. Here, we show that in western North America, Common Ravens (Corvus corax) have admixed mosaic genomes formed by the fusion of non-sister lineages ('California' and 'Holarctic') that diverged ~1.5 million years ago. Phylogenomic analyses and concordant patterns of geographic structuring in mtDNA, genome-wide SNPs and nuclear introns demonstrate long-term admixture and random interbreeding between the non-sister lineages. In contrast, our genomic data support reproductive isolation between Common Ravens and Chihuahuan Ravens (C. cryptoleucus) despite extensive geographic overlap and a sister relationship between Chihuahuan Ravens and the California lineage. These data suggest that the Common Raven genome was formed by secondary lineage fusion and most likely represents a case of ancient speciation reversal that occurred without anthropogenic causes.
Collapse
Affiliation(s)
- Anna M Kearns
- Natural History Museum, University of Oslo, P.O. Box 1172 Blindern, 0318, Oslo, Norway.
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA.
- Center for Conservation Genomics, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, 20013-7012, DC, USA.
| | - Marco Restani
- Department of Biological Sciences, St. Cloud State University, 720 Fourth Avenue, St. Cloud, MN, 56301-4498, USA
| | - Ildiko Szabo
- Cowan Tetrapod Collection, Beaty Biodiversity Museum, University of British Columbia, 2212 Main Mall, Vancouver, BC, V6T 1Z4, Canada
| | | | - Jin Ah Kim
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA
| | - Hayley M Richardson
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA
| | - John M Marzluff
- School of Environmental and Forest Sciences, University of Washington, Box 352100, Seattle, WA, 98195, USA
| | - Robert C Fleischer
- Center for Conservation Genomics, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, 20013-7012, DC, USA
| | - Arild Johnsen
- Natural History Museum, University of Oslo, P.O. Box 1172 Blindern, 0318, Oslo, Norway
| | - Kevin E Omland
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA
| |
Collapse
|
24
|
Adams RH, Schield DR, Card DC, Blackmon H, Castoe TA. GppFst: genomic posterior predictive simulations of FST and dXY for identifying outlier loci from population genomic data. Bioinformatics 2018; 33:1414-1415. [PMID: 28453670 DOI: 10.1093/bioinformatics/btw795] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2016] [Accepted: 12/13/2016] [Indexed: 11/12/2022] Open
Abstract
Summary We introduce GppFst, an open source R package that generates posterior predictive distributions of FST and dx under a neutral coalescent model to identify putative targets of selection from genomic data. Availability and Implementation GppFst is available at ( https://github.com/radamsRHA/GppFst ). Contact todd.castoe@uta.edu. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Richard H Adams
- Department of Biology, The University of Texas at Arlington, Arlington, TX 76019, USA
| | - Drew R Schield
- Department of Biology, The University of Texas at Arlington, Arlington, TX 76019, USA
| | - Daren C Card
- Department of Biology, The University of Texas at Arlington, Arlington, TX 76019, USA
| | - Heath Blackmon
- Department of Ecology, Evolution & Behavior, University of Minnesota, Saint Paul, MN 55108, USA
| | - Todd A Castoe
- Department of Biology, The University of Texas at Arlington, Arlington, TX 76019, USA
| |
Collapse
|
25
|
Comparative phylogeography clarifies the complexity and problems of continental distribution that drove A. R. Wallace to favor islands. Proc Natl Acad Sci U S A 2017; 113:7970-7. [PMID: 27432953 DOI: 10.1073/pnas.1601072113] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Deciphering the geographic context of diversification and distributional dynamics in continental biotas has long been an interest of biogeographers, ecologists, and evolutionary biologists. Thirty years ago, the approach now known as comparative phylogeography was introduced in a landmark study of a continental biota. Here, I use a set of 455 studies to explore the current scope of continental comparative phylogeography, including geographic, conceptual, temporal, ecological, and genomic attributes. Geographically, studies are more frequent in the northern hemisphere, but the south is catching up. Most studies focus on a Quaternary timeframe, but the Neogene is well represented. As such, explanations for geographic structure and history include geological and climatic events in Earth history, and responses include vicariance, dispersal, and range contraction-expansion into and out of refugia. Focal taxa are biased toward terrestrial or semiterrestrial vertebrates, although plants and invertebrates are well represented in some regions. The use of various kinds of nuclear DNA markers is increasing, as are multiple locus studies, but use of organelle DNA is not decreasing. Species distribution models are not yet widely incorporated into studies. In the future, continental comparative phylogeographers will continue to contribute to erosion of the simple vicariance vs. dispersal paradigm, including exposure of the widespread nature of temporal pseudocongruence and its implications for models of diversification; provide new templates for addressing a variety of ecological and evolutionary traits; and develop closer working relationships with earth scientists and biologists in a variety of disciplines.
Collapse
|
26
|
O'Connell KA, Streicher JW, Smith EN, Fujita MK. Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes. Mol Ecol 2017; 26:5729-5751. [PMID: 28802078 DOI: 10.1111/mec.14295] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2016] [Revised: 07/17/2017] [Accepted: 07/24/2017] [Indexed: 01/08/2023]
Abstract
Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.
Collapse
Affiliation(s)
- Kyle A O'Connell
- Department of Biology, The University of Texas at Arlington, Arlington, TX, USA.,The Amphibian and Reptile Diversity Research Center, University of Texas at Arlington, Arlington, TX, USA
| | | | - Eric N Smith
- Department of Biology, The University of Texas at Arlington, Arlington, TX, USA.,The Amphibian and Reptile Diversity Research Center, University of Texas at Arlington, Arlington, TX, USA
| | - Matthew K Fujita
- Department of Biology, The University of Texas at Arlington, Arlington, TX, USA.,The Amphibian and Reptile Diversity Research Center, University of Texas at Arlington, Arlington, TX, USA
| |
Collapse
|
27
|
Gamble T, Castoe TA, Nielsen SV, Banks JL, Card DC, Schield DR, Schuett GW, Booth W. The Discovery of XY Sex Chromosomes in a Boa and Python. Curr Biol 2017; 27:2148-2153.e4. [DOI: 10.1016/j.cub.2017.06.010] [Citation(s) in RCA: 63] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2016] [Revised: 04/21/2017] [Accepted: 06/02/2017] [Indexed: 10/19/2022]
|
28
|
Schield DR, Adams RH, Card DC, Perry BW, Pasquesi GM, Jezkova T, Portik DM, Andrew AL, Spencer CL, Sanchez EE, Fujita MK, Mackessy SP, Castoe TA. Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes. Ecol Evol 2017; 7:3951-3966. [PMID: 28616190 PMCID: PMC5468163 DOI: 10.1002/ece3.2996] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Accepted: 03/21/2017] [Indexed: 01/03/2023] Open
Abstract
Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.
Collapse
Affiliation(s)
- Drew R Schield
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Richard H Adams
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Daren C Card
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Blair W Perry
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Giulia M Pasquesi
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Tereza Jezkova
- Department of Ecology and Evolutionary Biology University of Arizona Tucson AZ USA
| | - Daniel M Portik
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Audra L Andrew
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Carol L Spencer
- Museum of Vertebrate Zoology University of California Berkeley CA USA
| | - Elda E Sanchez
- National Natural Toxins Research Center and Department of Chemistry Texas A&M University Kingsville Kingsville TX USA
| | - Matthew K Fujita
- Department of Biology The University of Texas at Arlington Arlington TX USA
| | - Stephen P Mackessy
- School of Biological Sciences University of Northern Colorado Greeley CO USA
| | - Todd A Castoe
- Department of Biology The University of Texas at Arlington Arlington TX USA
| |
Collapse
|
29
|
Shortt JA, Card DC, Schield DR, Liu Y, Zhong B, Castoe TA, Carlton EJ, Pollock DD. Whole Genome Amplification and Reduced-Representation Genome Sequencing of Schistosoma japonicum Miracidia. PLoS Negl Trop Dis 2017; 11:e0005292. [PMID: 28107347 PMCID: PMC5287463 DOI: 10.1371/journal.pntd.0005292] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2016] [Revised: 02/01/2017] [Accepted: 12/29/2016] [Indexed: 11/18/2022] Open
Abstract
Background In areas where schistosomiasis control programs have been implemented, morbidity and prevalence have been greatly reduced. However, to sustain these reductions and move towards interruption of transmission, new tools for disease surveillance are needed. Genomic methods have the potential to help trace the sources of new infections, and allow us to monitor drug resistance. Large-scale genotyping efforts for schistosome species have been hindered by cost, limited numbers of established target loci, and the small amount of DNA obtained from miracidia, the life stage most readily acquired from humans. Here, we present a method using next generation sequencing to provide high-resolution genomic data from S. japonicum for population-based studies. Methodology/Principal Findings We applied whole genome amplification followed by double digest restriction site associated DNA sequencing (ddRADseq) to individual S. japonicum miracidia preserved on Whatman FTA cards. We found that we could effectively and consistently survey hundreds of thousands of variants from 10,000 to 30,000 loci from archived miracidia as old as six years. An analysis of variation from eight miracidia obtained from three hosts in two villages in Sichuan showed clear population structuring by village and host even within this limited sample. Conclusions/Significance This high-resolution sequencing approach yields three orders of magnitude more information than microsatellite genotyping methods that have been employed over the last decade, creating the potential to answer detailed questions about the sources of human infections and to monitor drug resistance. Costs per sample range from $50-$200, depending on the amount of sequence information desired, and we expect these costs can be reduced further given continued reductions in sequencing costs, improvement of protocols, and parallelization. This approach provides new promise for using modern genome-scale sampling to S. japonicum surveillance, and could be applied to other schistosome species and other parasitic helminthes. Schistosomiasis is a devastating tropical disease that affects more than 200 million people worldwide. Over the past several decades, transmission control strategies implemented in China have reduced the prevalence and morbidity of Schistosoma japonicum in many areas. Infections still persist, however, and it is therefore of great interest to determine the sources of recurring infections. Genetic analysis is a promising means to achieve this. Towards this aim, we conducted a pilot study to assess the feasibility of using high-throughput sequencing to assess the geographic distribution of schistosome genetic variants. Because DNA yields from miracidia, the most easily accessible life stage, are insufficient for high throughput sequencing, we first employed whole genome amplification to obtain sufficient quantities of DNA. We then employed a technique that reproducibly sequences the same fraction of a genome across numerous samples. We successfully sequenced 6-year old S. japonicum samples from Sichuan Province, China, easily and economically identifying tens of thousands of variable loci, a sufficient number to discriminate fine-scale population structure. Further population sampling will help answer important questions concerning the persistence of infections, the sources of new infections, and whether parasite populations have undergone incipient evolution of drug resistance.
Collapse
Affiliation(s)
- Jonathan A. Shortt
- Department of Biochemistry & Molecular Genetics, University of Colorado School of Medicine, Aurora, CO, United States of America
| | - Daren C. Card
- Department of Biology, University of Texas Arlington, Arlington, TX, United States of America
| | - Drew R. Schield
- Department of Biology, University of Texas Arlington, Arlington, TX, United States of America
| | - Yang Liu
- Institute of Parasitic Disease, Sichuan Center for Disease Control and Prevention, Chengdu, The People’s Republic of China
| | - Bo Zhong
- Institute of Parasitic Disease, Sichuan Center for Disease Control and Prevention, Chengdu, The People’s Republic of China
| | - Todd A. Castoe
- Department of Biology, University of Texas Arlington, Arlington, TX, United States of America
| | - Elizabeth J. Carlton
- Department of Environmental and Occupational Health, University of Colorado, Colorado School of Public Health, Aurora, CO, United States of America
| | - David D. Pollock
- Department of Biochemistry & Molecular Genetics, University of Colorado School of Medicine, Aurora, CO, United States of America
- * E-mail:
| |
Collapse
|
30
|
Burns M, Starrett J, Derkarabetian S, Richart CH, Cabrero A, Hedin M. Comparative performance of double‐digest
RAD
sequencing across divergent arachnid lineages. Mol Ecol Resour 2016; 17:418-430. [DOI: 10.1111/1755-0998.12575] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2015] [Accepted: 06/23/2016] [Indexed: 11/28/2022]
Affiliation(s)
- Mercedes Burns
- Department of Biology San Diego State University 5500 Campanile Drive San Diego CA 92182 USA
| | - James Starrett
- Department of Biology San Diego State University 5500 Campanile Drive San Diego CA 92182 USA
| | - Shahan Derkarabetian
- Department of Biology San Diego State University 5500 Campanile Drive San Diego CA 92182 USA
- Department of Biology University of California 900 University Avenue Riverside CA 92521 USA
| | - Casey H. Richart
- Department of Biology San Diego State University 5500 Campanile Drive San Diego CA 92182 USA
- Department of Biology University of California 900 University Avenue Riverside CA 92521 USA
| | - Allan Cabrero
- Department of Biology San Diego State University 5500 Campanile Drive San Diego CA 92182 USA
| | - Marshal Hedin
- Department of Biology San Diego State University 5500 Campanile Drive San Diego CA 92182 USA
| |
Collapse
|
31
|
Malukiewicz J, Guschanski K, Grativol AD, Oliveira MAB, Ruiz-Miranda CR, Stone AC. Application of PE-RADSeq to the study of genomic diversity and divergence of two Brazilian marmoset species (Callithrix jacchusandC. penicillata). Am J Primatol 2016; 79:1-12. [DOI: 10.1002/ajp.22587] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2015] [Revised: 07/03/2016] [Accepted: 07/10/2016] [Indexed: 12/21/2022]
Affiliation(s)
- Joanna Malukiewicz
- Departamento de Bioquímica e Biologia Molecular; Universidade Federal de Viçosa; Viçosa MG Brazil
- School of Life Sciences; Arizona State University; Tempe Arizona
| | - Katerina Guschanski
- Department of Animal Ecology, Evolutionary Biology Centre; Uppsala University; Uppsala Sweden
| | - Adriana D. Grativol
- Laboratório de Ciências Ambientais, Centro de Biociências e Biotecnologia; Universidade Estadual do Norte Fluminense; Campos dos Goytacazes RJ Brazil
| | - Maria Adélia B. Oliveira
- Departamento de Morfologia e Fisiologia Animal; Universidade Federal Rural de Pernambuco; Recife PE Brazil
| | - Carlos R. Ruiz-Miranda
- Laboratório de Ciências Ambientais, Centro de Biociências e Biotecnologia; Universidade Estadual do Norte Fluminense; Campos dos Goytacazes RJ Brazil
| | - Anne C. Stone
- School of Human Evolution and Social Change; Arizona State University; Tempe Arizona
| |
Collapse
|
32
|
Card DC, Schield DR, Adams RH, Corbin AB, Perry BW, Andrew AL, Pasquesi GIM, Smith EN, Jezkova T, Boback SM, Booth W, Castoe TA. Phylogeographic and population genetic analyses reveal multiple species of Boa and independent origins of insular dwarfism. Mol Phylogenet Evol 2016; 102:104-16. [PMID: 27241629 DOI: 10.1016/j.ympev.2016.05.034] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2015] [Revised: 05/05/2016] [Accepted: 05/26/2016] [Indexed: 10/21/2022]
Abstract
Boa is a Neotropical genus of snakes historically recognized as monotypic despite its expansive distribution. The distinct morphological traits and color patterns exhibited by these snakes, together with the wide diversity of ecosystems they inhabit, collectively suggest that the genus may represent multiple species. Morphological variation within Boa also includes instances of dwarfism observed in multiple offshore island populations. Despite this substantial diversity, the systematics of the genus Boa has received little attention until very recently. In this study we examined the genetic structure and phylogenetic relationships of Boa populations using mitochondrial sequences and genome-wide SNP data obtained from RADseq. We analyzed these data at multiple geographic scales using a combination of phylogenetic inference (including coalescent-based species delimitation) and population genetic analyses. We identified extensive population structure across the range of the genus Boa and multiple lines of evidence for three widely-distributed clades roughly corresponding with the three primary land masses of the Western Hemisphere. We also find both mitochondrial and nuclear support for independent origins and parallel evolution of dwarfism on offshore island clusters in Belize and Cayos Cochinos Menor, Honduras.
Collapse
Affiliation(s)
- Daren C Card
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Drew R Schield
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Richard H Adams
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Andrew B Corbin
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Blair W Perry
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Audra L Andrew
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Giulia I M Pasquesi
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Eric N Smith
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA
| | - Tereza Jezkova
- Department of Ecology & Evolutionary Biology, University of Arizona, P.O. Box 210088, Tucson, AZ 85721, USA
| | - Scott M Boback
- Department of Biology, P.O. Box 1773, Dickinson College, Carlisle, PA 17013, USA
| | - Warren Booth
- Department of Biological Science, 800 South Tucker Drive, University of Tulsa, Tulsa, OK 74104, USA
| | - Todd A Castoe
- Department of Biology, 501 S. Nedderman Drive, University of Texas at Arlington, Arlington, TX 76019, USA.
| |
Collapse
|
33
|
Manthey JD, Robbins MB, Moyle RG. A genomic investigation of the putative contact zone between divergent Brown Creeper (Certhia americana) lineages: chromosomal patterns of genetic differentiation. Genome 2015; 59:115-25. [PMID: 26794151 DOI: 10.1139/gen-2015-0093] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Sky islands, or montane forest separated by different lowland habitats, are highly fragmented regions that potentially limit gene flow between isolated populations. In the sky islands of the Madrean Archipelago (Arizona, USA), various taxa display different phylogeographic patterns, from unrestricted gene flow among sky islands to complex patterns with multiple distinct lineages. Using genomic-level approaches allows the investigation of differential patterns of gene flow, selection, and genetic differentiation among chromosomes and specific genomic regions between sky island populations. Here, we used thousands of SNPs to investigate the putative contact zone of divergent Brown Creeper (Certhia americana) lineages in the Madrean Archipelago sky islands. We found the two lineages to be completely allopatric (during the breeding season) with a lack of hybridization and gene flow between lineages and no genetic structure among sky islands within lineages. Additionally, the two lineages inhabit different climatic and ecosystem conditions and have many local primary song dialects in the southern Arizona mountain ranges. We identified a positive relationship between genetic differentiation and chromosome size, but the sex chromosome (Z) was not found to be an outlier. Differential patterns of genetic differentiation per chromosome may be explained by genetic drift--possibly in conjunction with non-random mating and non-random gene flow--due to variance in recombination rates among chromosomes.
Collapse
Affiliation(s)
- Joseph D Manthey
- a Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA.,b Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA
| | - Mark B Robbins
- a Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA
| | - Robert G Moyle
- a Biodiversity Institute, University of Kansas, Lawrence, KS 66045, USA.,b Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66045, USA
| |
Collapse
|
34
|
Moseley MA, Cox CL, Streicher JW, Roelke CE, Chippindale PT. Phylogeography and lineage-specific patterns of genetic diversity and molecular evolution in a group of North American skinks. Biol J Linn Soc Lond 2015. [DOI: 10.1111/bij.12626] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Matthew A. Moseley
- Department of Biology; The University of Texas at Arlington; Arlington TX 76010 USA
| | - Christian L. Cox
- Department of Biology; The University of Texas at Arlington; Arlington TX 76010 USA
- Department of Biology; The University of Virginia; Charlottesville VA 22903 USA
- Department of Biology; Georgia Southern University; Statesboro GA USA
| | - Jeffrey W. Streicher
- Department of Biology; The University of Texas at Arlington; Arlington TX 76010 USA
- Department of Life Sciences; The Natural History Museum; London SW7 5BD UK
| | - Corey E. Roelke
- Department of Biology; The University of Texas at Arlington; Arlington TX 76010 USA
| | - Paul T. Chippindale
- Department of Biology; The University of Texas at Arlington; Arlington TX 76010 USA
| |
Collapse
|
35
|
Blair C, Campbell CR, Yoder AD. Assessing the utility of whole genome amplified DNA for next-generation molecular ecology. Mol Ecol Resour 2015; 15:1079-90. [PMID: 25619406 DOI: 10.1111/1755-0998.12376] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2014] [Revised: 01/16/2015] [Accepted: 01/19/2015] [Indexed: 01/23/2023]
Abstract
DNA quantity can be a hindrance in ecological and evolutionary research programmes due to a range of factors including endangered status of target organisms, available tissue type, and the impact of field conditions on preservation methods. A potential solution to low-quantity DNA lies in whole genome amplification (WGA) techniques that can substantially increase DNA yield. To date, few studies have rigorously examined sequence bias that might result from WGA and next-generation sequencing of nonmodel taxa. To address this knowledge deficit, we use multiple displacement amplification (MDA) and double-digest RAD sequencing on the grey mouse lemur (Microcebus murinus) to quantify bias in genome coverage and SNP calls when compared to raw genomic DNA (gDNA). We focus our efforts in providing baseline estimates of potential bias by following manufacturer's recommendations for starting DNA quantities (>100 ng). Our results are strongly suggestive that MDA enrichment does not introduce systematic bias to genome characterization. SNP calling between samples when genotyping both de-novo and with a reference genome are highly congruent (>98%) when specifying a minimum threshold of 20X stack depth to call genotypes. Relative genome coverage is also similar between MDA and gDNA, and allelic dropout is not observed. SNP concordance varies based on coverage threshold, with 95% concordance reached at ~12X coverage genotyping de-novo and ~7X coverage genotyping with the reference genome. These results suggest that MDA may be a suitable solution for next-generation molecular ecological studies when DNA quantity would otherwise be a limiting factor.
Collapse
Affiliation(s)
- Christopher Blair
- Department of Biology, Duke University, Box 90338, BioSci 130 Science Drive, Durham, NC, 27708, USA
| | - C Ryan Campbell
- Department of Biology, Duke University, Box 90338, BioSci 130 Science Drive, Durham, NC, 27708, USA
| | - Anne D Yoder
- Department of Biology, Duke University, Box 90338, BioSci 130 Science Drive, Durham, NC, 27708, USA
| |
Collapse
|