Climatic and geographic predictors of life history variation in Eastern Massasauga (Sistrurus catenatus): A range-wide synthesis.

Elucidating how life history traits vary geographically is important to understanding variation in population dynamics. Because many aspects of ectotherm life history are climate-dependent, geographic variation in climate is expected to have a large impact on population dynamics through effects on annual survival, body size, growth rate, age at first reproduction, size-fecundity relationship, and reproductive frequency. The Eastern Massasauga (Sistrurus catenatus) is a small, imperiled North American rattlesnake with a distribution centered on the Great Lakes region, where lake effects strongly influence local conditions. To address Eastern Massasauga life history data gaps, we compiled data from 47 study sites representing 38 counties across the range. We used multimodel inference and general linear models with geographic coordinates and annual climate normals as explanatory variables to clarify patterns of variation in life history traits. We found strong evidence for geographic variation in six of nine life history variables. Adult female snout-vent length and neonate mass increased with increasing mean annual precipitation. Litter size decreased with increasing mean temperature, and the size-fecundity relationship and growth prior to first hibernation both increased with increasing latitude. The proportion of gravid females also increased with increasing latitude, but this relationship may be the result of geographically varying detection bias. Our results provide insights into ectotherm life history variation and fill critical data gaps, which will inform Eastern Massasauga conservation efforts by improving biological realism for models of population viability and climate change.

Distinctiveness in the face of gene flow: hybridization between specialist and generalist gartersnakes.

Patterns of divergence and polymorphism across hybrid zones can provide important clues as to their origin and maintenance. Unimodal hybrid zones or hybrid swarms are composed predominantly of recombinant individuals whose genomes are patchworks of alleles derived from each parental lineage. In contrast, bimodal hybrid zones contain few identifiable hybrids; most individuals fall within distinct genetic clusters. Distinguishing between hybrid swarms and bimodal hybrid zones can be important for taxonomic and conservation decisions regarding the status and value of hybrid populations. In addition, the causes of bimodality are important in understanding the generation and maintenance of biological diversity. For example, are distinct clusters mostly reproductively isolated and co-adapted gene complexes, or can distinctiveness be maintained by a few 'genomic islands' despite rampant gene flow across much of the genome? Here we focus on three patterns of distinctiveness in the face of gene flow between gartersnake taxa in the Great Lakes region of North America. Bimodality, the persistence of distinct clusters of genotypes, requires strong barriers to gene flow and supports recognition of distinct specialist (Thamnophis butleri) and generalist (Thamnophis radix) taxa. Concordance of DNA-based clusters with morphometrics supports the hypothesis that trophic morphology is a key component of divergence. Finally, disparity in the level of differentiation across molecular markers (amplified fragment length polymorphisms) indicates that distinctiveness is maintained by strong selection on a few traits despite high gene flow currently or in the recent past.

Post-glacial recolonization of the Great Lakes region by the common gartersnake (Thamnophis sirtalis) inferred from mtDNA sequences.

Pleistocene events played an important role in the differentiation of North American vertebrate populations. Michigan, in particular, and the Great Lakes region, in general, were greatly influenced by the last glaciation. While several hypotheses regarding the recolonization of this region have been advanced, none have been strongly supported. We generated 148 complete ND2 mitochondrial DNA (mtDNA) sequences from common gartersnake (Thamnophis sirtalis) populations throughout the Great Lakes region to evaluate phylogeographic patterns and population structure and to determine whether the distribution of haplotypic variants is related to the post-Pleistocene retreat of the Wisconsinan glacier. The common gartersnake was utilized, as it is believed to have been one of the primary vertebrate invaders of the Great Lakes region following the most recent period of glacial retreat and because it has been a model species for a variety of evolutionary, ecological, behavioral, and physiological studies. Several genetically distinct evolutionary lineages were supported by both genealogical and molecular population genetic analyses, although to different degrees. The geographic distribution of the majority of these lineages is interpreted as reflecting post-glacial recolonization dynamics during the late Pleistocene. These findings generally support previous hypotheses of range expansion in this region.