Search

Supplementary Information: The online version contains supplementary material available at https://doi.org/ 10.1038/s41598-021-95996-3., Sequence data for UE2298/MS060 can be downloaded from the European Nucleotide Archive (accession number: PRJEB47085). Newly reported present-day mtDNA sequences are deposited into GenBank (MZ920249 - MZ921390). Additional requests should be addressed to: marina.silva@crick.ac.uk; gonzalo.oteo-garcia@hud. ac.uk; m.b.richards@hud.ac.uk., Historical records document medieval immigration from North Africa to Iberia to create Islamic al-Andalus. Here, we present a low-coverage genome of an eleventh century CE man buried in an Islamic necropolis in Segorbe, near Valencia, Spain. Uniparental lineages indicate North African ancestry, but at the autosomal level he displays a mosaic of North African and European-like ancestries, distinct from any present-day population. Altogether, the genome-wide evidence, stable isotope results and the age of the burial indicate that his ancestry was ultimately a result of admixture between recently arrived Amazigh people (Berbers) and the population inhabiting the Peninsula prior to the Islamic conquest. We detect differences between our sample and a previously published group of contemporary individuals from Valencia, exemplifying how detailed, small-scale aDNA studies can illuminate fine-grained regional and temporal differences. His genome demonstrates how ancient DNA studies can capture portraits of past genetic variation that have been erased by later demographic shifts-in this case, most likely the seventeenth century CE expulsion of formerly Islamic communities as tolerance dissipated following the Reconquista by the Catholic kingdoms of the north., We thank the Museo Municipal de Arqueologia y Etnologia de Segorbe for granting access to their collections and the Conselleria d'Educacio, Investigacio, Cultura i Esport de la Generalitat Valenciana for granting the permissions for the study. We thank Lara Cassidy, Valeria Mattiangeli and Dan Bradley for valuable advice and technical support. Part of this work was delivered via the BBSRC National Capability in Genomics and Single Cell Analysis (BBS/E/T/000PR9816) at Earlham Institute by members of the Genomics Pipelines and Core Bioinformatics Groups. We wish to acknowledge the use of the Orion High Performance Computing cluster at the School of Applied Sciences, University of Huddersfield. M.S., G.O.G., A.Fi., P.J., M.G.B.F., K.D., B.Y. were supported by a Leverhulme Doctoral Scholarship awarded to M.B.R. and M.P. P.S., M.B.R., and M.P. acknowledge FCT (Fundacao para a Ciencia e a Tecnologia) support through project PTDC/EPH-ARQ/4164/2014, partially funded by FEDER funds (COMPETE 2020 project 016899). P.S., A.Br., A.R. and T.R. acknowledge FCT support through project PTDC/SOC-ANT/30316/2017. P.S. acknowledges the "Contrato-Programa" UIDB/04050/2020 and contract CEECINST/0007772018 funded by FCT I.P. A.A., A.O., and A.T. acknowledge the support of the Italian Ministry of Education, University and Research for the projects "Dipartimenti di Eccellenza" Program (2018-2022) -Department of Biology and Biotechnology "L. Spallanzani," University of Pavia and PRIN2017 20174BTC4R. The KORA research platform (KORA, Cooperative Research in the Region of Augsburg) was initiated and financed by the Helmholtz Zentrum Munchen -German Research Center for Environmental Health, which is funded by the German Federal Ministry of Education and Research and by the State of Bavaria. Furthermore, KORA research was supported within the Munich Center of Health Sciences (MC Health), LudwigMaximilians-Universitat, as part of LMUinnovativ.

To define the matrilineal relationships between Bulgarians and other European populations, we have evaluated the mitochondrial DNA (mtDNA) variation in a sample of 855 Bulgarian subjects from the mtDNA perspective. The molecular survey was performed by sequencing ∼750 bp of the control region, which resulted in 557 different haplotypes, and by a subsequent restriction fragment length polymorphism analysis to confirm haplogroup/subhaplogroup affiliation. The classification was carried out according to the most updated criteria as reported by van Oven and Kayser (Hum Mutat 30:386-394, 2009), allowing the identification of 45 mitochondrial clades. The observed pattern of mtDNA variation indicates that the Bulgarian mitochondrial pool is geographically homogeneous across the country, and that is characterized by an overall extremely high frequency of western Eurasian lineages. In the principal component analysis, Bulgarians locate in an intermediate position between Eastern European and Mediterranean populations, which is in agreement with historical events. Thus, while the Mediterranean legacy could be attributed to the Thracians, indigenous people that firstly inhabited the Balkans, the Eastern contribution is likely due to the Proto-Bulgarians originating from the Middle East and to the Slavs migrating from northeast Europe.

Recently, Martorell et al1 have analyzed the sequence variation of complete mitochondrial genomes in six schizophrenic patients with an apparent maternal transmission of schizophrenia. According to the authors, the patients showed new mitochondrial DNA variants related to the disease outcome. We here show, by employing a phylogeny-based approach, that the published sequences and the resulting conclusions are affected by serious shortcomings.

Haplotype and phylogenetic analyses of "normal" mitochondrial DNAs (mtDNAs) have allowed a clarification of several controversial issues concerning the origin of humans, the time and colonization pattern of the various regions of the world, and the genetic relationships of modern human populations. More recently, the same type of analyses has also been applied to mtDNA disease studies. A review of these studies indicates that exhaustive screenings of "normal" mtDNA variation in all human populations associated with haplotype and phylogenetic analyses are essential if we are to understand the etiology of mitochondrial pathologies.

Background Human genetic diversity observed in Indian subcontinent is second only to that of Africa. This implies an early settlement and demographic growth soon after the first 'Out-of-Africa' dispersal of anatomically modern humans in Late Pleistocene. In contrast to this perspective, linguistic diversity in India has been thought to derive from more recent population movements and episodes of contact. With the exception of Dravidian, which origin and relatedness to other language phyla is obscure, all the language families in India can be linked to language families spoken in different regions of Eurasia. Mitochondrial DNA and Y chromosome evidence has supported largely local evolution of the genetic lineages of the majority of Dravidian and Indo-European speaking populations, but there is no consensus yet on the question of whether the Munda (Austro-Asiatic) speaking populations originated in India or derive from a relatively recent migration from further East. Results Here, we report the analysis of 35 novel complete mtDNA sequences from India which refine the structure of Indian-specific varieties of haplogroup R. Detailed analysis of haplogroup R7, coupled with a survey of ~12,000 mtDNAs from caste and tribal groups over the entire Indian subcontinent, reveals that one of its more recently derived branches (R7a1), is particularly frequent among Munda-speaking tribal groups. This branch is nested within diverse R7 lineages found among Dravidian and Indo-European speakers of India. We have inferred from this that a subset of Munda-speaking groups have acquired R7 relatively recently. Furthermore, we find that the distribution of R7a1 within the Munda-speakers is largely restricted to one of the sub-branches (Kherwari) of northern Munda languages. This evidence does not support the hypothesis that the Austro-Asiatic speakers are the primary source of the R7 variation. Statistical analyses suggest a significant correlation between genetic variation and geography, rather than between genes and languages. Conclusion Our high-resolution phylogeographic study, involving diverse linguistic groups in India, suggests that the high frequency of mtDNA haplogroup R7 among Munda speaking populations of India can be explained best by gene flow from linguistically different populations of Indian subcontinent. The conclusion is based on the observation that among Indo-Europeans, and particularly in Dravidians, the haplogroup is, despite its lower frequency, phylogenetically more divergent, while among the Munda speakers only one sub-clade of R7, i.e. R7a1, can be observed. It is noteworthy that though R7 is autochthonous to India, and arises from the root of hg R, its distribution and phylogeography in India is not uniform. This suggests the more ancient establishment of an autochthonous matrilineal genetic structure, and that isolation in the Pleistocene, lineage loss through drift, and endogamy of prehistoric and historic groups have greatly inhibited genetic homogenization and geographical uniformity.

Background Recent advances in the understanding of the maternal and paternal heritage of south and southwest Asian populations have highlighted their role in the colonization of Eurasia by anatomically modern humans. Further understanding requires a deeper insight into the topology of the branches of the Indian mtDNA phylogenetic tree, which should be contextualized within the phylogeography of the neighboring regional mtDNA variation. Accordingly, we have analyzed mtDNA control and coding region variation in 796 Indian (including both tribal and caste populations from different parts of India) and 436 Iranian mtDNAs. The results were integrated and analyzed together with published data from South, Southeast Asia and West Eurasia. Results Four new Indian-specific haplogroup M sub-clades were defined. These, in combination with two previously described haplogroups, encompass approximately one third of the haplogroup M mtDNAs in India. Their phylogeography and spread among different linguistic phyla and social strata was investigated in detail. Furthermore, the analysis of the Iranian mtDNA pool revealed patterns of limited reciprocal gene flow between Iran and the Indian sub-continent and allowed the identification of different assemblies of shared mtDNA sub-clades. Conclusions Since the initial peopling of South and West Asia by anatomically modern humans, when this region may well have provided the initial settlers who colonized much of the rest of Eurasia, the gene flow in and out of India of the maternally transmitted mtDNA has been surprisingly limited. Specifically, our analysis of the mtDNA haplogroups, which are shared between Indian and Iranian populations and exhibit coalescence ages corresponding to around the early Upper Paleolithic, indicates that they are present in India largely as Indian-specific sub-lineages. In contrast, other ancient Indian-specific variants of M and R are very rare outside the sub-continent.

The debate concerning the mechanisms underlying the prehistoric spread of farming to Southeast Europe is framed around the opposing roles of population movement and cultural diffusion. To investigate the possible involvement of local people during the transition of agriculture in the Balkans, we analysed patterns of Y-chromosome diversity in 1206 subjects from 17 population samples, mainly from Southeast Europe. Evidence from three Y-chromosome lineages, I-M423, E-V13 and J-M241, make it possible to distinguish between Holocene Mesolithic forager and subsequent Neolithic range expansions from the eastern Sahara and the Near East, respectively. In particular, whereas the Balkan microsatellite variation associated to J-M241 correlates with the Neolithic period, those related to E-V13 and I-M423 Balkan Y chromosomes are consistent with a late Mesolithic time frame. In addition, the low frequency and variance associated to I-M423 and E-V13 in Anatolia and the Middle East, support an European Mesolithic origin of these two clades. Thus, these Balkan Mesolithic foragers with their own autochthonous genetic signatures, were destined to become the earliest to adopt farming, when it was subsequently introduced by a cadre of migrating farmers from the Near East. These initial local converted farmers became the principal agents spreading this economy using maritime leapfrog colonization strategies in the Adriatic and transmitting the Neolithic cultural package to other adjacent Mesolithic populations. The ensuing range expansions of E-V13 and I-M423 parallel in space and time the diffusion of Neolithic Impressed Ware, thereby supporting a case of cultural diffusion using genetic evidence.

Recent advances in the understanding of the maternal and paternal heritage of south and southwest Asian populations have highlighted their role in the colonization of Eurasia by anatomically modern humans. Further understanding requires a deeper insight into the topology of the branches of the Indian mtDNA phylogenetic tree, which should be contextualized within the phylogeography of the neighboring regional mtDNA variation. Accordingly, we have analyzed mtDNA control and coding region variation in 796 Indian (including both tribal and caste populations from different parts of India) and 436 Iranian mtDNAs. The results were integrated and analyzed together with published data from South, Southeast Asia and West Eurasia. Four new Indian-specific haplogroup M sub-clades were defined. These, in combination with two previously described haplogroups, encompass approximately one third of the haplogroup M mtDNAs in India. Their phylogeography and spread among different linguistic phyla and social strata was investigated in detail. Furthermore, the analysis of the Iranian mtDNA pool revealed patterns of limited reciprocal gene flow between Iran and the Indian sub-continent and allowed the identification of different assemblies of shared mtDNA sub-clades. Since the initial peopling of South and West Asia by anatomically modern humans, when this region may well have provided the initial settlers who colonized much of the rest of Eurasia, the gene flow in and out of India of the maternally transmitted mtDNA has been surprisingly limited. Specifically, our analysis of the mtDNA haplogroups, which are shared between Indian and Iranian populations and exhibit coalescence ages corresponding to around the early Upper Paleolithic, indicates that they are present in India largely as Indian-specific sub-lineages. In contrast, other ancient Indian-specific variants of M and R are very rare outside the sub-continent.