Late‐glacial recolonization and phylogeography of European red deer (Cervus elaphus L.)

The Pleistocene was an epoch of extreme climatic and environmental changes. How individual species responded to the repeated cycles of warm and cold stages is a major topic of debate. For the European fauna and flora, an expansion–contraction model has been suggested, whereby temperate species were restricted to southern refugia during glacial times and expanded northwards during interglacials, including the present interglacial (Holocene). Here, we test this model on the red deer (Cervus elaphus) a large and highly mobile herbivore, using both modern and ancient mitochondrial DNA from the entire European range of the species over the last c. 40 000 years. Our results indicate that this species was sensitive to the effects of climate change. Prior to the Last Glacial Maximum (LGM) haplogroups restricted today to South‐East Europe and Western Asia reached as far west as the UK. During the LGM, red deer was mainly restricted to southern refugia, in Iberia, the Balkans and possibly in Italy and South‐Western Asia. At the end of the LGM, red deer expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, Scandinavia, Germany, Poland and Belarus. Ancient DNA data cannot rule out refugial survival of red deer in North‐West Europe through the LGM. Had such deer survived, though, they were replaced by deer migrating from Iberia at the end of the glacial. The Balkans served as a separate LGM refugium and were probably connected to Western Asia with genetic exchange between the two areas.     

Ancient DNA supports southern survival of Richardson's collared lemming (Dicrostonyx richardsoni) during the last glacial maximum

Collared lemmings (genus Dicrostonyx) are circumpolar Arctic arvicoline rodents associated with tundra. However, during the last glacial maximum (LGM), Dicrostonyx lived along the southern ice margin of the Laurentide ice sheet in communities comprising both temperate and boreal species. To better understand these communities and the fate of these southern individuals, we compare mitochondrial cytochrome b sequence data from three LGM‐age Dicrostonyx fossils from south of the Laurentide ice sheet to sequences from modern Dicrostonyx sampled from across their present‐day range. We test whether the Dicrostonyx populations from LGM‐age continental USA became extinct at the Pleistocene–Holocene transition ~11000 years ago or, alternatively, if they belong to an extant species whose habitat preferences can be used to infer the palaeoclimate along the glacial margin. Our results indicate that LGM‐age Dicrostonyx from Iowa and South Dakota belong to Dicrostonyx richardsoni, which currently lives in a temperate tundra environment west of Hudson Bay, Canada. This suggests a palaeoclimate south of the Laurentide ice sheet that contains elements similar to the more temperate shrub tundra characteristic of extant D. richardsoni habitat, rather than the very cold, dry tundra of the Northern Arctic. While more data are required to determine whether or not the LGM southern population is ancestral to extant D. richardsoni, it seems most probable that the species survived the LGM in a southern refugium.     

Colonization of the Scottish islands via long-distance Neolithic transport of red deer (Cervus elaphus)

Red deer (Cervus elaphus) have played a key role in human societies throughout history, with important cultural significance and as a source of food and materials. This relationship can be traced back to the earliest human cultures and continues to the present day. Humans are thought to be responsible for the movement of a considerable number of deer throughout history, although the majority of these movements are poorly described or understood. Studying such translocations allows us to better understand ancient human–wildlife interactions, and in the case of island colonizations, informs us about ancient human maritime practices. This study uses DNA sequences to characterise red deer genetic diversity across the Scottish islands (Inner and Outer Hebrides and Orkney) and mainland using ancient deer samples, and attempts to infer historical colonization events. We show that deer from the Outer Hebrides and Orkney are unlikely to have originated from mainland Scotland, implying that humans introduced red deer from a greater distance. Our results are also inconsistent with an origin from Ireland or Norway, suggesting long-distance maritime travel by Neolithic people to the outer Scottish Isles from an unknown source. Common haplotypes and low genetic differentiation between the Outer Hebrides and Orkney imply common ancestry and/or gene flow across these islands. Close genetic proximity between the Inner Hebrides and Ireland, however, corroborates previous studies identifying mainland Britain as a source for red deer introductions into Ireland. This study provides important information on the processes that led to the current distribution of the largest surviving indigenous land mammal in the British Isles.     

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

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

Validation of a single‐nucleotide polymorphism panel for parentage testing of farmed red deer

Red deer behaviour prevents the accurate physical matching of dams to fawns and, therefore, breeders have to rely on DNA‐based parentage testing. A panel consisting of 100 single nucleotide polymorphism markers, with an average minor allele frequency of 0.25, was able to resolve 92% of fawns to both parents. In comparison, an existing 12‐marker microsatellite panel was able to resolve 68% of fawns to both parents. When excluding dam DNA information, the single nucleotide polymorphism panel matched 81% of the fawns to their sires and the microsatellite panel 71%.     

The Tien Shan vole (Microtus ilaeus; Rodentia: Cricetidae) as a new species in the Late Pleistocene of Europe

Grey voles (subgenus Microtus) represent a complex of at least seven closely related and partly cryptic species. The range of these species extends from the Atlantic to the Altai Mountains, but most of them occur east of the Black Sea. Using ancient DNA analyses of the Late Pleistocene specimens, we identified a new mtDNA lineage of grey voles in Europe. Phylogenetic analysis of mitochondrial DNA cytochrome b sequences from 23 voles from three caves, namely, Emine–Bair–Khosar (Crimea, Ukraine), Cave 16 (Bulgaria), and Bacho Kiro (Bulgaria), showed that 14 specimens form a previously unrecognized lineage, sister to the Tien Shan vole. The average sequence divergence of this lineage and the extant Tien Shan vole was 4.8%, which is similar to the divergence of grey vole forms, which are considered distinct species or being on the verge of speciation; M. arvalis and M. obscurus or M. mystacinus and M. rossiaemeridionalis. We estimated the time to the most recent common ancestor of the grey voles to be 0.66 Ma, which is over twice the recent estimates, while the divergence of the extant Tien Shan vole and the new lineage to be 0.29 Ma. Our discovery suggests that grey voles may have been more diversified in the past and that their ranges may have differed substantially from current ones. It also underlines the utility of ancient DNA to decipher the evolutionary history of voles.     

Evolution of the Cyprinella lutrensis species group. III. Geographic variation in the mitochondrial DNA of Cyprinella lutrensis— the influence of Pleistocene glaciation on population dispersal and divergence

We employed restriction site variation in mitochondrial (mt)DNA to determine if significant phylogeographic structure occurs in the North American cyprinid fish Cyprinella lutrensis. Digestion patterns from 16 restriction endonucleases identified fifty mtDNA haplotypes among 127 individuals of Cyprinella lutrensis assayed from localities in the Gulf Coastal Plain, the Great Plains, and the Central Lowlands. Nucleotide sequence divergence among haplotypes was highly variable (mean ± SE: 2.87%± 0.08; range: 0.14–9.24%). Maximum-parsimony analysis and the neighbour joining method of tree construction revealed three major groupings (clades) of haplotypes that differed in geographic distribution. Divergence estimates between the basal clade, comprised of haplotypes primarily from the Brazos River in east Texas, and the remaining two clades, place C. lutrensis in the western Gulf Coastal Plain prior to Pleistocene glaciation. Nucleotide sequence divergence between the second clade, comprised of haplotypes from the Trinity and Calcasieu rivers in east Texas and southwestern Louisiana, respectively, and the third clade (comprised primarily of haplotypes from localities north of Texas and affected directly by Pleistocene glaciation), suggest that C. lutrensis colonized gladated regions to the north during the mid- to late Pleistocene. This hypothesis is supported by levels of intrapopulational nucleotide diversity in geographic localities outside of Texas and by geological evidence. Despite marked geographic variation in morphometries, meristics, and nuptial coloration, mtDNA variation in glaciated regions was not geographically structured, and subspecies of C. lutrensis were not identifiable by phylogenetic analysis of mtDNA.     

Assessing the phylogeographic history of the montane caddisfly Thremma gallicum using mitochondrial and restriction‐site‐associated DNA (RAD) markers

Repeated Quaternary glaciations have significantly shaped the present distribution and diversity of several European species in aquatic and terrestrial habitats. To study the phylogeography of freshwater invertebrates, patterns of intraspecific variation have been examined primarily using mitochondrial DNA markers that may yield results unrepresentative of the true species history. Here, population genetic parameters were inferred for a montane aquatic caddisfly, Thremma gallicum, by sequencing a 658‐bp fragment of the mitochondrial CO1 gene, and 12,514 nuclear RAD loci. T. gallicum has a highly disjunct distribution in southern and central Europe, with known populations in the Cantabrian Mountains, Pyrenees, Massif Central, and Black Forest. Both datasets represented rangewide sampling of T. gallicum. For the CO1 dataset, this included 352 specimens from 26 populations, and for the RAD dataset, 17 specimens from eight populations. We tested 20 competing phylogeographic scenarios using approximate Bayesian computation (ABC) and estimated genetic diversity patterns. Support for phylogeographic scenarios and diversity estimates differed between datasets with the RAD data favouring a southern origin of extant populations and indicating the Cantabrian Mountains and Massif Central populations to represent highly diverse populations as compared with the Pyrenees and Black Forest populations. The CO1 data supported a vicariance scenario (north–south) and yielded inconsistent diversity estimates. Permutation tests suggest that a few hundred polymorphic RAD SNPs are necessary for reliable parameter estimates. Our results highlight the potential of RAD and ABC‐based hypothesis testing to complement phylogeographic studies on non‐model species.     

Post‐glacial history and introgression in Abies (Pinaceae) species of the Russian Far East inferred from both nuclear and cytoplasmic markers

Aim The main aim of the present study is to infer the post‐glacial history of Abies species from north‐east Asia and to test the hypotheses that coastal Abies populations suffered less from climatic fluctuations during Pleistocene glacial periods than their more continental counterparts, and that Sakhalin was a major area of introgression. Location Natural ranges of the fir species Abies nephrolepis, Abies sachalinensis and Abies holophylla in the Russian Far East, and of Abies gracilis, which is endemic to the Kamchatka Peninsula. Methods Nineteen populations were sampled for allozyme analysis. Seventeen of these populations were also screened for variation at two paternally inherited chloroplast DNA microsatellite loci (cpSSR) and variation at one maternally inherited mitochondrial marker (nad4‐3/4). Finally a subset of 11 populations was analysed with amplified fragment length polymorphism (AFLP). Comparisons were made with already available Abies sibirica data. For all sets of markers, we estimated genetic diversity and differentiation using an analysis of molecular variance (AMOVA). Population clustering was assessed with a Bayesian approach implemented in structure v.2.3. Results Among the three major species, A. sibirica, A. nephrolepis and A. sachalinensis, A. sachalinensis demonstrated the highest cytoplasmic and nuclear diversity and the most continental species, A. sibirica, the lowest. Both nuclear and mitochondrial DNA markers revealed the presence of a transitional zone on Sakhalin Island between A. nephrolepis and A. sachalinensis of south Sakhalin. The structure analysis delivered very clear results confirming the admixed origin of A. sachalinensis, with a genetic contribution from A. nephrolepis. No variation in cytoplasmic markers was found in A. gracilis, suggesting the occurrence of a recent bottleneck. Main conclusions There is a clear reduction of genetic diversity in Abies species from the Pacific coast into the continent. The higher diversity in A. sachalinensis could have two causes: a larger effective population size in the islands due to relatively stable climatic conditions and consequently less pronounced demographic fluctuations in population size and/or hybridization with continental and Japanese populations. Sakhalin Island is a major transitional zone for conifer species. Finally, the fir from Kamchatka, A. gracilis, should be regarded as a separate species closely related to the A. nephrolepis–A. sachalinensis complex.     

Brief communication: Ancient nuclear DNA and kinship analysis: the case of a medieval burial in San Esteban Church in Cuellar (Segovia, Central Spain)

The aim of this work was to investigate a very common situation in the archaeological and anthropological context: the study of a burial site containing several individuals, probably related genetically, using ancient DNA techniques. We used available ancient DNA and forensic protocols to obtain reliable results on archaeological material. The results also enabled molecular sex determination to be compared with osteological data. Specifically, a modified ancient DNA extraction method combined with the amplification of nuclear markers with the AmpFlSTR®MiniFiler? kit(Applied Biosystems) was used. Seven medieval individuals buried in four niches dated in the 15th Century at San Esteban Church in Cuellar (Segovia, Central Spain) were analyzed by the proposed method, and four of seven provided complete autosomal short tandem repeat (STRs) profiles. Kinship analyses comprising paternity and sibship relations were carried out with pedigree‐specific software used in forensic casework. A 99.98% paternity probability was established between two individuals, although lower percentages (68%) were obtained in other cases, and some hypothetical kinship relations were excluded. The overall results could eventually provide evidence for reconstructing the historical record. Am J Phys Anthropol, 2011. © 2010 Wiley‐Liss, Inc.     

Phylogeography of the bigeye chub Hybopsis amblops (Teleostei: Cypriniformes): early Pleistocene diversification and post‐glacial range expansion

The bigeye chub, Hybopsis amblops, is a member of the Central Highlands ichthyofauna of eastern North America. Phylogenetic analyses of the H. amblops species group based on a 1059 bp fragment of the mitochondrial DNA cytochrome b gene did not recover a monophyletic group. The inclusion of Hybopsis hypsinotus in the species complex is questionable. Within H. amblops, five strongly supported clades were identified; two clades containing haplotypes from the Ozark Highlands and three clades containing haplotypes from the Eastern Highlands and previously glaciated regions of the Ohio and Wabash River drainages. Estimates of the timing of divergence indicated that prior to the onset of glaciation, vicariant events separated populations east and west of the Mississippi River. East of the Mississippi River glacial cycles associated with the blocking and rerouting of the Teays River system caused populations to be pushed southward into refugia of the upper Ohio River. Following the most recent Wisconsinan glaciation, populations expanded northward into previously glaciated regions and southward into the Cumberland River drainage. In the Ozarks, west of the Mississippi River, isolation of clades appears to be maintained by the lack of stream capture events between the upper Arkansas and the White River systems and a barrier formed by the Arkansas River.     

Mitochondrial DNA as a tool for reconstructing past life‐history traits in mammals

Reconstructing the ancestral characteristics of species is a major goal in evolutionary and comparative biology. Unfortunately, fossils are not always available and sufficiently informative, and phylogenetic methods based on models of character evolution can be unsatisfactory. Genomic data offer a new opportunity to estimate ancestral character states, through (i) the correlation between DNA evolutionary processes and species life‐history traits and (ii) available reliable methods for ancestral sequence inference. Here, we assess the relevance of mitochondrial DNA – the most popular molecular marker in animals – as a predictor of ancestral life‐history traits in mammals, using the order of Cetartiodactyla as a benchmark. Using the complete set of 13 mitochondrial protein‐coding genes, we show that the lineage‐specific nonsynonymous over synonymous substitution rate ratio (dN/dS) is closely correlated with the species body mass, longevity and age of sexual maturity in Cetartiodactyla and can be used as a marker of ancestral traits provided that the noise introduced by short branches is appropriately dealt with. Based on ancestral dN/dS estimates, we predict that the first cetartiodactyls were relatively small animals (around 20 kg). This finding is in accordance with Cope's rule and the fossil record but could not be recovered via continuous character evolution methods.     

Reconstructing the evolutionary history of China: a caveat about inferences drawn from ancient DNA

The decipherment of the meager information provided by short fragments of ancient mitochondrial DNA (mtDNA) is notoriously difficult but is regarded as a most promising way toward reconstructing the past from the genetic perspective. By haplogroup-specific hypervariable segment (HVS) motif search and matching or near-matching with available modern data sets, most of the ancient mtDNAs can be tentatively assigned to haplogroups, which are often subcontinent specific. Further typing for mtDNA haplogroup-diagnostic coding region polymorphisms, however, is indispensable for establishing the geographic/genetic affinities of ancient samples with less ambiguity. In the present study, we sequenced a fragment (approximately 982 bp) of the mtDNA control region in 76 Han individuals from Taian, Shandong, China, and we combined these data with previously reported samples from Zibo and Qingdao, Shandong. The reanalysis of two previously published ancient mtDNA population data sets from Linzi (same province) then indicates that the ancient populations had features in common with the modern populations from south China rather than any specific affinity to the European mtDNA pool. Our results highlight that ancient mtDNA data obtained under different sampling schemes and subject to potential contamination can easily create the impression of drastic spatiotemporal changes in the genetic structure of a regional population during the past few thousand years if inappropriate methods of data analysis are employed.     

Genetic continuity after the collapse of the Wari empire: Mitochondrial DNA profiles from Wari and post‐Wari populations in the ancient Andes

The Wari empire flourished in the central, highland Peruvian Andes from AD 600–1000, and although the events that led to its demise are unknown, archaeological evidence indicates that Wari control waned at the end of the first millennium. Here, we test the hypothesis that, despite the major shift in social and political organization at the fall of the Wari empire, the mitochondrial DNA (mtDNA) composition of populations from the Ayacucho Basin, the former imperial heartland of the empire, remained essentially unchanged. Results show that mtDNA haplogroup frequencies among the Wari and post‐Wari groups differ, but the difference is not statistically significant (χ² = 5.886, df = 3, P = 0.1172). This is the first study in the Andes to use haplotypic data to evaluate the observed genetic distance between two temporally distinct prehispanic populations (F_ST = 0.029) against modeled expectations of four possible evolutionary scenarios. None of these simulations allowed the rejection of continuity. In total, at both the haplogroup and haplotype levels these data do not allow us to reject the hypothesis that post‐Wari individuals sampled in this study are the maternal descendants of those sampled from the Wari era site of Conchopata. However, genetic homogeneity in the mitochondrial gene pool, as seen in the late prehispanic southern Andes, may also characterize our study region. But, prior to this research, this was unknown. If our new data show mtDNA homogeneity, then this could limit the detection of female migration if, in fact, it occurred. Nonetheless, the novel mtDNA data presented here currently do not support the hypothesis that there was an influx of genetically distinct females into the former Wari heartland after the Wari collapse. Am J Phys Anthropol, 2009. © 2009 Wiley‐Liss, Inc.