The first complete mitochondrial genome data of the Afghan pika Ochotona rufescens (Lagomorpha, Ochotonidae), near the type locality.

The Afghan pika Ochotona rufescens (Gray, 1842) is widely distributed across the mountains of Afghanistan, Iran, Pakistan, and southwestern Turkmenistan, most often at elevations between 2,000 and 3,000 m. Here we present, for the first time, the complete mitochondrial genomes of two specimens of the nominotypical subspecies Ochotona rufescens rufescens, de novo assembled from Illumina short reads of fragmented probe-enriched DNA. The lengths of the circular mitogenomes are 16,408 bp and 16,407 bp, respectively. Both mitogenomes contain 13 protein-coding genes (PCGs), two ribosomal RNAs (16S rRNA and 12S rRNA), 22 transfer RNA genes, and a control region. The gene NAD6 and the tRNA (Gln), tRNA (Ala), tRNA (Asn), tRNA (Cys), tRNA (Tyr), tRNA (Ser), tRNA (Glu), and tRNA (Pro) are encoded on the light strand while the rest are encoded on the heavy strand. The overall nucleotide composition was ∼30% for A, 25% for T, 15% for G, and 29% for C. The mitogenome data are available in the GenBank under the accession numbers ON859136 and ON859137.

Ancient DNA of narrow-headed vole reveal common features of the Late Pleistocene population dynamics in cold-adapted small mammals.

The narrow-headed vole, collared lemming and common vole were the most abundant small mammal species across the Eurasian Late Pleistocene steppe-tundra environment. Previous ancient DNA studies of the collared lemming and common vole have revealed dynamic population histories shaped by climatic fluctuations. To investigate the extent to which species with similar adaptations share common evolutionary histories, we generated a dataset comprised the mitochondrial genomes of 139 ancient and 6 modern narrow-headed voles from several sites across Europe and northwestern Asia covering approximately the last 100 thousand years (kyr). We inferred Bayesian time-aware phylogenies using 11 radiocarbon-dated samples to calibrate the molecular clock. Divergence of the main mtDNA lineages across the three species occurred during marine isotope stages (MIS) 7 and MIS 5, suggesting a common response of species adapted to open habitat during interglacials. We identified several time-structured mtDNA lineages in European narrow-headed vole, suggesting lineage turnover. The timing of some of these turnovers was synchronous across the three species, allowing us to identify the main drivers of the Late Pleistocene dynamics of steppe- and cold-adapted species.

Population dynamics and demographic history of Eurasian collared lemmings.

BACKGROUND: Ancient DNA studies suggest that Late Pleistocene climatic changes had a significant effect on population dynamics in Arctic species. The Eurasian collared lemming (Dicrostonyx torquatus) is a keystone species in the Arctic ecosystem. Earlier studies have indicated that past climatic fluctuations were important drivers of past population dynamics in this species. RESULTS: Here, we analysed 59 ancient and 54 modern mitogenomes from across Eurasia, along with one modern nuclear genome. Our results suggest population growth and genetic diversification during the early Late Pleistocene, implying that collared lemmings may have experienced a genetic bottleneck during the warm Eemian interglacial. Furthermore, we find multiple temporally structured mitogenome clades during the Late Pleistocene, consistent with earlier results suggesting a dynamic late glacial population history. Finally, we identify a population in northeastern Siberia that maintained genetic diversity and a constant population size at the end of the Pleistocene, suggesting suitable conditions for collared lemmings in this region during the increasing temperatures associated with the onset of the Holocene. CONCLUSIONS: This study highlights an influence of past warming, in particular the Eemian interglacial, on the evolutionary history of the collared lemming, along with spatiotemporal population structuring throughout the Late Pleistocene.

The first complete mitochondrial genome data of the pygmy rabbit Brachylagus idahoensis, the world's smallest leporid.

The pygmy rabbit Brachylagus idahoensis (Merriam, 1891) is the smallest extant leporid, which naturally occurs in the Great Basin and adjacent areas in western parts of the United States of America. Its distribution is strongly associated with the sagebrush (Artemisia ssp.) vegetation. Here we present, for the first time, the complete mitochondrial genome of Brachylagus idahoensis, de novo assembled from Illumina short reads of fragmented probe-enriched DNA. The circular mitogenome is 17,021 bp in length and contains 13 protein-coding genes (PCGs), two ribosomal RNAs (16S rRNA and 12S rRNA), 22 transfer RNA genes, and a control region. The gene NAD6 and the tRNA(Gln), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), tRNA(Ser), tRNA(Glu) and tRNA(Pro) are encoded on the light strand while the rest are encoded on the heavy strand. The overall nucleotide composition was 30.78% for A, 28.5% for T, 13.62% for G and 27.08% for C. The mitogenome data are available in the GenBank under the accession number OL436257.

Stable isotopes unveil one millennium of domestic cat paleoecology in Europe.

The domestic cat is the world's most popular pet and one of the most detrimental predators in terrestrial ecosystems. Effective protection of wildlife biodiversity demands detailed tracking of cat trophic ecology, and stable isotopes serve as a powerful proxy in dietary studies. However, a variable diet can make an isotopic pattern unreadable in opportunistic predators. To evaluate the usefulness of the isotopic method in cat ecology, we measured C and N isotope ratios in hundreds of archaeological cat bones. We determined trends in cat trophic paleoecology in northern Europe by exploiting population-scale patterns in animals from diverse locations. Our dataset shows a high variability of isotopic signals related to the socio-economic and/or geomorphological context. This points toward regularities in isotopic patterns across past cat populations. We provide a generalized guide to interpret the isotopic ecology of cats, emphasizing that regional isotopic baselines have a major impact on the isotopic signal.

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.

Ancient genomes reveal long-range influence of the pre-Columbian culture and site of Tiwanaku.

Tiwanaku civilization flourished in the Lake Titicaca basin between 500 and 1000 CE and at its apogee influenced wide areas across the southern Andes. Despite a considerable amount of archaeological data, little is known about the Tiwanaku population. We analyzed 17 low-coverage genomes from individuals dated between 300 and 1500 CE and demonstrated genetic continuity in the Lake Titicaca basin throughout this period, which indicates that the substantial cultural and political changes in the region were not accompanied by large-scale population movements. Conversely, the ritual center of Tiwanaku revealed high diversity, including individuals with primarily local genetic ancestry and those with foreign admixture or provenance from as far as the Amazon. Nonetheless, most human offerings associated with the Akapana platform exhibited pure Titicaca basin ancestry and dated to ca. 950 CE­the onset of Tiwanaku's decline as a sociopolitical center. Our results strengthen the view of Tiwanaku as a complex and far-reaching polity.

Ancestors of domestic cats in Neolithic Central Europe: Isotopic evidence of a synanthropic diet.

Cat remains from Poland dated to 4,200 to 2,300 y BCE are currently the earliest evidence for the migration of the Near Eastern cat (NE cat), the ancestor of domestic cats, into Central Europe. This early immigration preceded the known establishment of housecat populations in the region by around 3,000 y. One hypothesis assumed that NE cats followed the migration of early farmers as synanthropes. In this study, we analyze the stable isotopes in six samples of Late Neolithic NE cat bones and further 34 of the associated fauna, including the European wildcat. We approximate the diet and trophic ecology of Late Neolithic felids in a broad context of contemporary wild and domestic animals and humans. In addition, we compared the ecology of Late Neolithic NE cats with the earliest domestic cats known from the territory of Poland, dating to the Roman Period. Our results reveal that human agricultural activity during the Late Neolithic had already impacted the isotopic signature of rodents in the ecosystem. These synanthropic pests constituted a significant proportion of the NE cat's diet. Our interpretation is that Late Neolithic NE cats were opportunistic synanthropes, most probably free-living individuals (i.e., not directly relying on a human food supply). We explore niche partitioning between studied NE cats and the contemporary native European wildcats. We find only minor differences between the isotopic ecology of both these taxa. We conclude that, after the appearance of the NE cat, both felid taxa shared the ecological niches.

Highly divergent lineage of narrow-headed vole from the Late Pleistocene Europe.

During the Late Pleistocene, narrow-headed voles (Lasiopodomys gregalis) inhabited Eurasia's vast territories, frequently becoming the dominant small mammal species among steppe-tundra communities. We investigated the relationship between this species' European and Asiatic populations by sequencing the mtDNA genomes of two extant specimens from Russia and 10 individuals from five Central European sites, dated to the post-LGM period. Phylogenetic analyses based on a large portion of mtDNA genomes highly supported the positioning of L. gregalis within Arvicolinae. The phylogeny based on mtDNA cytochrome b sequences revealed a deep divergence of European narrow-headed voles from Asiatic ones and their sister position against the extant L. gregalis and L. raddei. The divergence of the European lineage was estimated to a minimum 230 thousand years ago. This suggest, contrary to the current biogeographic hypotheses, that during the interglacial periods narrow-headed vole did not retreat from Europe but survived the unfavourable conditions within the refugial areas. Based on this result, we propose to establish a cryptic species status for the Late Pleistocene European narrow-headed vole and to name this taxon Lasiopodomys anglicus.

Human-mediated dispersal of cats in the Neolithic Central Europe.

Archeological and genetic evidence suggest that all domestic cats derived from the Near Eastern wildcat (Felis silvestris lybica) and were first domesticated in the Near East around 10,000 years ago. The spread of the domesticated form in Europe occurred much later, primarily mediated by Greek and Phoenician traders and afterward by Romans who introduced cats to Western and Central Europe around 2000 years ago. We investigated mtDNA of Holocene Felis remains and provide evidence of an unexpectedly early presence of cats bearing the Near Eastern wildcat mtDNA haplotypes in Central Europe, being ahead of Roman period by over 2000 years. The appearance of the Near Eastern wildcats in Central Europe coincides with the peak of Neolithic settlement density, moreover most of those cats belonged to the same mtDNA lineages as those domesticated in the Near East. Thus, although we cannot fully exclude that the Near Eastern wildcats appeared in Central Europe as a result of introgression with European wildcat, our findings support the hypothesis that the Near Eastern wildcats spread across Europe together with the first farmers, perhaps as commensal animals. We also found that cats dated to the Neolithic period belonged to different mtDNA lineages than those brought to Central Europe in Roman times, this supports the hypothesis that the gene pool of contemporary European domestic cats might have been established from two different source populations that contributed in different periods.

Retreat and extinction of the Late Pleistocene cave bear (Ursus spelaeus sensu lato).

The cave bear (Ursus spelaeus sensu lato) is a typical representative of Pleistocene megafauna which became extinct at the end of the Last Glacial. Detailed knowledge of cave bear extinction could explain this spectacular ecological transformation. The paper provides a report on the youngest remains of the cave bear dated to 20,930 ± 140 14C years before present (BP). Ancient DNA analyses proved its affiliation to the Ursus ingressus haplotype. Using this record and 205 other dates, we determined, following eight approaches, the extinction time of this mammal at 26,100-24,300 cal. years BP. The time is only slightly earlier, i.e. 27,000-26,100 cal. years BP, when young dates without associated collagen data are excluded. The demise of cave bear falls within the coldest phase of the last glacial period, Greenland Stadial 3. This finding and the significant decrease in the cave bear records with cooling indicate that the drastic climatic changes were responsible for its extinction. Climate deterioration lowered vegetation productivity, on which the cave bear strongly depended as a strict herbivore. The distribution of the last cave bear records in Europe suggests that this animal was vanishing by fragmentation into subpopulations occupying small habitats. One of them was the Kraków-Czestochowa Upland in Poland, where we discovered the latest record of the cave bear and also two other, younger than 25,000 14C years BP. The relatively long survival of this bear in karst regions may result from suitable microclimate and continuous access to water provided by deep aquifers, indicating a refugial role of such regions in the Pleistocene for many species.

Synchronous genetic turnovers across Western Eurasia in Late Pleistocene collared lemmings.

Recent palaeogenetic studies indicate a highly dynamic history in collared lemmings (Dicrostonyx spp.), with several demographical changes linked to climatic fluctuations that took place during the last glaciation. At the western range margin of D. torquatus, these changes were characterized by a series of local extinctions and recolonizations. However, it is unclear whether this pattern represents a local phenomenon, possibly driven by ecological edge effects, or a global phenomenon that took place across large geographical scales. To address this, we explored the palaeogenetic history of the collared lemming using a next-generation sequencing approach for pooled mitochondrial DNA amplicons. Sequences were obtained from over 300 fossil remains sampled across Eurasia and two sites in North America. We identified five mitochondrial lineages of D. torquatus that succeeded each other through time across Europe and western Russia, indicating a history of repeated population extinctions and recolonizations, most likely from eastern Russia, during the last 50 000 years. The observation of repeated extinctions across such a vast geographical range indicates large-scale changes in the steppe-tundra environment in western Eurasia during the last glaciation. All Holocene samples, from across the species' entire range, belonged to only one of the five mitochondrial lineages. Thus, extant D. torquatus populations only harbour a small fraction of the total genetic diversity that existed across different stages of the Late Pleistocene. In North American samples, haplotypes belonging to both D. groenlandicus and D. richardsoni were recovered from a Late Pleistocene site in south-western Canada. This suggests that D. groenlandicus had a more southern and D. richardsoni a more northern glacial distribution than previously thought. This study provides significant insights into the population dynamics of a small mammal at a large geographical scale and reveals a rather complex demographical history, which could have had bottom-up effects in the Late Pleistocene steppe-tundra ecosystem.

Complete mitochondrial genome sequences of Atlantic sturgeon, Acipenser oxyrinchus oxyrinchus, Gulf sturgeon, A. o. desotoi and European sturgeon A. sturio (Acipenseriformes: Acipenseridae) obtained through next generation sequencing.

Complete mitochondrial genome sequences of European sturgeon and two subspecies of the North American, Atlantic and Gulf sturgeons were determined using MiSeq Illumina technology. All three genomes show typical vertebrate organization. They possess 22 tRNA genes, 13 protein-coding genes, 2 rRNA (ribosomal RNA) genes and a non-coding control region. Excluding ND6, all protein-coding genes are on the heavy strand. The whole mitogenome sequences have been deposited in GenBank under accession numbers KP997216-KP997218.

Locals, resettlers, and pilgrims: a genetic portrait of three pre-Columbian Andean populations.

The common practice of resettlement and the development of administrative and ceremonial systems shaped the population landscape of the Andean region under the Inca rule. The area surrounding Coropuna and Solimana volcanoes, in the Arequipa region (Peru), carried a high-density, multiethnic population. We studied the genetic variation among three pre-Columbian populations from three functionally diverse archaeological sites excavated in this region. By analyzing the genetic composition of a large ceremonial center (Acchaymarca), an isolated pastoral settlement (Tompullo 2), and an agricultural settlement characterized by architectural features rare in the region (Puca), we investigated the patterns of population movements and the distribution of genetic diversity. We obtained mitochondrial DNA sequences for 25 individuals and autosomal microsatellite profiles for 20 individuals from Acchaymarca and Puca sites. These were compared with previously published genetic data for Tompullo 2 and other pre-Columbian populations. We found differences among the genetic portraits of the three populations, congruent with the archaeologically described functions and characteristics of the sites. The Acchaymarca population had the highest genetic diversity and possessed the lowest number of unique mtDNA haplotypes. The Tompullo 2 population exhibited the lowest level of genetic diversity. The Puca population was distinct from the other two populations owing to a high frequency of haplogroup A haplotypes, what potentially explains the non-local character of the burial architecture. Our analyses of microsatellite data suggest that gene flow between sites was mostly mediated by females, which is consistent with ethnohistorical knowledge of the social organization of the pre-Columbian communities.

Ancient DNA reveals kinship burial patterns of a pre-Columbian Andean community.

BACKGROUND: A detailed genetic study of the pre-Columbian population inhabiting the Tompullo 2 archaeological site (department Arequipa, Peru) was undertaken to resolve the kin relationships between individuals buried in six different chullpas. Kin relationships were an important factor shaping the social organization in the pre-Columbian Andean communities, centering on the ayllu, a group of relatives that shared a common land and responsibilities. The aim of this study was to evaluate whether this Andean model of a social organization had an influence on mortuary practices, in particular to determine whether chullpas served as family graves. RESULTS: The remains of forty-one individuals were analyzed with both uniparental (mtDNA, Y-chromosome) and biparental (autosomal microsatellites) markers. Reproducible HVRI sequences, autosomal and Y chromosomal STR profiles were obtained for 24, 16 and 11 individuals, respectively. Mitochondrial DNA diversity was comparable to that of ancient and contemporary Andean populations. The Tompullo 2 population exhibited the closest relationship with the modern population from the same region. A kinship analysis revealed complex pattern of relations within and between the graves. However mean relatedness coefficients regarding the pairs of individuals buried in the same grave were significantly higher than those regarding pairs buried in different graves. The Y chromosome profiles of 11 males suggest that only members of one male line were buried in the same grave. CONCLUSIONS: Genetic investigation of the population that inhabited Tompullo 2 site shows continuity between pre-Columbian and modern Native Amerindian populations inhabiting the Arequipa region. This suggests that no major demographic processes have influenced the mitochondrial DNA diversity of these populations during the past five hundred years. The kinship analysis involving uni- and biparental markers suggests that the community that inhabited the Tompullo 2 site was organized into extended family groups that were buried in different graves. This finding is in congruence with known models of social organization of Andean communities.