skip to main content


Title: The influence of history, geography and environment on patterns of diversification in the western terrestrial garter snake
Abstract Aim

A central aim of biogeography is to understand how biodiversity is generated and maintained across landscapes. Here, we establish phylogenetic and population genetic patterns in a widespread reptile to quantify the influence of historical biogeography and current environmental variation on patterns of genetic diversity.

Location

Western North America.

Taxon

Western terrestrial garter snake,Thamnophis elegans.

Methods

We used double‐digest RADseq to estimate phylogenetic relationships and characterize population genetic structure across the three widespread subspecies ofTelegans:T. e. vagrans(wandering garter snake),Teelegans(mountain garter snake) andTeterrestris(coast garter snake). We assessed patterns of dispersal and vicariance across biogeographic regions using ancestral area reconstruction (AAR) and deviations from isolation‐by‐distance across the landscape using estimated effective migration surfaces (EEMS). We identified environmental variables potentially shaping local adaptation in regional lineages using genetic‐environment association (GEA) analyses.

Results

We recovered three well‐differentiated genetic groups that correspond to the three subspecies. AAR analyses inferred the eastern Cascade Range as the ancestral area, with dispersal to both the east and west across western North America. Populations ofT. e. elegansdisplayed a latitudinal gradient in genetic variation across the Sierra Nevada and northern California, while populations ofTeterrestrisshow discrete genetic breaks consistent with well‐known biogeographic barriers. Lastly, GEA analyses identified allele frequency shifts at loci associated with a common set of environmental variables in bothTeelegansandTeterrestris.

Main Conclusion

T. elegansis composed of distinct evolutionary lineages, each with its own geographic range and history of diversification.TeelegansandTeterrestrisshow unique patterns of diversification as populations dispersed from east to west and while adapting to the new environments they colonized. Historical events, landscape features and environmental variation have all contributed to patterns of differentiation inTelegans.

 
more » « less
NSF-PAR ID:
10247801
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Journal of Biogeography
Volume:
48
Issue:
9
ISSN:
0305-0270
Page Range / eLocation ID:
p. 2226-2245
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Aim

    Species with wide distributions spanning the African Guinean and Congolian rain forests are often composed of genetically distinct populations or cryptic species with geographic distributions that mirror the locations of the remaining forest habitats. We used phylogeographic inference and demographic model testing to evaluate diversification models in a widespread rain forest species, the African foam‐nest treefrogChiromantis rufescens.

    Location

    Guinean and Congolian rain forests, West and Central Africa.

    Taxon

    Chiromantis rufescens.

    Methods

    We collected mitochondrial DNA (mtDNA) and single‐nucleotide polymorphism (SNP) data for 130 samples ofC. rufescens. After estimating population structure and inferring species trees using coalescent methods, we tested demographic models to evaluate alternative population divergence histories that varied with respect to gene flow, population size change and periods of isolation and secondary contact. Species distribution models were used to identify the regions of climatic stability that could have served as forest refugia since the last interglacial.

    Results

    Population structure withinC. rufescensresembles the major biogeographic regions of the Guinean and Congolian forests. Coalescent‐based phylogenetic analyses provide strong support for an early divergence between the western Upper Guinean forest and the remaining populations. Demographic inferences support diversification models with gene flow and population size changes even in cases where contemporary populations are currently allopatric, which provides support for forest refugia and barrier models. Species distribution models suggest that forest refugia were available for each of the populations throughout the Pleistocene.

    Main conclusions

    Considering historical demography is essential for understanding population diversification, especially in complex landscapes such as those found in the Guineo–Congolian forest. Population demographic inferences help connect the patterns of genetic variation to diversification model predictions. The diversification history ofC. rufescenswas shaped by a variety of processes, including vicariance from river barriers, forest fragmentation and adaptive evolution along environmental gradients.

     
    more » « less
  2. Abstract Aim

    Current distributions of widespread North American (NA) species have been shaped by Pleistocene glacial cycles, latitudinal temperature gradients, sharp longitudinal habitat transitions and the vicariant effects of major mountain and river systems that subdivide the continent. Within these transcontinental species, genetic diversity patterns might not conform to established biogeographic breaks compared to more spatially restricted taxa due to intrinsic differences or spatiotemporal differences. In this study, we highlight the effects of these extrinsic variables on genetic structuring by investigating the phylogeographic history of a widespread generalist squamate found throughout NA.

    Location

    North America.

    Taxon

    Common gartersnake,Thamnophis sirtalis.

    Methods

    We evaluate the effects of major river basins and the forest‐grassland transition into the Interior Plains on genetic structure patterns using phylogenetic, spatially informed population structure and demographic analyses of single nucleotide polymorphism data and address range expansion history with ecological niche modelling using locality and historic climate data.

    Results

    We identify four phylogeographic lineages with varying degrees of connectivity between them. We find discordant population structure patterns between sex‐linked and autosomal loci with respect to the relationship between the central NA lineage relative to coastal lineages. We find support for southeast Pleistocene refugia where recent secondary contact occurred during the Last Glacial Maximum and evidence for both northern and southern refugia in western NA.

    Main Conclusion

    Our results provide strong evidence for a Pliocene origin forT. sirtalisin central‐southeastern NA preceding its rapid expansion across the continent prior to middle Pleistocene climate‐mediated lineage formation. We implicate major riverine networks within the Mississippi watershed in likely repeated westward expansion events across the Interior Plains. Finally, we corroborate prior conclusions that phenotypic differences between subspecies do not reflect shared evolutionary history and note that the degree of separation between inferred lineages warrants further investigation before any taxonomic revisions are proposed.

     
    more » « less
  3. Abstract Aim

    Refugial isolation during glaciation is an established driver of speciation; however, the opposing role of interglacial population expansion, secondary contact, and gene flow on the diversification process remains less understood. The consequences of glacial cycling on diversity are complex and especially so for archipelago species, which experience dramatic fluctuations in connectivity in response to both lower sea levels during glacial events and increased fragmentation during glacial recession. We test whether extended refugial isolation has led to the divergence of genetically and morphologically distinct species within Holarctic ermine (Mustela erminea), a small cosmopolitan carnivore species that harbours 34 extant subspecies, 14 of which are insular endemics.

    Location

    Holarctic.

    Methods

    We use genetic sequences (complete mitochondrial genomes, four nuclear genes) from >100 ermine (stoats) and geometric morphometric data for >200 individuals (27 of the 34 extant subspecies) from across their Holarctic range to provide an integrative perspective on diversification and endemism across this complex landscape. Multiple species delimitation methods (iBPP,bPTP) assessed congruence between morphometric and genetic data.

    Results

    Our results support the recognition of at least three species within theM. ermineacomplex, coincident with three of four genetic clades, tied to diversification in separate glacial refugia. We found substantial geographic variation within each species, with geometric morphometric results largely consistent with historical infraspecific taxonomy.

    Main conclusions

    Phylogeographic structure mirrors patterns of diversification in other Holarctic species, with a major Nearctic‐Palearctic split, but with greater intraspecific morphological diversity. Recognition of insular endemic speciesM. haidarumis consistent with a deep history of refugial persistence and highlights the urgency of mindful management of island populations along North America's North Pacific Coast. Significant environmental modification (e.g. industrial‐scale logging, mining) has been proposed for a number of these islands, which may elevate the risk of extinction of insular palaeoendemics.

     
    more » « less
  4. Abstract

    Reciprocal adaptation is the hallmark of arms race coevolution. Local coadaptation between natural enemies should generate a geographic mosaic pattern where both species have roughly matched abilities across their shared range. However, mosaic variation in ecologically relevant traits can also arise from processes unrelated to reciprocal selection, such as population structure or local environmental conditions. We tested whether these alternative processes can account for trait variation in the geographic mosaic of arms race coevolution between resistant garter snakes (Thamnophis sirtalis) and toxic newts (Taricha granulosa). We found that predator resistance and prey toxin levels are functionally matched in co-occurring populations, suggesting that mosaic variation in the armaments of both species results from the local pressures of reciprocal selection. By the same token, phenotypic and genetic variation in snake resistance deviates from neutral expectations of population genetic differentiation, showing a clear signature of adaptation to local toxin levels in newts. Contrastingly, newt toxin levels are best predicted by genetic differentiation among newt populations, and to a lesser extent, by the local environment and snake resistance. Exaggerated armaments suggest that coevolution occurs in certain hotspots, but prey population structure seems to be of particular influence on local phenotypic variation in both species throughout the geographic mosaic. Our results imply that processes other than reciprocal selection, like historical biogeography and environmental pressures, represent an important source of variation in the geographic mosaic of coevolution. Such a pattern supports the role of “trait remixing” in the geographic mosaic theory, the process by which non-adaptive forces dictate spatial variation in the interactions among species.

     
    more » « less
  5. Abstract Aim

    Islands provide opportunities for isolation and speciation. Many landmasses in the Indo‐Australian Archipelago (IAA) are oceanic islands, and founder‐event speciation is expected to be the predominant form of speciation of volant taxa on these islands. We studied the biogeographic history of flying foxes, a group with many endemic species and a predilection for islands, to test this hypothesis and infer the biogeographic origin of the group.

    Location

    Australasia, Indo‐Australian Archipelago, Madagascar, Pacific Islands.

    Taxon

    Pteropus(Pteropodidae).

    Methods

    To infer the biogeographic history ofPteropus, we sequenced up to 6,169 bp of genetic data from 10 markers and reconstructed a multilocus species tree of 34 currently recognizedPteropusspecies and subspecies with threeAcerodonoutgroups usingBEASTand subsequently estimated ancestral areas using models implemented inBioGeoBEARS.

    Results

    Species‐level resolution was occasionally low because of slow rates of molecular evolution and/or recent divergences. Older divergences, however, were more strongly supported and allow the evolutionary history of the group to be inferred. The genus diverged in Wallacea from its common ancestor withAcerodon; founder‐event speciation out of Wallacea was a common inference.Pteropusspecies in Micronesia and the western Indian Ocean were also inferred to result from founder‐event speciation.

    Main conclusions

    Dispersal between regions of the IAA and the islands found therein fostered diversification ofPteropusthroughout the IAA and beyond. Dispersal inPteropusis far higher than in most other volant taxa studied to date, highlighting the importance of inter‐island movement in the biogeographic history of this large clade of large bats.

     
    more » « less