Temperature-induced multi-species cohort effects in sympatric snakes.

In reptiles, reproductive maturity is often determined by size rather than age. Consequently, growth early in life may influence population dynamics through effects on generation time and survival to reproduction. Because reproductive phenology and pre- and post-natal growth are temperature dependent, environmental conditions may induce multi-species cohort effects on body size in sympatric reptiles. I present evidence of this using 10 years of neonatal size data for three sympatric viviparous snakes, Dekay's Brown snakes (Storeria dekayi), Red-bellied Snakes (S. occipitomaculata), and Common Garter snakes (Thamnophis sirtalis). End-of-season neonatal size varied in parallel across species such that snout-vent length was 36%-61% greater and mass was 65%-223% greater in years when gestating females could achieve higher April-May (vs. June-July or August-September) operative temperatures. Thus, temperature had a larger impact during follicular enlargement and ovulation than during gestation or post-natal growth. Multi-species cohort effects like these may affect population dynamics and the magnitude of these effects may increase with climate change.

Effects of management outweigh effects of plant diversity on restored animal communities in tallgrass prairies.

A primary goal of ecological restoration is to increase biodiversity in degraded ecosystems. However, the success of restoration ecology is often assessed by measuring the response of a single functional group or trophic level to restoration, without considering how restoration affects multitrophic interactions that shape biodiversity. An ecosystem-wide approach to restoration is therefore necessary to understand whether animal responses to restoration, such as changes in biodiversity, are facilitated by changes in plant communities (plant-driven effects) or disturbance and succession resulting from restoration activities (management-driven effects). Furthermore, most restoration ecology studies focus on how restoration alters taxonomic diversity, while less attention is paid to the response of functional and phylogenetic diversity in restored ecosystems. Here, we compared the strength of plant-driven and management-driven effects of restoration on four animal communities (ground beetles, dung beetles, snakes, and small mammals) in a chronosequence of restored tallgrass prairie, where sites varied in management history (prescribed fire and bison reintroduction). Our analyses indicate that management-driven effects on animal communities were six-times stronger than effects mediated through changes in plant biodiversity. Additionally, we demonstrate that restoration can simultaneously have positive and negative effects on biodiversity through different pathways, which may help reconcile variation in restoration outcomes. Furthermore, animal taxonomic and phylogenetic diversity responded differently to restoration, suggesting that restoration plans might benefit from considering multiple dimensions of animal biodiversity. We conclude that metrics of plant diversity alone may not be adequate to assess the success of restoration in reassembling functional ecosystems.

Sunning themselves in heaps, knots, and snarls: The extraordinary abundance and demography of island watersnakes.

Snakes represent a sizable fraction of vertebrate biodiversity, but until recently, data on their demography have been sparse. Consequently, generalizations regarding patterns of variation are weak and the potential for population projections is limited. We address this information gap through an analysis of spatial and temporal variation in demography (population size, annual survival, and realized population growth) of the Lake Erie Watersnake, Nerodia sipedon insularum, and a review of snake survival more generally. Our study spans a period during which the Lake Erie Watersnake was listed as threatened under the U.S. Endangered Species Act, recovered, and was delisted. We collected capture-mark-recapture data at 14 study sites over 20 years, accruing 20,000 captures of 13,800 individually marked adults. Lake Erie Watersnakes achieve extraordinary abundance, averaging 520 adults per km of shoreline (ca. 260 adult per ha) at our study sites (range = 160-1,600 adults per km; ca. 80-800 adults per ha) and surpassing population recovery and postdelisting monitoring criteria. Annual survival averages 0.68 among adult females and 0.76 among adult males, varies among sites, and is positively correlated with body size among study sites. Temporal process variance in annual survival is low, averaging 0.0011 or less than 4% of total variance; thus, stochasticity in annual survival may be of minor significance to snake extinction risk. Estimates of realized population growth indicate that population size has been stable or increasing over the course of our study. More generally, snake annual survival overlaps broadly across continents, climate zones, families, subfamilies, reproductive modes, body size categories, maturation categories, and parity categories. Differences in survival in relation to size, parity, and maturation are in the directions predicted by life history theory but are of small magnitude with much variation around median values. Overall, annual survival appears to be quite plastic, varying with food availability, habitat quality, and other ecological variables.

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.

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake, the Eastern Massasauga (Sistrurus catenatus).

Genetic diversity is fundamental to maintaining the long-term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.

Size Matters: Individual Variation in Ectotherm Growth and Asymptotic Size.

Body size, and, by extension, growth has impacts on physiology, survival, attainment of sexual maturity, fecundity, generation time, and population dynamics, especially in ectotherm animals that often exhibit extensive growth following attainment of sexual maturity. Frequently, growth is analyzed at the population level, providing useful population mean growth parameters but ignoring individual variation that is also of ecological and evolutionary significance. Our long-term study of Lake Erie Watersnakes, Nerodia sipedon insularum, provides data sufficient for a detailed analysis of population and individual growth. We describe population mean growth separately for males and females based on size of known age individuals (847 captures of 769 males, 748 captures of 684 females) and annual growth increments of individuals of unknown age (1,152 males, 730 females). We characterize individual variation in asymptotic size based on repeated measurements of 69 males and 71 females that were each captured in five to nine different years. The most striking result of our analyses is that asymptotic size varies dramatically among individuals, ranging from 631-820 mm snout-vent length in males and from 835-1125 mm in females. Because female fecundity increases with increasing body size, we explore the impact of individual variation in asymptotic size on lifetime reproductive success using a range of realistic estimates of annual survival. When all females commence reproduction at the same age, lifetime reproductive success is greatest for females with greater asymptotic size regardless of annual survival. But when reproduction is delayed in females with greater asymptotic size, lifetime reproductive success is greatest for females with lower asymptotic size when annual survival is low. Possible causes of individual variation in asymptotic size, including individual- and cohort-specific variation in size at birth and early growth, warrant further investigation.

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.

COLOR-PATTERN VARIATION IN LAKE ERIE WATER SNAKES: THE ROLE OF GENE FLOW.

In an effort to clarify the evolutionary processes influencing color-pattern variation in Lake Erie island water snake (Nerodia sipedon) populations, rates of gene flow among island and mainland populations were estimated from patterns of allozymic variation detected using electrophoresis. Rates of gene flow were high with Nm, the number of migrants per generation, averaging 25.5 among island sites, 9.2 between the Ontario mainland and the islands, and 3.6 between the Ohio mainland and the islands. Based on estimates of current population size from mark-recapture work and of past population size extrapolated from the extent of shoreline habitat, values of m between island and mainland populations ranged from 0.0008-0.01. Synthesis of estimates of the rate of gene flow with information on inheritance of color pattern, the strength of natural selection, and population history supports the hypothesis that color-pattern variation in island populations results from a balance between gene flow and natural selection. However, depending on the mode of inheritance of color pattern, stochastic processes such as drift may have been important in the initial stages of differentiation between island and mainland populations.

COLOR-PATTERN VARIATION IN LAKE ERIE WATER SNAKES: PREDICTION AND MEASUREMENT OF NATURAL SELECTION.

A classic example of natural selection, that of color-pattern variation in Lake Erie island populations of water snakes, was reexamined to overcome shortcomings resulting from classification of snakes into discrete color-pattern categories and use of cross-sectional data. Four continuously varying color-pattern components (DB, the number of dorsal blotches; LB, the number of lateral blotches; ROWS, the height of lateral blotches measured in scale rows; and VEXT, the extent of ventral pigmentation) were analyzed. Patterns of natural selection were predicted from the relationship between color-pattern scores and independent measures of relative crypsis. Tests for natural selection were carried out using longitudinal data on neonate to juvenile-aged snakes and cross-sectional data on juvenile to adult-aged snakes. As predicted, the form of selection differed between younger and older age classes of snakes: selection resulted in a reduction in DB and LB among neonate and juvenile snakes but had little influence on color-pattern components in older snakes. The correspondence between observed patterns of natural selection and predictions based on the relationship between color-pattern scores and relative crypsis supports the hypothesis that differential predation by visual predators on younger age classes of snakes is the mechanism of selection. Gene flow from mainland populations or the initial lack of an allele necessary for reduced pattern may explain why selection has not resulted in greater differentiation between island and mainland populations.

COLOR PATTERN POLYMORPHISM IN THE LAKE ERIE WATER SNAKE, NERODIA SIPEDON INSULARUM.

Populations of the water snake, Nerodia sipedon, on islands in western Lake Erie are polymorphic for color pattern. These populations include banded, intermediate, and unbanded morphs while surrounding mainland populations consist solely of the banded morph. The hypothesis that this polymorphism is maintained by strong selection and migration pressures is widely accepted. Unbanded morphs are apparently more cryptic along island shorelines while banded morphs are more cryptic on the mainland. Migration of banded morphs from the mainland explains their persistence in island populations. Data collected in a capture-mark-recapture program on six islands provide no evidence of differential selection among morphs; morph frequencies do not differ among age classes, between once-captured and multiply-captured snakes, or between scarred and unscarred snakes. Furthermore, herring gulls, the most common snake predators in the island area, appear to detect banded and unbanded model snakes with equal ease. High site fidelity of water snakes and the distribution of morphs among islands suggest that migration from the mainland is not common. However, islands close to each other are similar in morph frequency, and water snakes have colonized islands elsewhere in the Great Lakes, indicating that some migration does occur. Recently, the frequency of banded morphs has increased in island populations while adult population sizes have declined. This increase in banded morphs is interpreted as reflecting an increased impact of migration from the mainland into these reduced populations. One scenario for the evolution and maintenance of this polymorphism is that selection was important in establishing unbanded morphs in island populations as they became isolated from the mainland. As populations declined to their present size, the impact of migration from the mainland increased and is now swamping the effect of selection. Further declines in island population size may result in fixation of the banded morph.