High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end‐glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land‐bridge route deriving from a “Carpathian” glacial refugium and one via a north‐eastern route from an “Eastern” glacial refugium near the Ural Mountains. Clustering of genome‐wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red‐backed voles (Myodes rutilus) took place in Fennoscandia just after end‐glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.     

Niche differentiation in a postglacial colonizer,the bank vole Clethrionomys glareolus

Species‐level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species‐wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid‐Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.     

Beyond the Mediterranean peninsulas: evidence of central European glacial refugia for a temperate forest mammal species, the bank vole (Clethrionomys glareolus)

This study details the phylogeographic pattern of the bank vole, Clethrionomys glareolus, a European rodent species strongly associated with forest habitat. We used sequences of 1011 base pairs of the mitochondrial DNA cytochrome b gene from 207 bank voles collected in 62 localities spread throughout its distribution area. Our results reveal the presence of three Mediterranean (Spanish, Italian and Balkan) and three continental (western, eastern and 'Ural') phylogroups. The endemic Mediterranean phylogroups did not contribute to the post-glacial recolonization of much of the Palaearctic range of species. Instead, the major part of this region was apparently recolonized by bank voles that survived in glacial refugia in central Europe. Moreover, our phylogeographic analyses also reveal differentiated populations of bank voles in the Ural mountains and elsewhere, which carry the mitochondrial DNA of another related vole species, the ruddy vole (Clethrionomys rutilus). In conclusion, this study demonstrates a complex phylogeographic history for a forest species in Europe which is sufficiently adaptable that, facing climate change, survives in relict southern and northern habitats. The high level of genetic diversity characterizing vole populations from parts of central Europe also highlights the importance of such regions as a source of intraspecific genetic biodiversity.     

Rivers,not refugia,drove diversification in arboreal,sub‐Saharan African snakes

The relative roles of rivers versus refugia in shaping the high levels of species diversity in tropical rainforests have been widely debated for decades. Only recently has it become possible to take an integrative approach to test predictions derived from these hypotheses using genomic sequencing and paleo‐species distribution modeling. Herein, we tested the predictions of the classic river, refuge, and river‐refuge hypotheses on diversification in the arboreal sub‐Saharan African snake genus Toxicodryas. We used dated phylogeographic inferences, population clustering analyses, demographic model selection, and paleo‐distribution modeling to conduct a phylogenomic and historical demographic analysis of this genus. Our results revealed significant population genetic structure within both Toxicodryas species, corresponding geographically to river barriers and divergence times from the mid‐Miocene to Pliocene. Our demographic analyses supported the interpretation that rivers are indications of strong barriers to gene flow among populations since their divergence. Additionally, we found no support for a major contraction of suitable habitat during the last glacial maximum, allowing us to reject both the refuge and river‐refuge hypotheses in favor of the river‐barrier hypothesis. Based on conservative interpretations of our species delimitation analyses with the Sanger and ddRAD data sets, two new cryptic species are identified from east‐central Africa. This study highlights the complexity of diversification dynamics in the African tropics and the advantages of integrative approaches to studying speciation in tropical regions.     

The role of anthropogenic dispersal in shaping the distribution and genetic composition of a widespread North American tree species

Dispersal and colonization are among the most important ecological processes for species persistence as they allow species to track changing environmental conditions. During the last glacial maximum (LGM), many cold‐intolerant Northern Hemisphere plants retreated to southern glacial refugia. During subsequent warming periods, these species expanded their ranges northward. Interestingly, some tree species with limited seed dispersal migrated considerable distances after the LGM ~19,000 years before present (YBP). It has been hypothesized that indigenous peoples may have dispersed valued species, in some cases beyond the southern limits of the Laurentide Ice Sheet. To investigate this question, we employed a molecular genetics approach on a widespread North American understory tree species whose fruit was valued by indigenous peoples. Twenty putative anthropogenic (near pre‐Columbian habitations) and 62 wild populations of Asimina triloba (pawpaw), which produces the largest edible fruit of any North American tree, were genetically assayed with nine microsatellite loci. Putative anthropogenic populations were characterized by reduced genetic diversity and greater excess heterozygosity relative to wild populations. Anthropogenic populations in regions that were glaciated during the LGM had profiles consistent with founder effects and reduced gene flow, and shared rare alleles with wild populations hundreds of kilometers away (mean = 723 km). Some of the most compelling evidence for human‐mediated dispersal is that putative anthropogenic and wild populations sharing rare alleles were separated by significantly greater distances (mean = 695 km) than wild populations sharing rare alleles (mean = 607 km; p = .014). Collectively, the genetic data suggest that long‐distance dispersal played an important role in the distribution of pawpaw and is consistent with the hypothesized role of indigenous peoples.     

Post‐glacial colonization of eastern Europe from the Carpathian refugium: evidence from mitochondrial DNA of the common vole Microtus arvalis

There is now considerable evidence for the survival of temperate species within glacial refugia that were situated at relatively high latitudes, notably the Carpathian Basin and Dordogne region in Europe. However, the prevalence of fossil remains in such locations is rarely matched by molecular evidence for their contribution to subsequent geographical and demographic expansion of the species in question. One obstacle to this has been insufficient analysis of modern samples from the relevant areas, in particular the parts of eastern Europe that surround the Carpathian refugium. In the present study, we examine the patterns of variation in mitochondrial DNA of the common vole (Microtus arvalis), obtained from existing museum specimens and from newly‐collected samples obtained in this area. We show that common voles from one of six extant mitochondrial DNA lineages have colonized most of the species' range in eastern Europe. We contend that the post‐glacial dispersal of this lineage most likely originated from the Carpathian refugium, adding support to the argument that such northern refugia made an important contribution to existing genetic diversity in Europe. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, ●● , ●●–●●.     

Multiple introductions and overwintering shape the progressive invasion of Aedes albopictus beyond the Alps

Aedes albopictus originates from Southeast Asia and is considered one of the most invasive species globally. This mosquito is a nuisance and a disease vector of significant public health relevance. In Europe, Ae. albopictus is firmly established and widespread south of the Alps, a mountain range that forms a formidable biogeographic barrier to many organisms. Recent reports of Ae. albopictus north of the Alps raise questions of (1) the origins of its recent invasion, and (2) if this mosquito has established overwintering populations north of the Alps. To answer these questions, we analyzed population genomic data from >4000 genome‐wide SNPs obtained through double‐digest restriction site‐associated DNA sequencing. We collected SNP data from specimens from six sites in Switzerland, north and south of the Alps, and analyzed them together with specimens from other 33 European sites, five from the Americas, and five from its Asian native range. At a global level, we detected four genetic clusters with specimens from Indonesia, Brazil, and Japan as the most differentiated, whereas specimens from Europe, Hong Kong, and USA largely overlapped. Across the Alps, we detected a weak genetic structure and high levels of genetic admixture, supporting a scenario of rapid and human‐aided dispersal along transportation routes. While the genetic pattern suggests frequent re‐introductions into Switzerland from Italian sources, the recovery of a pair of full siblings in two consecutive years in Strasbourg, France, suggests the presence of an overwintering population north of the Alps. The suggestion of overwintering populations of Ae. albopictus north of the Alps and the expansion patterns identified points to an increased risk of further northward expansion and the need for increased surveillance of mosquito populations in Northern Europe.     

Adaptive phylogeography: functional divergence between haemoglobins derived from different glacial refugia in the bank vole

Over the years, researchers have used presumptively neutral molecular variation to infer the origins of current species'' distributions in northern latitudes (especially Europe). However, several reported examples of genic and chromosomal replacements suggest that end-glacial colonizations of particular northern areas may have involved genetic input from different source populations at different times, coupled with competition and selection. We investigate the functional consequences of differences between two bank vole (Clethrionomys glareolus) haemoglobins deriving from different glacial refugia, one of which partially replaced the other in Britain during end-glacial climate warming. This allows us to examine their adaptive divergence and hence a possible role of selection in the replacement. We determine the amino acid substitution Ser52Cys in the major expressed β-globin gene as the allelic difference. We use structural modelling to reveal that the protein environment renders the 52Cys thiol a highly reactive functional group and we show its reactivity in vitro. We demonstrate that possessing the reactive thiol in haemoglobin increases the resistance of bank vole erythrocytes to oxidative stress. Our study thus provides striking evidence for physiological differences between products of genic variants that spread at the expense of one another during colonization of an area from different glacial refugia.     

Genetic structure and colonization processes in European populations of the common vole, Microtus arvalis

The level of genetic differentiation within and between evolutionary lineages of the common vole (Microtus arvalis) in Europe was examined by analyzing mitochondrial sequences from the control region (mtDNA) and 12 nuclear microsatellite loci (nucDNA) for 338 voles from 18 populations. The distribution of evolutionary lineages and the affinity of populations to lineages were determined with additional sequence data from the mitochondrial cytochrome b gene. Our analyses demonstrated very high levels of differentiation between populations (overall FST: mtDNA 70%; nucDNA 17%). The affinity of populations to evolutionary lineages was strongly reflected in mtDNA but not in nucDNA variation. Patterns of genetic structure for both markers visualized in synthetic genetic maps suggest a postglacial range expansion of the species into the Alps, as well as a potentially more ancient colonization from the northeast to the southwest of Europe. This expansion is supported by estimates for the divergence times between evolutionary lineages and within the western European lineage, which predate the last glacial maximum (LGM). Furthermore, all measures of genetic diversity within populations increased significantly with longitude and showed a trend toward increase with latitude. We conclude that the detected patterns are difficult to explain only by range expansions from separate LGM refugia close to the Mediterranean. This suggests that some M. arvalis populations persisted during the LGM in suitable habitat further north and that the gradients in genetic diversity may represent traces of a more ancient colonization of Europe by the species.     

Female differential allocation in response to extrapair offspring and social mate attractiveness

Renewed debate over what benefits females might gain from producing extra‐pair offspring emphasizes the possibility that apparent differences in quality between within‐pair and extra‐pair offspring are confounded by greater maternal investment in extra‐pair offspring. Moreover, the attractiveness of a female''s social mate can also influence contributions of both partners to a reproductive attempt. Here, we explore the complexities involved in parental investment decisions in response to extra‐pair offspring and mate attractiveness with a focus on the female point of view. Adult zebra finches paired and reproduced in a colony setting. A male''s early‐life diet quality and his extra‐pair reproductive success were used as metrics of his mating attractiveness. Females paired with males that achieved extra‐pair success laid heavier eggs than other females and spent less time attending their nests than their mates or other females. Extra‐pair nestlings were fed more protein‐rich hen''s egg than within‐pair nestlings. Females producing extra‐pair offspring had more surviving sons than females producing only within‐pair offspring. Collectively, results show that females differentially allocate resources in response to offspring extra‐pair status and their social mate''s attractiveness. Females may also obtain fitness benefits through the production of extra‐pair offspring.     

Signature of mid‐Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within‐ and among‐population genetic diversity using phylogeographic reconstructions, tests of isolation‐by‐distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation‐by‐distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine–Corsica–Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid‐Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east–west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.     

Genetically diverse but with surprisingly little geographical structure: the complex history of the widespread herb Carex nigra (Cyperaceae)

Aim Based on extensive range‐wide sampling, we address the phylogeographical history of one of the most widespread and taxonomically complex sedges in Europe, Carex nigra s. lat. We compare the genetic structure of the recently colonized northern areas (front edge) and the long‐standing southern areas (rear edge), and assess the potential genetic basis of suggested taxonomic divisions at the rank of species and below. Location Amphi‐Atlantic, central and northern Europe, circum‐Mediterranean mountain ranges, central Siberia, Himalayas. Methods A total of 469 individuals sampled from 83 populations, covering most of the species’ range, were analysed with amplified fragment length polymorphism (AFLP) and chloroplast DNA (cpDNA) markers. Bayesian clustering, principal coordinates analysis, and estimates of diversity and differentiation were used for the analysis of AFLP data. CpDNA data were analysed with statistical parsimony networks and maximum parsimony and Bayesian inference of phylogenetic trees. Results Overall genetic diversity was high, but differentiation among populations was limited. Major glacial refugia were inferred in the Mediterranean Basin and in western Russia; in addition, there may have been minor refugia in the North Atlantic region. In the southern part of the range, we found high levels, but geographically quite poorly structured genetic diversity, whereas the levels of genetic diversity varied among different areas in the north. North American populations were genetically very similar to the European populations. Main conclusions The data are consistent with extensive gene flow, which has obscured the recent history of the taxon. The limited differentiation in the south probably results from the mixing of lineages expanding from several local refugia. Northward post‐glacial colonization resulted in a leading‐edge pattern of low diversity in the Netherlands, Belgium, Scotland and Iceland, whereas the observed high diversity levels in Fennoscandia suggest broad‐fronted colonization from the south as well as from the east. The patterns found in the American populations are consistent with post‐glacial colonization, possibly even with anthropogenic introduction from Europe. Our data also suggest that the tussock‐forming populations of C. nigra, often referred to as a distinct species (Carex juncella), represent an ecotype that has originated repeatedly from different populations with creeping rhizomes.     

Niche partitioning of avian predators in northern grasslands amended by biosolids

Many food webs are affected by bottom‐up nutrient addition, as additional biomass or productivity at a given trophic level can support more consumers. In turn, when prey are abundant, predators may converge on the same diets rather than partitioning food resources. Here, we examine the diets and habitat use of predatory and omnivorous birds in response to biosolids amendment of northern grasslands used as grazing range for cattle in British Columbia, Canada. From an ecosystem management perspective, we test whether dietary convergence occurred and whether birds preferentially used the pastures with biosolids. Biosolids treatments increased Orthoptera densities and our work occurred during a vole (Microtus spp.) population peak, so both types of prey were abundant. American Kestrels (Falco sparverius) consumed both small mammals and Orthoptera. Short‐eared Owls (Asio flammeus) and Long‐eared owls (Asio otus) primarily ate voles (>97% of biomass consumed) as did Northern Harriers (Circus hudsonius, 88% vole biomass). Despite high dietary overlap, these species had minimal spatial overlap, and Short‐eared Owls strongly preferred pastures amended with biosolids. Common Ravens (Corvus corax), Black‐billed Magpies (Pica hudsonia), and American Crows (Corvus brachyrhynchos) consumed Orthoptera, Coleoptera, vegetation, and only a few small mammals; crows avoided pastures with biosolids. Thus, when both insect and mammalian prey were abundant, corvids maintained omnivorous diets, whereas owls and Harriers specialized on voles. Spatial patterns were more complex, as birds were likely responding to prey abundance, vegetation structure, and other birds in this consumer guild.     

A complex scenario of glacial survival in Mediterranean and continental refugia of a temperate continental vole species (Microtus arvalis) in Europe

The role of glacial refugia in shaping contemporary species distribution is a long-standing question in phylogeography and evolutionary ecology. Recent studies are questioning previous paradigms on glacial refugia and postglacial recolonization pathways in Europe, and more flexible phylogeographic scenarios have been proposed. We used the widespread common vole Microtus arvalis as a model to investigate the origin, locations of glacial refugia, and dispersal pathways, in the group of “Continental” species in Europe. We used a Bayesian spatiotemporal diffusion analysis (relaxed random walk model) of cytochrome b sequences across the species range, including newly collected individuals from 10 Iberian localities and published sequences from 68 localities across 22 European countries. Our data suggest that the species originated in Central Europe, and we revealed the location of multiple refugia (in both southern peninsulas and continental regions) for this continental model species. Our results confirm the monophyly of Iberian voles and the pre-LGM divergence between Iberian and European voles. We found evidence of restricted postglacial dispersal from refugia in Mediterranean peninsulas. We inferred a complex evolutionary and demographic history of M. arvalis in Europe over the last 50,000 years that does not adequately fit previous glacial refugial scenarios. The phylogeography of M. arvalis provides a paradigm of ice-age survival of a temperate continental species in western and eastern Mediterranean peninsulas (sources of endemism) and multiple continental regions (sources of postglacial spread). Our findings also provide support for a major role of large European river systems in shaping geographic boundaries of M. arvalis in Europe.     

Female extra‐pair behavior is not associated with reduced paternal care in Thorn‐tailed Rayadito

Extra‐pair behavior is present in 76% of socially monogamous bird species with biparental care. This behavior may produce costs to females related to a reduction in paternal care. We estimated the percentage of extra‐pair offspring and quantified paternal care in 44 nests of Thorn‐tailed Rayadito (Aphrastura spinicauda) to assess whether males reduce their parental care when females obtain extra‐pair fertilizations. We used data from a sub‐Antarctic population of Rayadito located on Navarino Island (55°4′S, 67°40′W), southern Chile. We found no statistical support for a relationship between variation in paternal care and the percentage of extra‐pair offspring. We discuss how the inability of breeding males to assess their genetic paternity and potential restrictions on behavioral flexibility may explain this result. Additionally, if paternal care is subjected to sexual selection, this could limit a facultative response to female extra‐pair behavior by males. Finally, it is possible that a reduction in paternal care might not have evolved in this particular locality given the low frequency of extra‐pair paternity in our study population.     

Increased genetic structuring of isolated Salamandra salamandra populations (Caudata: Salamandridae) at the margins of the Carpathian Mountains

The traditional southern Pleistocene refugia hypothesis in Europe has lately been challenged for several animal and plant species. The Carpathian Basin, especially at the marginal regions, is one of the recently recognized biodiversity hotspots in Europe. Marginal populations are prone to have lower genetic diversity and higher genetic differentiation than central populations. Here, we examined one mitochondrial DNA fragment (D‐loop) and nine nuclear (microsatellite) loci to describe the genetic diversity and phylogeographical pattern of fire salamander (Salamandra salamandra) populations in the Carpathian Basin with focusing on the southern margins of the Western Carpathians, where isolated populations of this species are present. Analyses of microsatellites indicated reduced genetic diversity for most of the isolated populations. Based on the mitochondrial DNA, only two haplotypes were found, whereas the analyses with the nuclear markers revealed a more recent genetic split between Western (Alpine) and Eastern (Carpathian) populations, and separated the Apuseni Mountains population (part of the Western Carpathians). Using approximate Bayesian computation analyses, we identified the most probable colonization scenario for the isolated North Hungarian Carpathian Basin populations. The split between isolated salamander populations from the central populations in the Carpathian Mountains dates back to the beginning of the Late Pleistocene, while the split between most of the Hungarian populations can be associated with the Last Glacial Maximum. We found evidence for long‐time isolation between the marginal Carpathian Basin and central populations. Our results also show that S. salamandra survived glacial periods in the temperate forests of north‐east Pannonia (North Hungarian Mountains), confirming that the Carpathian Basin served as important northerly refugia during the Pleistocene climatic oscillations.     

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five‐prime region of cytochrome c oxidase subunit I (COI‐5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI‐5P lineages. The distribution of COI‐5P haplotypes may not demonstrate significant phylogeographic structure (N _ST > G _ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI‐5P haplotypes. For both COI‐5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica''s three major mitochondrial COI‐5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015–0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species’ distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.     

Genomic association with pathogen carriage in bighorn sheep (Ovis canadensis)

Genetic composition can influence host susceptibility to, and transmission of, pathogens, with potential population‐level consequences. In bighorn sheep (Ovis canadensis), pneumonia epidemics caused by Mycoplasma ovipneumoniae have been associated with severe population declines and limited recovery across North America. Adult survivors either clear the infection or act as carriers that continually shed M. ovipneumoniae and expose their susceptible offspring, resulting in high rates of lamb mortality for years following the outbreak event. Here, we investigated the influence of genomic composition on persistent carriage of M. ovipneumoniae in a well‐studied bighorn sheep herd in the Wallowa Mountains of Oregon, USA. Using 10,605 SNPs generated using RADseq technology for 25 female bighorn sheep, we assessed genomic diversity metrics and employed family‐based genome‐wide association methodologies to understand variant association and genetic architecture underlying chronic carriage. We observed no differences among genome‐wide diversity metrics (heterozygosity and allelic richness) between groups. However, we identified two variant loci of interest and seven associated candidate genes, which may influence carriage status. Further, we found that the SNP panel explained ~55% of the phenotypic variance (SNP‐based heritability) for M. ovipneumoniae carriage, though there was considerable uncertainty in these estimates. While small sample sizes limit conclusions drawn here, our study represents one of the first to assess the genomic factors influencing chronic carriage of a pathogen in a wild population and lays a foundation for understanding genomic influence on pathogen persistence in bighorn sheep and other wildlife populations. Future research should incorporate additional individuals as well as distinct herds to further explore the genomic basis of chronic carriage.     

Life in the desert: The impact of geographic and environmental gradients on genetic diversity and population structure of Ivesia webberi

For range‐restricted species with disjunct populations, it is critical to characterize population genetic structure, gene flow, and factors that influence functional connectivity among populations in order to design effective conservation programs. In this study, we genotyped 314 individuals from 16 extant populations of Ivesia webberi, a United States federally threatened Great Basin Desert using six microsatellite loci. We assessed the effects of Euclidean distance, landscape features, and ecological dissimilarity on the pairwise genetic distance of the sampled populations, while also testing for a potential relationship between I. webberi genetic diversity and diversity in the vegetative communities. The results show low levels of genetic diversity overall (H _e = 0.200–0.441; H _o = 0.192–0.605) and high genetic differentiation among populations. Genetic diversity was structured along a geographic gradient, congruent with patterns of isolation by distance. Populations near the species’ range core have relatively high genetic diversity, supporting in part a central‐marginal pattern, while also showing some evidence for a metapopulation dynamic. Peripheral populations have lower genetic diversity, significantly higher genetic distances, and higher relatedness. Genotype cluster admixture results suggest a complex dispersal pattern among populations with dispersal direction and distance varying on the landscape. Pairwise genetic distance strongly correlates with elevation, actual evapotranspiration, and summer seasonal precipitation, indicating a role for isolation by environment, which the observed phenological mismatches among the populations also support. The significant correlation between pairwise genetic distance and floristic dissimilarity in the germinated soil seed bank suggests that annual regeneration in the plant communities contribute to the maintenance of genetic diversity in I. webberi.     

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.