HOT SPOTS OF GENETIC DIVERSITY DESCENDED FROM MULTIPLE PLEISTOCENE REFUGIA IN AN ALPINE UNGULATE

Snow avalanches are a primary climate-linked driver of mountain ungulate populations

Snow is a major, climate-sensitive feature of the Earth's surface and catalyst of fundamentally important ecosystem processes. Understanding how snow influences sentinel species in rapidly changing mountain ecosystems is particularly critical. Whereas effects of snow on food availability, energy expenditure, and predation are well documented, we report how avalanches exert major impacts on an ecologically significant mountain ungulate - the coastal Alaskan mountain goat (Oreamnos americanus). Using long-term GPS data and field observations across four populations (421 individuals over 17 years), we show that avalanches caused 23-65% of all mortality, depending on area. Deaths varied seasonally and were directly linked to spatial movement patterns and avalanche terrain use. Population-level avalanche mortality, 61% of which comprised reproductively important prime-aged individuals, averaged 8% annually and exceeded 22% when avalanche conditions were severe. Our findings reveal a widespread but previously undescribed pathway by which snow can elicit major population-level impacts and shape demographic characteristics of slow-growing populations of mountain-adapted animals.

Potential and progress of studying mountain biodiversity by means of butterfly genetics and genomics

Mountains are rich in biodiversity, and butterflies are species-rich and have a good ecological and evolutionary research foundation. This review addresses the potential and progress of studying mountain biodiversity using butterflies as a model. We discuss the uniqueness of mountain ecosystems, factors influencing the distribution of mountain butterflies, representative genetic and evolutionary models in butterfly research, and evolutionary studies of mountain biodiversity involving butterfly genetics and genomics. Finally, we demonstrate the necessity of studying mountain butterflies and propose future perspectives. This review provides insights for studying the biodiversity of mountain butterflies as well as a summary of research methods for reference.

New age constraints for human entry into the Americas on the north Pacific coast

The timing of the initial peopling of the Americas is unresolved. Because the archaeological record necessitates discussion of human entry from Beringia into southern North America during the last glaciation, addressing this problem routinely involves evaluating environmental parameters then targeting areas suitable for human settlement. Vertebrate remains indicate landscape quality and are a key dataset for assessing coastal migration theories and the viability of coastal routes. Here, radiocarbon dates on vertebrate specimens and archaeological sites are calibrated to document species occurrences and the ages of human settlements across the western expansion and decay of the Cordilleran Ice Sheet (CIS) during the Late Wisconsin Fraser Glaciation in four subregions of the north Pacific coast of North America. The results show archaeological sites occur after glacial maxima and are generally consistent with the age of other securely dated earliest sites in southern North America. They also highlight gaps in the vertebrate chronologies around CIS maxima in each of the subregions that point to species redistributions and extirpations and signal times of low potential for human settlement and subsistence in a key portion of the proposed coastal migration route. This study, therefore, defines new age constraints for human coastal migration theories in the peopling of the Americas debate.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Variation in age of primiparity in mountain goats (Oreamnos americanus) estimated from horn growth increments

Variation in age of primiparity is important for population dynamics and wildlife management because it can affect population growth. Using a novel technique based on the trade-off between annual horn growth and reproduction, we estimated the age of primiparity for 2274 female mountain goats (Oreamnos americanus (Blainville, 1816)) harvested across British Columbia, Canada, from 1976 to 2019. We then investigated spatio-temporal variation in the probability that harvested females were primiparous when aged 3, 4, or 5 years and older using Bayesian ordinal regressions. We found that the probability of primiparity at 3 years decreased over time in nearly all mountain ranges. In the Coastal Mountain range, however, the probability of primiparity at age 3 significantly increased. These results suggest that the large coastal populations of mountain goats could be more resilient to harvest than other populations in British Columbia, which may be experiencing environmental effects promoting later primiparity. Models predicting age of primiparity from annual growth measures are a valuable tool for wildlife management and could help conservation of many species.

Area Not Geographic Isolation Mediates Biodiversity Responses of Alpine Refugia to Climate Change

Climate refugia, where local populations of species can persist through periods of unfavorable regional climate, play a key role in the maintenance of regional biodiversity during times of environmental change. However, the ability of refugia to buffer biodiversity change may be mediated by the landscape context of refugial habitats. Here, we examined how plant communities restricted to refugial sky islands of alpine tundra in the Colorado Rockies are changing in response to rapid climate change in the region (increased temperature, declining snowpack, and earlier snow melt-out) and if these biodiversity changes are mediated by the area or geographic isolation of the sky island. We resampled plant communities in 153 plots at seven sky islands distributed across the Colorado Rockies at two time points separated by 12 years (2007/2008–2019/2020) and found changes in taxonomic, phylogenetic, and functional diversity over time. Specifically, we found an increase in species richness, a trend toward increased phylogenetic diversity, a shift toward leaf traits associated with the stress-tolerant end of leaf economics spectrum (e.g., lower specific leaf area, higher leaf dry matter content), and a decrease in the functional dispersion of specific leaf area. Importantly, these changes were partially mediated by refugial area but not by geographic isolation, suggesting that dispersal from nearby areas of tundra does not play a strong role in mediating these changes, while site characteristics associated with a larger area (e.g., environmental heterogeneity, larger community size) may be relatively more important. Taken together, these results suggest that considering the landscape context (area and geographic isolation) of refugia may be critical for prioritizing the conservation of specific refugial sites that provide the most conservation value.

Cryptic genetic diversity and cytonuclear discordance characterize contact among Canada jay (Perisoreus canadensis) morphotypes in western North America

Three distinct Canada jay (Perisoreus canadensis) morphotypes with easily recognizable plumage traits come into contact in western North America. Recent work demonstrated high genetic structure across the species’ range; however, patterns of genetic variation in these contact zones remain unknown. We categorized 605 individuals into one of three morphotypes (Pacific, Rocky Mountain, and Boreal) based on plumage, and genotyped individuals at the mtDNA control region and 12 microsatellite loci to assess the extent of hybridization between morphotypes. Our data showed cryptic genetic diversity and high cytonuclear discordance among morphotypes within contact zones, which is likely the result of recent and historical admixture. The distributions of the Boreal and Pacific morphotypes each showed a strong association with a single, distinct genetic group, whereas the Rocky Mountain morphotype exhibited higher genetic diversity and was associated with multiple genotypes. Our analyses show the importance of considering both plumage and genetic traits when examining contact zones between closely related taxa. Finally the data presented in this study reaffirm that the Pacific morphotype is distinct from the Boreal and Rocky Mountain morphotypes based on genetic, phenotypic and ecological data, indicating that the Pacific morphotype should be re-elevated to a full species.

Predicting the spread-risk potential of chronic wasting disease to sympatric ungulate species

Wildlife disease incidence is increasing, resulting in negative impacts on the economy, biodiversity, and potentially human health. Chronic wasting disease (CWD) is a fatal, transmissible spongiform encephalopathy of cervids (wild and captive) which continues to spread geographically resulting in exposure to potential new host species. The disease agent (PrPCWD) is a misfolded conformer of the cellular prion protein (PrPC). In Canada, the disease is endemic in Alberta and Saskatchewan, affecting mule and white-tail deer, with lesser impact on elk and moose. As the disease continues to expand, additional wild ungulate species including bison, bighorn sheep, mountain goat, and pronghorn antelope may be exposed. To better understand the species-barrier, we reviewed the current literature on taxa naturally or experimentally exposed to CWD to identify susceptible and resistant species. We created a phylogeny of these taxa using cytochrome B and found that CWD susceptibility followed the species phylogeny. Using this phylogeny we estimated the probability of CWD susceptibility for wild ungulate species. We then compared PrPC amino acid polymorphisms among these species to identify which sites segregated between susceptible and resistant species. We identified sites that were significantly associated with susceptibility, but they were not fully discriminating. Finally, we sequenced Prnp from 578 wild ungulates to further evaluate their potential susceptibility. Together, these data suggest the host-range for CWD will potentially include pronghorn, mountain goat and bighorn sheep, but bison are likely to be more resistant. These findings highlight the need for monitoring potentially susceptible species as CWD continues to expand.

Appendicular skeletal morphology of North American Martes reflect independent modes of evolution in conjunction with Pleistocene glacial cycles

Pleistocene glacial cycles are thought to have driven ecological niche shifts, including novel niche formation. North American pine martens, Martes americana and M. caurina, are exemplar taxa thought to have diverged molecularly and morphologically during Pleistocene glaciation. Previous research found correlations between Martes limb morphology with biome and climate, suggesting that appendicular evolution may have occurred via adaptation to selective pressures imposed by novel and shifting habitats. Such variation can also be achieved through non-adaptive means such as genetic drift. Here, we evaluate whether regional genetic differences reflect limb morphology differences among populations of M. americana and M. caurina by analyzing evolutionary tempo and mode of six limb elements. Our comparative phylogenetic models indicate that genetic structure predicts limb shape better than size. Marten limb size has low phylogenetic signal, and the best supported model of evolution is punctuational (kappa), with morphological and genetic divergence occurring simultaneously. Disparity through time analysis suggests that the tempo of limb evolution in Martes tracks Pleistocene glacial cycles, such that limb size may be responding to shifting climates rather than population genetic structure. Contrarily, we find that limb shape is strongly tied to genetic relationships, with high phylogenetic signal and a lambda mode of evolution. Overall, this pattern of limb size and shape variation may be the result of geographic isolation during Pleistocene glacial advance, while declines in disparity suggest hybridization during interglacial periods. Future inclusion of extinct populations of Martes, which were more morphologically and ecologically diverse, may further clarify Martes phenotypic evolution.

Spatiotemporal analyses suggest the role of glacial history and the ice-free corridor in shaping American badger population genetic variation

Recurring glacial cycles through the Quaternary period drastically altered the size and distribution of natural populations of North American flora and fauna. The "southerly refugia model" has been the longstanding framework for testing the effects of glaciation on contemporary genetic patterns; however, insights from ancient DNA have contributed to the reconstruction of more complex histories for some species. The American badger, Taxidea taxus, provides an interesting species for exploring the genetic legacy of glacial history, having been hypothesized to have postglacially emerged from a single, southerly refugium to recolonize northern latitudes. However, previous studies have lacked genetic sampling from areas where distinct glacial refugia have been hypothesized, including the Pacific Northwest and American Far North (Yukon, Alaska). In order to further investigate the phylogeographic history of American badgers, we collected mitochondrial DNA sequence data from ancient subfossil material collected within the historical range (Alaska, Yukon) and combined them with new and previously published data from across the species' contemporary distribution (n = 1,207). We reconstructed a mostly unresolved phylogenetic tree and star-like haplotype network indicative of emergence from a largely panmictic glacial refugium and recent population expansion, the latter further punctuated by significantly negative Tajima's D and Fu's Fs values. Although directionality of migration cannot be unequivocally inferred, the moderate to high levels of genetic variation exhibited by American badgers, alongside the low frequency of haplotypes with indels in the Midwest, suggest a potential recolonization into central North America after the hypothesized ice-free corridor reopened ~13,000 years ago. Overall, the expanded reconstruction of phylogeographic history of American badgers offers a broader understanding of contemporary range-wide patterns and identifies unique genetic units that can likely be used to inform conservation of at-risk populations at the northern periphery.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

Genome Assembly and Analysis of the North American Mountain Goat (Oreamnos americanus) Reveals Species-Level Responses to Extreme Environments

The North American mountain goat (Oreamnos americanus) is an iconic alpine species that faces stressors from climate change, industrial development, and recreational activities. This species’ phylogenetic position within the Caprinae lineage has not been resolved and their phylogeographic history is dynamic and controversial. Genomic data could be used to address these questions and provide valuable insights to conservation and management initiatives. We sequenced short-read genomic libraries constructed from a DNA sample of a 2.5-year-old female mountain goat at 80X coverage. We improved the short-read assembly by generating Chicago library data and scaffolding using the HiRise approach. The final assembly was 2,506 Mbp in length with an N50 of 66.6 Mbp, which is within the length range and in the upper quartile for N50 published ungulate genome assemblies. Comparative analysis identified 84 gene families unique to the mountain goat. The species demographic history in terms of effective population size generally mirrored climatic trends over the past one hundred thousand years and showed a sharp decline during the last glacial maximum. This genome assembly will provide a reference basis for future population and comparative genomic analyses.

Northern areas as refugia for temperate species under current climate warming: Atlantic salmon (Salmo salar L.) as a model in northern Europe

In this work, patterns of geographical genetic diversity in Atlantic salmon Salmo salar were studied across the whole Atlantic Arc; whether these patterns (and thus genetic population structure) were affected by water temperatures was also evaluated. Salmo salar populations were characterized using microsatellite loci and then analysed with reference to ocean surface temperature data from across the region. Analysis showed the presence of a latitudinal cline of genetic variability (higher in northern areas) and water temperatures (sea surface temperatures) determining genetic population structure (the latter in combination with genetic drift in southern populations). Under the current global change scenario, northern areas of Europe would constitute refugia for diversity in the future. This is effectively the inverse of what appears to have happened in glacial refugia during the last glacial maximum. From this perspective, the still abundant and large northern populations S. salar should be considered as precious as the small almost relict southern ones and given appropriate protection. Careful management of the species, coordinated across countries and latitudes, is needed in order to avoid its extinction in Europe.

Quantitatively evaluate the environmental impact factors of the life expectancy in Tibet, China

Life expectancy is influenced by both innate human self-factors and external environmental factors. Due to the individual difference of self-factors, the external environment impacts on life expectancy are especially important under the geoscience perspective. However, few studies discussed impacts of various socioeconomic and environmental factors comprehensively, and quantitative analysis of multi-factors is lacking. We chose 13 potential factors of socio-economy, ecological environment, geological environment and geographical environment to quantitatively analyze their impacts and their interactive influences on Tibet's life expectancy in 2010 by using geographical detector (Geodetector) algorithm, and figure out their suitable ranges or types. Our results indicated that the high life expectancy at birth in Tibet are distributed like a strip; the main control factor of life expectancy in Tibet is socioeconomy, followed by geological environment; furthermore, the highest range of socioeconomic factors and the igneous rock exposed areas correspond to the maximum value of life expectancy at birth; factors of ecological environment have a certain impact on the life expectancy of Tibetan residents, while factors of geographical environment have few effects on the life expectancy. Lower radiation, moderate temperature and rainfall may have favorable effects on the increase in the life expectancy. The life expectancy of Tibet is mainly affected by socioeconomy, geological and ecological environment; interactions between these factors can increase the impact on the life expectancy of Tibet. Consequently, this study can better understand the impact factors of the life expectancy in Tibet and could provide a reference basis for local government's policy making of relocation and population management.

The genetic legacy of 50 years of desert bighorn sheep translocations

Conservation biologists have increasingly used translocations to mitigate population declines and restore locally extirpated populations. Genetic data can guide the selection of source populations for translocations and help evaluate restoration success. Bighorn sheep (Ovis canadensis) are a managed big game species that suffered widespread population extirpations across western North America throughout the early 1900s. Subsequent translocation programs have successfully re-established many formally extirpated bighorn herds, but most of these programs pre-date genetically informed management practices. The state of Nevada presents a particularly well-documented case of decline followed by restoration of extirpated herds. Desert bighorn sheep (O. c. nelsoni) populations declined to less than 3,000 individuals restricted to remnant herds in the Mojave Desert and a few locations in the Great Basin Desert. Beginning in 1968, the Nevada Department of Wildlife translocated ~2,000 individuals from remnant populations to restore previously extirpated areas, possibly establishing herds with mixed ancestries. Here, we examined genetic diversity and structure among remnant herds and the genetic consequences of translocation from these herds using a genotyping-by-sequencing approach to genotype 17,095 loci in 303 desert bighorn sheep. We found a signal of population genetic structure among remnant Mojave Desert populations, even across geographically proximate mountain ranges. Further, we found evidence of a genetically distinct, potential relict herd from a previously hypothesized Great Basin lineage of desert bighorn sheep. The genetic structure of source herds was clearly reflected in translocated populations. In most cases, herds retained genetic evidence of multiple translocation events and subsequent admixture when founded from multiple remnant source herds. Our results add to a growing literature on how population genomic data can be used to guide and monitor restoration programs.

Protected areas' effectiveness under climate change: a latitudinal distribution projection of an endangered mountain ungulate along the Andes Range

BACKGROUND

Climate change is one of the greatest threats to biodiversity, pushing species to shift their distribution ranges and making existing protected areas inadequate. Estimating species distribution and potential modifications under climate change are then necessary for adjusting conservation and management plans; this is especially true for endangered species. An example of this issue is the huemul (Hippocamelus bisulcus), an endemic endangered deer from the southern Andes Range, with less than 2,000 individuals. It is distributed in fragmented populations along a 2,000 km latitudinal gradient, in Chile and Argentina. Several threats have reduced its distribution to <50% of its former range.

METHODS

To estimate its potential distribution and protected areas effectiveness, we constructed a species distribution model using 2,813 huemul presence points throughout its whole distribution range, together with 19 bioclimatic layers and altitude information from Worldclim. Its current distribution was projected for years 2050 and 2070 using five different Global Climate Models estimated for scenarios representing two carbon Representative Concentration Routes (RCP)-RCP4.5 and RCP6.0.

RESULTS

Based on current huemul habitat variables, we estimated 91,617 km2 of suitable habitat. In future scenarios of climate change, there was a loss of suitable habitat due to altitudinal and latitudinal variation. Future projections showed a decrease of 59.86-60.26% for the year 2050 and 58.57-64.34% for the year 2070 according to RCP4.5 and RCP6.0, respectively. Protected areas only covered only 36.18% of the present distribution, 38.57-34.94% for the year 2050 and 30.79-31.94% for 2070 under climate change scenarios.

DISCUSSION

Modeling current and future huemul distributions should allow the establishment of priority conservation areas in which to focus efforts and funds, especially areas without official protection. In this way, we can improve management in areas heavily affected by climate change to help ensure the persistence of this deer and other species under similar circumstances worldwide.

Multilocus genetic analyses and spatial modeling reveal complex population structure and history in a widespread resident North American passerine (Perisoreus canadensis)

An increasing body of studies of widely distributed, high latitude species shows a variety of refugial locations and population genetic patterns. We examined the effects of glaciations and dispersal barriers on the population genetic patterns of a widely distributed, high latitude, resident corvid, the gray jay (Perisoreus canadensis), using the highly variable mitochondrial DNA (mtDNA) control region and microsatellite markers combined with species distribution modeling. We sequenced 914 bp of mtDNA control region for 375 individuals from 37 populations and screened seven loci for 402 individuals from 27 populations across the gray jay range. We used species distribution modeling and a range of phylogeographic analyses (haplotype diversity, Φ_ST, SAMOVA, F_ST, Bayesian clustering analyses) to examine evolutionary history and population genetic structure. MtDNA and microsatellite markers revealed significant genetic differentiation among populations with high concordance between markers. Paleodistribution models supported at least five potential areas of suitable gray jay habitat during the last glacial maximum and revealed distributions similar to the gray jay's contemporary during the last interglacial. Colonization from and prolonged isolation in multiple refugia is evident. Historical climatic fluctuations, the presence of multiple dispersal barriers, and highly restricted gene flow appear to be responsible for strong genetic diversification and differentiation in gray jays.

The death of the circulatory system diseases in China: provincial socioeconomic and environmental perspective

Few studies have explored the association between circulatory system diseases (CSDs) and provincial socioeconomic and environmental factors from spatial perspective, although large literature have focused on CSD. The numbers of death of hypertension disease (HD), ischemic heart disease (IHD), and cerebrovascular disease (CVD) are investigated, and 14 representative socioeconomic and environmental factors are collected. Stepwise regression model (SRM) and geographically weighted regression model (GWRM) are applied to determine the spatial correlation between the number of death of those diseases and selected factors. The results are the following: (1) diseases exhibit a pattern of zonal distribution. Higher HD is mostly distributed in south district, whereas higher IHD and CVD are observed in the north area. (2) SO₂ emission amount (SO₂ EA) is significantly positively related with HD, while coal consumption (CC) and PM_2.5 are notably positively correlated with IHD and CVD. (3) A 10,000 tons increase in SO₂ EA results in three increases in the numbers of death of HD. For every 100 ten thousand tons (TTTs) increase in CC, the death of IHD and CVD increases by 11.1 and 15.7, while for every 1 μg/m³ increase in ambient PM_2.5 concentration, the numbers of death of IHD and CVD increase by 34.773 and 43.222, respectively. (4) Our findings show that there exist spatial differences for SO₂ EA, CC, and PM_2.5 influencing HD, IHD, and CVD. This study is expected to provide a reference for HD, IHD, and CVD control in different regions.

Population genetic isolation and limited connectivity in the purple finch (Haemorhous purpureus)

Using a combination of mitochondrial and z-linked sequences, microsatellite data, and spatio-geographic modeling, we examined historical and contemporary factors influencing the population genetic structure of the purple finch (Haemorhous purpureus). Mitochondrial DNA data show the presence of two distinct groups corresponding to the two subspecies, H. p. purpureus and H. p. californicus. The two subspecies likely survived in separate refugia during the last glacial maximum, one on the Pacific Coast and one east of the Rocky Mountains, and now remain distinct lineages with little evidence of gene flow between them. Southwestern British Columbia is a notable exception, as subspecies mixing between central British Columbia and Vancouver Island populations suggests a possible contact zone in this region. Z-linked data support two mitochondrial groups; however, Coastal Oregon and central British Columbia sites show evidence of mixing. Contemporary population structure based on microsatellite data identified at least six genetic clusters: three H. p. purpureus clusters, two H. p. californicus clusters, and one mixed cluster, which likely resulted from high site fidelity and isolation by distance, combined with sexual selection on morphological characters reinforcing subspecies differences.

Temporal and spatial mosaics: deep host association and shallow geographic drivers shape genetic structure in a widespread pinworm, Rauschtineria eutamii

Climate and host demographic cycling often shape both parasite genetic diversity and host distributions, processes that transcend a history of strict host-parasite association. We explored host associations and histories based on an evaluation of mitochondrial and nuclear sequences to reveal the underlying history and genetic structure of a pinworm, Rauschtineria eutamii, infecting 10 species of western North American chipmunks (Rodentia:Tamias, subgenus Neotamias). Rauschtineria eutamii contains divergent lineages influenced by the diversity of hosts and variation across the complex topography of western North America. We recovered six reciprocally monophyletic R. eutamii mitochondrial clades, largely supported by nuclear gene trees, exhibiting divergence levels comparable to intraspecific variation reported for other nematodes. Phylogenetic relationships among pinworm clades suggest that R. eutamii colonized an ancestral lineage of western chipmunks and lineages persisted during historical isolation in diverging Neotamias species or species groups. Pinworm diversification, however, is incongruent and asynchronous relative to host diversification. Secondarily, patterns of shallow divergence were shaped by geography through events of episodic colonization reflecting an interaction of taxon pulses and ecological fitting among assemblages in recurrent sympatry. Pinworms occasionally infect geographically proximal host species; however, host switching may be unstable or ephemeral, as there is no signal of host switching in the deeper history of R. eutamii.

Phylogeographic History of Atraphaxis Plants in Arid Northern China and the Origin of A. bracteata in the Loess Plateau

In China, species of Atraphaxis (Polygonaceae) primarily inhabit arid zones across temperate steppe and desert regions. The complex geologic history (e.g., expansion of deserts) and extreme climate shifts of the region appear to have played an important role in shaping the phylogeography of Atraphaxis. The present study focuses on species-level phylogeographic patterns of Atraphaxis in China, with the goal of determining the impact of past environmental changes, in northern China, on the evolutionary history of the genus. Five hundred and sixty-four individuals distributed among 71 populations of 11 species of Atraphaxis from across the geographic range of the genus were studied using sequence data from two plastid spacers, psbK-psbI and psbB-psbH. The results demonstrate that most chloroplast haplotypes are species-specific, except for some present among widespread species. The phylogeny of Atraphaxis was well structured, and molecular dating analyses suggest that the main divergence events occurred during the late Pliocene and Pleistocene (5.73–0.03 million years ago). The statistical dispersal-vicariance analysis (S-DIVA) results provide evidence that phylogeographic patterns for the genus were characterized by both vicariance events and regional dispersal. The presented data suggest that the rapid expansion of deserts and climatic changes in northern China during the late Pliocene and Pleistocene have driven the diversification and spread of Atraphaxis in the region. The expansion of the Tengger Desert provided appropriate conditions for the origin of A. bracteata. Additionally, a contact zone in the north of the Hexi Corridor was identified as having played a significant role as a migratory route for species in adjacent areas.

Influence of landscape features on the microgeographic genetic structure of a resident songbird

Landscape features influence individual dispersal and as a result can affect both gene flow and genetic variation within and between populations. The landscape of British Columbia, Canada, is already highly heterogeneous because of natural ecological and geological transitions, but disturbance from human-mediated processes has further fragmented continuous habitat, particularly in the central plateau region. In this study, we evaluated the effects of landscape heterogeneity on the genetic structure of a common resident songbird, the black-capped chickadee (Poecile atricapillus). Previous work revealed significant population structuring in British Columbia that could not be explained by physical barriers, so our aim was to assess the pattern of genetic structure at a microgeographic scale and determine the effect of different landscape features on genetic differentiation. A total of 399 individuals from 15 populations were genotyped for fourteen microsatellite loci revealing significant population structuring in this species. Individual- and population-based analyses revealed as many as nine genetic clusters with isolation in the north, the central plateau and the south. Moreover, a mixed modelling approach that accounted for non-independence of pairwise distance values revealed a significant effect of land cover and elevation resistance on genetic differentiation. These results suggest that barriers in the landscape influence dispersal which has led to the unexpectedly high levels of population isolation. Our study demonstrates the importance of incorporating landscape features when interpreting patterns of population differentiation. Despite taking a microgeographic approach, our results have opened up additional questions concerning the processes influencing dispersal and gene flow at the local scale.

Beringian sub-refugia revealed in blackfish (Dallia): implications for understanding the effects of Pleistocene glaciations on Beringian taxa and other Arctic aquatic fauna

Background

Pleistocene climatic instability had profound and diverse effects on the distribution and abundance of Arctic organisms revealed by variation in phylogeographic patterns documented in extant Arctic populations. To better understand the effects of geography and paleoclimate on Beringian freshwater fishes, we examined genetic variability in the genus Dallia (blackfish: Esociformes: Esocidae). The genus Dallia groups between one and three nominal species of small, cold- and hypoxia-tolerant freshwater fishes restricted entirely in distribution to Beringia from the Yukon River basin near Fairbanks, Alaska westward including the Kuskokwim River basin and low-lying areas of Western Alaska to the Amguema River on the north side of the Chukotka Peninsula and Mechigmen Bay on the south side of the Chukotka Peninsula. The genus has a non-continuous distribution divided by the Bering Strait and the Brooks Range. We examined the distribution of genetic variation across this range to determine the number and location of potential sub-refugia within the greater Beringian refugium as well as the roles of the Bering land bridge, Brooks Range, and large rivers within Beringia in shaping the current distribution of populations of Dallia. Our analyses were based on DNA sequence data from two nuclear gene introns (S7 and RAG1) and two mitochondrial genome fragments from nineteen sampling locations. These data were examined under genetic clustering and coalescent frameworks to identify sub-refugia within the greater Beringia refugium and to infer the demographic history of different populations of Dallia.

Results

We identified up to five distinct genetic clusters of Dallia. Four of these genetic clusters are present in Alaska: (1) Arctic Coastal Plain genetic cluster found north of the Brooks Range, (2) interior Alaska genetic cluster placed in upstream locations in the Kuskokwim and Yukon river basins, (3) a genetic cluster found on the Seward Peninsula, and (4) a coastal Alaska genetic cluster encompassing downstream Kuskokwim River and Yukon River basin sample locations and samples from Southwest Alaska not in either of these drainages. The Chukotka samples are assigned to their own genetic cluster (5) similar to the coastal Alaska genetic cluster. The clustering and ordination analyses implemented in Discriminant Analysis of Principal Components (DAPC) and STRUCTURE showed mostly concordant groupings and a high degree of differentiation among groups. The groups of sampling locations identified as genetic clusters correspond to geographic areas divided by likely biogeographic barriers including the Brooks Range and the Bering Strait. Estimates of sequence diversity (θ) are highest in the Yukon River and Kuskokwim River drainages near the Bering Sea. We also infer asymmetric migration rates between genetic clusters. The isolation of Dallia on the Arctic Coastal Plain of Alaska is associated with very low estimated migration rates between the coastal Alaska genetic cluster and the Arctic Coastal Plain genetic cluster.

Conclusions

Our results support a scenario with multiple aquatic sub-refugia in Beringia during the Pleistocene and the preservation of that structure in extant populations of Dallia. An inferred historical presence of Dallia across the Bering land bridge explains the similarities in the genetic composition of Dallia in West Beringia and western coastal Alaska. In contrast, historic and contemporary isolation across the Brooks Range shaped the distinctiveness of present day Arctic Coastal Plain Dallia. Overall this study uncovered a high degree of genetic structuring among populations of Dallia supporting the idea of multiple Beringian sub-refugia during the Pleistocene and which appears to be maintained to the present due to the strictly freshwater nature and low dispersal ability of this genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0413-2) contains supplementary material, which is available to authorized users.

Testing Taxon Tenacity of Tortoises: evidence for a geographical selection gradient at a secondary contact zone

We examined a secondary contact zone between two species of desert tortoise, Gopherus agassizii and G. morafkai. The taxa were isolated from a common ancestor during the formation of the Colorado River (4–8 mya) and are a classic example of allopatric speciation. However, an anomalous population of G. agassizii comes into secondary contact with G. morafkai east of the Colorado River in the Black Mountains of Arizona and provides an opportunity to examine reinforcement of species' boundaries under natural conditions. We sampled 234 tortoises representing G. agassizii in California (n - 103), G. morafkai in Arizona (n - 78), and 53 individuals of undetermined assignment in the contact zone including and surrounding the Black Mountains. We genotyped individuals for 25 STR loci and determined maternal lineage using mtDNA sequence data. We performed multilocus genetic clustering analyses and used multiple statistical methods to detect levels of hybridization. We tested hypotheses about habitat use between G. agassizii and G. morafkai in the region where they co-occur using habitat suitability models. Gopherus agassizii and G. morafkai maintain independent taxonomic identities likely due to ecological niche partitioning, and the maintenance of the hybrid zone is best described by a geographical selection gradient model.

Comparative phylogeography highlights the double-edged sword of climate change faced by arctic- and alpine-adapted mammals

Recent studies suggest that alpine and arctic organisms may have distinctly different phylogeographic histories from temperate or tropical taxa, with recent range contraction into interglacial refugia as opposed to post-glacial expansion out of refugia. We use a combination of phylogeographic inference, demographic reconstructions, and hierarchical Approximate Bayesian Computation to test for phylodemographic concordance among five species of alpine-adapted small mammals in eastern Beringia. These species (Collared Pikas, Hoary Marmots, Brown Lemmings, Arctic Ground Squirrels, and Singing Voles) vary in specificity to alpine and boreal-tundra habitat but share commonalities (e.g., cold tolerance and nunatak survival) that might result in concordant responses to Pleistocene glaciations. All five species contain a similar phylogeographic disjunction separating eastern and Beringian lineages, which we show to be the result of simultaneous divergence. Genetic diversity is similar within each haplogroup for each species, and there is no support for a post-Pleistocene population expansion in eastern lineages relative to those from Beringia. Bayesian skyline plots for four of the five species do not support Pleistocene population contraction. Brown Lemmings show evidence of late Quaternary demographic expansion without subsequent population decline. The Wrangell-St. Elias region of eastern Alaska appears to be an important zone of recent secondary contact for nearctic alpine mammals. Despite differences in natural history and ecology, similar phylogeographic histories are supported for all species, suggesting that these, and likely other, alpine- and arctic-adapted taxa are already experiencing population and/or range declines that are likely to synergistically accelerate in the face of rapid climate change. Climate change may therefore be acting as a double-edged sword that erodes genetic diversity within populations but promotes divergence and the generation of biodiversity.

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range-wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid-latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black-tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.

When east meets west: population structure of a high-latitude resident species, the boreal chickadee (Poecile hudsonicus)

The population genetic structure of northern boreal species has been strongly influenced both by the Quaternary glaciations and the presence of contemporary barriers, such as mountain ranges and rivers. We used a combination of mitochondrial DNA (mtDNA), nuclear microsatellites and spatial distribution modelling to study the population genetic structure of the boreal chickadee (Poecile hudsonicus), a resident passerine, and to investigate whether historical or contemporary barriers have influenced this northern species. MtDNA data showed evidence of eastern and western groups, with secondary admixture occurring in central Canada. This suggests that the boreal chickadee probably persisted in multiple glacial refugia, one in Beringia and at least one in the east. Palaeo-distribution modelling identified suitable habitat in Beringia (Alaska), Atlantic Canada and the southern United States, and correspond to divergence dates of 60-96 kya. Pairwise FST values for both mtDNA and microsatellites were significant for all comparisons involving Newfoundland, though mtDNA data suggest a more recent separation. Furthermore, unlike mtDNA data, nuclear data support population connectivity among the continental populations, possibly due to male-biased dispersal. Although both are significant, the isolation-by-distance signal is much stronger for mtDNA (r(2)=0.51) than for microsatellites (r(2)=0.05), supporting the hypothesis of male-biased dispersal. The population structure of the boreal chickadee was influenced by isolation in multiple refugia and contemporary barriers. In addition to geographical distance, physical barriers such as the Strait of Belle Isle and northern mountains in Alaska are restricting gene flow, whereas the Rocky Mountains in the west are a porous barrier.

Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus

Background

Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia.

Results

To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers.

Conclusions

The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.

Polyploidization in Heuchera cylindrica (Saxifragaceae) did not result in a shift in climatic requirements

PREMISE OF THE STUDY

Polyploidization is a key factor involved in the diversification of plants. Although polyploids are commonly found, there remains controversy on the mechanisms that lead to their successful establishment. One major problem that has been identified is that newly formed polyploids lack mates of the appropriate ploidy level and may experience severely reduced fertility due to nonproductive intercytotype crosses. Niche differentiation has been proposed as a primary mechanism that can alleviate this reproductive disadvantage and facilitate polyploid establishment. Here we test whether the establishment of tetraploid cytotypes of Heuchera cylindrica (Saxifragaceae) is consistent with climatic niche differentiation. •

METHODS

We use a combination of field surveys, flow cytometry and species distribution models to: (1) examine the distribution of diploid and tetraploid cytotypes; and (2) determine whether tetraploid Heuchera cylindrica occupy climates that differ from those of its diploid progenitors. •

KEY RESULTS

The geographic distributions of diploid and tetraploid cytotypes are largely allopatric as an extensive survey of 636 plants from 43 locations failed to detect any populations with both cytotypes. Although diploids and tetraploids occur in different geographic areas, polyploid Heuchera cylindrica occur almost exclusively in environments that are predicted to be suitable to diploid populations. •

CONCLUSIONS

Climatic niche differentiation does not explain the geographic distribution of tetraploid Heuchera cylindrica. We propose instead that tetraploid lineages were able to establish by taking advantage of glacial retreat and expanding into previously unoccupied sites.

Parasites in ungulates of Arctic North America and Greenland: a view of contemporary diversity, ecology, and impact in a world under change

Parasites play an important role in the structure and function of arctic ecosystems, systems that are currently experiencing an unprecedented rate of change due to various anthropogenic perturbations, including climate change. Ungulates such as muskoxen, caribou, moose and Dall's sheep are also important components of northern ecosystems and are a source of food and income, as well as a focus for maintenance of cultural traditions, for northerners. Parasites of ungulates can influence host health, population dynamics and the quality, quantity and safety of meat and other products of animal origin consumed by people. In this article, we provide a contemporary view of the diversity of nematode, cestode, trematode, protozoan and arthropod parasites of ungulates in arctic and subarctic North America and Greenland. We explore the intricate associations among host and parasite assemblages and identify key issues and gaps in knowledge that emerge in a regime of accelerating environmental transition.

The impact of Pleistocene climate change on an ancient arctic-alpine plant: multiple lineages of disparate history in Oxyria digyna

The ranges of arctic-alpine species have shifted extensively with Pleistocene climate changes and glaciations. Using sequence data from the trnH-psbA and trnT-trnL chloroplast DNA spacer regions, we investigated the phylogeography of the widespread, ancient (>3 million years) arctic-alpine plant Oxyria digyna (Polygonaceae). We identified 45 haplotypes and six highly divergent major lineages; estimated ages of these lineages (time to most recent common ancestor, T(MRCA)) ranged from ∼0.5 to 2.5 million years. One lineage is widespread in the arctic, a second is restricted to the southern Rocky Mountains of the western United States, and a third was found only in the Himalayan and Altai regions of Asia. Three other lineages are widespread in western North America, where they overlap extensively. The high genetic diversity and the presence of divergent major cpDNA lineages within Oxyria digyna reflect its age and suggest that it was widespread during much of its history. The distributions of individual lineages indicate repeated spread of Oxyria digyna through North America over multiple glacial cycles. During the Last Glacial Maximum it persisted in multiple refugia in western North America, including Beringia, south of the continental ice, and within the northern limits of the Cordilleran ice sheet. Our data contribute to a growing body of evidence that arctic-alpine species have migrated from different source regions over multiple glacial cycles and that cryptic refugia contributed to persistence through the Last Glacial Maximum.

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.

Discovery and description of a new trichostrongyloid species (Nematoda: Ostertagiinae), abomasal parasites in mountain goat, Oreamnos americanus, from the Western Cordillera of North America

Marshallagia lichtenfelsi sp. n. is a dimorphic ostertagiine nematode occurring in the abomasum of mountain goats, Oreamnos americanus, from the Western Cordillera of North America. Major and minor morphotype males and females are characterized and distinguished relative to the morphologically similar Marshallagia marshalli / Marshallagia occidentalis from North America and Marshallagia dentispicularis, along with other congeners, from the Palearctic region. The configuration of the convoluted and irregular synlophe in the cervical region of males and females of M. lichtenfelsi is apparently unique, contrasting with a continuous and parallel system of ridges among those species of Marshallagia, including M. marshalli/M. occidentalis, which have been evaluated. Specimens of M. lichtenfelsi are further defined by the rectangular form of the accessory bursal membrane (width > length) in the major morphotype and by the trapezoidal Sjöberg's organ in the minor morphotype, in addition to specific attributes of the spicules and spicule tips. We regard 12 species, including the proposed new taxon, to be valid. Primary diagnostic characters are reviewed for Marshallagia and a framework is presented for standardization of future descriptions incorporating the synlophe in males and females and the structure of the spicules and genital cone in major and minor morphotype males. The center of diversity for species of Marshallagia is the mountain-steppe region of central Eurasia where 11 species (including the Holarctic M. marshalli) are recognized in association with Caprini, Rupicaprini, and Antelopinae; only 2 species occur in the Nearctic. In this assemblage, M. lichtenfelsi is endemic to North America and limited in host distribution to mountain goats. An intricate history for refugial isolation and population fragmentation demonstrated for mountain goats and wild sheep indicates the potential for considerable cryptic diversity for Marshallagia and other nematodes. Shifting patterns of contact and sympatry among assemblages of ungulates during the Pleistocene are consistent with geographic and host colonization as a process involved in diversification of these parasites.

Northern host-parasite assemblages: history and biogeography on the borderlands of episodic climate and environmental transition

Diversity among assemblages of mammalian hosts and parasites in northern terrestrial ecosystems was structured by a deep history of biotic and abiotic change that overlies a complex geographic arena. Since the Pliocene, Holarctic ecosystems assembled in response to shifting climates (glacial and interglacial stages). Cycles of episodic dispersal/isolation and diversification defined northern diversity on landscape to regional scales. Episodes of geographic expansion and colonisation linked Eurasia and North America across Beringia and drove macroevolutionary structure of host and parasite associations. Asynchronous dispersal from centres of origin in Eurasia into the Nearctic resulted in gradients in parasite diversity in the carnivoran, lagomorph, rodent and artiodactyl assemblages we reviewed. Recurrent faunal interchange and isolation in conjunction with episodes of host colonisation have produced a mosaic structure for parasite faunas and considerable cryptic diversity among nematodes and cestodes. Mechanisms of invasion and geographic colonisation leading to the establishment of complex faunal assemblages are equivalent in evolutionary and ecological time, as demonstrated by various explorations of diversity in these high-latitude systems. Our ability to determine historical responses to episodic shifts in global climate may provide a framework for predicting the cascading effects of contemporary environmental change.