Extreme environmental perturbations offer opportunities to observe the effects of natural selection in wild populations. During the winter of 2013-2014, the southeastern United States endured an extreme cold event. We used thermal performance, transcriptomics, and genome scans to measure responses of lizard populations to storm-induced selection. We found significant increases in cold tolerance at the species' southern limit. Gene expression in southern survivors shifted toward patterns characteristic of northern populations. Comparing samples before and after the extreme winter, 14 genomic regions were differentiated in the surviving southern population; four also exhibited signatures of local adaptation across the latitudinal gradient and implicate genes involved in nervous system function. Together, our results suggest that extreme winter events can rapidly produce strong selection on natural populations at multiple biological levels that recapitulate geographic patterns of local adaptation.

No abstract

Anolis carolinensis is an emerging model species and the sole member of its genus native to the United States. Considerable morphological and physiological variation has been described in the species, and the recent sequencing of its genome makes it an attractive system for studies of genome variation. To inform future studies of molecular and phenotypic variation within A. carolinensis, a rigorous account of intraspecific population structure and relatedness is needed. Here, we present the most extensive phylogeographic study of this species to date. Phylogenetic analyses of mitochondrial DNA sequence data support the previous hypothesis of a western Cuban origin of the species. We found five well-supported, geographically distinct mitochondrial haplotype clades throughout the southeastern United States. Most Florida populations fall into one of three divergent clades, whereas the vast majority of populations outside Florida belong to a single, shallowly diverged clade. Genetic boundaries do not correspond to major rivers, but may reflect effects of Pleistocene glaciation events and the Appalachian Mountains on migration and expansion of the species. Phylogeographic signal should be examined using nuclear loci to complement these findings.

Understanding the mechanisms that produce variation in thermal performance is a key component to investigating climatic effects on evolution and adaptation. However, disentangling the effects of local adaptation and phenotypic plasticity in shaping patterns of geographic variation in natural populations can prove challenging. Additionally, the physiological mechanisms that cause organismal dysfunction at extreme temperatures are still largely under debate. Using the green anole, Anolis carolinensis, we integrate measures of cold tolerance (CT ), standard metabolic rate, heart size, blood lactate concentration and RNAseq data from liver tissue to investigate geographic variation in cold tolerance and its underlying mechanisms along a latitudinal cline. We found significant effects of thermal acclimation and latitude of origin on variation in cold tolerance. Increased cold tolerance correlates with decreased rates of oxygen consumption and blood lactate concentration (a proxy for oxygen limitation), suggesting elevated performance is associated with improved oxygen economy during cold exposure. Consistent with these results, co-expression modules associated with blood lactate concentration are enriched for functions associated with blood circulation, coagulation and clotting. Expression of these modules correlates with thermal acclimation and latitude of origin. Our findings support the oxygen and capacity-limited thermal tolerance hypothesis as a potential contributor to variation in reptilian cold tolerance. Moreover, differences in gene expression suggest regulation of the blood coagulation cascade may play an important role in reptilian cold tolerance and may be the target of natural selection in populations inhabiting colder environments.