Trauma and violence in the Wari empire of the Peruvian Andes: warfare, raids, and ritual fights

This study examines bioarchaeological evidence for violence during the period of Wari imperialism in the Peruvian Andes through analysis of skeletal trauma from three populations dating to AD 650-800. The samples are from contemporaneous archaeological sites: Conchopata, a Wari heartland site in central highland Peru; Beringa, a community of commoners in the Majes valley of the southern Wari hinterland; and La Real, a high status mortuary site, also in the Majes valley. Given the expansionist nature of Wari and its military-related iconography and weaponry, it is hypothesized that Wari imperialism was concomitant with greater levels of violence relative to other prehispanic groups in the Andes. It is also hypothesized that differential articulation with the Wari empire (e.g., heartland vs. hinterland groups) affected the frequency and patterning of trauma. Results show that cranial trauma frequency of the three Wari era samples is significantly greater than several other Andean skeletal populations. This suggests that Wari rule was associated with high levels of violence, though it may not have always been related to militarism. The three adult samples show similar frequencies of cranial trauma (Conchopata = 26%; Beringa = 33%; La Real = 31%). This may suggest that differential positioning in the Wari empire had little effect on exposure to violence. Sex-based differences in cranial trauma frequencies are present only at La Real, but wound patterning differs between the sexes: females display more wounds on the posterior of the cranium, while males show more on the anterior. These data suggest that Wari rule may have contributed to violence.     

Investigating regional mobility in the southern hinterland of the Wari Empire: biogeochemistry at the site of Beringa, Peru

Empires have transformed political, social, and environmental landscapes in the past and present. Although much research on archaeological empires focuses on large-scale imperial processes, we use biogeochemistry and bioarchaeology to investigate how imperialism may have reshaped regional political organization and regional migration patterns in the Wari Empire of the Andean Middle Horizon (ca. AD 600-1000). Radiogenic strontium isotope analysis of human remains from the site of Beringa in the Majes Valley of southern Peru identified the geographic origins of individuals impacted by the Wari Empire. At Beringa, the combined archaeological human enamel and bone values range from (87)Sr/(86)Sr = 0.70802 - 0.70960, with a mean (87)Sr/(86)Sr = 0.70842 ± 0.00027 (1σ, n = 52). These data are consistent with radiogenic strontium isotope data from the local fauna in the Majes Valley and imply that most individuals were local inhabitants, rather than migrants from the Wari heartland or some other locale. There were two outliers at Beringa, and these "non-local" individuals may have derived from other parts of the South Central Andes. This is consistent with our understanding of expansive trade networks and population movement in the Andean Middle Horizon, likely influenced by the policies of the Wari Empire. Although not a Wari colony, the incorporation of small sites like Beringa into the vast social and political networks of the Middle Horizon resulted in small numbers of migrants at Beringa.     

Phylogeography of the brown hare (Lepus europaeus) in Europe: a legacy of south‐eastern Mediterranean refugia?

Aim We analysed the population genetics of the brown hare (Lepus europaeus) in order to test the hypothesis that this species migrated into central Europe from a number of late glacial refugia, including some in Asia Minor. Location Thirty‐three localities in Greece, Bulgaria, Italy, Croatia, Serbia, Poland, Switzerland, Austria, France, Germany, the Netherlands, Spain, the United Kingdom, Turkey and Israel. Methods In total, 926 brown hares were analysed for mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism (RFLP) performed on polymerase chain reaction‐amplified products spanning cytochrome b (cyt b)/control region (CR), cytochrome oxidase I (COI) and 12S–16S rRNA. In addition, sequence analysis of the mtDNA CR‐I region was performed on 69 individuals, and the data were compared with 137 mtDNA CR‐I sequences retrieved from GenBank. Results The 112 haplotypes detected were partitioned into five phylogeographically well‐defined major haplogroups, namely the ‘south‐eastern European type haplogroup’ (SEEh), ‘Anatolian/Middle Eastern type haplogroup’ (AMh), ‘European type haplogroup, subgroup A’ (EUh‐A), ‘European type haplogroup, subgroup B’ (EUh‐B) and ‘Intermediate haplogroup’ (INTERh). Sequence data retrieved from GenBank were consistent with the haplogroups determined in this study. In Bulgaria and north‐eastern Greece numerous haplotypes of all five haplogroups were present, forming a large overlap zone. Main conclusions The mtDNA results allow us to infer post‐glacial colonization of large parts of Europe from a late glacial/early Holocene source population in the central or south‐central Balkans. The presence of Anatolian/Middle Eastern haplotypes in the large overlap zone in Bulgaria and north‐eastern Greece reveals gene flow from Anatolia to Europe across the late Pleistocene Bosporus land‐bridge. Although various restocking operations could be partly responsible for the presence of unexpected haplotypes in certain areas, we nevertheless trace a strong phylogeographic signal throughout all regions under study. Throughout Europe, mtDNA results indicate that brown hares are not separated into discernable phyletic groups.     

NATURALLY OCCURRING MITOCHONDRIAL DNA HAPLOTYPES EXHIBIT METABOLIC DIFFERENCES: INSIGHT INTO FUNCTIONAL PROPERTIES OF MITOCHONDRIA

Linking the mitochondrial genotype and the organismal phenotype is of paramount importance in evolution of mitochondria. In this study, we determined the differences in catalytic properties of mitochondria dictated by divergences in the siII and siIII haplogroups of Drosophila simulans using introgressions of siII mtDNA type into the siIII nuclear background. We used a novel in situ method (permeabilized fibers) that allowed us to accurately measure the consumption of oxygen by mitochondria in constructed siII‐introgressed flies and in siIII‐control flies. Our results showed that the catalytic capacity of the electron transport system is not impaired by introgressions, suggesting that the functional properties of mitochondria are tightly related to the mtDNA haplogroup and not to the nuclear DNA or to the mito‐nuclear interactions. This is the first study, to our knowledge, that demonstrates a naturally occurring haplogroup can confer specific functional differences in aspects of mitochondrial metabolism. This study illustrates the importance of mtDNA changes on organelle evolution and highlights the potential bioenergetic and metabolic impacts that divergent mitochondrial haplogroups may have upon a wide variety of species including humans.     

Dynamics of haplogroup frequencies and survival rates in a contact zone of two mtDNA lineages of the lizard Lacerta vivipara

The analysis of contact zones between lineages that were previously isolated in allopatry can lead to important insights on evolutionary processes such as selection and adaptation. In this paper we conducted a comparative demographic study of two mitochondrial DNA (mtDNA) lineages of the lizard Lacerta vivipara in the western Pyrénées to provide detail on the dynamics of their contact zone. By surveying haplogroup frequency across the contact area, we revealed the existence of a stable and very narrow contact zone between two parapatric lineages, which we infer to demonstrate a role for selection in the maintenance of this contact zone. We suggest these two lineages evolved in allopatry after retreating to different refugia during the Pleistocene glaciations, and subsequently came into secondary contact after the last glacial maximum. Although haplogroup frequencies were stable over time, we found significant age and environment (temperature) dependent survival differences between mtDNA haplogroups in one contact population sampled yearly from 2002 to 2009. Therefore, temperature‐induced demographic differences between the two mtDNA lineages may be responsible for the stability of this narrow contact zone. This is one of the first demographic studies conducted under natural conditions indicating the possibility of selection on mtDNA.     

Do Haplogroups H and U Act to Increase the Penetrance of Alzheimer’s Disease?

1. Alzheimer’s disease (AD) is the most common form of dementia in the elderly in which interplay between genes and the environment is supposed to be involved. Mitochondrial DNA (mtDNA) has the only noncoding regions at the displacement loop (D-loop) region that contains two hypervariable segments (HVS-I and HVS-II) with high polymorphism. mtDNA has already been fully sequenced and many subsequent publications have shown polymorphic sites, haplogroups, and haplotypes. Haplogroups could have important implications to understand the association between mutability of the mitochondrial genome and the disease. 2. To assess the relationship between mtDNA haplogroup and AD, we sequenced the mtDNA HVS-I in 30 AD patients and 100 control subjects. We could find that haplogroups H and U are significantly more abundant in AD patients (P = 0.016 for haplogroup H and P = 0.0003 for haplogroup U), Thus, these two haplogroups might act synergistically to increase the penetrance of AD disease.     

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill‐coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function‐altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal‐level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.     

Multiple Pleistocene refugia and post‐glacial colonization in the European chub (Squalius cephalus) revealed by combined use of nuclear and mitochondrial markers

Aim To analyse patterns of nuclear and mitochondrial genetic variation in the European chub, Squalius cephalus (Linnaeus, 1758), in order to understand the evolutionary history of this species and to test biogeographical hypotheses for the existence of co‐distributed European freshwater fish species. Location Rivers in Europe (Finland, Poland, Czech Republic, France, Bulgaria, Spain, Italy). Methods We genotyped 12 polymorphic microsatellite markers derived from 310 individuals collected from across the distribution of S. cephalus in Europe (including a total of 15 populations) and sequenced mitochondrial DNA (mtDNA) from a subset of 75 individuals. Sequences of mtDNA cytochrome b were analysed using both phylogenetic (median‐joining networks) and population genetic methods (tests for demographic history, mismatch distributions, Bayesian coalescent analysis). Geographical structure in microsatellite loci was examined using a distance method (F_ST), factorial correspondence analysis (FCA) and a Bayesian clustering method (structure ). Results The mtDNA network showed a clear split into four different haplogroup lineages: Western (separated into Atlantic and Danubian sublineages), Eastern, Aegean (occurring in two distinct sublineages in the Balkans and in Spain) and Adriatic. Our results indicate recent population expansion in the Eastern and Western Atlantic lineages and the admixture of two previously separate sublineages (Atlantic and Danubian) in the Western lineage. Bayesian structure analysis as well as FCA results roughly corresponded to the mtDNA‐based structure, separating the sampled individuals into almost non‐overlapping groups. Main conclusions Our results support hypotheses suggesting origins of extant lineages of freshwater fishes in multiple refugia and the subsequent post‐glacial colonization of Europe via different routes. We confirmed the previously proposed two‐step expansion scenario from the Danube refuge, the existence of a secondary (Atlantic) refuge during the last glaciation (probably in the Rhone River) and population expansion of this lineage. Conspicuous divergences among Mediterranean populations reflect their different origin, as well as their low contribution to the recent genetic pool of chub in central Europe.     

Analysis of Mitochondrial DNA Lineages in Yakuts

To study the mitochondrial gene pool structure in Yakuts, polymorphism of mtDNA hypervariable segment I (16,024–16,390) was analyzed in 191 people sampled from the indigenous population of the Sakha Republic. In total, 67 haplotypes of 14 haplogroups were detected. Most (91.6%) haplotypes belonged to haplogroups A, B, C, D, F, G, M*, and Y, which are specific for East Eurasian ethnic groups; 8.4% haplotypes represented Caucasian haplogroups H, HV1, J, T, U, and W. A high frequency of mtDNA types belonging to Asian supercluster M was peculiar for Yakuts: mtDNA types belonging to haplogroup C, D, or G and undifferentiated mtDNA types of haplogroup M (M*) accounted for 81% of all haplotypes. The highest diversity was observed for haplogroups C and D, which comprised respectively 22 (44%) and 18 (30%) haplotypes. Yakuts showed the lowest genetic diversity (H = 0.964) among all Turkic ethnic groups. Phylogenetic analysis testified to common genetic substrate of Yakuts, Mongols, and Central Asian (Kazakh, Kyrgyz, Uighur) populations. Yakuts proved to share 21 (55.5%) mtDNA haplotypes with the Central Asian ethnic groups and Mongols. Comparisons with modern Paleoasian populations (Chukcha, Itelmen, Koryaks) revealed three (8.9%) haplotypes common for Yakuts and Koryaks. The results of mtDNA analysis disagree with the hypothesis of an appreciable Paleoasian contribution to the modern Yakut gene pool.     

Y‐chromosome and mtDNA variation confirms independent domestications and directional hybridization in South American camelids

Investigations of genetic diversity and domestication in South American camelids (SAC) have relied on autosomal microsatellite and maternally‐inherited mitochondrial data. We present the first integrated analysis of domestic and wild SAC combining male and female sex‐specific markers (male specific Y‐chromosome and female‐specific mtDNA sequence variation) to assess: (i) hypotheses about the origin of domestic camelids, (ii) directionality of introgression among domestic and/or wild taxa as evidence of hybridization and (iii) currently recognized subspecies patterns. Three male‐specific Y‐chromosome markers and control region sequences of mitochondrial DNA are studied here. Although no sequence variation was found in SRY and ZFY, there were seven variable sites in DBY generating five haplotypes on the Y‐chromosome. The haplotype network showed clear separation between haplogroups of guanaco–llama and vicuña–alpaca, indicating two genetically distinct patrilineages with near absence of shared haplotypes between guanacos and vicuñas. Although we document some examples of directional hybridization, the patterns strongly support the hypothesis that llama (Lama glama) is derived from guanaco (Lama guanicoe) and the alpaca (Vicugna pacos) from vicuña (Vicugna vicugna). Within male guanacos we identified a haplogroup formed by three haplotypes with different geographical distributions, the northernmost of which (Peru and northern Chile) was also observed in llamas, supporting the commonly held hypothesis that llamas were domesticated from the northernmost populations of guanacos (L. g. cacilensis). Southern guanacos shared the other two haplotypes. A second haplogroup, consisting of two haplotypes, was mostly present in vicuñas and alpacas. However, Y‐chromosome variation did not distinguish the two subspecies of vicuñas.     

Phylogeography and ecological niche modeling unravel the evolutionary history of the African green toad,Bufotes boulengeri boulengeri (Amphibia: Bufonidae), through the Quaternary

Recent integration of ecological niche models in phylogeographic studies is improving our understanding of the processes structuring genetic variation across landscapes. Previous studies on the amphibian Bufotes boulengeri boulengeri uncovered a surprisingly weak intraspecific differentiation across the Maghreb region. We widely sampled this species from Morocco to Egypt and sequenced one nuclear and three mitochondrial (mtDNA) genes to determine the level of genetic variability across its geographic range. We evaluated these data with ecological niche modeling to reveal its evolutionary history in response to climate change during the Quaternary. Our results highlight some mtDNA phylogeographic structure within this species, with one haplogroup endemic to coastal Morocco, and one haplogroup widely distributed throughout North Africa. No or little genetic differentiation is observed between isolated populations from the Hoggar Mountains, the Sabha district and the islands of Kerkennah and Lampedusa, compared to others populations. This can be explained by the expansion of the distribution range of B. b. boulengeri during glacial periods. This might have facilitated the species’ dispersal and subsequent gene flow between most North African localities.     

Population Structure of the Giant Garter Snake, Thamnophis gigas

The giant garter snake, Thamnophis gigas, is a threatened species endemic to California’s Central Valley. We tested the hypothesis that current watershed boundaries have caused genetic differentiation among populations of T. gigas. We sampled 14 populations throughout the current geographic range of T. gigas and amplified 859 bp from the mitochondrial gene ND4 and one nuclear microsatellite locus. DNA sequence variation from the mitochondrial gene indicates there is some genetic structuring of the populations, with high F_ST values and unique haplotypes occurring at high frequency in several populations. We found that clustering populations by watershed boundary results in significant between-region genetic variance for mtDNA. However, analysis of allele frequencies at the microsatellite locus NSU3 reveals very low F_ST values and little between-region variation in allele frequencies. The discordance found between mitochondrial and microsatellite data may be explained by aspects of molecular evolution and/or T. gigas life history characteristics. Differences in effective population size between mitochondrial and nuclear DNA, or male-biased gene flow, result in a lower migration rate of mitochondrial haplotypes relative to nuclear alleles. However, we cannot exclude homoplasy as one explanation for homogeneity found for the single microsatellite locus. The mitochondrial nucleotide sequence data supports conservation practices that identify separate management units for T. gigas.     

Mitochondrial D-Loop Variation in Persian Multiple Sclerosis Patients: K and A Haplogroups as a Risk Factor!!

Multiple Sclerosis (MS) is a multifocal demyelinating central nervous system disorder in which interplay between genes and the environment are supposed to be involved. Mitochondrial DNA (mtDNA) has the only non-coding regions at the displacement loop (D-loop) region that contains two hypervariable segments (HVS-I and HVS-II) with high polymorphism. mtDNA has already been fully sequenced and many subsequent publications have showed polymorphic sites, haplogroups and haplotypes. Haplogroups could have important implications to understand association between mutability of the mitochondrial genome and disease. To assess relationship between mtDNA haplogroups and MS, we have sequenced the mtDNA HVS-I in 54 MS patients and 100 control subjects. We have found that haplogroups A and K are significantly more abundant in MS patients (P=0.042 for haplogroup A and P=0.0005 for haplogroup K). Thus, these two haplogroups might act synergistically to increase the penetrance of MS disease.     

E(y)2/Sus1 is required for blocking PRE silencing by the Wari insulator in Drosophila melanogaster

Chromatin insulators affect interactions between promoters and enhancers/silencers and function as barriers to the spread of repressive chromatin. Recently, we have found an insulator, named Wari, located on the 3′ side of the white gene. Here, we show that the previously identified 368-bp core of this insulator is sufficient for blocking Polycomb response element-mediated silencing. Although Wari does not contain binding sites for known insulator proteins, the E(y)2 and CP190 proteins bind to Wari as well as to the Su(Hw)-containing insulators in vivo. It may well be that these proteins are recruited to the insulator by as yet unidentified DNA-binding protein. Partial inactivation of E(y)2 in a weak e(y)2 ^u1 mutation impairs only the anti-silencing but not the enhancer-blocking activity of the Wari insulator. Thus, the E(y)2 protein in different Drosophila insulators serves to protect gene expression from silencing.     

Patterns of genetic isolation in a widely distributed pelagic fish,the narrow‐barred Spanish mackerel (Scomberomorus commerson)

Although migratory pelagic fishes generally exhibit little geographic differentiation across oceans, as expected from their life history (broadcast spawning, pelagic larval life, swimming ability of adults) and the assumed homogeneity of the pelagic habitat, exceptions to the rule deserve scrutiny. One such exception is the narrow‐barred Spanish mackerel (Scomberomorus commerson Lacepède, 1800), where strong genetic heterogeneity at the regional scale has been previously reported. We investigated the genetic composition of S. commerson across the Indo‐West Pacific range using control‐region sequences (including previously published data sets), cytochrome b gene partial sequences, and eight microsatellite loci, to further explore its phylogeographic structure. All haplotypes sampled from the Indo‐Malay‐Papua archipelago (IMPA) and the south‐western Pacific coalesced into a clade (clade II) that was deeply separated (14.5% nucleotide divergence) from a clade grouping all haplotypes from the Persian Gulf and Oman Sea (clade I). Such a high level of genetic divergence suggested the occurrence of two sister species. Further phylogeographic partition was evident between the western IMPA and the regions sampled east and south of it, i.e. northern Australia, West Papua, and the Coral Sea. Strong allele‐frequency differences were found between local populations in the south‐western Pacific, both at the mitochondrial locus (Φ_st = 0.282–0.609) and at microsatellite loci ( = 0.202–0.313). Clade II consisted of four deeply divergent subclades (9.0–11.8% nucleotide divergence for the control region; 0.3–2.5% divergence at the cytochrome b locus). Mitochondrial subclades within clade II generally had narrow geographic distribution, demonstrating further genetic isolation. However, one particular haplogroup within clade II was present throughout the central Indo‐West Pacific: this haplogroup was found to be the sister group to a haplogroup restricted to West Papua and the Coral Sea, yielding evidence of recent secondary westward colonization. Such a complex structure is in sharp contrast with the generally weak phylogeographic patterns uncovered to date in other widely distributed, large pelagic fishes with pelagic eggs and larvae. We hypothesize that in S. commerson and possibly other Scomberomorus species, philopatric migration may play a role in maintaining the geographic isolation of populations by annihilating the potential consequences of passive dispersal. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 886–902.     

Mitochondrial DNA Variation in Russian Populations of Stavropol Krai, Orel and Saratov Oblasts

Mitochondrial DNA (mtDNA) polymorphism was examined in three Russian populations from the European part of Russia (Stavropol krai, Orel oblast, and Saratov oblast). This analysis showed that mitochondrial gene pool of Russians was represented by the mtDNA types belonging to haplogroups H, V, HV*, J, T, U, K, I, W, and X. A mongoloid admixture (1.5%) was revealed in the form of mtDNA types of macrohaplogroup M. Comparative analysis of the mtDNA haplogroup frequency distribution patterns in six Russian populations from the European part of Russia indicated the absence of substantial genetic differences between them. However, in Russian populations from the southern and central regions the frequency of haplogroup V (average frequency 8%) was higher than in the populations from more northern regions. Based on the data on mtDNA HVS1 sequence variation, it was shown that the diversity of haplogroup V in Russians (h= 0.742) corresponded to the highest h values observed in Europe. The reasons for genetic differentiation of the Russian population (historical, ecological, and adaptive) are discussed.     

Mitochondrial DNA phylogeography of least cisco Coregonus sardinella in Alaska

This study presents the first detailed analysis of the mitochondrial DNA diversity of least cisco Coregonus sardinella in Alaska using a 678 bp segment of the control region (D‐loop) of the mitochondrial genome. Findings suggest that the history of C. sardinella in Alaska differs from that of other species of Coregonus present in the state and surrounding regions. The examined populations of C. sardinella are genetically diverse across Alaska. Sixty‐eight distinct mitochondrial haplotypes were identified among 305 individuals sampled from nine locations. The haplotype minimum spanning network and phylogeny showed a modest level of geographic segregation among haplotypes, suggesting high levels of on‐going or recent connectivity among distant populations. Observed Φ_ST values and the results of homogeneity and AMOVAs indicate incipient genetic differentiation between aggregations in three broad regional groups. Sites north of the Brooks Range formed one group, sites in the Yukon and Selawik Rivers formed a second group and sites south of the Yukon drainage formed the third group. Overall, the sequence data showed that a large proportion of mtDNA genetic variation in C. sardinella is shared across Alaska, but this variation is not homogeneously distributed across all regions and for all haplotype groups.     

Mitochondrial DNA haplogrouping of the brown bear,Ursus arctos (Carnivora: Ursidae) in Asia,based on a newly developed APLP analysis

Sequence analyses of the complete brown bear, Ursus arctos, mitochondrial DNA (mtDNA) genome have detected scattered single nucleotide polymorphisms (SNPs) that define distinct mtDNA haplogroups in phylogeographical studies. The degraded DNA in historical samples, such as stuffed or excavated specimens, however, is often not suitable for sequence analyses. To address this problem, we developed an amplified product length polymorphism (APLP) analysis for mtDNA‐haplogrouping U. arctos specimens by detecting haplogroup‐specific SNPs. We verified the validity and utility of this method by analysing up to 170‐year‐old skin samples from U. arctos specimens collected widely across continental Eurasia. We detected some of the same haplogroups as those occurring in eastern Hokkaido (Japan) and eastern Alaska in continental Eurasia (the Altai and the Caucasus). Our results show that U. arctos in eastern Hokkaido and eastern Alaska descended from a common ancestor in continental Eurasia, and suggest that U. arctos occupied several refugia in southern Asia during the Last Glacial Maximum. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 627–635.     

Applicability of mitochondrial DNA for the identification of Arvicolid species from faecal samples: a case study from the threatened Cabrera's vole

Arvicolid mitochondrial genomes evolve faster than in any other mammalian lineage. The genetic diversity exhibited by these rodents contrasts sharply with their phenotypic homogeneity. Furthermore, faecal droppings from Arvicolid rodents of similar body size are almost undistinguishable on the basis of pellet morphology and content. In this study, we advantaged from their high genetic diversity vs. phenotypic homogeneity to document the applicability of mtDNA extraction from vole droppings for latter identification of such via a rapid and efficient nested PCR‐based technique using the threatened Microtus cabrerae as a model species. We sequenced the mitochondrial control region from 75 individuals belonging to 11 species of Arvicolinae from Spain, Portugal, Greece and Italy, and an additional 19 sequences from ten Microtus species from other countries were downloaded from Genbank. Based on these control region sequences, we successfully designed and applied a nested PCR for M. cabrerae‐specific and arvicolid‐generic mtDNA markers to differentiate Cabrera’s vole faecal samples among other species of the Arvicolinae subfamily. Although this study used Cabrera’s vole as a model species, similar techniques based on mtDNA sequences may find a broader applicability for noninvasive genetic conservation of vole species and their populations.