Testing hypotheses driving genetic structure in the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla

Habitat fragmentation has often been implicated in the decline of many species. For habitat specialists and/or sedentary species, loss of habitat can result in population isolation and lead to negative genetic effects. However, factors other than fragmentation can often be important and also need to be considered when assessing the genetic structure of a species. We genotyped individuals from 13 populations of the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla in Florida to test three alternative hypotheses regarding the effects that habitat fragmentation might have on genetic structure. A map of potential habitat developed from recent satellite imagery suggested that Brown‐headed Nuthatch populations in southern Florida occupied smaller and more isolated habitat patches (i.e. were more fragmented) than populations in northern Florida. We also genotyped individuals from a small, isolated Brown‐headed Nuthatch population on Grand Bahama Island. We found that populations associated with more fragmented habitat in southern Florida had lower allelic richness than populations in northern Florida (P = 0.02), although there were no differences in heterozygosity. Although pairwise estimates of F_ST were low overall, values among southern populations were generally higher than northern populations. Population assignment tests identified K = 3 clusters corresponding to a northern cluster, a southern cluster and a unique population in southeast Florida; using sampling localities as prior information revealed K = 7 clusters, with greater structure only among southern Florida populations. The Bahamas population showed moderate to high differentiation compared with Florida populations. Overall, our results suggest that fragmentation could affect gene flow in Brown‐headed Nuthatch populations and is likely to become more pronounced over time.     

Genetic effects of habitat fragmentation and population isolation on Etheostoma raneyi (Percidae)

The use of genetic methods to quantify the effects of anthropogenic habitat fragmentation on population structure has become increasingly common. However, in today’s highly fragmented habitats, researchers have sometimes concluded that populations are currently genetically isolated due to habitat fragmentation without testing the possibility that populations were genetically isolated before European settlement. Etheostoma raneyi is a benthic headwater fish restricted to river drainages in northern Mississippi, USA, that has a suite of adaptive traits that correlate with poor dispersal ability. Aquatic habitat within this area has been extensively modified, primarily by flood-control projects, and populations in headwater streams have possibly become genetically isolated from one another. We used microsatellite markers to quantify genetic structure as well as contemporary and historical gene flow across the range of the species. Results indicated that genetically distinct populations exist in each headwater stream analyzed, current gene flow rates are lower than historical rates, most genetic variation is partitioned among populations, and populations in the Yocona River drainage show lower levels of genetic diversity than populations in the Tallahatchie River drainage and other Etheostoma species. All populations have negative F_IS scores, of which roughly half are significant relative to Hardy–Weinberg expectations, perhaps due to small population sizes. We conclude that anthropogenic habitat alteration and fragmentation has had a profoundly negative impact on the species by isolating E. raneyi within headwater stream reaches. Further research is needed to inform conservation strategies, but populations in the Yocona River drainage are in dire need of management action. Carefully planned human-mediated dispersal and habitat restoration should be explored as management options across the range of the species.     

The effects of agriculture and the availability of edge habitat on populations of Helmeted Guineafowl Numida meleagris and on the diversity and composition of associated bird assemblages in KwaZulu-Natal province, South Africa

We investigated the effects of agriculture and the availability of edge habitat on populations of Helmeted Guineafowl Numida meleagris and associated avian diversity and species composition in woodland and grassland biomes in the Midlands of KwaZulu-Natal province, South Africa. Study sites within woodland biome had greater species diversity than those in grassland, whereas adjacent, high-quality, protected habitat in grassland sites, enhanced diversity within this biome. Both guineafowl populations and overall avian diversity declined with increasingly intensive agriculture and disappearance of edge habitat and associated, optimally fragmented habitat mosaic. Furthermore, traditional agriculture in the form of contouring in a pesticide-free environment, resulted in extensive edge habitat that appeared to provide additional food and cover for birds. This, in turn, caused an increase in overall bird diversity, and in guineafowl populations in particular. The widespread decline in Helmeted Guineafowl populations in the Midlands that started in the 1980s, and possibly the decline in species associated with this variegated landscape, was therefore caused by the loss of the habitat mosaic to intensive, modern, monoculture, crop agriculture. Maintaining species diversity and healthy guineafowl populations within these habitats requires the persistence or re-creation of a habitat mosaic and the resulting edge habitat on a landscape scale.     

Does fragmentation of Urtica habitats affect phytophagous and predatory insects differentially?

Effects of habitat fragmentation on the insect community of stinging nettle (Urtica dioica L.) were studied, using 32 natural nettle patches of different area and degree of isolation in an agricultural landscape. Habitat fragmentation reduced the species richness of Heteroptera, Auchenorrhyncha, and Coleoptera, and the abundance of populations. Habitat isolation and area reduction did not affect all insect species equally. Monophagous herbivores had a higher probability of absence from small patches than all (monophagous and polyphagous) herbivore species, and the percentage of monophagous herbivores increased with habitat area. Abundance and population variability of species were negatively correlated and could both be used as a predictor of the percentage of occupied habitats. Species richness of herbivores correlated (positively) with habitat area, while species richness of predators correlated (negatively) with habitat isolation. In logistic regressions, the probability of absence of monophagous herbivores from habitat patches could only be explained by habitat area (in 4 out of 10 species) and predator absence probability only by habitat isolation (in 3 out of 14 species). Presumably because of the instability of higher-trophic-level populations and dispersal limitation, predators were more affected by habitat isolation than herbivores, while they did not differ from herbivore populations with respect to abundance or variability. Thus increasing habitat connectivity in the agricultural landscape should primarily promote predator populations. Received: 4 February 1998 / Accepted: 4 May 1998     

Habitat use and ecological specialization within lake Daphnia populations

Many species of planktonic cladocerans display substantial variation in habitat use (mean depth and diel vertical migration), both among and within populations. We examined whether clonal segregation and specialization contributes to such behavioral variation within several lake populations of the cladoceran, Daphnia pulicaria. Electrophoretic and quantitative genetic analysis of clonal lines isolated from different depths at night revealed that clonal habitat specialization was common. Clones that utilized shallow water at night were genetically smaller at maturity and lower fecundity under standard laboratory conditions than the deep-water clones. The magnitude of this clonal habitat specialization varied among lakes: populations displaying broad use of depth habitats contained greater genetic variance than populations with more constrained habitat use. These results are consistent with known differences in selective factors in different depth habitats and suggest that substantial clonal specialization can occur within single populations. Since previous work has discovered a heritable basis to habitat selection in several Daphnia species, including D. pulicaria in our study lakes, it is likely that clonal/depth specialization is an important factor affecting the trophic ecology of Daphnia. Received: 18 April 1996 / Accepted: 25 September 1996     

Genetic structure,diversity and distribution of a threatened lizard affected by widespread habitat fragmentation

Understanding the demographic consequences of habitat loss on populations is essential for the conservation of threatened species. The threatened swamp skink (Lissolepis coventryi) is restricted to fragmented wetland habitats in Victoria and southeast South Australia. It has experienced significant habitat loss in the last 150 years, particularly around the Melbourne metropolitan area, where several small and isolated populations remain. Using mtDNA and nuDNA SNPs, we examined distribution patterns and population structure to infer evolutionary history and genetic distinctiveness of populations throughout the species’ range. For populations in the Melbourne metropolitan area, we examined genetic diversity. We found the species to be highly divergent, separating into two distinct lineages to the east and west of Melbourne, likely due to geological and climate influences causing isolation of populations. Species’ detectability was low, particularly in the far east despite relatively intact habitat and presumed higher abundance. Melbourne populations showed signs of limited genetic diversity. We suggest that translocations to promote gene diversity amongst these populations, together with habitat restoration and protection, present an important management strategy for L. coventryi.

     

Forest fragmentation and plant reproductive success: a case study in four perennial herbs

We assessed the effects of habitat fragmentation on reproductive success in natural populations of four forest herbs with differing life-history traits and whose distribution patterns appeared to be negatively affected by decreased habitat size and/or increased isolation: Carex sylvatica, Galium odoratum, Sanicula europaea and Veronica montana. Our aims were to test (1) whether habitat size and isolation are positively correlated with population size and isolation, respectively, (2) whether plant reproductive success, a major component of plant fitness, is reduced in small and/or isolated populations when also accounting for differences in habitat quality (edaphic conditions, light intensity) and the effects of plant size, and (3) whether species with different life histories are affected differently. There were significant positive relationships between habitat and population size and between habitat and population isolation in some, but not all of the species. We mostly found no negative effects of small population size or isolation on reproduction. However, reproductive success was reduced in small populations of Sanicula, and this effect was independent of differences in plant size and environmental conditions. The reduced fecundity in small populations may be a consequence of the Allee-effect, a possible mechanism being pollen limitation. Furthermore, the proportion of flowering ramets was reduced in small and isolated populations of Galium, which may have been caused by changes in population structure. Lastly, we found some evidence for largely outcrossing, non-clonal species to be more sensitive to reductions in population size, at least in terms of their reproductive success.     

GENETIC VARIATION IN CACTOPHILIC DROSOPHILA FOR OVIPOSITION ON NATURAL YEAST SUBSTRATES

Theory predicts that environmental heterogeneity in space or in time can maintain genetic polymorphism. Stable polymorphisms are expected to be more readily maintained if there are genotype specific habitat preferences. Genotype specific preferences for oviposition sites in Drosophila could be a major factor promoting habitat selection, and thus the maintenance of genetic variation. This hypothesis is being tested using the cactophilic species, D. buzzatii and D. aldrichi, where available evidence indicates a potential for such habitat selection, the habitats (oviposition sites) being yeast species found in the natural environment of these flies (cactus rots). Genetic variation for oviposition preferences was tested using isofemale lines—for D. buzzatii, a total of 60 lines from seven localities widely distributed through the species range in Australia, and for D. aldrichi, 21 lines from three of these localities. Females were given a choice of five yeast species as oviposition sites. Genetic variation for oviposition preferences on these natural substrates was demonstrated. There was significant variation among isofemale lines within populations in their patterns of preferences for oviposition on the five yeast species. However, analyses of preferences for each yeast species separately showed that the genetic variation for preferences relates to only three of the five species. Heritabilities of individual female preferences for these three species were low, ranging up to 9%. Little geographic differentiation was apparent among populations, most likely due to similar selection regimes within each population. Within populations, this kind of habitat selection could act to maintain polymorphisms, both at loci determining the habitat preferences and at other loci in linkage disequilibrium with them.     

Population genetic structure of three pond-inhabiting Daphnia species on a regional scale (Flanders, Belgium)

1. Only a few studies have compared patterns of genetic variation among populations of different Daphnia species on a regional scale. The present study addresses this gap and examines the relationship between diversity as revealed by allozyme variation and habitat size for populations of Daphnia pulex, D. obtusa and D. curvirostris in Flanders (Belgium). In addition, we examined whether patterns of isolation‐by‐distance could be observed in each of these three Daphnia species. 2. The relationship between genetic diversity and habitat size varied among Daphnia species that occur in the same region. In D. pulex and D. obtusa populations, a positive relationship between local genetic diversity and habitat size was found, whereas the relationship was negative in D. curvirostris populations. 3. Regional genetic diversity was lower than expected from patterns of local genetic diversity in D. pulex and D. obtusa populations in Flanders. This suggests that the subdivision of local Daphnia populations in a region did not obviously increase genetic diversity. 4. Genetic differentiation among populations of these three species in Flanders was moderate and comparable with values observed in other Daphnia species. Patterns of isolation‐by‐distance could be observed, but the scatter was high (D. pulex) or the slope was very low (D. obtusa).     

Contrasting genetic responses to population fragmentation in a coevolving fig and fig wasp across a mainland–island archipelago

Interacting species of pollinator–host systems, especially the obligate ones, are sensitive to habitat fragmentation, due to the nature of mutual dependence. Comparative studies of genetic structure can provide insights into how habitat fragmentation contributes to patterns of genetic divergence among populations of the interacting species. In this study, we used microsatellites to analyse genetic variation in Chinese populations of a typical mutualistic system – Ficus pumila and its obligate pollinator Wiebesia sp. 1 – in a naturally fragmented landscape. The plants and wasps showed discordant patterns of genetic variation and geographical divergence. There was no significant positive relationship in genetic diversity between the two species. Significant isolation‐by‐distance (IBD) patterns occurred across the populations of F. pumila and Wiebesia sp. 1 as whole, and IBD also occurred among island populations of the wasps, but not the plants. However, there was no significant positive relationship in genetic differentiation between them. The pollinator populations had significantly lower genetic variation in small habitat patches than in larger patches, and three island pollinator populations showed evidence of a recent bottleneck event. No effects of patch size or genetic bottlenecks were evident in the plant populations. Collectively, the results indicate that, in more fragmented habitats, the pollinators, but not the plants, have experienced reduced genetic variation. The contrasting patterns have multiple potential causes, including differences in longevity and hence number of generations experiencing fragmentation; different dispersal patterns, with the host's genes dispersed as seeds as well as a result of pollen dispersal via the pollinator; asymmetrical responses to fluctuations in partner populations; and co‐existence of a rare second pollinating wasp on some islands. These results indicate that strongly interdependent species may respond in markedly different ways to habitat fragmentation.     

The role of demography,intra‐species variation,and species distribution models in species' projections under climate change

Organisms are projected to shift their distribution ranges under climate change. The typical way to assess range shifts is by species distribution models (SDMs), which predict species’ responses to climate based solely on projected climatic suitability. However, life history traits can impact species’ responses to shifting habitat suitability. Additionally, it remains unclear if differences in vital rates across populations within a species can offset or exacerbate the effects of predicted changes in climatic suitability on population viability. In order to obtain a fuller understanding of the response of one species to projected climatic changes, we coupled demographic processes with predicted changes in suitable habitat for the monocarpic thistle Carlina vulgaris across northern Europe. We first developed a life history model with species‐specific average fecundity and survival rates and linked it to a SDM that predicted changes in habitat suitability through time with changes in climatic variables. We then varied the demographic parameters based upon observed vital rates of local populations from a translocation experiment. Despite the fact that the SDM alone predicted C. vulgaris to be a climate ‘winner’ overall, coupling the model with changes in demography and small‐scale habitat suitability resulted in a matrix of stable, declining, and increasing patches. For populations predicted to experience declines or increases in abundance due to changes in habitat suitability, altered fecundity and survival rates can reverse projected population trends.     

Genetic variability and the effect of habitat fragmentation in spotted suslik Spermophilus suslicus populations from two different regions

Endangered species worldwide exist in remnant populations, often within fragmented landscapes. Although assessment of genetic diversity in fragmented habitats is very important for conservation purposes, it is usually impossible to evaluate the amount of diversity that has actually been lost. Here, we compared population structure and levels of genetic diversity within populations of spotted suslik Spermophilus suslicus, inhabiting two different parts of the species range characterized by different levels of habitat connectivity. We used microsatellites to analyze 10 critically endangered populations located at the western part of the range, where suslik habitat have been severely devastated due to agriculture industrialization. Their genetic composition was compared with four populations from the eastern part of the range where the species still occupies habitat with reasonable levels of connectivity. In the western region, we detected extreme population structure (F _ST = 0.20) and levels of genetic diversity (Allelic richness ranged from 1.45 to 3.07) characteristic for highly endangered populations. Alternatively, in the eastern region we found significantly higher allelic richness (from 5.09 to 5.81) and insignificant population structure (F _ST = 0.03). As we identified a strong correlation between genetic and geographic distance and a lack of private alleles in the western region, we conclude that extreme population structure and lower genetic diversity is due to recent habitat loss. Results from this study provide guidelines for conservation and management of this highly endangered species.     

Strongly diverging population genetic patterns of three skipper species: the role of habitat fragmentation and dispersal ability

The fragmentation of landscapes has an important impact on the conservation of biodiversity, and the genetic diversity is an important factor for a populations viability, influenced by the landscape structure. However, different species with differing ecological demands react rather different on the same landscape pattern. To address this feature, we studied three skipper species with differing habitat requirements (Lulworth Skipper Thymelicus acteon: a habitat specialist with low dispersal ability, Small Skipper Thymelicus sylvestris: a habitat generalist with low dispersal ability, Essex Skipper Thymelicus lineola: a habitat generalist with higher dispersal ability). We analysed 18 allozyme loci for 1,063 individuals in our western German study region with adjoining areas in Luxembourg and north-eastern France. The genetic diversity of all three species were intermediate in comparison with other butterfly species. The F _ST was relatively high for T. acteon (5.1%), low for T. sylvestris (1.6%) and not significant for T. lineola. Isolation by distance analyses revealed a significant correlation for T. sylvestris explaining 20.3% of its differentiation, but no such structure was found for the two other species. Most likely, the high dispersal ability of T. lineola in comparison with T. sylvestris leads to a more or less panmictic structure and hence impedes isolation by distance. On the other hand, the isolation of the populations of T. acteon seems to be so strict that the populations develop independently. Although no general genetic impoverishing was observed for the endangered T. acteon, small populations had significantly lower genetic diversities than big populations, and therefore the high degree of isolation among populations might threaten its local and regional survival.     

Adaptive radiation versus intraspecific differentiation: morphological variation in Caribbean Anolis lizards

Lizards in the genus Anolis have experienced adaptive radiation in the Greater Antilles, producing a suite of species morphologically adapted to use different parts of the environment. In the Lesser Antilles, adaptive radiation has not occurred, but on some islands, interpopulational variation is high and represents adaptation to different habitats. We compared the extent of morphological differentiation among Greater Antillean habitat specialists with that exhibited among populations of two species, Anolis marmoratus and A. oculatus, from the Lesser Antillean islands of Guadeloupe and Dominica. Although extensive, intraspecific divergence in the Lesser Antilles is substantially less in magnitude than the differences among habitat specialists in the Greater Antilles. All populations of A. marmoratus are most similar to Greater Antillean trunk‐crown habitat specialists, but populations of A. oculatus differ in their affinities: some are similar to trunk‐crown anoles, but others are more similar to trunk‐ground habitat specialists.     

Impact of urban fragmentation on the genetic structure of the eastern red-backed salamander

Urban development is a major cause of habitat loss and fragmentation. Few studies, however, have dealt with fragmentation in an urban landscape. In this paper, we examine the genetic structure of isolated populations of the eastern red-backed salamander (Plethodon cinereus) in a metropolitan area. We sampled four populations located on a mountain in the heart of Montréal (Québec, Canada), which presents a mosaic of forested patches isolated by roads, graveyards and buildings. We assessed the genetic structure of these populations using microsatellite loci and compared it to the genetic structure of four populations located in a continuous habitat in southern Québec. Our results indicate that allelic richness and heterozygosity are lower in the urban populations. Exact differentiation tests and pairwise F _ST also show that the populations found in the fragmented habitat are genetically differentiated, whereas populations located in the continuous habitat are genetically homogeneous. These results raise conservation concerns for these populations as well as for rare or threatened species inhabiting urban landscapes.     

WEAK CROSSABILITY BARRIER BUT STRONG JUVENILE SELECTION SUPPORTS ECOLOGICAL SPECIATION OF THE HYBRID PINE PINUS DENSATA ON THE TIBETAN PLATEAU

Determining how a new hybrid lineage can achieve reproductive isolation is a key to understanding the process and mechanisms of homoploid hybrid speciation. Here, we evaluated the degree and nature of reproductive isolation between the ecologically successful hybrid species Pinus densata and its parental species P. tabuliformis and P. yunnanensis. We performed interspecific crosses among the three species to assess their crossability. We then conducted reciprocal transplantation experiments to evaluate their fitness differentiation, and to examine how natural populations representing different directions of introgression differ in adaptation. The crossing experiments revealed weak genetic barriers among the species. The transplantation trials showed manifest evidence of local adaptation as the three species all performed best in their native habitats. Pinus densata populations from the western edge of its distribution have evolved a strong local adaptation to the specific habitat in that range; populations representing different directions of introgressants with the two parental species all showed fitness disadvantages in this P. densata habitat. These observations illustrate that premating isolation through selection against immigrants from other habitat types or postzygotic isolation through selection against backcrosses between the three species is strong. Thus, ecological selection in combination with endogenous components and geographic isolation has likely played a significant role in the speciation of P. densata.     

Genetic variation in time and space: the use of herbarium specimens to reconstruct patterns of genetic variation in the endangered orchid Anacamptis palustris

Habitat alteration, fragmentation and destruction as a consequence of human impact are a global phenomenon. Here we document changes in genetic variation in the marsh orchid Anacamptis palustris as a consequence of such habitat changes. We examined historical specimens that were collected during the nineteenth and early twentieth centuries, prior to the most recent massive habitat changes affecting this species. Sequences of a hypervariable region in the plastid DNA, located in the tRNA^LEU intron, from herbarium vouchers were compared with those from a near-exhaustive survey of the extant A. palustris populations on the Italian peninsula. It was found that private haplotypes and alleles found in small, extant populations were once widespread and more common in historic populations and that alleles, once present in historic populations, have gone extinct. In addition, genetic differentiation among populations has increased over time and haplotype frequencies significantly differ among historic and extant populations. These results document that human induced habitat changes have reduced genetic diversity and increased the importance of random genetic drift in this species. It is concluded that the analysis of herbarium specimens may provide important insights into changes of genetic diversity over time and may be critical for correct inference of the evolutionary history of rare and endangered species.     

Habitat preferences and habitat restoration options for small‐bodied and juvenile fish species in the northern Murray–Darling Basin

Degradation of instream habitats in the northern Murray–Darling Basin has occurred through numerous stressors, including siltation, clearing of bankside vegetation, intrusion of livestock and impacts of pest species. A better understanding of habitat preferences of native fish species could help guide future instream habitat restoration actions. The habitat choices of seven native fish species, juvenile Murray Cod (Maccullochella peelii), juvenile Golden Perch (Macquaria ambigua ambigua), juvenile Silver Perch (Bidyanus bidyanus), adult Murray–Darling Rainbowfish (Melanotaenia fluviatilis), adult Olive Perchlet (Ambassis agassizii), adult Un‐specked Hardyhead (Craterocephalus stercusmuscarum fulvus) and adult carp gudgeons (Hypseleotris spp.) were tested in preference troughs to help inform potential habitat restoration actions in the Condamine catchment. Each species was given a choice between pair combinations of open sandy habitat, submerged macrophytes, emergent plants and rocky rubble. Habitat preferences varied between species. Murray Cod, Golden Perch, carp gudgeons and Olive Perchlets preferred structure over open sandy habitat, whilst juvenile Silver Perch, Un‐specked Hardyhead and Murray–Darling Rainbowfish did not avoid open sandy habitats. Juvenile Murray Cod preferred rocky rubble habitat over all other habitat choices. Use of complex rock piles to provide nursery habitat for Murray Cod populations is a potential restoration option. Introduction of rock could also benefit Golden Perch and carp gudgeons. Use of emergent plants, submerged macrophytes and rocky rubble for habitat restoration all appear to have merit for one or more species of small‐bodied fishes or juvenile stages of larger sized fishes. Rocky rubble or floating attached macrophytes could be viable restoration options in areas too turbid to establish submerged macrophytes. These habitat interventions would complement existing actions such as re‐snagging and provision of fish passage to assist with sustainable management of native fish populations.     

Survey of Helminth Parasites in Populations of <Emphasis Type="Italic">Alouatta palliata mexicana</Emphasis> and <Emphasis Type="Italic">A. pigra</Emphasis> in Continuous and in Fragmented Habitat in Southern Mexico

The 2 howler species that occur in southern Mexico, Alouatta palliata mexicana and Alouatta pigra are endangered, mainly as a result of habitat loss and fragmentation from human activity. Little is known about the gastrointestinal parasite communities affecting their populations, and lack of baseline information for populations of howler species in continuous forest habitats, makes evaluations of gastrointestinal parasite prevalence in populations in fragmented landscapes difficult. We report the results of a one-time broad survey of gastrointestinal parasites in fecal samples of individuals from several demographically stable populations of Alouatta palliata mexicana and A. pigra existing in continuous and/or protected forests. We further report similar data for populations of both species in human-fragmented landscapes. We detected 6 parasites for each howler monkey species, but only 3 of them (Trematode I, Controrchis biliophilus, Trypanoxyuris sp.) were common to both species. While parasitic prevalence in populations of both howler species was, in general, higher in the fragmented habitat than in continuous and/or protected forests. The difference is only marginally significant in Alouatta pigra. Some parasites (Coccidia and Strongylid) only appeared in populations in fragmented landscapes. Preliminary data suggest that adult males tended to have higher parasite prevalence values than those of adult females in both howler species. Parasite prevalence is associated to average group size, but not to population density in Alouatta pigra.     

Habitat selection by grassland snakes at northern range limits: Implications for conservation

Animal populations at northern range limits may use habitat differently from those at range cores, requiring distinct conservation plans. Snakes are ectotherms that often have very specific requirements, but few studies have focused on the effect of northern latitudes on habitat selection by grassland snakes. We studied movement and habitat selection of 2 sympatric snake species at their northern range limits on the North American Great Plains: the eastern yellow-bellied racer (Coluber constrictor flaviventris, hereafter racer), a Threatened species in Canada, and the bullsnake (Pituophis catenifer sayi), which is listed as Data Deficient. Both of these species are potentially vulnerable to extinction in Canada because of habitat loss. Snakes from our study populations traveled up to 10-times farther from winter dens and occupied home ranges 3–104 times larger than populations further south. Both snake species moved from winter dens in the slopes of a major river valley to habitat in adjacent lowlands, including riparian zones (racers) and hilly areas with native grass species (bullsnakes). Multivariate modeling revealed that proximity to retreat sites was a significant predictor of snake site use for both species. Considering the need for winter dens and summering areas, our data suggest that snakes in northern latitudes should ideally have much larger protected areas compared to snakes near the core of their range. An alternative strategy is to conserve corridors linking wintering dens and summer habitats. Retreat sites such as burrows and shrubs are critical components of local habitat and should be included in conservation plans. © 2011 The Wildlife Society.