1
|
Holding ML, Sovic MG, Colston TJ, Gibbs HL. The scales of coevolution: comparative phylogeography and genetic demography of a locally adapted venomous predator and its prey. Biol J Linn Soc Lond 2020. [DOI: 10.1093/biolinnean/blaa192] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Abstract
Coevolutionary theory predicts that differences in the genetic demography of interacting species can influence patterns of local adaptation by affecting the potential of local populations to respond to selection. We conducted a comparative phylogeographical study of venomous rattlesnakes and their venom-resistant ground squirrel prey across California, and assessed how effective population size (Ne) estimates correspond with a previously documented pattern of rattlesnake local adaptation. Using RAD sequencing markers, we detected lineage relationships among both the rattlesnakes (Crotalus oreganus ssp.) and ground squirrels (Otospermophilus sp.) that are incongruent with previous phylogenetic hypotheses. Both rattlesnakes and squirrels share a deep divergence at the Sacramento–San Joaquin River Delta. At this broad phylogeographical scale, we found that the locally adapted rattlesnakes had higher Ne than squirrels. At the population scale, snakes also had larger Ne accompanied by larger values of several metrics of population genetic diversity. However, the specific magnitude of local adaptation of venom activity to ground squirrel venom resistance was not significantly correlated with local differences in Ne or other diversity statistics between predator and prey populations, suggesting that other factors in the geographic mosaic of coevolution contribute to the specific local-scale outcomes of this interaction. These results suggest an evolutionary mechanism that may explain some (but clearly not all) of rattlesnake local adaptation in this coevolutionary interaction – larger population sizes raise the adaptive potential of rattlesnakes compared to ground squirrels.
Collapse
Affiliation(s)
- Matthew L Holding
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
- Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - Michael G Sovic
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| | - Timothy J Colston
- Department of Biological Science, Florida State University, Tallahassee, FL, USA
| | - H Lisle Gibbs
- Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, USA
| |
Collapse
|
2
|
Genome-wide genetic variation coupled with demographic and ecological niche modeling of the dusky-footed woodrat (Neotoma fuscipes) reveal patterns of deep divergence and widespread Holocene expansion across northern California. Heredity (Edinb) 2020; 126:521-536. [PMID: 33323954 DOI: 10.1038/s41437-020-00393-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2019] [Revised: 10/28/2020] [Accepted: 11/18/2020] [Indexed: 01/25/2023] Open
Abstract
Understanding how species have responded to past climate change may help refine projections of how species and biotic communities will respond to future change. Here, we integrate estimates of genome-wide genetic variation with demographic and niche modeling to investigate the historical biogeography of an important ecological engineer: the dusky-footed woodrat, Neotoma fuscipes. We use RADseq to generate a genome-wide dataset for 71 individuals from across the geographic distribution of the species in California. We estimate population structure using several model-based methods and infer the demographic history of regional populations using a site frequency spectrum-based approach. Additionally, we use ecological niche modeling to infer current and past (Last Glacial Maximum) environmental suitability across the species' distribution. Finally, we estimate the directionality and possible spatial origins of regional population expansions. Our analyses indicate this species is subdivided into three regionally distinct populations, with the deepest divergence occurring ~1.7 million years ago across the modern-day San Francisco-Bay Delta region; a common biogeographic barrier for the flora and fauna of California. Our models of environmental suitability through time coincide with our estimates of population expansion, with relative long-term stability in the southern portion of the range, and more recent expansion into the northern end of the range. Our study illustrates how the integration of genome-wide data with spatial and demographic modeling can reveal the timing and spatial extent of historic events that determine patterns of biotic diversity and may help predict biotic response to future change.
Collapse
|
3
|
Card DC, Adams RH, Schield DR, Perry BW, Corbin AB, Pasquesi GIM, Row K, Van Kleeck MJ, Daza JM, Booth W, Montgomery CE, Boback SM, Castoe TA. Genomic Basis of Convergent Island Phenotypes in Boa Constrictors. Genome Biol Evol 2020; 11:3123-3143. [PMID: 31642474 PMCID: PMC6836717 DOI: 10.1093/gbe/evz226] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/10/2019] [Indexed: 12/14/2022] Open
Abstract
Convergent evolution is often documented in organisms inhabiting isolated environments with distinct ecological conditions and similar selective regimes. Several Central America islands harbor dwarf Boa populations that are characterized by distinct differences in growth, mass, and craniofacial morphology, which are linked to the shared arboreal and feast-famine ecology of these island populations. Using high-density RADseq data, we inferred three dwarf island populations with independent origins and demonstrate that selection, along with genetic drift, has produced both divergent and convergent molecular evolution across island populations. Leveraging whole-genome resequencing data for 20 individuals and a newly annotated Boa genome, we identify four genes with evidence of phenotypically relevant protein-coding variation that differentiate island and mainland populations. The known roles of these genes involved in body growth (PTPRS, DMGDH, and ARSB), circulating fat and cholesterol levels (MYLIP), and craniofacial development (DMGDH and ARSB) in mammals link patterns of molecular evolution with the unique phenotypes of these island forms. Our results provide an important genome-wide example for quantifying expectations of selection and convergence in closely related populations. We also find evidence at several genomic loci that selection may be a prominent force of evolutionary change—even for small island populations for which drift is predicted to dominate. Overall, while phenotypically convergent island populations show relatively few loci under strong selection, infrequent patterns of molecular convergence are still apparent and implicate genes with strong connections to convergent phenotypes.
Collapse
Affiliation(s)
- Daren C Card
- Department of Biology, University of Texas Arlington.,Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts.,Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts
| | | | | | - Blair W Perry
- Department of Biology, University of Texas Arlington
| | | | | | | | | | - Juan M Daza
- Grupo Herpetológico de Antioquia, Instituto de Biología, Universidad de Antioquia, Medellín, Colombia
| | - Warren Booth
- Department of Biological Science, University of Tulsa, Oklahoma
| | | | - Scott M Boback
- Department of Biology, Dickinson College, Carlisle, Pennsylvania
| | - Todd A Castoe
- Department of Biology, University of Texas Arlington
| |
Collapse
|
4
|
Renaud G, Hanghøj K, Korneliussen TS, Willerslev E, Orlando L. Joint Estimates of Heterozygosity and Runs of Homozygosity for Modern and Ancient Samples. Genetics 2019; 212:587-614. [PMID: 31088861 PMCID: PMC6614887 DOI: 10.1534/genetics.119.302057] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2019] [Accepted: 05/01/2019] [Indexed: 11/18/2022] Open
Abstract
Both the total amount and the distribution of heterozygous sites within individual genomes are informative about the genetic diversity of the population they belong to. Detecting true heterozygous sites in ancient genomes is complicated by the generally limited coverage achieved and the presence of post-mortem damage inflating sequencing errors. Additionally, large runs of homozygosity found in the genomes of particularly inbred individuals and of domestic animals can skew estimates of genome-wide heterozygosity rates. Current computational tools aimed at estimating runs of homozygosity and genome-wide heterozygosity levels are generally sensitive to such limitations. Here, we introduce ROHan, a probabilistic method which substantially improves the estimate of heterozygosity rates both genome-wide and for genomic local windows. It combines a local Bayesian model and a Hidden Markov Model at the genome-wide level and can work both on modern and ancient samples. We show that our algorithm outperforms currently available methods for predicting heterozygosity rates for ancient samples. Specifically, ROHan can delineate large runs of homozygosity (at megabase scales) and produce a reliable confidence interval for the genome-wide rate of heterozygosity outside of such regions from modern genomes with a depth of coverage as low as 5-6× and down to 7-8× for ancient samples showing moderate DNA damage. We apply ROHan to a series of modern and ancient genomes previously published and revise available estimates of heterozygosity for humans, chimpanzees and horses.
Collapse
Affiliation(s)
- Gabriel Renaud
- Lundbeck Foundation GeoGenetics Center, Globe Institute, University of Copenhagen, 1350K, Denmark
| | - Kristian Hanghøj
- Lundbeck Foundation GeoGenetics Center, Globe Institute, University of Copenhagen, 1350K, Denmark
- Laboratoire d'Anthropobiologie Moléculaire et d'Imagerie de Synthèse, CNRS UMR 5288, Université de Toulouse, Université Paul Sabatier, 31000, France
| | | | - Eske Willerslev
- Lundbeck Foundation GeoGenetics Center, Globe Institute, University of Copenhagen, 1350K, Denmark
- Department of Zoology, University of Cambridge, CB2 3EJ, UK
- The Wellcome Trust Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SA, UK
- The Danish Institute for Advanced Study at The University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Ludovic Orlando
- Lundbeck Foundation GeoGenetics Center, Globe Institute, University of Copenhagen, 1350K, Denmark
- Laboratoire d'Anthropobiologie Moléculaire et d'Imagerie de Synthèse, CNRS UMR 5288, Université de Toulouse, Université Paul Sabatier, 31000, France
| |
Collapse
|
5
|
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
|
6
|
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]
|
7
|
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]
|