Biogeographic barriers and historic climate shape the phylogeography and demography of the common gartersnake

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.North America.Common gartersnake, Thamnophis sirtalis.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.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.Our results provide strong evidence for a Pliocene origin for T. sirtalis in 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.

[1]  D. Weisrock,et al.  The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini). , 2023, Molecular phylogenetics and evolution.

[2]  T. Parchman,et al.  Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution. , 2021, Molecular phylogenetics and evolution.

[3]  J. Maldonado,et al.  Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone , 2021, Frontiers in Ecology and Evolution.

[4]  T. Parchman,et al.  The influence of history, geography and environment on patterns of diversification in the western terrestrial garter snake , 2021, Journal of Biogeography.

[5]  S. Edwards,et al.  Evidence for the Pleistocene Arc Hypothesis from genome‐wide SNPs in a Neotropical dry forest specialist, the Rufous‐fronted Thornbird (Furnariidae: Phacellodomus rufifrons) , 2020, Molecular ecology.

[6]  E. Brodie,et al.  The geographic mosaic of arms race coevolution is closely matched to prey population structure , 2020, Evolution letters.

[7]  E. A. Myers,et al.  Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex) , 2020, bioRxiv.

[8]  W. Funk,et al.  Coalescent-based species delimitation is sensitive to geographic sampling and isolation by distance , 2020 .

[9]  E. A. Myers,et al.  Biogeographic barriers, Pleistocene refugia, and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex) , 2020, Molecular ecology.

[10]  C. J. Battey Evidence of linked selection on the Z chromosome of hybridizing hummingbirds * , 2019, Evolution; international journal of organic evolution.

[11]  B. Rannala,et al.  A Bayesian Implementation of the Multispecies Coalescent Model with Introgression for Phylogenomic Analysis , 2019, Molecular biology and evolution.

[12]  Olga Chernomor,et al.  IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era , 2019, bioRxiv.

[13]  Sami R. Cherikh,et al.  A global catalog of primary reptile type specimens. , 2019, Zootaxa.

[14]  Ren A. Thompson,et al.  Middle Pleistocene formation of the Rio Grande Gorge, San Luis Valley, south-central Colorado and north-central New Mexico, USA: Process, timing, and downstream implications , 2019, Quaternary Science Reviews.

[15]  Nicolas C. Rochette,et al.  Stacks 2: Analytical methods for paired‐end sequencing improve RADseq‐based population genomics , 2019, Molecular ecology.

[16]  A. Haywood,et al.  PaleoClim, high spatial resolution paleoclimate surfaces for global land areas , 2018, Scientific Data.

[17]  Jeffery P. Demuth,et al.  Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome , 2018, Genome biology and evolution.

[18]  M. Suchard,et al.  Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7 , 2018, Systematic biology.

[19]  Yun S. Song,et al.  Efficiently inferring the demographic history of many populations with allele count data , 2018, bioRxiv.

[20]  Jeremy M. Brown,et al.  Impact of Model Violations on the Inference of Species Boundaries Under the Multispecies Coalescent , 2018, Systematic biology.

[21]  W. Salzburger,et al.  Bayesian Divergence-Time Estimation with Genome-Wide Single-Nucleotide Polymorphism Data of Sea Catfishes (Ariidae) Supports Miocene Closure of the Panamanian Isthmus , 2018, Systematic biology.

[22]  Daniel L. Ayres,et al.  Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10 , 2018, Virus evolution.

[23]  Simon H. Martin,et al.  Patterns of Z chromosome divergence among Heliconius species highlight the importance of historical demography , 2017, bioRxiv.

[24]  Stephen E. Fick,et al.  WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas , 2017 .

[25]  A. Leaché,et al.  Discordance between genomic divergence and phenotypic variation in a rapidly evolving avian genus (Motacilla). , 2017, Molecular phylogenetics and evolution.

[26]  A. von Haeseler,et al.  UFBoot2: Improving the Ultrafast Bootstrap Approximation , 2017, bioRxiv.

[27]  C. Feldman,et al.  Convergent adaptation to dangerous prey proceeds through the same first‐step mutation in the garter snake Thamnophis sirtalis , 2017, Evolution; international journal of organic evolution.

[28]  K. Karlstrom,et al.  Birth and evolution of the Rio Grande fluvial system in the past 8 Ma: Progressive downward integration and the influence of tectonics, volcanism, and climate , 2017 .

[29]  Thomas K. F. Wong,et al.  ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates , 2017, Nature Methods.

[30]  E. A. Myers,et al.  Asynchronous demographic responses to Pleistocene climate change in Eastern Nearctic vertebrates. , 2016, Ecology letters.

[31]  Philipp W. Messer,et al.  Plumage Genes and Little Else Distinguish the Genomes of Hybridizing Warblers , 2016, Current Biology.

[32]  Scott V Edwards,et al.  Reticulation, divergence, and the phylogeography–phylogenetics continuum , 2016, Proceedings of the National Academy of Sciences.

[33]  B. Riddle Comparative phylogeography clarifies the complexity and problems of continental distribution that drove A. R. Wallace to favor islands , 2016, Proceedings of the National Academy of Sciences.

[34]  D. J. Hill The non-analogue nature of Pliocene temperature gradients , 2015 .

[35]  John D. McVay,et al.  Diversification of North American natricine snakes , 2015 .

[36]  F. Burbrink,et al.  Considering gene flow when using coalescent methods to delimit lineages of North American pitvipers of the genus Agkistrodon , 2015 .

[37]  Matthew Stephens,et al.  Visualizing spatial population structure with estimated effective migration surfaces , 2014, Nature Genetics.

[38]  Dong Xie,et al.  BEAST 2: A Software Platform for Bayesian Evolutionary Analysis , 2014, PLoS Comput. Biol..

[39]  G. Schuurman,et al.  Evidence of Hybridization between Common Gartersnakes (Thamnophis sirtalis) and Butler's Gartersnakes (Thamnophis butleri) in Wisconsin, USA , 2013 .

[40]  J. J. Emerson,et al.  Comparative Sex Chromosome Genomics in Snakes: Differentiation, Evolutionary Strata, and Lack of Global Dosage Compensation , 2013, PLoS biology.

[41]  Danielle Venton,et al.  Phylogeography , 2018, The International Encyclopedia of Biological Anthropology.

[42]  Matthew J. Smith,et al.  The Effects of Sampling Bias and Model Complexity on the Predictive Performance of MaxEnt Species Distribution Models , 2013, PloS one.

[43]  Peter L. Ralph,et al.  DISENTANGLING THE EFFECTS OF GEOGRAPHIC AND ECOLOGICAL ISOLATION ON GENETIC DIFFERENTIATION , 2013, Evolution; international journal of organic evolution.

[44]  Craig Moritz,et al.  Coalescent-based species delimitation in an integrative taxonomy. , 2012, Trends in ecology & evolution.

[45]  M. Kearney,et al.  Integrating phylogeography and physiology reveals divergence of thermal traits between central and peripheral lineages of tropical rainforest lizards , 2012, Philosophical Transactions of the Royal Society B: Biological Sciences.

[46]  R. A. Pyron,et al.  Out of Asia: Natricine snakes support the Cenozoic Beringian Dispersal Hypothesis , 2012 .

[47]  H. Hoekstra,et al.  Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species , 2012, PloS one.

[48]  Yuchi Zheng,et al.  Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: a case study of salamanders (order Caudata). , 2011, Molecular biology and evolution.

[49]  M. Pfrender,et al.  Genetic architecture of a feeding adaptation: garter snake (Thamnophis) resistance to tetrodotoxin bearing prey , 2010, Proceedings of the Royal Society B: Biological Sciences.

[50]  Remco R. Bouckaert,et al.  DensiTree: making sense of sets of phylogenetic trees , 2010, Bioinform..

[51]  David Bryant,et al.  Next-generation sequencing reveals phylogeographic structure and a species tree for recent bird divergences. , 2009, Molecular phylogenetics and evolution.

[52]  A. Leaché,et al.  Quantifying ecological, morphological, and genetic variation to delimit species in the coast horned lizard species complex (Phrynosoma) , 2009, Proceedings of the National Academy of Sciences.

[53]  Thibaut Jombart,et al.  adegenet: a R package for the multivariate analysis of genetic markers , 2008, Bioinform..

[54]  Paul R. Martin,et al.  Impacts of climate warming on terrestrial ectotherms across latitude , 2008, Proceedings of the National Academy of Sciences.

[55]  F. Burbrink,et al.  Phylogeography across a continent: the evolutionary and demographic history of the North American racer (Serpentes: Colubridae: Coluber constrictor). , 2008, Molecular phylogenetics and evolution.

[56]  M. Siddall,et al.  Phylogeography of Diadophis punctatus: extensive lineage diversity and repeated patterns of historical demography in a trans-continental snake. , 2008, Molecular phylogenetics and evolution.

[57]  L. Gustafsson,et al.  Sex Chromosome-Linked Species Recognition and Evolution of Reproductive Isolation in Flycatchers , 2007, Science.

[58]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[59]  W. Jetz,et al.  Environmental and historical constraints on global patterns of amphibian richness , 2007, Proceedings of the Royal Society B: Biological Sciences.

[60]  G. Burghardt,et al.  Post-glacial recolonization of the Great Lakes region by the common gartersnake (Thamnophis sirtalis) inferred from mtDNA sequences. , 2007, Molecular phylogenetics and evolution.

[61]  J. Overpeck,et al.  Simulating Arctic Climate Warmth and Icefield Retreat in the Last Interglaciation , 2006, Science.

[62]  Robert P. Anderson,et al.  Maximum entropy modeling of species geographic distributions , 2006 .

[63]  N. Swenson,et al.  Clustering of Contact Zones, Hybrid Zones, and Phylogeographic Breaks in North America , 2005, The American Naturalist.

[64]  P. Boag,et al.  Discordant temporal and geographic patterns in maternal lineages of eastern North American frogs, Rana catesbeiana (Ranidae) and Pseudacris crucifer (Hylidae). , 2004, Molecular phylogenetics and evolution.

[65]  John Wakeley,et al.  Estimating Divergence Times from Molecular Data on Phylogenetic and Population Genetic Timescales , 2002 .

[66]  R. Zink Methods in Comparative Phylogeography, and Their Application to Studying Evolution in the North American Aridlands1 , 2002, Integrative and comparative biology.

[67]  R. Durrett Probability Models for DNA Sequence Evolution , 2002 .

[68]  R. Lawson,et al.  Phylogenetic relationships of North American garter snakes (Thamnophis) based on four mitochondrial genes: how much DNA sequence is enough? , 2002, Molecular phylogenetics and evolution.

[69]  S. J. Arnold,et al.  Molecular systematics and evolution of Regina and the thamnophiine snakes. , 2001, Molecular phylogenetics and evolution.

[70]  F. Burbrink,et al.  MITOCHONDRIAL DNA PHYLOGEOGRAPHY OF THE POLYTYPIC NORTH AMERICAN RAT SNAKE (ELAPHE OBSOLETA): A CRITIQUE OF THE SUBSPECIES CONCEPT , 2000, Evolution; international journal of organic evolution.

[71]  G. Hewitt Some genetic consequences of ice ages, and their role in divergence and speciation , 1996 .

[72]  J. Falls,et al.  Thermal biology of the common garter snake Thamnophis sirtalis (L.) , 1979, Oecologia.

[73]  H. Akaike Maximum likelihood identification of Gaussian autoregressive moving average models , 1973 .

[74]  J. Bull,et al.  Karyotypic Studies of Thirty-Eight Species of North American Snakes , 1972 .

[75]  T. Schoener The Anolis Lizards of Bimini: Resource Partitioning in a Complex Fauna , 1968 .

[76]  E. Brodie Investigations on the Skin Toxin of the Adult Rough-Skinned Newt, Taricha granulosa , 1968 .

[77]  S. Wright,et al.  Isolation by Distance. , 1943, Genetics.

[78]  OUP accepted manuscript , 2021, Systematic Biology.

[79]  OUP accepted manuscript , 2021, Genome Biology and Evolution.

[80]  D. Doak,et al.  Incorporating intraspecific variation into species distribution models improves distribution predictions, but cannot predict species traits for a wide‐spread plant species , 2019, Ecography.

[81]  F. Burbrink,et al.  Ecological divergence in the yellow-bellied kingsnake (Lampropeltis calligaster) at two North American biodiversity hotspots. , 2017, Molecular phylogenetics and evolution.

[82]  A. Qvarnström,et al.  Speciation through evolution of sex-linked genes , 2009, Heredity.

[83]  A. Leaché,et al.  Molecular systematics of the Eastern Fence Lizard (Sceloporus undulatus): a comparison of Parsimony, Likelihood, and Bayesian approaches. , 2002, Systematic biology.

[84]  D. Soltis,et al.  Comparative phylogeography of north-western North America : a synthesis , 2001 .

[85]  R. Mason Chemical ecology of the red-sided garter snake, Thamnophis sirtalis parietalis. , 1993, Brain, behavior and evolution.