Digest: How the snowshoe hare got its brown coat: Convergent evolution or gene flow?*

The winter-brown phenotype of snowshoe hares in the Pacific Northwest was acquired through hybridization with black-tailed jackrabbits. Some snowshoe hares in more northern boreal populations exhibit the same phenotype, but how did they acquire it? Jones and colleagues show that the phenotype in the boreal populations is the outcome of convergent evolution, highlighting the importance of understanding the genetic basis of a trait in reconstructing its evolution.     

The hidden history of the snowshoe hare,Lepus americanus: extensive mitochondrial DNA introgression inferred from multilocus genetic variation

Hybridization drives the evolutionary trajectory of many species or local populations, and assessing the geographic extent and genetic impact of interspecific gene flow may provide invaluable clues to understand population divergence or the adaptive relevance of admixture. In North America, hares (Lepus spp.) are key species for ecosystem dynamics and their evolutionary history may have been affected by hybridization. Here we reconstructed the speciation history of the three most widespread hares in North America – the snowshoe hare (Lepus americanus), the white‐tailed jackrabbit (L. townsendii) and the black‐tailed jackrabbit (L. californicus) – by analysing sequence variation at eight nuclear markers and one mitochondrial DNA (mtDNA) locus (6240 bp; 94 specimens). A multilocus–multispecies coalescent‐based phylogeny suggests that L. americanus diverged ~2.7 Ma and that L. californicus and L. townsendii split more recently (~1.2 Ma). Within L. americanus, a deep history of cryptic divergence (~2.0 Ma) was inferred, which coincides with major speciation events in other North American species. While the isolation‐with‐migration model suggested that nuclear gene flow was generally rare or absent among species or major genetic groups, coalescent simulations of mtDNA divergence revealed historical mtDNA introgression from L. californicus into the Pacific Northwest populations of L. americanus. This finding marks a history of past reticulation between these species, which may have affected other parts of the genome and influence the adaptive potential of hares during climate change.     

Genetic structure among greater white‐fronted goose populations of the Pacific Flyway

An understanding of the genetic structure of populations in the wild is essential for long‐term conservation and stewardship in the face of environmental change. Knowledge of the present‐day distribution of genetic lineages (phylogeography) of a species is especially important for organisms that are exploited or utilize habitats that may be jeopardized by human intervention, including climate change. Here, we describe mitochondrial (mtDNA) and nuclear genetic (microsatellite) diversity among three populations of a migratory bird, the greater white‐fronted goose (Anser albifrons), which breeds discontinuously in western and southwestern Alaska and winters in the Pacific Flyway of North America. Significant genetic structure was evident at both marker types. All three populations were differentiated for mtDNA, whereas microsatellite analysis only differentiated geese from the Cook Inlet Basin. In sexual reproducing species, nonrandom mate selection, when occurring in concert with fine‐scale resource partitioning, can lead to phenotypic and genetic divergence as we observed in our study. If mate selection does not occur at the time of reproduction, which is not uncommon in long‐lived organisms, then mechanisms influencing the true availability of potential mates may be obscured, and the degree of genetic and phenotypic diversity may appear incongruous with presumed patterns of gene flow. Previous investigations revealed population‐specific behavioral, temporal, and spatial mechanisms that likely influence the amount of gene flow measured among greater white‐fronted goose populations. The degree of observed genetic structuring aligns well with our current understanding of population differences pertaining to seasonal movements, social structure, pairing behavior, and resource partitioning.     

Population genetic structure of the cyclic snowshoe hare (Lepus americanus) in southwestern Yukon,Canada

Spatial population structure has important ecological and evolutionary consequences. Little is known about the population structure of snowshoe hares (Lepus americanus), despite their ecological importance in North American boreal forests. We used seven variable microsatellite DNA loci to determine the spatial genetic structure of snowshoe hares near Kluane Lake, Yukon during a cyclic population peak. We sampled 317 hares at 12 sites separated by distances ranging from 3 to 140 km, and used 46 additional samples from Alaska and Montana. The level of genetic variation was high (13.4 alleles/locus, 0.67 expected heterozygosity) and the distribution of alleles and genotypes was not homogeneous across the sites. The degree of differentiation was low among Yukon sites (FST = 0.015) and between Yukon and Alaska (FST = 0.012), but the Montana site was highly differentiated (FST = 0.20). A weak pattern of isolation by distance was found over the Yukon study area, with an indication that local genetic drift may be important in shaping the regional genetic structure. Landscape barriers expected to influence gene flow did not consistently affect genetic structure, although there was evidence for a partial barrier effect of Kluane Lake. The high level of inferred gene flow confirms that snowshoe hare dispersal is widespread, successful and equal between the sexes. A stepping-stone model of gene flow, potentially influenced by the synchronous density cycle, appears to best explain the observed genetic structure. Our results suggest that despite their dramatic fluctuations in density, snowshoe hares in the northern boreal forest have a large evolutionary effective population size and are not strongly subdivided by either physical or social barriers to gene flow.     

Assessment of Canada Lynx Research and Conservation Needs in the Southern Range: Another Kick at the Cat

Abstract The ecology of Canada lynx (Lynx canadensis) and their main prey, snowshoe hares (Lepus americanus), is poorly understood in southern Canada and the contiguous United States compared to the boreal forest of Canada and Alaska, USA, where both species are well studied. However, given recent listing of lynx under the Endangered Species Act, accurate understanding of lynx and snowshoe hare ecology and conservation requirements in the United States is a high priority. We critically examined unchallenged perceptions and important research needs related to lynx and hare ecology and conservation at the southern extent of their range. Contrary to popular dogma, lynx do not require old-growth forest for denning, but further research on lynx and hare use of fragmented landscapes at lower latitudes is required. The contention that southern lynx are subject to higher interference or exploitative competition compared to their northern counterparts remains without strong empirical support. Lynx rely more on red squirrels (Tamiasciurus hudsonicus) and possibly other alternate prey at lower latitudes, but hares are the predominant food type for lynx across their range. Southern lynx and hare populations do not exhibit periodic cyclicity, but harvest statistics suggest that lynx abundance in the southern range is highly variable, implying that numerical fluctuations likely are fueled by immigration from Canada. Southern lynx population viability in the absence of ingress is suspect and thus maintaining connectivity with northern areas of occupancy should be a priority. Successful conservation of lynx populations in the contiguous United States will require 1) improved understanding of lynx population and habitat ecology at lower latitudes, 2) protection and management of large tracts of lynx and snowshoe hare habitat, and 3) ensured connectivity between lynx populations at the core and periphery of the species' range. However, in light of the numerous challenges facing conservation of populations of many species at their southern distributional limit, the long-term prognosis for lynx in the southern range currently is uncertain.     

Biochemical genetic variability in brown hares (Lepus europaeus) from Greece

Allozyme variability of 91 brown hares (Lepus europaeus) from seven regions in Greece was compared to existing data of Bulgarian populations to test the hypothesis of the occurrence of specific alleles in Greece, likely stemming from an isolated Late Pleistocene refugial population in the southern Balkans. This hypothesis is particularly suggested by some subfossil Late Pleistocene hare remains in Greece and the reported high mtDNA diversity in Greek hares. Allozymic diversity could be higher in Greek hares than in hares from neighboring regions as a result of the accumulation of variants in a long-lasting Pleistocene refugium. Conversely, Greek hares could exhibit reduced genetic diversity because of long-lasting low effective population sizes during the Late Glacial Maximum and a lower chance of postglacial gene flow from other populations into this rather marginal part in the southern Balkans. Horizontal starch gel electrophoresis of proteins from 35~loci revealed three alleles (Es-1 ^–162, Pep-2 ¹¹⁴, Mpi ⁸⁸) at low frequencies, which were not found in Bulgarian or any other brown hare population. In contrast, some alleles from the populations from Bulgaria and other regions of Europe were absent in the Greek samples. Population genetic statistics indicated only a slight tendency of increased gene pool diversity in Greek hares, little substructuring in Greek and Bulgarian populations, respectively, as well as an only slightly lower level of gene flow between the two neighboring regions, as compared to the gene flow within each region. The results conform to the hypothesis of a Late Pleistocene refugial population in the southern Balkans, with some few specific nuclear gene pool characteristics, but little effect on the overall genetic differentiation between Greek and Bulgarian hares.     

Unexpectedly complex gradation of coral population structure in the Nansei Islands,Japan

To establish effective locations and sizes of potential protected areas for reef ecosystems, detailed information about source and sink relationships between populations is critical, especially in archipelagic regions. Therefore, we assessed population structure and genetic diversity of Acropora tenuis, one of the dominant stony coral species in the Pacific, using 13 microsatellite markers to investigate 298 colonies from 15 locations across the Nansei Islands in southwestern Japan. Genetic diversity was not significant among sampling locations, even in possibly peripheral locations. In addition, our results showed that there are at least two populations of A. tenuis in the study area. The level of genetic differentiation between these populations was relatively low, but significant between many pairs of sampling locations. Directions of gene flow, which were estimated using a coalescence‐based approach, suggest that gene flow not only occurs from south to north, but also from north to south in various locations. Consequently, the Yaeyama Islands and the Amami Islands are potential northern and southern sources of corals. On the other hand, the Miyako Islands and west central Okinawa Island are potential sink populations. The Kerama Islands and the vicinity of Taketomi Island are potential contact points of genetic subdivision of coral populations in the Nansei Islands. We found that genetic population structure of A. tenuis in the Nansei Islands is more complex than previously thought. These cryptic populations are very important for preserving genetic diversity and should be maintained.     

High mtDNA haplotype diversity among introduced Swedish brown haresLepus europaeus

The brown hareLepus europaeus Pallas, 1778 occurs naturally in central Eurasia, but has been introduced to parts of northern Europe, South- and North America, Australia and New Zealand. Brown hares were introduced to Sweden from central Europe for hunting purposes during the 19th century. We investigated how the human--mediated brown hare colonisation of Sweden is reflected in the amount of genetic variation present by assessing variation and composition of mitochondrial DNA (mtDNA) lineages among Swedish brown hares. MtDNA from a total of 40 brown hare specimens from 15 localities were analysed for Restriction Fragment Length Polymorphisms. The haplotype diversity is surprisingly high (0.893 ± 0.002) when compared to the mtDNA diversity among brown hares on the European continent as well as to other mammalian species. Admixture of haplotypes from different source populations combined with a reduced effect of random genetic drift and a relaxed selection pressure due to rapid population growth after introduction are mechanisms that are likely to account for the observed high mtDNA haplotype diversity.     

Multiscale habitat relationships of snowshoe hares (Lepus americanus) in the mixed conifer landscape of the Northern Rockies,USA: Cross‐scale effects of horizontal cover with implications for forest management

Snowshoe hares (Lepus americanus) are an ecologically important herbivore because they modify vegetation through browsing and serve as a prey resource for multiple predators. We implemented a multiscale approach to characterize habitat relationships for snowshoe hares across the mixed conifer landscape of the northern Rocky Mountains, USA. Our objectives were to (1) assess the relationship between horizontal cover and snowshoe hares, (2) estimate how forest metrics vary across the gradient of snowshoe hare use and horizontal cover, and (3) model and map snowshoe hare occupancy and intensity of use. Results indicated that both occupancy and intensity of use by snowshoe hares increased with horizontal cover and that the effect became stronger as intensity of use increased. This underscores the importance of dense horizontal cover to achieve high use, and likely density, of snowshoe hares. Forest structure in areas with high snowshoe hare use and horizontal cover was characterized as multistoried with dense canopy cover and medium‐sized trees (e.g., 12.7–24.4 cm). The abundance of lodgepole pine (Pinus contorta) was associated with snowshoe hare use within a mixed conifer context, and the only species to increase in abundance with horizontal cover was Engelmann spruce (Picea engelmannii) and subalpine fir (Abies lasiocarpa). Our landscape‐level modeling produced similar patterns in that we observed a positive effect of lodgepole pine and horizontal cover on both occupancy and use by snowshoe hares, but we also observed a positive yet parabolic effect of snow depth on snowshoe hare occupancy. This work is among the first to characterize the multiscale habitat relationships of snowshoe hares across a mixed conifer landscape as well as to map their occupancy and intensity of use. Moreover, our results provide stand‐ and landscape‐level insights that directly relate to management agencies, which aids in conservation efforts of snowshoe hares and their associated predators.     

Population Estimates of Snowshoe Hares in the Southern Rocky Mountains

Abstract: In 1999 Canada lynx (Lynx canadensis) were reintroduced to the southern Rocky Mountains and in 2000 the species was listed as threatened under the Endangered Species Act in the contiguous United States (Colorado Division of Wildlife 2000, U.S. Fish and Wildlife Service 2000). To better evaluate the progress of this reintroduction, we conducted field studies to estimate population densities of snowshoe hares (Lepus americanus), the primary prey of lynx in Colorado, USA. We conducted our field studies in southwestern Colorado in winters 2002 and 2003. We estimated population densities in forested stands of mature Engelmann spruce (Picea engelmannii)-subalpine fir (Abies lasiocarpa) and mature lodgepole pine (Pinus contorta) using mark-recapture data and 3 methods for estimating effective area trapped: half trap interval, mean maximum distance moved (MMDM), and half MMDM. In Engelmann spruce-subalpine fir, we found density estimates ranged from 0.08 ± 0.03 (SE) hares/ha to 1.32 ± 0.15 hares/ha and in lodgepole pine, density estimates ranged from 0.06 ± 0.01 hares/ha to 0.34 ± 0.06 hares/ha, depending on year and method used for estimating effective area trapped. Our density estimates are similar to those reported at the low phase of the hare cycle in populations to the north (<0.1–1.1 hares/ha), where Canada lynx persist (Hodges 2000a). Although density estimates are a useful comparative tool, they depend upon methods used to estimate effective area trapped. Therefore, we urge caution in comparing our density estimates with those from other areas, which may have used dissimilar methods. We also examined effects of temperature and moon phase on capture success of snowshoe hares; extremely low temperatures affected capture success but moon phase did not. Capture success can be improved by trapping snowshoe hares at temperatures above their lower critical temperature (T_lc). If abundance estimates are derived from mark-recapture studies then effects of temperature should be included when modeling capture probabilities.     

Genetic diversity loss in a biodiversity hotspot: ancient DNA quantifies genetic decline and former connectivity in a critically endangered marsupial

The extent of genetic diversity loss and former connectivity between fragmented populations are often unknown factors when studying endangered species. While genetic techniques are commonly applied in extant populations to assess temporal and spatial demographic changes, it is no substitute for directly measuring past diversity using ancient DNA (aDNA). We analysed both mitochondrial DNA (mtDNA) and nuclear microsatellite loci from 64 historical fossil and skin samples of the critically endangered Western Australian woylie (Bettongia penicillata ogilbyi), and compared them with 231 (n = 152 for mtDNA) modern samples. In modern woylie populations 15 mitochondrial control region (CR) haplotypes were identified. Interestingly, mtDNA CR data from only 29 historical samples demonstrated 15 previously unknown haplotypes and detected an extinct divergent clade. Through modelling, we estimated the loss of CR mtDNA diversity to be between 46% and 91% and estimated this to have occurred in the past 2000–4000 years in association with a dramatic population decline. In addition, we obtained near‐complete 11‐loci microsatellite profiles from 21 historical samples. In agreement with the mtDNA data, a number of ‘new’ microsatellite alleles was only detected in the historical populations despite extensive modern sampling, indicating a nuclear genetic diversity loss >20%. Calculations of genetic diversity (heterozygosity and allelic rarefaction) showed that these were significantly higher in the past and that there was a high degree of gene flow across the woylie's historical range. These findings have an immediate impact on how the extant populations are managed and we recommend the implementation of an assisted migration programme to prevent further loss of genetic diversity. Our study demonstrates the value of integrating aDNA data into current‐day conservation strategies.     

Home-loving boreal hare mitochondria survived several invasions in Iberia: the relative roles of recurrent hybridisation and allele surfing

Genetic introgression from a resident species into an invading close relative can result from repeated hybridisation along the invasion front and/or allele surfing on the expansion wave. Cases where the phenomenon is massive and systematic, such as for hares (genus Lepus) in Iberia, would be best explained by recurrent hybridisation but this is difficult to prove because the donor populations are generally extinct. In the Pyrenean foothills, Lepus europaeus presumably replaced Lepus granatensis recently and the present species border is parallel to the direction of invasion, so that populations of L. granatensis in the contact zone represent proxies of existing variation before the invasion. Among three pairs of populations sampled across this border, we find less differentiation of mitochondrial DNA (mtDNA) across than along it, as predicted under recurrent hybridisation at the invasion front. Using autosomal microsatellite loci and X- and Y-linked diagnostic loci, we show that admixture across the border is quasi-absent, making it unlikely that lack of interspecific mtDNA differentiation results from ongoing gene flow. Furthermore, we find that the local species ranges are climatically contrasted, making it also unlikely that ongoing ecology-driven movement of the contact account for mtDNA introgression. The lack of mtDNA differentiation across the boundary is mostly due to sharing of mtDNA from a boreal species currently extinct in Iberia (Lepus timidus) whose mitochondria have thus remained in place since the last deglaciation despite successive invasions by two other species. Home-loving mitochondria thus witness past species distribution rather than ongoing exchanges across stabilised contact zones.     

Genetic differentiation of the pine processionary moth at the southern edge of its range: contrasting patterns between mitochondrial and nuclear markers

The pine processionary moth (Thaumetopoea pityocampa) is an important pest of coniferous forests at the southern edge of its range in Maghreb. Based on mitochondrial markers, a strong genetic differentiation was previously found in this species between western (pityocampa clade) and eastern Maghreb populations (ENA clade), with the contact zone between the clades located in Algeria. We focused on the moth range in Algeria, using both mitochondrial (a 648 bp fragment of the tRNA‐cox2) and nuclear (11 microsatellite loci) markers. A further analysis using a shorter mtDNA fragment and the same microsatellite loci was carried out on a transect in the contact zone between the mitochondrial clades. Mitochondrial diversity showed a strong geographical structure and a well‐defined contact zone between the two clades. In particular, in the pityocampa clade, two inner subclades were found whereas ENA did not show any further structure. Microsatellite analysis outlined a different pattern of differentiation, with two main groups not overlapping with the mitochondrial clades. The inconsistency between mitochondrial and nuclear markers is probably explained by sex‐biased dispersal and recent afforestation efforts that have bridged isolated populations.     

Survival of dispersing versus philopatric juvenile snowshoe hares: do dispersers die?

We used radio-telemetry to monitor the survival of dispersing and philopatric juvenile snowshoe hares ( Lepus americanus ) in southwestern Yukon Territory, Canada, during a cyclic population increase. Neither 28-d survival nor the proportion of hares surviving to breed differed significantly between juvenile hares that dispersed and those that did not, nor was there a significant relationship between dispersal distance and fate (dead or alive). Our results indicate that the overall survival cost associated with natal dispersal is low for snowshoe hares during the early increase of the hare cycle.     

High latitudes and high genetic diversity: phylogeography of a widespread boreal bird, the gray jay (Perisoreus canadensis)

We describe range-wide phylogeographic variation in gray jays (Perisoreus canadensis), a boreal Nearctic corvid that occurs today primarily in recently glaciated regions. Phylogenetic analysis of mitochondrial DNA (1041 base pairs ND2 gene; N=205, 50 localities) revealed four reciprocally monophyletic groups. One widespread clade occurs across the North American boreal zone, from Newfoundland to Alaska and southwest into Utah. Three other clades occur at lower latitudes in the montane West in Colorado, the northern Rocky Mountains, and the Pacific Northwest respectively. The geographic distribution of clades in gray jays corresponds with a general pattern that is emerging for boreal taxa, having one widespread northern clade and one or more geographically restricted southwestern clades. Population genetic analyses indicate that the larger boreal clade is genetically structured and harbors significantly more genetic diversity than those clades occurring at lower latitudes. Species distribution modeling (SDM) revealed multiple putative Pleistocene refugia including several occurring at higher latitudes. We suggest that multiple post-glacial colonization routes, some of which originate from these northern refugia, are responsible for the relatively high genetic diversity at high latitudes. Conversely, lower latitude clades show little variation, probably as a result of historical restriction to smaller geographical areas with smaller long-term population sizes. This 'upside-down' pattern of genetic diversity contrasts with the conventional view that populations of north-temperate species occupying previously glaciated habitats should possess lower levels of diversity than their southern counterparts.     

Is Size of Fecal Pellets a Reliable Indicator of Species of Leporids in the Southern Rocky Mountains?

Abstract: In the southern Rocky Mountains, USA, snowshoe hares (Lepus americanus) and mountain cottontails (Sylvilagus nuttallii) are syntopic. Previous researchers used size of fecal pellets to identify leporid species in the Southern Rockies, a common criterion for identifying species of cervids. We measured 655 fecal pellets from 10 mountain cottontails and 2,374 fecal pellets from 23 snowshoe hares. We found no relationship between the body weight of mountain cottontails and the size of their fecal pellets (r = 0.04, F = 0.01, P = 0.91) but found one for snowshoe hares (r = 0.48, F = 9.3, P = 0.005). Although the 2 species differed in the size of their fecal pellets, the difference between means (1.2 mm) was sufficiently small to require measuring individual pellets and is only applicable to individuals of adult size. Although fecal pellet counts may be used to estimate presence and relative abundance of snowshoe hares in the absence of syntopic leporids, where multiple species of leporids are syntopic this method may yield misleading results.     

Intraspecific genetic variation in the common midwife toad (Alytes obstetricans): subspecies assignment using mitochondrial and microsatellite markers

The common midwife toad (Alytes obstetricans), widely distributed in the northern half of the Iberian Peninsula and part of Western Europe, is currently subdivided into four subspecies: A. o. obstetricans, A. o. boscai, A. o. pertinax and A. o. almogavarii. However, the delimitation of these subspecies and their ranges are still under discussion because strong discordances have been found between morphological and molecular data, and especially among different genetic markers. Here, we screen a set of novel microsatellite loci and mtDNA sequences of A. obstetricans populations representative of all currently recognized subspecies to investigate the correspondence between genetic groupings inferred from clustering analysis of microsatellite genotypes and the described subspecies and test whether patterns of mtDNA variation are concordant with those genetic clusters. Our results confirm previous expectations of extremely high intraspecific diversity in A. obstetricans in Iberia. Analyses of microsatellite and mtDNA data were concordant in recovering five well‐defined groups, of which three correspond to previously defined subspecies, while the two additional clusters correspond to populations of subspecies A. o. boscai separated by the Douro River. Our results suggest the occurrence of two distinct genetic units within A. o. boscai that likely result from a long independent evolutionary history, thus deserving special attention from a conservation point of view.     

Influence of northern limit range on genetic diversity and structure in a widespread North American tree,sugar maple (Acer saccharum Marshall)

Due to climate change, the ranges of many North American tree species are expected to shift northward. Sugar maple (Acer saccharum Marshall) reaches its northern continuous distributional limit in northeastern North America at the transition between boreal mixed‐wood and temperate deciduous forests. We hypothesized that marginal fragmented northern populations from the boreal mixed wood would have a distinct pattern of genetic structure and diversity. We analyzed variation at 18 microsatellite loci from 23 populations distributed along three latitudinal transects (west, central, and east) that encompass the continuous–discontinuous species range. Each transect was divided into two zones, continuous (temperate deciduous) and discontinuous (boreal mixed wood), based on sugar maple stand abundance. Respective positive and negative relationships were found between the distance of each population to the northern limit (D_north), and allelic richness (A_R) and population differentiation (F_ST). These relations were tested for each transect separately; the pattern (discontinuous–continuous) remained significant only for the western transect. structure analysis revealed the presence of four clusters. The most northern populations of each transect were assigned to a distinct group. Asymmetrical gene flow occurred from the southern into the four northernmost populations. Southern populations in Québec may have originated from two different postglacial migration routes. No evidence was found to validate the hypothesis that northern populations were remnants of a larger population that had migrated further north of the species range after the retreat of the ice sheet. The northernmost sugar maple populations possibly originated from long‐distance dispersal.     

Disrupted phylogeographical microsatellite and chloroplast DNA patterns indicate a vicariance rather than long‐distance dispersal origin for the disjunct distribution of the Chilean endemic Dioscorea biloba (Dioscoreaceae) around the Atacama Desert

Aim The Chilean endemic Dioscorea biloba (Dioscoreaceae) is a dioecious geophyte that shows a remarkable 600 km north–south disjunction in the peripheral arid area of the Atacama Desert. Its restricted present‐day distribution and probable Neogene origin indicate that its populations have a history linked to that of the Atacama Desert, making this an ideal model species with which to investigate the biogeography of the region. Location Chile, Atacama Desert and peripheral arid area. Methods Two hundred and seventy‐five individuals from nine populations were genotyped for seven nuclear microsatellite loci, and plastid trnL–F and trnT–L sequences were obtained for a representative subset of these. Analyses included the estimation of genetic diversity and population structure through clustering, Bayesian and analysis of molecular variance analyses, and statistical parsimony networks of chloroplast haplotypes. Isolation by distance was tested against alternative dispersal hypotheses. Results Microsatellite markers revealed moderate to high levels of genetic diversity within populations, with those from the southern Limarí Valley showing the highest values and northern populations showing less exclusive alleles. Bayesian analysis of microsatellite data identified three genetic groups that corresponded to geographical ranges. Chloroplast phylogeography revealed no haplotypes shared between northern and southern ranges, and little haplotype sharing between the two neighbouring southern valleys. Dispersal models suggested the presence of extinct hypothetical populations between the southern and northern ranges. Main conclusions Our results are consistent with prolonged isolation of the northern and southern groups, mediated by the life‐history traits of the species. Significant isolation was revealed at both large and moderate distances as gene flow was not evident even between neighbouring valleys. Bayesian analyses of microsatellite and chloroplast haplotype diversity identified the southern area of Limarí as the probable area of origin of the species. Our data do not support recent dispersal of D. biloba from the southern range into Antofagasta, but indicate the fragmentation of an earlier wider range, concomitant with the Pliocene–Pleistocene climatic oscillations, with subsequent extinctions of the Atacama Desert populations and the divergence of the peripheral ones as a consequence of genetic drift.     

Novel wildlife in the Arctic: the influence of changing riparian ecosystems and shrub habitat expansion on snowshoe hares

Warming during the 20th century has changed the arctic landscape, including aspects of the hydrology, vegetation, permafrost, and glaciers, but effects on wildlife have been difficult to detect. The primary aim of this study is to examine the physical and biological processes contributing to the expanded riparian habitat and range of snowshoe hares (Lepus americanus) in northern Alaska. We explore linkages between components of the riparian ecosystem in Arctic Alaska since the 1960s, including seasonality of stream flow, air temperature, floodplain shrub habitat, and snowshoe hare distributions. Our analyses show that the peak discharge during spring snowmelt has occurred on average 3.4 days per decade earlier over the last 30 years and has contributed to a longer growing season in floodplain ecosystems. We use empirical correlations between cumulative summer warmth and riparian shrub height to reconstruct annual changes in shrub height from the 1960s to the present. The effects of longer and warmer growing seasons are estimated to have stimulated a 78% increase in the height of riparian shrubs. Earlier spring discharge and the estimated increase in riparian shrub height are consistent with observed riparian shrub expansion in the region. Our browsing measurements show that snowshoe hares require a mean riparian shrub height of at least 1.24–1.36 m, a threshold which our hindcasting indicates was met between 1964 and 1989. This generally coincides with observational evidence we present suggesting that snowshoe hares became established in 1977 or 1978. Warming and expanded shrub habitat is the most plausible reason for recent snowshoe hare establishment in Arctic Alaska. The establishment of snowshoe hares and other shrub herbivores in the Arctic in response to increasing shrub habitat is a contrasting terrestrial counterpart to the decline in marine mammals reliant on decreasing sea ice.