Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia
Total Page:16
File Type:pdf, Size:1020Kb
RESEARCH ARTICLE Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia Sara De Fanti1☯, Chiara Barbieri1☯¤*, Stefania Sarno1, Federica Sevini2,3, Dario Vianello2,3, Erika Tamm4,5, Ene Metspalu4,5, Mannis van Oven4,6, Alexander Hübner7, Marco Sazzini1, Claudio Franceschi2,3,8,9, Davide Pettener1, Donata Luiselli1* 1 Department of Biological, Geological and Environmental Sciences, University of Bologna, Bologna, Italy, 2 Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, Bologna, Italy, 3 C.I.G. Interdepartmental Centre L. Galvani for Integrated Studies on Bioinformatics, Biophysics and Biocomplexity, University of Bologna, Bologna, Italy, 4 Estonian Biocentre, Evolutionary Biology group, Tartu, Estonia, 5 Department of Evolutionary Biology, University of Tartu, Tartu, Estonia, 6 Department of Forensic Molecular Biology, Erasmus MC - University Medical Center Rotterdam, Rotterdam, The Netherlands, 7 Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany, 8 IRCCS, Institute of Neurological Sciences of Bologna, Ospedale Bellaria, Bologna, Italy, 9 CNR, Institute of Organic Synthesis and Photoreactivity (ISOF), Bologna, Italy OPEN ACCESS ☯ These authors contributed equally to this work. ¤ Current Address: Department of Linguistic and Cultural Evolution, Max Planck Institute for the Science of Citation: De Fanti S, Barbieri C, Sarno S, Sevini F, Human History, Jena, Germany Vianello D, Tamm E, et al. (2015) Fine Dissection of * [email protected] (CB); [email protected] (DL) Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia. PLoS ONE 10(12): e0144391. doi:10.1371/ journal.pone.0144391 Abstract Editor: Ron Pinhasi, University College Dublin, Genetic signatures from the Paleolithic inhabitants of Eurasia can be traced from the early IRELAND divergent mitochondrial DNA lineages still present in contemporary human populations. Received: April 15, 2015 Previous studies already suggested a pre-Neolithic diffusion of mitochondrial haplogroup Accepted: November 17, 2015 HV*(xH,V) lineages, a relatively rare class of mtDNA types that includes parallel branches Published: December 7, 2015 mainly distributed across Europe and West Asia with a certain degree of structure. Up till now, variation within haplogroup HV was addressed mainly by analyzing sequence data Copyright: © 2015 De Fanti et al. This is an open access article distributed under the terms of the from the mtDNA control region, except for specific sub-branches, such as HV4 or the widely Creative Commons Attribution License, which permits distributed haplogroups H and V. In this study, we present a revised HV topology based on unrestricted use, distribution, and reproduction in any full mtDNA genome data, and we include a comprehensive dataset consisting of 316 com- medium, provided the original author and source are credited. plete mtDNA sequences including 60 new samples from the Italian peninsula, a previously underrepresented geographic area. We highlight points of instability in the particular topol- Data Availability Statement: All the newly generated 55 sequences are available in GenBank ogy of this haplogroup, reconstructed with BEAST-generated trees and networks. We also (http://www.ncbi.nlm.nih.gov/genbank/) with confirm a major lineage expansion that probably followed the Late Glacial Maximum and accession codes KP340126-KP340180. preceded Neolithic population movements. We finally observe that Italy harbors a reservoir Funding: This work was supported by the PRIN of mtDNA diversity, with deep-rooting HV lineages often related to sequences present in the Progetti di Ricerca di Interesse Nazionale Caucasus and the Middle East. The resulting hypothesis of a glacial refugium in Southern 2010EL8TXP_006 grant to DL, http://prin.miur.it/; the Italy has implications for the understanding of late Paleolithic population movements and is European Research Council ERC-2011-AdG 295733 (Langelin) to DP and CB, http://erc.europa.eu/; the discussed within the archaeological cultural shifts occurred over the entire continent. Institutional Research Funding from the Estonian Research Council [IUT24-1], http://www.etag.ee/en/ estonian-research-council/; and by the European PLOS ONE | DOI:10.1371/journal.pone.0144391 December 7, 2015 1/19 Fine Dissection of Human Mitochondrial DNA Haplogroup HV Regional Development Fund (European Union), Introduction http://ec.europa.eu/regional_policy/index.cfm/en/ funding/erdf/, through the Centre of Excellence in In the past years, an increasing number of studies reported whole mitochondrial genome data Genomics to Estonian Biocentre and University of from various human populations from all the continents, providing new insights into the Tartu, http://genomics.ebc.ee/. MvO was supported in shape of the matrilineal phylogeny and the distribution of its clades. These fine-grained studies part by an Estonian Biocentre EBC ECOGENE helped to both understand demic migration flows and disentangle the phylogeographic distri- student fellowship, https://edukad.etag.ee/project/ bution and the time of divergence of characteristic mitochondrial DNA (mtDNA) lineages [1– 3010. The funders had no role in study design, data collection and analysis, decision to publish, or 8]. Uniparental genetic markers (i.e. mtDNA and the non-recombinant portion of the Y chro- preparation of the manuscript. mosome) have been widely studied for reconstructing the prehistory of European and Mediter- ranean populations from living individuals [9–13], as well as from ancient DNA (aDNA) Competing Interests: The authors have declared that no competing interests exist. samples [13]. Accordingly, the available published data from these regions provide an extensive geographical coverage, even though at a shallow level of phylogenetic resolution. In fact, until a couple of years ago, population-level analyses of mtDNA focused almost exclusively on the hypervariable segment(s) (HVS) located within the D-loop or control region. Only recently, this kind of studies began to include full mtDNA genome data to explore the variability of spe- cific lineages at a higher resolution and to get deeper insights into the matrilineal prehistory of a given population [5, 14–16]. A recent comprehensive survey of the Italian uniparental genetic landscape, comprising almost a thousand of individuals, revealed that mtDNA haplogroup HV is the lineage in Italy with the oldest coalescence, and highlighted a divergent structure between the Northern and the Southern regions of the peninsula [17]. This early signal can be paralleled with noticeable findings obtained by aDNA analysis on remains from a Mesolithic site in Sicily (Favignana, 14 kilo years ago, kya) [18]), which included a specimen assigned to haplogroup HV1. A Paleo- lithic site in Puglia (Paglicci, 28 kya) also reported data from an individual assigned to HVÃ (or pre-HV/R0) [19], but such finding could not be replicated in a more recent study [20]. While haplogroup HV was recognized as a crucial component of early human dispersal in Eurasia [9,21], patterns of its internal variability have been poorly investigated in previous works, which only focused on the major clades HV4 [4], V and H [21–27]. Haplogroup HV is a major mitochondrial clade within haplogroup R0 [9] characterized by the T14766C mutation and comprising at least 18 recognized basal subclades. Seventeen of these are designated with a consecutive numerical label, i.e. HV0-HV17 (the HV3 label is not used) (http://www.phylotree.org/, Build 16) [28]. Additionally, it includes haplogroup H as one of its direct subclades, while haplogroup V is a nested subhaplogroup within the HV0 clade. The basal structure of the HV phylogeny is partially characterized by mutations located within the hypervariable segments. Some of these branch-defining variants occur for example at nucleotide positions 72, 73, 152, 195, 16311, which are recognized as recurrent sites throughout the whole mtDNA phylogeny—the latter three for example, appears more than 80 times in the tree [29]—and could therefore obscure the topology of the reconstructed HV tree. The HV clade as a whole, including its H subclade, encompasses Eurasian haplogroups that likely arose between Western and Central Asia [30]. Haplogroup HVÃ(xH,V) (i.e. the whole of HV excluding its subclades H and V) is not particularly common in Europe, with frequencies that range between 0% and 10%. In particular, high HVÃ occurrences are observed in Southern Europe (e.g., Italy and Spain) and exceptionally high frequencies are found in Iran (19–24%) [31], although somewhat lower values (9–14%) are reported for the same populations in Farjadian et al. [32]. An average frequency of 4.05% is observed in the Italian peninsula [17]. On the other hand, the major HV subclade, haplogroup H, experienced a vast diffusion during the late Neolithic time, becoming the most common haplogroup in Europe today [27,33].Otherlineageswithinhap- logroup HV were likely spread in more ancient times, dating back to more than 15 kya and having independent diffusions within Europe of Early, Middle or Late Upper Paleolithic origins [9,21]. PLOS ONE | DOI:10.1371/journal.pone.0144391 December 7, 2015 2/19 Fine Dissection of Human Mitochondrial DNA Haplogroup HV Several scholars associated the diffusion of Paleolithic characteristic lineages with the re- peopling