Evidence for cryptic glacial refugia from North American mountain sheep mitochondrial DNA

Local thermal adaptation and local temperature regimes drive the performance of a parasitic helminth under climate change: The case of Marshallagia marshalli from wild ungulates

Across a species' range, populations are exposed to their local thermal environments, which on an evolutionary scale, may cause adaptative differences among populations. Helminths often have broad geographic ranges and temperature-sensitive life stages but little is known about whether and how local thermal adaptation can influence their response to climate change. We studied the thermal responses of the free-living stages of Marshallagia marshalli, a parasitic nematode of wild ungulates, along a latitudinal gradient. We first determine its distribution in wild sheep species in North America. Then we cultured M. marshalli eggs from different locations at temperatures from 5 to 38°C. We fit performance curves based on the metabolic theory of ecology to determine whether development and mortality showed evidence of local thermal adaptation. We used parameter estimates in life-cycle-based host-parasite models to understand how local thermal responses may influence parasite performance under general and location-specific climate-change projections. We found that M. marshalli has a wide latitudinal and host range, infecting wild sheep species from New Mexico to Yukon. Increases in mortality and development time at higher temperatures were most evident for isolates from northern locations. Accounting for location-specific parasite parameters primarily influenced the magnitude of climate change parasite performance, while accounting for location-specific climates primarily influenced the phenology of parasite performance. Despite differences in development and mortality among M. marshalli populations, when using site-specific climate change projections, there was a similar magnitude of impact on the relative performance of M. marshalli among populations. Climate change is predicted to decrease the expected lifetime reproductive output of M. marshalli in all populations while delaying its seasonal peak by approximately 1 month. Our research suggests that accurate projections of the impacts of climate change on broadly distributed species need to consider local adaptations of organisms together with local temperature profiles and climate projections.

Geography vs. past climate: the drivers of population genetic structure of the Himalayan langur

Background

Contemporary species distribution, genetic diversity and evolutionary history in many taxa are shaped by both historical and current climate as well as topography. The Himalayas show a huge variation in topography and climatic conditions across its entire range, and have experienced major climatic fluctuations in the past. However, very little is known regarding how this heterogenous landscape has moulded the distribution of Himalayan fauna. A recent study examined the effect of these historical events on the genetic diversity of the Himalayan langurs in Nepal Himalaya. However, this study did not include the samples from the Indian Himalayan region (IHR). Therefore, here we revisit the questions addressed in the previous study with a near complete sampling from the IHR, along with the samples from the Nepal Himalaya. We used the mitochondrial Cytochrome-b (Cyt-b, 746 bp) region combined with multiple phylogeographic analyses and palaeodistribution modelling.

Results

Our dataset contained 144 sequences from the IHR as well as the Nepal Himalaya. Phylogenetic analysis showed a low divergent western clade nested within high divergent group of eastern lineages and in the network analysis we identified 22 haplotypes over the entire distribution range of the Himalayan langurs. Samples from the Nepal Himalaya showed geographically structured haplotypes corresponding to different river barriers, whereas samples from IHR showed star-like topology with no structure. Our statistical phylogeography analysis using diyABC supported the model of east to west colonisation of these langurs with founder event during colonisation. Analysis of demographic history showed that the effective population size of the Himalayan langurs decreased at the onset of last glacial maximum (LGM) and started increasing post LGM. The palaeodistribution modelling showed that the extent of suitable habitat shifted from low elevation central Nepal, and adjoining parts of north India, during LGM to the western Himalaya at present.

Conclusion

The current genetic diversity and distribution of Himalayan langurs in the Nepal Himalaya has been shaped by river barriers, whereas the rivers in the IHR had relatively less time to act as a strong genetic barrier after the recent colonisation event. Further, the post LGM expansion could have had confounding effect on Himalayan langur population structure in both Nepal Himalaya and IHR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02054-1.

Phylogeography and population genetics of a widespread cold-adapted ant, Prenolepis imparis

As species arise, evolve, and diverge, they are shaped by forces that unfold across short and long time scales and at both local and vast geographic scales. It is rare, however, to be able document this history across broad sweeps of time and space in a single species. Here, we report the results of a continental-scale phylogenomic analysis across the entire range of a widespread species. We analyzed sequences of 1,402 orthologous Ultraconserved Element (UCE) loci from 75 individuals to identify population genetic structure and historical demographic patterns across the continent-wide range of a cold-adapted ant, the winter ant, Prenolepis imparis. We recovered five well-supported, genetically isolated clades representing lineages that diverged from 8.2-2.2 million years ago. These include: 1) an early diverging lineage located in Florida, 2) a lineage that spans the southern United States, 3) populations that extend across the midwestern and northeastern United States, 4) populations from the western United States, and 5) populations in southwestern Arizona and Mexico. Population genetic analyses revealed little or no gene flow among these lineages, but patterns consistent with more recent gene flow among populations within lineages, and localized structure with migration in the western United States. High support for five major geographic lineages and lack of evidence of contemporary gene flow indicate in situ diversification across the species' range, producing relatively ancient lineages that persisted through subsequent climate change and glaciation during the Quaternary.

Neutral and Selective Processes Shape MHC Diversity in Roe Deer in Slovenia

Simple Summary

Disease prevention and appropriate wildlife management are among the major challenges in wildlife conservation. In the present study, we made a first assessment of the variability of major histocompatibility complex (MHC) genes in roe deer in Slovenia and evaluated local population adaptation by comparing MHC variability with neutral microsatellites. We discovered three new MHC DRB exon 2 alleles in addition to seven previously described in the literature. Moreover, we found evidence of historical positive selection, as selection analysis indicated that approx. 10% of the encoded amino acids were subjected to episodic positive selection. This study provides the basis for further research on immunogenetic variation in roe deer and highlights opportunities to incorporate genetic data into science-based population management.

Abstract

Disease control and containment in free-ranging populations is one of the greatest challenges in wildlife management. Despite the importance of major histocompatibility complex (MHC) genes for immune response, an assessment of the diversity and occurrence of these genes is still rare in European roe deer, the most abundant and widespread large mammal in Europe. Therefore, we examined immunogenetic variation in roe deer in Slovenia to identify species adaptation by comparing the genetic diversity of the MHC genes with the data on neutral microsatellites. We found ten MHC DRB alleles, three of which are novel. Evidence for historical positive selection on the MHC was found using the maximum likelihood codon method. Patterns of MHC allelic distribution were not congruent with neutral population genetic findings. The lack of population genetic differentiation in MHC genes compared to existing structure in neutral markers suggests that MHC polymorphism was influenced primarily by balancing selection and, to a lesser extent, by neutral processes such as genetic drift, with no clear evidence of local adaptation. Selection analyses indicated that approx. 10% of amino acids encoded under episodic positive selection. This study represents one of the first steps towards establishing an immunogenetic map of roe deer populations across Europe, aiming to better support science-based management of this important game species.

Ancient hybridization patterns between bighorn and thinhorn sheep

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

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.

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.

Global phylogeography of toxic cyanobacteria Moorea producens reveals distinct genetic partitioning influenced by Proterozoic glacial cycles

Lyngbya majuscula is a marine filamentous cyanobacteria belonging to the family Oscillatoriaceae. Recent phylogenetic analyses of L. majuscula have reclassified a subset of this species into various genera such as Moorea, Okeania and Dapis. From the genus Moorea, Moorea producens is a toxic invasive cyanobacterium that produces bioactive secondary metabolites that can cause severe inflammation and blistering. Despite the global distribution of M. producens, little information is available on their origin, patterns of dispersal and population structure. In this study, the spatial population structure of M. producens was investigated using near-complete 16S rRNA sequences. Analysis of the global population of M. producens by Isolation by Distance and STRUCTURE revealed two significantly distinct cosmopolitan populations that were separated by a genetic break. Lineage-specific divergence estimates of 147 cyanobacterial taxa, based on a relaxed molecular clock indicated the first global emergence of M. producens during the Mesoarchean and a subsequent global recolonization during the Mesoproterozoic period. We conclude that the genetic discontinuity between both cosmopolitan populations is attributed to refugia associated with Proterozoic glacial cycles.

Current evidence allows multiple models for the peopling of the Americas

Some recent academic and popular literature implies that the problem of the colonization of the Americas has been largely resolved in favor of one specific model: a Pacific coastal migration, dependent on high marine productivity, from the Bering Strait to South America, thousands of years before Clovis, the earliest widespread cultural manifestation south of the glacial ice. Speculations on maritime adaptations and typological links (stemmed points) across thousands of kilometers have also been advanced. A review of the current genetic, archeological, and paleoecological evidence indicates that ancestral Native American population expansion occurred after 16,000 years ago, consistent with the archeological record, particularly with the earliest securely dated sites after ~15,000 years ago. These data are largely consistent with either an inland (ice-free corridor) or Pacific coastal routes (or both), but neither can be rejected at present. Systematic archeological and paleoecological investigations, informed by geomorphology, are required to test each hypothesis.

Dental enamel hypoplasias in Holocene bighorn sheep (Oviscanadensis) in eastern Washington state, USA

Dental enamel hypoplasias have been documented in extant and fossil mammal species and attributed to several kinds of physiological stress. They have not previously been reported among bighorn sheep (Ovis canadensis Shaw, 1804). Forty-six (36.8%) of 125 mandibular molars (m1, m2, m3) of bighorn recovered from disturbed Holocene archaeological deposits in eastern Washington state display several kinds of hypoplasias. The exact ontogenetic age of the individual animals when hypoplasias formed cannot be determined. The majority of the hypoplasias occur near the root–enamel junction of the m3, suggesting that most individuals were young adults when the defect formed. Physiological stress associated with reproductive costs, winter nutritional deficits, or both seems likely.

Cattle phenotypes can disguise their maternal ancestry

Background

Cattle are bred for, amongst other factors, specific traits, including parasite resistance and adaptation to climate. However, the influence and inheritance of mitochondrial DNA (mtDNA) are not usually considered in breeding programmes. In this study, we analysed the mtDNA profiles of cattle from Victoria (VIC), southern Australia, which is a temperate climate, and the Northern Territory (NT), the northern part of Australia, which has a tropical climate, to determine if the mtDNA profiles of these cattle are indicative of breed and phenotype, and whether these profiles are appropriate for their environments.

Results

A phylogenetic tree of the full mtDNA sequences of different breeds of cattle, which were obtained from the NCBI database, showed that the mtDNA profiles of cattle do not always reflect their phenotype as some cattle with Bos taurus phenotypes had Bos indicus mtDNA, whilst some cattle with Bos indicus phenotypes had Bos taurus mtDNA. Using D-loop sequencing, we were able to contrast the phenotypes and mtDNA profiles from different species of cattle from the 2 distinct cattle breeding regions of Australia. We found that 67 of the 121 cattle with Bos indicus phenotypes from NT (55.4%) had Bos taurus mtDNA. In VIC, 92 of the 225 cattle with Bos taurus phenotypes (40.9%) possessed Bos indicus mtDNA. When focusing on oocytes from cattle with the Bos taurus phenotype in VIC, their respective oocytes with Bos indicus mtDNA had significantly lower levels of mtDNA copy number compared with oocytes possessing Bos taurus mtDNA (P < 0.01). However, embryos derived from oocytes with Bos indicus mtDNA had the same ability to develop to the blastocyst stage and the levels of mtDNA copy number in their blastocysts were similar to blastocysts derived from oocytes harbouring Bos taurus mtDNA. Nevertheless, oocytes originating from the Bos indicus phenotype exhibited lower developmental potential due to low mtDNA copy number when compared with oocytes from cattle with a Bos taurus phenotype.

Conclusions

The phenotype of cattle is not always related to their mtDNA profiles. MtDNA profiles should be considered for breeding programmes as they also influence phenotypic traits and reproductive capacity in terms of oocyte quality.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0523-5) contains supplementary material, which is available to authorized users.

Elevation Shift in Abies Mill. (Pinaceae) of Subtropical and Temperate China and Vietnam-Corroborative Evidence from Cytoplasmic DNA and Ecological Niche Modeling

The "elevational shift" scenario has been proposed as a model to explain the response of cold-adapted organisms to Quaternary climatic fluctuations in Europe and North America. However, the elevational shift model has not been well-explored in eastern Asia, which is more topographically complex than the other Northern Hemisphere biogeographic regions. Here, we evaluated the role of elevational shift in the closely related firs, or Abies Mill., of subtropical and temperate China. These firs are typical alpine trees with sensitivity to climate change. We tested the elevational shift hypothesis in firs of China using phylogeographic methods and ecological niche models. Our phylogeographic analyses comprised mitochondrial and chloroplast polymorphisms surveyed across 479 individuals from 43 populations representing 11 species. M1 of the 11 mitotypes and C1 of the 25 chlorotypes were inferred as the ancestral haplotype, and they had the widest distribution. The results of our phylogeographic survey revealed multiple centers of genetic diversity in distinct geographic regions and no latitudinal trend. Moreover, our results showed range expansions for seven taxa during the last glacial (64.9-18.2 or 32.5-9.1 kya), and this was consistent with the Quaternary fossil record of Abies in China. Taken together, our findings support a historical biogeographic pattern in firs of glacial expansions, probably through corridors at lower elevation, and interglacial fragmentations, through isolation at higher elevation peaks. Therefore, Abies in China probably undergoes elevational shift in response to climate change. Facing the forecasting global warming, the risk of several critically endangered firs was further enhanced as these species would have little escape space in situ to higher altitudes. According to our ENMs, we proposed an ex situ conservation strategy in the southern Hengduan Mountains region of south western China.

Phylogeographic structure in long-tailed voles (Rodentia: Arvicolinae) belies the complex Pleistocene history of isolation, divergence, and recolonization of Northwest North America's fauna

Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole, Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.

Fine-scale assessment of genetic diversity of trembling aspen in northwestern North America

BACKGROUND

In North America, the last ice age is the most recent event with severe consequences on boreal species' ranges. Phylogeographic patterns of range expansion in trembling aspen (Populus tremuloides) suggested that Beringia is likely to be a refugium and the "ice-free corridor" in Alberta may represent a region where small populations persisted during the last glacial maximum (LGM). The purpose of this study was to ascertain whether the origins of trembling aspen in western North America are reflected in the patterns of neutral genetic diversity and population structure. A total of 28 sites were sampled covering the northwestern part of aspen's distribution, from Saskatchewan to Alaska. Twelve microsatellite markers were used to describe patterns of genetic diversity. The genetic structure of trembling aspen populations was assessed by using multivariate analyses, Mantel correlograms, neighbor-joining trees and Bayesian analysis.

RESULTS

Microsatellite markers revealed little to no neutral genetic structure of P. tremuloides populations in northwestern North America. Low differentiation among populations and small isolation by distance (IBD) were observed. The most probable number of clusters detected by STRUCTURE was K = 3 (∆K = 5.9). The individuals in the populations of the 3 clusters share a common gene pool and showed a high level of admixture. No evidence was found that either Beringia or the "ice-free corridor" were refugia. Highest allelic richness (AR) and lowest heterozygosity (H_o) were observed in Alberta foothills of the Rocky Mountains.

CONCLUSIONS

Contrary to our hypothesis, our results showed that microsatellite markers revealed little to no genetic structure in P. tremuloides populations. Consequently, no divergent populations were observed near supposed refugia. The lack of detectable refugia in Beringia and in the "ice-free corridor" was due to high levels of gene flow between trembling apsen populations. More favorable environmental conditions for sexual reproduction and successful trembling aspen seedling establishment may have contributed to increase allelic richness through recombination in populations from the Albertan foothills of the Rocky Mountains.

Genetic Diversity and Phylogenetic Evolution of Tibetan Sheep Based on mtDNA D-Loop Sequences

The molecular and population genetic evidence of the phylogenetic status of the Tibetan sheep (Ovis aries) is not well understood, and little is known about this species' genetic diversity. This knowledge gap is partly due to the difficulty of sample collection. This is the first work to address this question. Here, the genetic diversity and phylogenetic relationship of 636 individual Tibetan sheep from fifteen populations were assessed using 642 complete sequences of the mitochondrial DNA D-loop. Samples were collected from the Qinghai-Tibetan Plateau area in China, and reference data were obtained from the six reference breed sequences available in GenBank. The length of the sequences varied considerably, between 1031 and 1259 bp. The haplotype diversity and nucleotide diversity were 0.992±0.010 and 0.019±0.001, respectively. The average number of nucleotide differences was 19.635. The mean nucleotide composition of the 350 haplotypes was 32.961% A, 29.708% T, 22.892% C, 14.439% G, 62.669% A+T, and 37.331% G+C. Phylogenetic analysis showed that all four previously defined haplogroups (A, B, C, and D) were found in the 636 individuals of the fifteen Tibetan sheep populations but that only the D haplogroup was found in Linzhou sheep. Further, the clustering analysis divided the fifteen Tibetan sheep populations into at least two clusters. The estimation of the demographic parameters from the mismatch analyses showed that haplogroups A, B, and C had at least one demographic expansion in Tibetan sheep. These results contribute to the knowledge of Tibetan sheep populations and will help inform future conservation programs about the Tibetan sheep native to the Qinghai-Tibetan Plateau.

Phylogeographic and population genetic structure of bighorn sheep ( Ovis canadensis ) in North American deserts

Fossil data are ambiguous regarding the evolutionary origin of contemporary desert bighorn sheep ( Ovis canadensis subspecies). To address this uncertainty, we conducted phylogeographic and population genetic analyses on bighorn sheep subspecies found in southwestern North America. We analyzed 515 base pairs of mtDNA control region sequence and 39 microsatellites in 804 individuals from 58 locations. Phylogenetic analyses revealed 2 highly divergent clades concordant with Sierra Nevada ( O. c. sierrae ) and Rocky Mountain ( O. c. canadensis ) bighorn and showed that these 2 subspecies both diverged from desert bighorn prior to or during the Illinoian glaciation (~315-94 thousand years ago [kya]). Desert bighorn comprised several more recently diverged haplogroups concordant with the putative Nelson ( O. c. nelsoni ), Mexican ( O. c. mexicana ), and Peninsular ( O. c. cremnobates ) subspecies. Corresponding estimates of effective splitting times (~17-3 kya), and haplogroup ages (~85-72 kya) placed the most likely timeframe for divergence among desert bighorn subspecies somewhere within the last glacial maximum. Median-joining haplotype network and Bayesian skyline analyses both indicated that desert bighorn collectively comprised a historically large and haplotype-diverse population, which subsequently lost much of its diversity through demographic decline. Using microsatellite data, discriminant analysis of principle components (DAPC) and Bayesian clustering analyses both indicated genetic structure concordant with the geographic distribution of 3 desert subspecies. Likewise, microsatellite and mitochondrial-based FST comparisons revealed significant fixation indices among the desert bighorn genetic clusters. We conclude these desert subspecies represent ancient lineages likely descended from separate Pleistocene refugial populations and should therefore be managed as distinct taxa to preserve maximal biodiversity. Los datos de fósiles sobre el origen evolutivo de las ovejas del desierto ( Ovis canadensis subespecies) contemporáneas son ambiguos. Para dilucidar esta incertidumbre, llevamos a cabo análisis filogeográficos y de genética de poblaciones entre cinco subespecies de ovejas del suroccidente de Norteamérica. Analizamos 515 pb de secuencia de la región control del ADN mitocondrial y 39 microsatélites en 804 ovejas de 58 localidades. Los análisis filogenéticos revelaron 2 clados altamente divergentes concordantes con ovejas de la Sierra Nevada ( O. c. sierrae ) y de las Montañas Rocosas ( O. c. canadensis ), y demostraron que estas dos subespecies divergieron antes o durante la glaciación de Illinois (315,000-94,000 años). Las ovejas del desierto formaron varios haplogrupos recientemente derivados concordantes con las subespecies de Nelson ( O. c. nelsoni ), México ( O. c. mexicana ) y peninsular ( O. c. cremnobates ). Las estimaciones correspondientes al tiempo de separación efectiva (17,000-3,000 años) y edades de haplogrupos (85,000-72,000 años) son los plazos más probables para las divergencias entre subespecies de ovejas del desierto dentro de la última glaciación máxima. Análisis de redes de haplotipos de unión de medias y análisis bayesianos de líneas de horizonte indicaron que las ovejas del desierto formaron una población históricamente grande y diversa en términos de haplotipos, que luego perdieron gran parte de su diversidad a través de un descenso demográfico. Utilizando datos de microsatélites los análisis DAPC y TESS indicaron agrupamiento genético concordante con la distribución geográfica actual de las tres subespecies. Asimismo, comparaciones de FST con datos de microsatélites y mitocondriales revelaron índices de fijación significativos entre los grupos genéticos de ovejas del desierto. Concluimos que estas subespecies de ovejas del desierto representan linajes antiguos que probablemente descienden de poblaciones de distintos refugios del Pleistoceno, y que por lo tanto deben ser manejadas como taxones distintos para preservar su biodiversidad máxima.

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.

Pleistocene glaciations, demographic expansion and subsequent isolation promoted morphological heterogeneity: A phylogeographic study of the alpine Rosa sericea complex (Rosaceae)

While most temperate plants probably underwent glacial constriction to refugia and interglacial expansion, another type of interglacial refugia might have existed to maintain alpine plants during warm periods. To test this hypothesis, we applied phylogeographic methods to 763 individuals (62 populations) which belong to 7 taxonomically difficult species of the Rosa sericea complex distributed in alpine regions of the temperate and subtropical zones in eastern Asia. We used three chloroplast (cp) DNA fragments (trnL-trnF, ndhF-rpl32 and ndhJ-trnF) approximately 3,100 bp and nuclear microsatellite (nSSR) on eight sites to determine whether cold tolerant plants experienced expansion during the Pleistocene. The neutral test and mismatch distribution analysis (MDA) indicated that whole populations and major lineages of the Qinghai-Tibet Plateau (QTP) underwent expansion during the middle to late Pleistocene. Environmental niche modeling (ENM) indicates more suitable habitats during the Last Glacial Maximum (LGM) than at present. We concluded that the demographic history of R. sericea, which diverged in the middle Pleistocene, was mostly affected by climatic oscillations instead of by geographical barriers. The low genetic divergence, as well as the weak phylogenetic structure in the R. sericea complex both support treating this complex as a single taxon.

Phylogeography and post-glacial recolonization in wolverines (Gulo gulo) from across their circumpolar distribution

Interglacial-glacial cycles of the Quaternary are widely recognized in shaping phylogeographic structure. Patterns from cold adapted species can be especially informative - in particular, uncovering additional glacial refugia, identifying likely recolonization patterns, and increasing our understanding of species' responses to climate change. We investigated phylogenetic structure of the wolverine, a wide-ranging cold adapted carnivore, using a 318 bp of the mitochondrial DNA control region for 983 wolverines (n=209 this study, n=774 from GenBank) from across their full Holarctic distribution. Bayesian phylogenetic tree reconstruction and the distribution of observed pairwise haplotype differences (mismatch distribution) provided evidence of a single rapid population expansion across the wolverine's Holarctic range. Even though molecular evidence corroborated a single refugium, significant subdivisions of population genetic structure (0.01< ΦST <0.99, P<0.05) were detected. Pairwise ΦST estimates separated Scandinavia from Russia and Mongolia, and identified five main divisions within North America - the Central Arctic, a western region, an eastern region consisting of Ontario and Quebec/Labrador, Manitoba, and California. These data are in contrast to the nearly panmictic structure observed in northwestern North America using nuclear microsatellites, but largely support the nuclear DNA separation of contemporary Manitoba and Ontario wolverines from northern populations. Historic samples (c. 1900) from the functionally extirpated eastern population of Quebec/Labrador displayed genetic similarities to contemporary Ontario wolverines. To understand these divergence patterns, four hypotheses were tested using Approximate Bayesian Computation (ABC). The most supported hypothesis was a single Beringia incursion during the last glacial maximum that established the northwestern population, followed by a west-to-east colonization during the Holocene. This pattern is suggestive of colonization occurring in accordance with glacial retreat, and supports expansion from a single refugium. These data are significant relative to current discussions on the conservation status of this species across its range.

Genetic variation and phylogenetic relationship between three serow species of the genus Capricornis based on the complete mitochondrial DNA control region sequences

The molecular evidence of phylogenetic status regarding the Formosan serow (Capricornis swinhoei) is not robust and little is known about the genetic diversity of the Sumatran serow (Capricornis sumatraensis), which partly is due to the hardness in sample collection. Here we determined the sequences of the complete mitochondrial DNA control region (1,014 bp) of 19 Sumatran-serow individuals. Nine new haplotypes were defined based on 78 variable sites. Combined analysis with other 32 haplotypes downloaded from the public database, including 1 Sumatran-serow, 11 Formosan-serow and 20 Japanese-serow (Capricornis crispus) haplotypes, a relatively high level of nucleotide diversity was first observed in Sumatran serow (π = 0.0249). By comparative analysis with structural consensus sequences from other mammals, we have identified central, left and right domains and depicted the putative functional structure, including extend termination associated sequences and conserve sequence blocks, in mtDNA control region. The alignment of mtDNA control region revealed that both Sumatran and Japanese serow have two tandem repeats (TRs), but three TRs in Formosan serow. Phylogenetic analyses revealed that the Formosan serow is distinct species with the Japanese serow, but a sister group with the Sumatran serow. The divergence time estimated among three serow species revealed that Pleistocene climate changes and the uplift of Qinghai-Tibetan plateau might play an important role in the genetic differentiation of the serows. These results mainly provide the convinced evidence on the genetic relationship between three serow species.

Phylogenetic relationships of three Italian merino-derived sheep breeds evaluated through a complete mitogenome analysis

In Italy, the crisis of the wool industry triggered the necessity to reconvert the two traditional Merino-derived breeds, Gentile di Puglia and Sopravissana, to meat production, by creating the Merinizzata Italiana. The aim of the present study was to assess the genetic diversity of these three Italian Merino-derived (IMd) breeds by examining the molecular information encoded in the maternally-inherited mitochondrial DNA (mtDNA). A parallel molecular investigation was performed on the putative paternal and maternal breeds, the Merino from Spain and the Appenninica from Italy, respectively, as well as on three unrelated dairy breeds (Sarda and Comisana from Italy, and Lacaune from France). Firstly, the mtDNA control region of 291 samples was analyzed. When comparing the overall genetic distances among the eight stocks, the three IMd breeds clustered together close to the Appenninica, thus confirming its parental role. Among the 90 IMd samples, 82 different haplotypes were observed, almost all belonging to haplogroup B, and only one to A. For 23 mtDNAs, including nine IMd, the analysis was then brought to the level of entire mitogenomes. Three distinct sub-haplogroups within B were found to encompass the IMd samples, with one clade (B1a2a1) apparently restricted to those sheep. Thus, despite experiencing a drastic reduction in number (mainly due to changes in breeding practices driven by the economy), the IMd breeds still represent a reservoir of distinctive mitochondrial variants, which could potentially contribute to the development of conservation and management programs of Italian sheep breeds.

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.

Evaluation of experimental genetic management in reintroduced bighorn sheep

Positive demographic responses have been reported in several species where the immigration or supplementation of genetically distinct individuals into wild populations has resulted in a genetic rescue effect. However, rarely have researchers incorporated what could be considerable risk of outbreeding depression into planning for genetic management programs. We assess the genetic effects of an experiment in genetic management involving replicate populations of California bighorn sheep (Ovis canadensis californiana) in Oregon, USA, which previously experienced poor productivity and numerical declines. In the experiment, two declining populations were supplemented with ewes from a more genetically diverse population of California bighorn sheep in Nevada. We incorporated analysis of genetic samples representing both experimental populations prior to supplementation, samples from the supplemented individuals, and samples collected from both experimental populations approximately one generation after supplementation. We used genetic analyses to assess the integration of supplemented and resident populations by identifying interpopulation hybrids. Further, we incorporated demographic simulations to assess the risk of outbreeding depression as a result of the experimental augmentation. Finally, we used data from microsatellites and mitochondrial sequences to determine if genetic management increased genetic diversity in the experimental populations. Our analyses demonstrated the success of genetic management by documenting interpopulation hybrids, identifying no evidence for outbreeding depression as a result of contact between the genetically distinct supplemented and resident populations, and by identifying increased population-level metrics of genetic diversity in postsupplementation populations compared with presupplementation levels.

Hybridization among Arctic white-headed gulls (Larus spp.) obscures the genetic legacy of the Pleistocene

We studied the influence of glacial oscillations on the genetic structure of seven species of white-headed gull that breed at high latitudes (Larus argentatus, L. canus, L. glaucescens, L. glaucoides, L. hyperboreus, L. schistisagus, and L. thayeri). We evaluated localities hypothesized as ice-free areas or glacial refugia in other Arctic vertebrates using molecular data from 11 microsatellite loci, mitochondrial DNA (mtDNA) control region, and six nuclear introns for 32 populations across the Holarctic. Moderate levels of genetic structure were observed for microsatellites (F_ST= 0.129), introns (Φ_ST= 0.185), and mtDNA control region (Φ_ST= 0.461), with among-group variation maximized when populations were grouped based on subspecific classification. Two haplotype and at least two allele groups were observed across all loci. However, no haplotype/allele group was composed solely of individuals of a single species, a pattern consistent with recent divergence. Furthermore, northernmost populations were not well differentiated and among-group variation was maximized when L. argentatus and L. hyberboreus populations were grouped by locality rather than species, indicating recent hybridization. Four populations are located in putative Pleistocene glacial refugia and had larger τ estimates than the other 28 populations. However, we were unable to substantiate these putative refugia using coalescent theory, as all populations had genetic signatures of stability based on mtDNA. The extent of haplotype and allele sharing among Arctic white-headed gull species is noteworthy. Studies of other Arctic taxa have generally revealed species-specific clusters as well as genetic structure within species, usually correlated with geography. Aspects of white-headed gull behavioral biology, such as colonization ability and propensity to hybridize, as well as their recent evolutionary history, have likely played a large role in the limited genetic structure observed.

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.

Inheritance of mitochondrial DNA in serially recloned pigs by somatic cell nuclear transfer (SCNT)

Somatic cell nuclear transfer (SCNT) has been established for the transmission of specific nuclear DNA. However, the fate of donor mitochondrial DNA (mtDNA) remains unclear. Here, we examined the fate of donor mtDNA in recloned pigs through third generations. Fibroblasts of recloned pigs were obtained from offspring of each generation produced by fusion of cultured fibroblasts from a Minnesota miniature pig (MMP) into enucleated oocytes of a Landrace pig. The D-loop regions from the mtDNA of donor and recipient differ at nucleotide sequence positions 16050 (A→T), 16062 (T→C), and 16135 (G→A). In order to determine the fate of donor mtDNA in recloned pigs, we analyzed the D-loop region of the donor's mtDNA by allele-specific PCR (AS-PCR) and real-time PCR. Donor mtDNA was successfully detected in all recloned offspring (F1, F2, and F3). These results indicate that heteroplasmy that originate from donor and recipient mtDNA is maintained in recloned pigs, resulting from SCNT, unlike natural reproduction.

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.

Large scale mitochondrial sequencing in Mexican Americans suggests a reappraisal of Native American origins

Background

The Asian origin of Native Americans is largely accepted. However uncertainties persist regarding the source population(s) within Asia, the divergence and arrival time(s) of the founder groups, the number of expansion events, and migration routes into the New World. mtDNA data, presented over the past two decades, have been used to suggest a single-migration model for which the Beringian land mass plays an important role.

Results

In our analysis of 568 mitochondrial genomes, the coalescent age estimates of shared roots between Native American and Siberian-Asian lineages, calculated using two different mutation rates, are A4 (27.5 ± 6.8 kya/22.7 ± 7.4 kya), C1 (21.4 ± 2.7 kya/16.4 ± 1.5 kya), C4 (21.0 ± 4.6 kya/20.0 ± 6.4 kya), and D4e1 (24.1 ± 9.0 kya/17.9 ± 10.0 kya). The coalescent age estimates of pan-American haplogroups calculated using the same two mutation rates (A2:19.5 ± 1.3 kya/16.1 ± 1.5 kya, B2:20.8 ± 2.0 kya/18.1 ± 2.4 kya, C1:21.4 ± 2.7 kya/16.4 ± 1.5 kya and D1:17.2 ± 2.0 kya/14.9 ± 2.2 kya) and estimates of population expansions within America (~21-16 kya), support the pre-Clovis occupation of the New World. The phylogeography of sublineages within American haplogroups A2, B2, D1 and the C1b, C1c andC1d subhaplogroups of C1 are complex and largely specific to geographical North, Central and South America. However some sub-branches (B2b, C1b, C1c, C1d and D1f) already existed in American founder haplogroups before expansion into the America.

Conclusions

Our results suggest that Native American founders diverged from their Siberian-Asian progenitors sometime during the last glacial maximum (LGM) and expanded into America soon after the LGM peak (~20-16 kya). The phylogeography of haplogroup C1 suggest that this American founder haplogroup differentiated in Siberia-Asia. The situation is less clear for haplogroup B2, however haplogroups A2 and D1 may have differentiated soon after the Native American founders divergence. A moderate population bottle neck in American founder populations just before the expansion most plausibly resulted in few founder types in America. The similar estimates of the diversity indices and Bayesian skyline analysis in North America, Central America and South America suggest almost simultaneous (~ 2.0 ky from South to North America) colonization of these geographical regions with rapid population expansion differentiating into more or less regional branches across the pan-American haplogroups.

Diversification at high latitudes: speciation of buntings in the genus Plectrophenax inferred from mitochondrial and nuclear markers

High-latitude diversification is a process characterized by speciation and extinction due to climatically driven vicariance and dispersal events. McKay's buntings (Plectrophenax hyperboreus) are high-latitude island endemic songbirds, and their global range is restricted to Beringia. Snow buntings (P. nivalis), their closest relatives, are distributed throughout the Holarctic, breeding in available habitat surrounding the island range of McKay's buntings. We sequenced 1123 base pairs of mitochondrial DNA for 40 individuals of each species and analysed a total of 913 AFLPs for 57 individuals. Both marker types suggested weak but significant genetic differentiation. Analysis of sequence data indicated divergence occurring when the current breeding range of McKay's buntings was a hill on the Beringian steppe (approximately 18,400 to approximately 73,700 years before present), suggesting that snow buntings were restricted to lower latitudes by ice sheets. Ancestral effective population size estimates indicate a founder event in McKay's buntings followed by an expansion and then a reduction in effective size. Rising sea levels and asymmetric hybridization from McKay's buntings into the postglacially-colonizing population of snow buntings could account for this reduction. Reproductive isolation is likely maintained through differential arrival dates on breeding grounds and the high breeding density of McKay's buntings. This recent, high-latitude divergence best fits a model of founder event speciation driven by vicariance and oscillations in habitat due to climate change.

Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae)

New insights for the systematic and evolution of the wild sheep are provided by molecular phylogenies inferred from Maximum parsimony, Bayesian, Maximum likelihood, and Neighbor-Joining methods. The phylogeny of the wild sheep was based on cytochrome b sequences of 290 samples representative of most of the sub-species described in the genus Ovis. The result was confirmed by a combined tree based on cytochrome b and nuclear sequences for 79 Ovis samples representative of the robust clades established with mitochondrial data. Urial and mouflon, which are either considered as a single or two separate species, form two monophyletic groups (O. orientalis and O. vignei). Their hybrids appear in one or the other group, independently from their geographic origin. The European mouflon O. musimon is clearly in the O. orientalis clade. The others species, O. dalli, O. canadensis, O. nivicola, and O. ammon are monophyletic. The results support an Asiatic origin of the genus Ovis, followed by a migration to North America through North-Eastern Asia and the Bering Strait and a diversification of the genus in Eurasia less than 3 million years ago. Our results show that the evolution of the genus Ovis is a striking example of successive speciation events occurring along the migration routes propagating from the ancestral area.

Refugia revisited: individualistic responses of species in space and time

Climate change in the past has led to significant changes in species' distributions. However, how individual species respond to climate change depends largely on their adaptations and environmental tolerances. In the Quaternary, temperate-adapted taxa are in general confined to refugia during glacials while cold-adapted taxa are in refugia during interglacials. In the Northern Hemisphere, evidence appears to be mounting that in addition to traditional southern refugia for temperate species, cryptic refugia existed in the North during glacials. Equivalent cryptic southern refugia, to the south of the more conventional high-latitude polar refugia, exist in montane areas during periods of warm climate, such as the current interglacial. There is also a continental/oceanic longitudinal gradient, which should be included in a more complete consideration of the interaction between species ranges and climates. Overall, it seems clear that there is large variation in both the size of refugia and the duration during which species are confined to them. This has implications for the role of refugia in the evolution of species and their genetic diversity.

The genetic impact of translocations and habitat fragmentation in chamois (Rupicapra) spp

The chamois is a useful species with which to investigate the combined genetic impact of habitat fragmentation, over hunting, and translocations. Genetic variation within and between chamois (genus Rupicapra) populations was analyzed in 259 individuals from 16 sampling sites located in Italy, Spain, Slovakia, and the Czech Republic. Two mitochondrial DNA markers (control region and cytochrome b) and 11 nuclear microsatellites were typed. The principal results of this study can be summarized as follows: 1) high and significant differentiation between almost all chamois populations is observed even on a microgeographical scale, probably caused by the patchy distribution of this species, sharp geographical barriers to gene flow, and drift effects related to recent bottlenecks; 2) historical translocation events have left a clear genetic signature, including interspecific hybridization in some Alpine localities; 3) the Apennine subspecies of chamois, Rupicapra pyrenaica ornata, shows a high and similar level of divergence (about 1.5 My) from the Pyrenean (Rupicapra pyrenaica pyrenaica) and the Alpine (Rupicapra rupicapra) chamois; therefore, the specific status of these taxa should be revised. These results confirm the potential of population genetic analyses to dissect and interpret complex patterns of diversity in order to define factors important to conservation and management.

Historic hybridization and persistence of a novel mito-nuclear combination in red-backed voles (genus Myodes)

BACKGROUND

The role of hybridization in generating diversity in animals is an active area of discovery and debate. We assess hybridization across a contact zone of northern (Myodes rutilus) and southern (M. gapperi) red-backed voles using variation in skeletal features and both mitochondrial and nuclear loci. This transect extends approximately 550 km along the North Pacific Coast of North America and encompasses 26 populations (n = 485). We establish the history, geographic extent and directionality of hybridization, determine whether hybridization is ongoing, and assess the evolutionary stability of novel genomic combinations.

RESULTS

Identification of M. rutilus and M. gapperi based on the degree of closure of the post-palatal bridge was concordant with the distribution of diagnostic nuclear MYH6 alleles; however, an 80 km zone of introgressed populations was identified. The introgressant form is characterized by having mitochondrial haplotypes closely related to the northern M. rutilus on a nuclear background and morphological characteristics of southern M. gapperi.

CONCLUSION

Introgression appears to have been historic as pure populations of M. rutilus are now isolated to the north from introgressants or pure M. gapperi by the LeConte Glacier. As we do not find pure M. rutilus or M. gapperi individuals throughout the distribution of the introgressant form, it appears that the introgressants are a self-sustaining entity not requiring continued hybridization between pure parental forms to generate this novel combination of characters.

Distinctive Paleo-Indian migration routes from Beringia marked by two rare mtDNA haplogroups

BACKGROUND

It is widely accepted that the ancestors of Native Americans arrived in the New World via Beringia approximately 10 to 30 thousand years ago (kya). However, the arrival time(s), number of expansion events, and migration routes into the Western Hemisphere remain controversial because linguistic, archaeological, and genetic evidence have not yet provided coherent answers. Notably, most of the genetic evidence has been acquired from the analysis of the common pan-American mitochondrial DNA (mtDNA) haplogroups. In this study, we have instead identified and analyzed mtDNAs belonging to two rare Native American haplogroups named D4h3 and X2a.

RESULTS

Phylogeographic analyses at the highest level of molecular resolution (69 entire mitochondrial genomes) reveal that two almost concomitant paths of migration from Beringia led to the Paleo-Indian dispersal approximately 15-17 kya. Haplogroup D4h3 spread into the Americas along the Pacific coast, whereas X2a entered through the ice-free corridor between the Laurentide and Cordilleran ice sheets. The examination of an additional 276 entire mtDNA sequences provides similar entry times for all common Native American haplogroups, thus indicating at least a dual origin for Paleo- Indians.

CONCLUSIONS

A dual origin for the first Americans is a striking novelty from the genetic point of view, and it makes plausible a scenario positing that within a rather short period of time, there may have been several entries into the Americas from a dynamically changing Beringian source. Moreover, this implies that most probably more than one language family was carried along with the Paleo-Indians.

Species-wide phylogeography of North American mule deer (Odocoileus hemionus): cryptic glacial refugia and postglacial recolonization

Quaternary climatic oscillations greatly influenced the present-day population genetic structure of animals and plants. For species with high dispersal and reproductive potential, phylogeographic patterns resulting from historical processes can be cryptic, overshadowed by contemporary processes. Here we report a study of the phylogeography of Odocoileus hemionus, a large, vagile ungulate common throughout western North America. We examined sequence variation of mitochondrial DNA (control region and cytochrome b) within and among 70 natural populations across the entire range of the species. Among the 1766 individual animals surveyed, we recovered 496 haplotypes. Although fine-scale phylogenetic structure was weakly resolved using phylogenetic methods, network analysis clearly revealed the presence of 12 distinct haplogroups. The spatial distribution of haplogroups showed a strong genetic discontinuity between the two morphological types of O. hemionus, mule deer and black-tailed deer, east and west of the Cascade Mountains in the Pacific Northwest. Within the mule deer lineage, we identified several haplogroups that expanded before or during the Last Glacial Maximum, suggesting that mule deer persisted in multiple refugia south of the ice sheets. Patterns of genetic diversity within the black-tailed deer lineage suggest a single refugium along the Pacific Northwest coast, and refute the hypothesis that black-tailed deer persisted in one or more northern refugia. Our data suggest that black-tailed deer recolonized areas in accordance with the pattern of glacial retreat, with initial recolonization northward along a coastal route and secondary recolonization inland.

Bilateral perirenal hemorrhage in two Stone's sheep (Ovis dalli stonei): a possible manifestation of malignant catarrhal fever

Malignant catarrhal fever (MCF)-like disease was diagnosed at postmortem in 2 Stone's sheep (Ovis dalli stonei). On gross examination, the predominant abnormality in both sheep was severe perirenal hemorrhage and multiple renal infarcts. Microscopically, there was severe, multisystemic lymphocytic arteritis. Both sheep were positive for Ovine herpesvirus 2 (OvHV-2) on polymerase chain reaction, and partial sequencing of the viral DNA polymerase genes from the 2 sheep revealed >99% homology, with 96% similarity to the reference GenBank OvHV-2 viral sequence. Based on the histological lesions, polymerase chain reaction results, and viral DNA polymerase gene sequencing, a diagnosis of OvHV-2-mediated MCF was made. Massive perirenal hemorrhage has not been described previously as a manifestation of MCF.