Mitogenomes from Two Uncommon Haplogroups Mark Late Glacial/Postglacial Expansions from the Near East and Neolithic Dispersals within Europe

The current human mitochondrial (mtDNA) phylogeny does not equally represent all human populations but is biased in favour of representatives originally from north and central Europe. This especially affects the phylogeny of some uncommon West Eurasian haplogroups, including I and W, whose southern European and Near Eastern components are very poorly represented, suggesting that extensive hidden phylogenetic substructure remains to be uncovered. This study expanded and re-analysed the available datasets of I and W complete mtDNA genomes, reaching a comprehensive 419 mitogenomes, and searched for precise correlations between the ages and geographical distributions of their numerous newly identified subclades with events of human dispersal which contributed to the genetic formation of modern Europeans. Our results showed that haplogroups I (within N1a1b) and W originated in the Near East during the Last Glacial Maximum or pre-warming period (the period of gradual warming between the end of the LGM, ∼19 ky ago, and the beginning of the first main warming phase, ∼15 ky ago) and, like the much more common haplogroups J and T, may have been involved in Late Glacial expansions starting from the Near East. Thus our data contribute to a better definition of the Late and postglacial re-peopling of Europe, providing further evidence for the scenario that major population expansions started after the Last Glacial Maximum but before Neolithic times, but also evidencing traces of diffusion events in several I and W subclades dating to the European Neolithic and restricted to Europe.     

Ancient DNA suggests modern wolves trace their origin to a Late Pleistocene expansion from Beringia

Grey wolves (Canis lupus) are one of the few large terrestrial carnivores that have maintained a wide geographical distribution across the Northern Hemisphere throughout the Pleistocene and Holocene. Recent genetic studies have suggested that, despite this continuous presence, major demographic changes occurred in wolf populations between the Late Pleistocene and early Holocene, and that extant wolves trace their ancestry to a single Late Pleistocene population. Both the geographical origin of this ancestral population and how it became widespread remain unknown. Here, we used a spatially and temporally explicit modelling framework to analyse a data set of 90 modern and 45 ancient mitochondrial wolf genomes from across the Northern Hemisphere, spanning the last 50,000 years. Our results suggest that contemporary wolf populations trace their ancestry to an expansion from Beringia at the end of the Last Glacial Maximum, and that this process was most likely driven by Late Pleistocene ecological fluctuations that occurred across the Northern Hemisphere. This study provides direct ancient genetic evidence that long‐range migration has played an important role in the population history of a large carnivore, and provides insight into how wolves survived the wave of megafaunal extinctions at the end of the last glaciation. Moreover, because Late Pleistocene grey wolves were the likely source from which all modern dogs trace their origins, the demographic history described in this study has fundamental implications for understanding the geographical origin of the dog.     

Were inefficient mitochondrial haplogroups selected during migrations of modern humans? A test using modular kinetic analysis of coupling in mitochondria from cybrid cell lines

We introduce a general test of the bioenergetic importance of mtDNA (mitochondrial DNA) variants: modular kinetic analysis of oxidative phosphorylation in mitochondria from cybrid cells with constant nuclear DNA but different mtDNA. We have applied this test to the hypothesis [Ruiz-Pesini, Mishmar, Brandon, Procaccio and Wallace (2004) Science 303, 223-226] that particular mtDNA haplogroups (specific combinations of polymorphisms) that cause lowered coupling efficiency, leading to generation of less ATP and more heat, were positively selected during radiations of modern humans into colder climates. Contrary to the predictions of this hypothesis, mitochondria from Arctic haplogroups had similar or even greater coupling efficiency than mitochondria from tropical haplogroups.     

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.     

Complete mitochondrial genomes reveal neolithic expansion into Europe

The Neolithic transition from hunting and gathering to farming and cattle breeding marks one of the most drastic cultural changes in European prehistory. Short stretches of ancient mitochondrial DNA (mtDNA) from skeletons of pre-Neolithic hunter-gatherers as well as early Neolithic farmers support the demic diffusion model where a migration of early farmers from the Near East and a replacement of pre-Neolithic hunter-gatherers are largely responsible for cultural innovation and changes in subsistence strategies during the Neolithic revolution in Europe. In order to test if a signal of population expansion is still present in modern European mitochondrial DNA, we analyzed a comprehensive dataset of 1,151 complete mtDNAs from present-day Europeans. Relying upon ancient DNA data from previous investigations, we identified mtDNA haplogroups that are typical for early farmers and hunter-gatherers, namely H and U respectively. Bayesian skyline coalescence estimates were then used on subsets of complete mtDNAs from modern populations to look for signals of past population expansions. Our analyses revealed a population expansion between 15,000 and 10,000 years before present (YBP) in mtDNAs typical for hunters and gatherers, with a decline between 10,000 and 5,000 YBP. These corresponded to an analogous population increase approximately 9,000 YBP for mtDNAs typical of early farmers. The observed changes over time suggest that the spread of agriculture in Europe involved the expansion of farming populations into Europe followed by the eventual assimilation of resident hunter-gatherers. Our data show that contemporary mtDNA datasets can be used to study ancient population history if only limited ancient genetic data is available.     

The western and eastern roots of the Saami--the story of genetic "outliers" told by mitochondrial DNA and Y chromosomes

The Saami are regarded as extreme genetic outliers among European populations. In this study, a high-resolution phylogenetic analysis of Saami genetic heritage was undertaken in a comprehensive context, through use of maternally inherited mitochondrial DNA (mtDNA) and paternally inherited Y-chromosomal variation. DNA variants present in the Saami were compared with those found in Europe and Siberia, through use of both new and previously published data from 445 Saami and 17,096 western Eurasian and Siberian mtDNA samples, as well as 127 Saami and 2,840 western Eurasian and Siberian Y-chromosome samples. It was shown that the “Saami motif” variant of mtDNA haplogroup U5b is present in a large area outside Scandinavia. A detailed phylogeographic analysis of one of the predominant Saami mtDNA haplogroups, U5b1b, which also includes the lineages of the “Saami motif,” was undertaken in 31 populations. The results indicate that the origin of U5b1b, as for the other predominant Saami haplogroup, V, is most likely in western, rather than eastern, Europe. Furthermore, an additional haplogroup (H1) spread among the Saami was virtually absent in 781 Samoyed and Ob-Ugric Siberians but was present in western and central European populations. The Y-chromosomal variety in the Saami is also consistent with their European ancestry. It suggests that the large genetic separation of the Saami from other Europeans is best explained by assuming that the Saami are descendants of a narrow, distinctive subset of Europeans. In particular, no evidence of a significant directional gene flow from extant aboriginal Siberian populations into the haploid gene pools of the Saami was found.     

Origin and post-glacial dispersal of mitochondrial DNA haplogroups C and D in northern Asia

More than a half of the northern Asian pool of human mitochondrial DNA (mtDNA) is fragmented into a number of subclades of haplogroups C and D, two of the most frequent haplogroups throughout northern, eastern, central Asia and America. While there has been considerable recent progress in studying mitochondrial variation in eastern Asia and America at the complete genome resolution, little comparable data is available for regions such as southern Siberia--the area where most of northern Asian haplogroups, including C and D, likely diversified. This gap in our knowledge causes a serious barrier for progress in understanding the demographic pre-history of northern Eurasia in general. Here we describe the phylogeography of haplogroups C and D in the populations of northern and eastern Asia. We have analyzed 770 samples from haplogroups C and D (174 and 596, respectively) at high resolution, including 182 novel complete mtDNA sequences representing haplogroups C and D (83 and 99, respectively). The present-day variation of haplogroups C and D suggests that these mtDNA clades expanded before the Last Glacial Maximum (LGM), with their oldest lineages being present in the eastern Asia. Unlike in eastern Asia, most of the northern Asian variants of haplogroups C and D began the expansion after the LGM, thus pointing to post-glacial re-colonization of northern Asia. Our results show that both haplogroups were involved in migrations, from eastern Asia and southern Siberia to eastern and northeastern Europe, likely during the middle Holocene.     

Mitochondrial DNA Variation in Two Russian Populations from Novgorod Oblast

Mitochondrial DNA (mtDNA) polymorphism was examined in two Russian populations of Novgorod oblast, from the city of Velikii Novgorod (n = 81), and the settlement of Volot (n = 79). This analysis showed that the mitochondrial gene pool of Russians examined was represented by the mtDNA types belonging to 20 haplogroups and subhaplogroups distributed predominantly among the European populations. Haplogroups typical of the indigenous populations of Asia were found in the population sample from Velikii Novgorod with the average frequency of 3.7% (haplogroups A, Z, and D5), and with the frequency of 6.3% (haplogroups Z, D, and M*) in the Volot population. It was demonstrated that the frequency of the mitochondrial lineages combination, D5, Z, U5b-16144, and U8, typical of the Finnish-speaking populations of Northeastern Europe, was somewhat higher in the urban population (7.4%) compared to rural one (3.8%). The problem of genetic differentiation of Russians from Eastern Europe inferred from mtDNA data, is discussed.     

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.     

Comprehensive association testing of common mitochondrial DNA variation in metabolic disease

Many lines of evidence implicate mitochondria in phenotypic variation: (a) rare mutations in mitochondrial proteins cause metabolic, neurological, and muscular disorders; (b) alterations in oxidative phosphorylation are characteristic of type 2 diabetes, Parkinson disease, Huntington disease, and other diseases; and (c) common missense variants in the mitochondrial genome (mtDNA) have been implicated as having been subject to natural selection for adaptation to cold climates and contributing to "energy deficiency" diseases today. To test the hypothesis that common mtDNA variation influences human physiology and disease, we identified all 144 variants with frequency >1% in Europeans from >900 publicly available European mtDNA sequences and selected 64 tagging single-nucleotide polymorphisms that efficiently capture all common variation (except the hypervariable D-loop). Next, we evaluated the complete set of common mtDNA variants for association with type 2 diabetes in a sample of 3,304 diabetics and 3,304 matched nondiabetic individuals. Association of mtDNA variants with other metabolic traits (body mass index, measures of insulin secretion and action, blood pressure, and cholesterol) was also tested in subsets of this sample. We did not find a significant association of common mtDNA variants with these metabolic phenotypes. Moreover, we failed to identify any physiological effect of alleles that were previously proposed to have been adaptive for energy metabolism in human evolution. More generally, this comprehensive association-testing framework can readily be applied to other diseases for which mitochondrial dysfunction has been implicated.     

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.     

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.     

Phylogeography of Asian sockeye salmon (Oncorhynchus nerka) based on analysis of mtDNA control region polymorphism

Analysis of the contemporary geographic distributions of mtDNA haplotypic diversity of sockeye salmon (Oncorhynchus nerka Walbaum) on the Asian part of its range and consideration of possible scenarios for microevolution of the species were performed. In total, 172 specimens of sockeye salmon from different rivers of Kamchatka, Chukotka, north coast of the Sea of Okhotsk, and from the Kuril Islands were analyzed using the entire mtDNA control region (1,010 bp) sequencing. After multiple alignments of the sequences 16 polymorphic sites were identified. Most of substitutions were detected in the first and the last thirds of the D-loop sequence, as well as one insertion was registered in the poly-T region. Two mass haplotypes (Hap_1_10T and Hap_2_13T) were revealed; they were present in most samples and differed by three substitutions and two indels. All the sequence variants were distributed among two haplogroups in compliance with the mass haplotypes. The star-shaped topology of the MST-tree points to a rapid increase in the diversity of groups, indicating a rapid expansion of the species after its number reducing in the past. The diversity of the first haplogroup was higher, but its main fraction fell on the samples from the middle reach of Kamchatka River, which seems to be a refugium during the Late Pleistocene glaciation. Presence of both haplogroups in most populations in approximately equal proportions indicates that the entire Asian part of the sockeye salmon range is a zone of secondary contact. We hypothesize that Hap_2_13T had more ancient origin and the Hap_1_10T was formed somewhat later, apparently in Asia, because almost all transitional forms between the haplogroups were found in Iturup Island populations.     

Diversity of Mitochondrial DNA Haplogroups in Ethnic Populations of the Volga–Ural Region

The mtDNA polymorphism was analyzed in eight ethnic groups (N = 979) of the Volga–Ural region. Most mtDNA variants belonged to haplogroups H, U, T, J, W, I, R, and N1 characteristic of West Eurasian populations. The most frequent were haplogroups H (12–42%) and U (18–44%). East Eurasian mtDNA types (A, B, Y, F, M, N9) were also observed. Genetic diversity was higher in Turkic than in Finno-Ugric populations. The frequency of mtDNA types characteristic of Siberian and Central Asian populations substantially increased in the ethnic groups living closer to the Urals, a boundary between Europe and Asia. Geographic distances, rather than linguistic barriers, were assumed to play the major role in distribution of mtDNA types in the Volga–Ural region. Thus, as concerns the maternal lineage, the Finno-Ugric populations of the region proved to be more similar to their Turkic neighbors rather than to linguistically related Balto-Finnish ethnic groups.     

Mitochondrial Haplogroup U5b3: A Distant Echo of the Epipaleolithic in Italy and the Legacy of the Early Sardinians

There are extensive data indicating that some glacial refuge zones of southern Europe (Franco-Cantabria, Balkans, and Ukraine) were major genetic sources for the human recolonization of the continent at the beginning of the Holocene. Intriguingly, there is no genetic evidence that the refuge area located in the Italian Peninsula contributed to this process. Here we show, through phylogeographic analyses of mitochondrial DNA (mtDNA) variation performed at the highest level of molecular resolution (52 entire mitochondrial genomes), that the most likely homeland for U5b3—a haplogroup present at a very low frequency across Europe—was the Italian Peninsula. In contrast to mtDNA haplogroups that expanded from other refugia, the Holocene expansion of haplogroup U5b3 toward the North was restricted by the Alps and occurred only along the Mediterranean coasts, mainly toward nearby Provence (southern France). From there, ∼7,000–9,000 years ago, a subclade of this haplogroup moved to Sardinia, possibly as a result of the obsidian trade that linked the two regions, leaving a distinctive signature in the modern people of the island. This scenario strikingly matches the age, distribution, and postulated geographic source of a Sardinian Y chromosome haplogroup (I2a2-M26), a paradigmatic case in the European context of a founder event marking both female and male lineages.     

Late Quaternary loss of genetic diversity in muskox (<Emphasis Type="Italic">Ovibos</Emphasis>)

Background  

The modern wildherd of the tundra muskox (Ovibos moschatus) is native only to the New World (northern North America and Greenland), and its genetic diversity is notably low. However, like several other megafaunal mammals, muskoxen enjoyed a holarctic distribution during the late Pleistocene. To investigate whether collapse in range and loss of diversity might be correlated, we collected mitochondrial sequence data (hypervariable region and cytochrome b) from muskox fossil material recovered from localities in northeastern Asia and the Arctic Archipelago of northern North America, dating from late Pleistocene to late Holocene, and compared our results to existing databases for modern muskoxen.     

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male‐biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal‐biased X‐chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.     

Holarctic genetic structure and range dynamics in the woolly mammoth

Ancient DNA analyses have provided enhanced resolution of population histories in many Pleistocene taxa. However, most studies are spatially restricted, making inference of species-level biogeographic histories difficult. Here, we analyse mitochondrial DNA (mtDNA) variation in the woolly mammoth from across its Holarctic range to reconstruct its history over the last 200 thousand years (kyr). We identify a previously undocumented major mtDNA lineage in Europe, which was replaced by another major mtDNA lineage 32–34 kyr before present (BP). Coalescent simulations provide support for demographic expansions at approximately 121 kyr BP, suggesting that the previous interglacial was an important driver for demography and intraspecific genetic divergence. Furthermore, our results suggest an expansion into Eurasia from America around 66 kyr BP, coinciding with the first exposure of the Bering Land Bridge during the Late Pleistocene. Bayesian inference indicates Late Pleistocene demographic stability until 20–15 kyr BP, when a severe population size decline occurred.     

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.     

Phylogeography and demographic history of Shaw's Jird (Meriones shawii complex) in North Africa

Palaeoenvironmental data and climatic reconstructions show that the Mediterranean ecoregion of North Africa underwent drastic ecological changes during the Pleistocene. Given its rich palaeontological record, North Africa is a pertinent region for documenting the role of climate change and human mediated‐habitat changes on the demography and genetic structure of faunal species. In the present study, we collected data from this species in Morocco, Algeria, and Tunisia, and we combined molecular (mitochondrial and nuclear DNA sequences, microsatellites), fossil, palaeoenvironmental, and human context data to propose an explanation for the fluctuations of populations belonging to the Meriones shawii complex in the past. Genetic and fossil data both indicate a strong bottleneck in Moroccan populations at the Middle Holocene (last interglacial optimum) compared to the Late Pleistocene. Our mitochondrial DNA data suggest a diversification event within Morocco corresponding to the 130–125 kya interglacial optimum. Given that (1) major demographic changes in the M. shawii complex coincide with the interglacial optimums, and (2) the impact of human activities on the landscape and faunal communities was moderate during the Middle Holocene (beginnings of the Neolithic culture), our results demonstrate that climate, rather than anthropogenic influences, likely explains the M. shawii complex population decline in the Holocene.