Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya origin in the
conquering Hungarians
plos.org
PUBLISH
ABOUT
BROWSE
create account
sign in
SEARCH
advanced search
OPEN ACCESS
PEER-REVIEWED
RESEARCH ARTICLE
Mitogenomic data indicate admixture components of
Central-Inner Asian and Srubnaya origin in the conquering
Hungarians
Endre Neparáczki
, Zoltán Maróti
Loading metrics
, Tibor Kalmár, Klaudia Kocsy, Kitti Maár, Péter Bihari, István Nagy, Erzsébet Fóthi,
Ildikó Pap, Ágnes Kustár, György Pálfi, István Raskó, Albert Zink, Tibor Török
Published: October 18, 2018
Article
https://doi.org/10.1371/journal.pone.0205920
Authors
Metrics
Comments
Media Coverage
Print
Reader Comments (0)
Media Coverage
Figures
Download PDF
Share
Correction
29 Nov 2018: Neparáczki E, Maróti Z, Kalmár T, Kocsy K, Maár K, et al. (2018)
Correction: Mitogenomic data indicate admixture components of Central-Inner Asian and
Srubnaya origin in the conquering Hungarians. PLOS ONE 13(11): e0208295.
https://doi.org/10.1371/journal.pone.0208295 | View correction
Figures
Related PLOS Articles
Correction: Mitogenomic data
indicate admixture
components of Central-Inner
Asian and Srubnaya origin in
the conquering Hungarians
ADVERTISEMENT
Abstract
It has been widely accepted that the Finno-Ugric Hungarian language, originated from proto
Uralic people, was brought into the Carpathian Basin by the conquering Hungarians. From the
middle of the 19th century this view prevailed against the deep-rooted Hungarian Hun tradition,
maintained in folk memory as well as in Hungarian and foreign written medieval sources, which
claimed that Hungarians were kinsfolk of the Huns. In order to shed light on the genetic origin of
the Conquerors we sequenced 102 mitogenomes from early Conqueror cemeteries and
compared them to sequences of all available databases. We applied novel population genetic
algorithms, named Shared Haplogroup Distance and MITOMIX, to reveal past admixture of
maternal lineages. Our results show that the Conquerors assembled from various nomadic
groups of the Eurasian steppe. Population genetic results indicate that they had closest
connection to the Onogur-Bulgar ancestors of Volga Tatars. Phylogenetic results reveal that
more than one third of the Conqueror maternal lineages were derived from Central-Inner Asia
and their most probable ultimate sources were the Asian Scythians and Asian Huns, giving
support to the Hungarian Hun tradition. The rest of the lineages most likely originated from the
Bronze Age Potapovka-Poltavka-Srubnaya cultures of the Pontic-Caspian steppe. Available data
imply that the Conquerors did not have a major contribution to the gene pool of the Carpathian
Basin.
Citation: Neparáczki E, Maróti Z, Kalmár T, Kocsy K, Maár K, Bihari P, et al. (2018)
Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya
origin in the conquering Hungarians. PLoS ONE 13(10): e0205920.
https://doi.org/10.1371/journal.pone.0205920
Editor: David Caramelli, University of Florence, ITALY
Received: May 9, 2018; Accepted: October 2, 2018; Published: October 18, 2018
Copyright: © 2018 Neparáczki et al. This is an open access article distributed under the
terms of the Creative Commons Attribution License, which permits unrestricted use,
distribution, and reproduction in any medium, provided the original author and source are
credited.
Data Availability: Data are available at the European Nucleotide Archive
(http://www.ebi.ac.uk/ena) under accession number PRJEB21279.
Funding: This work was supported by grants from the National Research, Development
and Innovation Office (K-124350 to TT), Bolyai Research Scholarship of the Hungarian
Academy of Sciences (to IN), and the Avicennna Middle East Research Institute
(GF/JSZF, Award Number: 814/9/2015 to IR). The Bolyai Research Scholarship provided
support in the form of salaries for one of the authors (IN). The specific roles of this author
are articulated in the ‘author contributions’ section. The funders had no role in study
design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Dr. Bihari and Dr. Nagy had consulting positions at SeqOmics
Biotechnology Ltd. at the time the study was conceived. SeqOmics Biotechnology Ltd.
was not directly involved in the design and execution of the experiments or in the writing
of the manuscript. This affiliation does not alter our adherence to PLOS ONE policies on
sharing data and materials.
Introduction
Foundation of the Hungarian state is connected to the conquering Hungarians, which arrived
from the Pontic steppes and occupied the Carpathian Basin at 895–905 AD as a confederation of
seven tribes under the leadership of prince Árpád. Modern Hungarians are generally identified as
successors of the conquering Hungarians (hence shortened as Conquerors). Until the middle of
the 19th century it was generally accepted that Hungarians were kinsfolk of the Huns and
Scythians, besides Árpád was a direct descendant of the great Hun leader Attila. Hun-Hungarian
affinity was declared in Hungarian and foreign written sources and has been maintained in
Hungarian folk memory [1–3]. In the second half of the 19th century the Hungarian language was
reclassified as belonging to the Uralic branch of the Finno-Ugric language family [4]. Philological
arguments launched a reevaluation of previous assumptions and as a result, the credibility of
medieval historical sources, including Hun-Hungarian relations, has been questioned. In
following decades the conquering Hungarians were deemed descendants of hypothetic proto
Uralic people, the putative common ancestors of people belonging to this language family.
Lately most philologists proclaim separability of linguistic and genetic relations, but appearance
of the Hungarian language in the Carpathian Basin is explicitly linked to the Conquerors [5].
The possible genetic relation of modern Hungarians to Finno-Ugric groups was tested in several
studies [6–8], however all these found Hungarians being genetically unrelated to Uralic people.
One of the latest studies [9] reported that a Y-chromosome haplogroup ( N-L1034) is shared
between 4% of the Hungarian Seklers (Hungarian-speaking ethnic group living in Transylvania)
and 15% of the closest language relatives the Mansis, though the same marker is also present in
Central Asian Uzbeks and has been detected just in one Hungarian [10]. These results indicated
that Uralic genetic links hardly exist in modern Hungarians.
The genetic composition of the Conquerors was also analyzed in several ancient DNA (aDNA)
studies [11–13] and indeed, all these detected significant presence of east Eurasian major
mtDNA haplogroups (Hg-s), which are rare in modern Hungarians but are found in Uralic people.
Another study [14] showed the presence of N-Tat (M46) Y-chromosome marker (a major clade
of the above mentioned N-L1034) in two of the Conqueror samples and one living Sekler, which
was interpreted as a Finno-Ugric link. It is notable that in the latest studies [12,13] population
genetic analysis also indicated considerable Central Asian affinity of the Conquerors. However in
these studies just hypervariable regions (HVR) of the mtDNA were analyzed, and more reliable
Next Generation Sequencing (NGS) data [15] have not been available from the Conquerors yet.
Entire mitochondrial genome sequences enable a much higher resolution analysis, as most
Subject Areas
Europe
Population genetics
Haplogroups
Mitochondrial DNA
Paleogenetics
Asia
Sequence databases
Sequence alignment
?
variable sites of mtDNA are located outside HVR [16].
In order to elicit the genetic origin and relationships of the Conquerors, we set out to assemble a
full length mtDNA sequence database from the earliest Conqueror cemeteries. Full length
mitogenomes are the most informative source of maternal population histories, as some of the
subclades have very distinctive geographic distribution [17,18], reviewed in [ 19]. Thus the
availability of ancient mitogenomes obtained with NGS greatly enhanced the resolution of the
phylogeographic approach, making it possible to refine the view of peopling of the Americas [20]
and Europe [21–23]. We also made use of this approach by comparing the mtDNA genomes of
102 Conqueror individuals to available public databases. Applying phylogenetic analysis we
could allocate the presumptive geographical origin of individual Conqueror Hg lineages to distant
regions of East and West Eurasia, while population genetic results pointed at source populations
in Volga district, today’s Belarus, Tuva and Central Asia, providing new information about the
origin of the Conquerors which is reconcilable with historical sources.
Materials and methods
Archaeological background
In the 10 th century a uniform well distinguishable new archaeological culture appeared in the
Carpathian Basin which can be connected to the historical record of the conquering Hungarians.
We extracted ancient DNA from 102 Conqueror individuals, derived from 8 different cemeteries
(Fig 1).
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Fig 1. Location of the Hungarian Conqueror cemeteries.
Red dots indicate cemeteries reported in this study, blue dots indicate cemeteries from
which HVR sequences were reported in [11,12]. Numbers indicate the following sites: 1.
Karos-Eperjesszög, 2. Kenézlő-Fazekaszug-II, 3. Harta-Freifelt, 4. Magyarhomoróg, 5.
Orosháza-Görbicstanya, 6. Szabadkígyós-Pálliget, 7. Sárrétudvari-Hízóföld, 8. SzegvárOromdűlő, 9. Balatonújlak-Erdődűlő, 10. Levice-Géňa, 11. Kiskundorozsma-Hosszúhát,
12. Baks-Iskola, 13. Szeged-Öthalom, 14. M43 no. 25 site Makó-Igási járandó, 15.
Szentes-Derekegyháza, 16. Nyíregyháza-Oross Megapark, 17. Kiszombor, 18. IzsákBalázspuszta, 19. Aldebrő-Mocsáros, 20. Besenyőtelek-Szőrhát, 21. EgerSzépasszonyvölgy, 22. Fadd-Jegeshegy, 23. Mözs-Szárazdomb, 24. Örménykút, 25.
Zalavár-Kápolna, 26. Lébény-Kaszás. Map was created with the maps package of R [24].
https://doi.org/10.1371/journal.pone.0205920.g001
As one of our purposes was to characterize the entire population from a few early Conqueror
cemeteries, the majority of samples came from three cemeteries of Karos-Eperjesszög,
representing the earliest Conquerors in the Carpathian Basin. These three cemeteries are
located in the upper Tisza river region on neighboring sand dunes a few 100 meters from each
other, with the richest archaeological findings of the period, and were probably used by
contemporary neighboring communities from the last years of the ninth century to the middle of
the tenth century, based on dating with coins and comparative analysis of archaeological
findings [25]. Basic archaeological description of the cemeteries were given in [ 13], further details
are provided in S1 Appendix and S1A Table.
Specimen numbers are the same as complete repository information, which correspond to the
sample names provided in the paper, summarized in S1A Table. Sample geographic locations
are provided in S1 Appendix and Fig 1. Name of permit issuing authority for this study:
Department of Anthropology; Hungarian Natural History Museum; Budapest, H-1083; Hungary
Permit number: M-2017-006.
NGS sequencing
Details of the aDNA purification, hybridization capture, sequencing and sequence analysis
methods are given in [15] and were deposited to http://dx.doi.org/10.17504/protocols.io.rmvd466.
In order to authenticate the results, we considered the latest recommendations of [26] throughout
of the experiments. We tried to apply the modifications recommended by [27] on a few samples
(Kenézlő-Fazekaszug/ 1027, 1044, 1045 and 10936, Sárrétudvari-Hízóföld/ 66 and 103), but in
our hands this method gave rather varying coverage. In some of the samples (Karos2/2, 17, 18,
33, 44, 67, Karos3/7, 9, 11, 13, 17, 18, Sárrétudvari-Hízóföld/9-anc11, Sárrétudvari-H/81, 136
and Kenézlő-F/1025, 1031, 1036, 1041, 1042) we decreased the recommended USER and UGI
concentrations of [28] to half (0.03 U/μL) and at the same time increased the incubation time
from 30 to 40 minutes. This modification removed uracils with comparable efficiency to the
original method.
Details of NGS data are shown in S2A and S2B Table . Most genomes had satisfactory coverage,
but we also included several low coverage sequences, whose Hg-s could be unmistakably
classified, as these revealed meaningful maternal relationships within and between cemeteries.
Contamination was estimated with two methods; a) using the Schmutzi algorithm and b)
calculating the proportion of reads which did not correspond to the consensus sequence in
diagnostic positions as in [15], the two methods gave consistent results. The raw nucleotide
sequence data of the 102 samples were deposited to the European Nucleotide Archive
(http://www.ebi.ac.uk/ena) under accession number PRJEB21279.
Phylogenetic study
We have downloaded all available modern (n = 32683) and ancient (n = 564, S3C and S3D
Table) complete mtDNA genome sequences from the NCBI and European Nucleotide Archive
databases or requested them from the authors. This database was also augmented with 314
mitogenomes including 272 new Hungarian ones described in [29]. Then we determined the
haplogroups of all sequences with the HaploFind program [30], and arranged them according to
haplogroups. Next we selected each subset of sequences (28-180/Hg) corresponding to the Hg
of individual Conqueror samples. Selected sequence subsets were aligned with MAFFT version
7 [31,32] using progressive G-INS-1 setting. Aligned multifasta groups were converted into
Nexus file with MEGA [33], then Median-Joining networks [ 34] were drawn with PopART [ 35].
Finally phylogeographic connections were inferred by looking up the geographic origin of the
closest matching samples from the literature (S1 Fig).
Population genetic study
We have created an Eurasian population database by grouping those mtDNA genomes
according to their geographic origin, for which this information was available (S3A Table). Our
population database contains 12224 modern samples from 62 Eurasian populations (S3B
Table), not considering India and Southeast Asia. In cases when populations were
underrepresented we grouped related neighboring groups, like Mansis with Khantys, Belgians
with Dutch etc., as listed in S3B Table. We also created a similar mitogenomic population
database from 25 ancient Eurasian populations including 496 sequences, though most of these
contain low number of samples (S3C Table).
We compared the genetic similarity of populations with two independent methods. We applied
the traditional sequence based method calculating pair-wise population differentiation values
(Fst) with Arlequin 3.5.2.2 [36] from entire mtDNA genomes ( S4A Table) assuming a Tamura &
Nei substitution model (Tamura and Nei, 1993) with a gamma value of 0.325. Significant
variations in Fst values were tested by 10,000 permutations between populations. As individual
insertions and deletions make the alignment of multiple mtDNA genomes troublesome, only
variable positions were aligned, and insertions and deletions were recoded to SNP-s as follows.
Whole mtDNA genome fasta files were aligned to the NC_012920 human mtDNA reference
sequence by an IUPAC code aware in-house aligner using the Needleman–Wunsch algorithm
with weight parameters: match 6, IUPAC2match (R, Y, M, W, S, K) 3, IUPAC3match (B, D, H, V)
2, IUPAC4match (N) 1, mismatch -12, gap open -24, gap extend -6. Modern sequences with
more than 500 missing or uncertain nucleotides (nt.) were excluded from further analysis. Then
all nt. positions where any variation was detected were outputted to VCF files. Since Arlequin
cannot manage VCF files SNPs, deletions and insertions were recoded by the following rules: nts with no variation at the given position were coded as the reference nt.; SNPs with variation
were coded as the alternate allele; all insertions were coded as additional nt. letters, C for
samples with reference sequence and T for samples containing the insertion; all deletions were
also coded as additional nt. letters, T for samples with reference sequence and C for samples
containing the deletion. Then Arlequin input files (arp) were generated from the recoded DNA
sequences.
Multidimensional scaling (MDS) was applied on the matrix of linearized Slatkin Fst values [ 37]
and visualized in the two-dimensional space using the cmdscale function implemented in R 3.0.3
[38].
In a second novel approach we also calculated so called Shared Haplogroup Distance (SHD)
values between populations [29]. This method considers that all individuals within the same subHg were descended from a single foremother, therefore their maternal lineages are more closely
related to each-other than to individuals of neighboring sub-Hg-s. Thus presence of identical
terminal subgroups in two populations testifies shared ancestry or past admixture. While Fst
based calculations are best suited for measuring evolutionary distances between not admixing
populations, we demonstrate that SHD based distance reveals recent admixtures more
accurately. We also show that SHD (and MITOMIX see below) results are in accord with Fst
calculations. SHD calculations give a distance value between 0–1, which is minimum between
populations containing the same sub-haplogroups with identical frequencies, and maximum
between populations with no sub-Hg overlap. We used corrected SHD vales, which also takes
into account the mutation and fixation rate on the mtDNA genome, thereby allows some
connection between progenitor and progeny Hg lineages [29]. Pair-wise SHD distances were
calculated between all 87 ancient and modern populations from the frequency of 1942 sub-Hg-s
occurring in any of them (S4B Table).
As an additional benefit, SHD enables a hypothesis independent computation to reveal plausible
past admixture events. Thus we have also introduced another novel algorithm called MITOMIX,
which computes all possible combinations and proportions of K populations to find the best
fitting admixtures with the smallest SHD values from a test population [29]. In our experience
most test populations are adequately admixed from 3–6 other populations, as K values greater
than 6 do not significantly improve the result. Theoretically MITOMIX can accurately reconstruct
past population admixtures if representative data are available from all periods and locations, but
it allows meaningful insights even from limited data [29]. With this method we have calculated the
best population admixtures giving the most similar mitogenome composition to that of the
Conquerors, as well as for their possible source populations (S5 Table).
Craniofacial reconstructions
The sculpting craniofacial reconstructions of three skulls from the Karos cemeteries ( Fig 2) were
carried out by Gyula Skultéty in cooperation with the Hungarian Natural History Museum,
Department of Anthropology [39,40]. During the facial reconstruction, soft tissue layers were
grafted back onto the plaster copy of the skulls carefully following the bone conformation to
accurately recreate the facial features according to the published guidelines [41–43]. Facial
reconstruction was performed by traditional sculpting anatomy, that is plasticine muscles were
attached in their anatomically correct position [40,44]. The width of a muscle was determined by
the ruggedness of the bone surface by means of a table compiled from measurements taken
from 45 different points of the skull. These data have been collected by scientific methods [45].
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Fig 2. Skulls and sculpting craniofacial reconstructions of Hungarian Conqueror
individuals.
A: Karos2/52 mature aged leader with Europid anthropological features. B: Karos2/60
senile aged man with Europo-Mongoloid features. C: Karos2/47 adult woman with
Europo-Mongoloid features.
https://doi.org/10.1371/journal.pone.0205920.g002
Results
Phylogenetic study
Using the NGS sequencing method combined with target enrichment, we could obtain 102
ancient mitogenome sequences, 78 of which are first reported in this paper, while 24 had been
reported in [15]. The 102 sequences belong to 67 sub-Hg-s, and first we elucidated the
phylogenetic relations of each Hg-s using M-J Networks as shown in S1 Fig. The closest
sequence matches pointed at a well-defined geographical region in most cases, which is
indicated next to the phylogenetic trees and is summarized on Fig 3.
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Fig 3. Phylogeographic origin of the 102 Conqueror maternal lineages.
Data are summarized from S1 Fig. Origin of modern individuals with closest matches to
Conqueror sequences are listed next to the indicated regions, ordered according to the
frequency of appearances.
https://doi.org/10.1371/journal.pone.0205920.g003
Phylogenetic trees revealed that the Conqueror maternal lineages originated from two distant
geographical regions; 31 were unequivocally derived from East Eurasia, while 60 from West
Eurasia. The remaining 11 Conqueror Hg-s are ubiquitous in Eurasia. Out of the 60 west
Eurasian lineages 13 are characteristic for modern Northwestern Europeans, while 7 have
primarily Caucasus-Middle-East distribution.
As high similarity of mitogenomes infer recent common maternal ancestor, sequence correlation
levels provide important phylogeographic information. Origin of modern individuals with closest
matches to Conqueror sequences are listed on Fig 3. We detected a very prominent frequency
of Hg N1a1a1a1a, represented by 7 Conqueror samples, while two more samples belonged to
the progenitor N1a1a1a1 lineage (S1 Fig; Network 36). N1a1a1a1a has Central Asian origin, as
its current distribution is restricted to Kazakhstan, Altai, Buryat Republic and Russia, attesting
that these areas were the center of expansion [46]. This Hg was detected in a Bronze Age
Sintashta sample from Kazakhstan [47], an Iron Age Pazyryk Scythian [ 48] and an early
Sarmatian sample [49], while its progenitor Hg N1a1a1a1 has a wide Eurasian distribution [ 46].
Our phylogeographic data imply a probable expansion of N1a1a1a1 from the European Pontic
Steppe to Central Asia around the Bronze Age and its sub-clade N1a1a1a1a from Central Asia
both to Inner Asia and back to Europe from the Iron Age. Besides many of the sequences
showed close matches with samples from ancient cultures (S1 Fig), the most recurrent being the
Bronze Age Srubnaya (Timber-grave) culture with 6 closely related Conqueror sequences. More
than one related sequences were found to samples from Neolith-Bronze Age Hungary,
Yamnaya-Eneolith Samara, Armenian Neolith-Bronze Age.
Distribution of the closest east Eurasian sequence matches outlines a well-defined geographic
region (Fig 4, red heat map) centering around modern Buryatia-Northern Mongolia, with some
extension through Tuva into Central Asia, an area well corresponding to the center and range of
the ancient Asian Hun (Xiongnu) Empire especially considering that Yakuts, Evenks and Evens
lived more south in the past [50].
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Fig 4. The most feasible origin and migration route of different components of the
Hungarian Conquerors based on this study.
Red heat map displays the geographic distribution of closest East Eurasian sequence
matches to individual Conqueror samples. Stars denote geolocations of East Eurasian
ethnic groups listed on S1 Fig (summarized on S1A Table), map was drawn from their
frequency of occurence. Heat map designate the area from which the East Eurasian
lineages most likely originated, well corresponding to the range of the ancient Xiongnu
Empire outlined by dashed line. Areas where Asian and European Scythian remains were
found are labeled green. Asian Scythians around Tuva correspond to the most probable
sources of Eurasian lineages. Pink label shows the presumptive range of the Srubnaya
culture, from where European lineages were most likely derived. Bluish line frames the
Eurasian steppe zone, within which all presumptive ancestors of the Conquerors were
found. The map was created using QGIS 2.18.4[51].
https://doi.org/10.1371/journal.pone.0205920.g004
Even though phylogenetic analysis may indicate the ultimate source region of individual maternal
lineages, but these together do not necessarily correspond to actual populations, which has to
be studied by population genetic methods.
Population genetic study
As the studied samples apparently represent real Conqueror populations we measured their
genetic distances from all recent and ancient populations. For increasing the resolution of the
method we compared mitogenomes of populations, albeit this inherently reduces sample
representativeness. Besides the traditional Fst distance calculations we used a novel approach
[29], which calculates so called Shared Haplogroup Distances (SHD). The simple logic behind
SHD is that sub-Hg-s originated from a single most recent common ancestor, thus presence of
identical subgroups links population histories in an extent of sharing, which is proportional to the
SHD value. Both pair-wise distance matrices are shown in S4A and S4B Table . The Fst and
SHD methods gave comparable results (S4C Table), thus close distance values measured with
both methods can be considered very plausible relationships. Latter populations are
summarized on Table 1 and the MDS plot from linearized Slatkin Fst values of this subset is
displayed on Fig 5.
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Fig 5. MDS plot from linearized Slatkin Fst values of S4A Table.
Only populations from Table 1 were depicted, which showed close Fst and SHD distance
values to the Conquerors. Abbreviations of population names are given in S3B Table.
https://doi.org/10.1371/journal.pone.0205920.g005
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Table 1. Fst and SHD distances of modern (rec) and ancient (arch) populations
measured from different Conqueror subpopulations.
https://doi.org/10.1371/journal.pone.0205920.t001
From the Conqueror population Volga Tatars have the smallest overall distance with both
methods (Table 1, all 102 samples), and accordingly they are positioned very close on the MDS
plot implying closest genetic relation at the population level.
As for the further analysis it is important to point at important differences between the Fst and
SHD methods, which is illuminated by their different patterns on Table 1 and S4C Table. Fst is
best suited to measure genetic distance between isolated populations where the effect of
evolutionary sequence divergence is determining [29] and in case of population admixture it
shall identify populations with similar admixture patterns or at best the predominant source. In
contrast SHD is capable of indicating admixing sources and its value is proportional to admixing
ratios [29]. The presence of 30% east Eurasian and 60% west Eurasian Hg-s in the Conquerors
is a clear indication of past admixtures whose sources should be identified. Besides Volga Tatars
the Fst similarity list includes exclusively west Eurasian populations (Table 1) most likely related
to the majority admixture components. In contrast SHD clearly identifies potential east Eurasian
admixture sources; Buryats (Bur), Central Asians (CenA), Mongolians (Mon) and Tuvans (Tuv).
The efficacy of the SHD method can be demonstrated by an artificial partition of the Conqueror
population into east and west Eurasian components (S1B Table) and performing the same
analysis with each component. Though this grouping correspond just to speculative populations
not real ones, the results clearly confirm our above claims. The entire Conqueror population
shows nearly identical Fst distance patterns to that of its dominant European subset (Table 1),
irrespective of the presence of the 11 Eurasian lineages and accordingly the European and
European+Eurasian Conqueror subpopulations map very close on the MDS plot (Fig 5). On the
other hand now both methods identify the same east Eurasian populations close to the
Conqueror east Eurasian subsets which had been identified within the entire population just by
SHD, moreover nearly the same SHD distance values are measured in the subsets as in the
total Conqueror population (Table 1). Thus SHD is much more sensitive in correctly identifying
genetic similarity levels to multiple potential source populations, however it does not necessarily
inform about real admixing sources, as admixtures could have happened multiple times.
Our another novel algorithm MITOMIX [ 29], can reveal more details of admixture histories as it
performs a hypothesis independent search to find the best admixture combinations from
available populations giving the smallest SHD value from a test population. MITOMIX indicates
that if all modern and ancient populations are considered as potential source, the Conquerors
are best admixed from 26–38% modern Belarusians, 19–34% Tuvans, 18% ancient Baltic Late
Bronze Age and 13% Srubnaya populations (S5A Table). Other possible admix components
may include 9–26% Volga Tatars, Poltavka-Potapovka, Sintastha and Combed Ware
populations. Thus MITOMIX principally derives East Eurasian Conqueror lineages from Tuvans,
Belarusians and Volga Tatars, though latter two are located now in Europe. Belarusians
comprise 22% Lipka Tatars in our dataset [52], who arrived to Europe after the Conquerors’ era,
but seemingly with similar Hg-s. Belarusians are best admixed from Russians, Romanians and
Central Asians (S5B Table), while Tuvans are best admixed from Central Asians and
Mongolians with some Yakut and European elements (S5C Table). Main admixture components
of Volga Tatars are 24–46% Conquerors, 20–50% Russians, 16–30% Mansis+Khantys (Yug)
and 9–18% Norvegians (S5F Table).
Though these results should be interpreted with caution due to imperfect population data, they
sketch the most feasible population processes; ancestors of modern Tuvans could be the
nearest sources of east Eurasian Conqueror lineages and accordingly they map very close to
the eastern Conqueror subset (Fig 5). In turn Tuvans originated mainly from Inner and Central
Asian ancestors, and all these are detected in the entire Conqueror population by SHD (Table
1). Belarusians come to play just because of their similar European and Lipka Tatar
components, while Volga Tatars still harbor a salient Conqueror like subpopulation consisting of
both eastern and western lineages (S5F Table, Table 1). Thus MITOMIX confirms the direct
genetic relation to Volga Tatars, detected by both Fst and SHD, and this relation must be closest
in time.
MITOMIX derives west Eurasian Conqueror lineages by augmenting the European components
of above populations, with admixtures from Baltic Bronze Age (BalBA), Srubnaya (Sru) and
Poltavka-Potapovka (PoBA) populations (S5A Table). This is again in line with Fst data, as
Srubnaya maps closest to the west Eurasian Conqueror subset while PoBA maps very close to
the entire Conqueror group (Fig 5). Although we have shown that limited sample size may give
meaningful results [29], these inferences vindicate caution as we have just 8 mitogenomes from
PoBA, 14 from Srubnaya and no sequences from BalBA.
When only ancient populations are considered as a source, the best admix includes 36–44%
Poltavka-Potapovka, 18–20% Baltic Bronze Age, 11–29% Combed Ware, 14–18% Sintashta
and 14% Srubnaya components (S5A Table), all of which are comprised of solely west Eurasian
Hg-s. However ancient MITOMIX gives significantly higher SHD distances signifying that our
ancient database lacks important east Eurasian components.
Discussion
The most plausible interpretation of the phylogenetic and population genetic results is that the
majority of eastern lineages were ultimately derived from Inner Asia which then migrated to
Central Asia where they admixed with Eurasian lineages before moving to Europe, where they in
turn incorporated west Eurasian elements. As the Conquer population was apparently
assembled from multiple sources this raises the questions as to when did the admixtures
happen, which ancient populations could have been the source and how can our results be
reconciled with historical, archaeological, anthropological and other genetic data.
Relation to Volga Tatars
Our data testify closest genetic relation to this modern population. Volga Tatars incorporate three
main ethnic components [53]; the Volga Bulgars, which arrived in the 8 th century, and
intermingled with local Scythian and Finno-Ugric populations, then in the 13th century Kipchak
Tatars of the Golden Horde brought a final Central-Inner Asian genetic layer and their language
to the region. MITOMIX seems to identify these historical components, as Finno-Ugric Mansis
and Khantys (Yug) comprise a major component of Volga Tatars besides Russians, while
Scythians also appear among their potential sources (S5F Table). Our remarkable result is that
the Conquerors seem to provide a predominant (26–41%) component of Volga Tatars (S5F
Table), while the opposite value is significantly lower (9–26%; S5A Table). This asymmetry is
due to the absence of some Tatar components, like Finno-Ugric ones, from the Conquerors.
Thus our data indicate that rather Volga Tatars harbor a “Conqueror like” genetic component
than the opposite, which may be linked historically to the Volga Bulgars.
This assumption is well supported by archaeological, anthropological and historical sources;
Volga Bulgars were one of the few groups which had the same partial horse burial customs [54]
and similar grave goods as the Conquerors. Both groups are characterized by similar
anthropological types [55], and practiced identical symbolic trepanation customs [ 56] which is
documented with such a high frequency just among the related Danube Bulgars [57]. Historical
data link both groups to the Onogurs [58], the Conquerors must have belonged to the Onogur
tribal union, as the name “Hungarian” is derived from “Onogur” [5,59]. Historical sources imply
that ruling dynasties of both groups might be traced back to the Hun ruling dynasty [60]. Taken
together the direct genetic relation of the Conquerors to Onogur-Bulgar ancestors of Volga
Tatars is very feasible.
East Eurasian relations
Identifying admixture sources further back in time is more precarious, but ancient DNA and
historical data allow drawing some inferences. We may rely on the better grounded Bulgar
prehistory and the Tuvan genetic affinity of the Conquerors, which define a time window, a
geographic region and a migration route through Central Asia to the Pontic steppes.
Both anthropological [ 61] and genetic data [ 47,62] indicate that until the Bronze Age Asia was
populated mainly by Europid Sintashta-Andronovo people west of the Altai, while populations
with Mongoloid traits and genes were confined east of the Altai. The first eastern Hg lineages
appeared in West Siberia at the beginning of Bronze Age [63], in the Altai at the Middle Bronze
Age [64], while in Central Asia just around the 6 th century BC corresponding to the Xiongnu
invasions [65].
In the Iron Age the Tuva region was inhabited by Scytho-Siberians, which were already an
admixed population of east and west Eurasians [49]. During the Iron Age Scytho-Siberians
further admixed with European Scythians in both directions, giving rise to 18–26% eastern
lineages in European Scythians by the 2nd century BC [ 49,66]. Before 200 AD Tuva became part
of the Asian Hun (Xiongnu) empire and Hun migration from Mongolia to west through Altai and
Tuva lead to a significant increase of Mongoloid anthropological components in Central Asia
between the 3 rd century BC and 2 nd century AD [ 61,67]. Thus western (Eurasian) lineages in the
Tuva region can be attributed to Andronovo and Scythian periods, while appearance of east
Eurasian lineages to Asian Scythian and Xiongnu periods. Genetic similarity between Xiongnu
and modern Turcic and Mongolian speaking groups indicate that the Xiongnu period played a
determining role in shaping the genetic profile of Eastern and Central Asia [68], supporting our
phylogeographic implications (Fig 4) that Xiongnus could be among the ancestors of the
Conquerors. A HVR based population genetic study [69] has indeed shown similarity between
Xiongnus and among others Conquerors, as well as Volga Tatars. At any case, the eastern Hg
lineages must have been brought to Europe by nomadic groups originating from this region.
During the first centuries AD Northern Xiongnus were expelled from Inner Asia and escaped
westward [70], leading to another major wave of east Eurasian gene flow into Central Asia, then
further to the Pontic steppes. According to some archaeologists traces of European Huns can be
detected on the Pontic steppe already in the 2nd century AD [ 71], but European Huns entered
history just from the middle of the 4th century as an empire. The Xiongnu origin of European
Huns has been accepted by most historians [72–74], but evidences are scarce.
A decade after the fall of the European Hun empire (472 AD) another grouping of Turkic tribes,
the Ogurs appeared on the Pontic steppe from Central Asia. The Onogurs are the first nomadic
groups from the east, which are reliably connected by historical sources to the later appearing
Bulgars, and less reliably to the Conquerors [58]. Onogurs had been part of the Hunnic people,
and after the death of Attila’s son Irnik, European Hun remains fused with the Onοgurs [58]. The
ensuing Avar invasion brought Onogur groups to the Carpathian Basin, others became part of
the later Danube Bulgar and Volga Bulgar states.
The succeeding group arriving from East Eurasia to the Pontic steppes in the middle of 6th
century were the Avars, who established an empire in the Carpathian Basin lasting for three
centuries [75]. It is relevant to note that none of the Hungarian medieval sources know about
Avars, presumably because they were not distinguished from the Huns [2], as many foreign
medieval sources also identified Avars with the Huns [3].
Subsequent east-west migrations are connected to Göktürk, Kipchak and Mongolian groups, but
these could have minor effect on the Conquerors as mostly arrived after the 10th century,
moreover most Turkic loanwords in Hungarian originate from West Old Turkic [76], the Oghur
Turkic branch associated with previous Turkic speaking groups as Onogurs, Bulgars, Khazars
and maybe the Avars.
Taken together genetic and historical data refer to four major groups delivering significant east
Eurasian lineages to Europe which could be connected to the Conquerors; Asian Scythians,
Huns, Onogurs and Avars. Of these groups we have mitogenome sequences just from
European Scythians [66]. Despite the presence of eastern lineages in European Scythians, they
rather resemble to the European component of the Conquerors (Table 1) suggesting that
eastern Conqueror lineages arrived with later invasions. Thus our genetic data are in line with
historical sources which indicate that Onogurs could have been a major source of the Conqueror
population, nevertheless it is obvious that Hun, Avar and Onogur waves intermingled with each
other and local populations.
West Eurasian relations
According to our data the best fitting sources of the west Eurasian lineages are the Late Bronze
Age Srubnaya (Timber-grave) culture (~1,850–1,200 BC) and its ancestors the Potapovka
(~2,500–1,900 BC) and Poltavka (~2,900–2,200 BC) cultures (Table 1). The Srubnaya was a
nomadic culture on the Pontic-Caspian steppe, both their genetic composition and life style
being closely related to the partly contemporary eastern Andronovo and Sintashta cultures,
together constituting the steppe Middle-Late Bronze Age (MLBA) population. Latter was
descended from the genetically tightly clustering steppe Early-Middle Bronze Age (EMBA)
Yamnaya-Afanasievo-Poltavka cultures with the addition of an European Neolithic farmer
genetic layer [77,78]. As a result, the steppe MLBA population very much resembled genetically
to the European Late Neolithic/Bronze Age (EULNBA) populations [77], providing an explanation
for the similarity of the Conquerors to EULNBA populations (Table 1, Fig 5), the appearance of a
considerable number of modern European and Northwestern European maternal lineages close
to the Conquerors (S1 Fig) and the presence of European Y-chromosomal Hg-s R1b-M269 and
I2a in the Conquerors, reported in our previous study [ 13].
The Armenian Bronze Age (ArmBA) population also appears very close to the Conquerors
(Table 1, Fig 5), that may be explained by the 48–58% Armenian-like Near East ancestry of the
steppe EMBA populations [77], which was ultimately derived from early Iranian farmers [ 78]. This
genetic layer may also explain the appearance of modern populations from the Caucasus region
(Cau, Adg, Arm) close to the Conquerors both in population genetic (Table 1) and phylogenetic
analysis, (S1 Fig). Nevertheless a more recent admixture from this region is also plausible, as all
presumptive carriers of the east Eurasian lineages contacted the Caucasus region during their
westward migrations.
Finno-Ugric relations
Surprisingly we did not find significant genetic relations to Finno-Ugric groups. Though
population genetic analysis indicates some connection of the European Conqueror component
to modern Finnish (Fin) and Baltic (Balt) people, but no relation to Saamis (Sam), Mansis and
Kanthys (Yug) (S4C and S5A Tables). The Baltic relation of the European component seems to
appear already in the Baltic Late Bronze Age (BalBA, 1000–230 BC), [79] measured with the
SHD method (Table 1). BalBA genomes cluster with modern Lithuanians and Estonians, and
lack eastern mtDNA Hg-s and Y-chromosomal haplogroup N-tat, (new name N1a1) which is
typical for Uralic speaking groups, thus Estonians must have received their east Asian-Siberian
components after the BalBA period, from a different source [79]. According to our data BalBA is
best admixed from the closely related Scandinavian Neolith-Bronze Age (NNBA), Afanasevo and
European Neolithic populations (S5D Table), so it is unlikely connected to Finno-Ugric groups.
As only 7 Estonian mitogenomes are available, they were grouped with other modern Baltic
populations (Balt; S3B Table), so the similarity of these to the Conquerors probably derives from
BalBA heritage. The connection to modern Finnish population can also be explained from BalBA
and steppe MLBA components which are present in modern Scandinavians, as Finnish
sequence matches regularly appear together with Danish ones on our phylogenetic trees (S1
Fig, Networks; 14, 15, 19, 25, 27, 30, 35, 40, 42, 43, 49, 52, 56).
Moreover, Y, B and N1a1a1a1a Hg-s have not been detected in Finno-Ugric populations
[80–84], implying that the east Eurasian component of the Conquerors and Finno-Ugric people
are probably not directly related. The same inference can be drawn from phylogenetic data, as
only two Mansi samples appeared in our phylogenetic trees on the side branches (S1 Fig,
Networks; 1, 4) suggesting that ancestors of the Mansis separated from Asian ancestors of the
Conquerors a long time ago. This inference is also supported by genomic Admixture analysis of
Siberian and Northeastern European populations [85], which revealed that Mansis received their
eastern Siberian genetic component approximately 5–7 thousand years ago from ancestors of
modern Even and Evenki people. Most likely the same explanation applies to the Y-chromosome
N-Tat marker which originated from China [86,87] and its subclades are now widespread
between various language groups of North Asia and Eastern Europe [88].
It must be emphasized that Finno-Ugric groups are underrepresented in our population
database, as we have no mitogenomic data from Komis, Maris, Mordvins and Udmurts and only
limited samples from Mansis, Kanthys, Saamis and Estonians. Therefore appearance of FinnoUgric matches from a more representative dataset cannot be excluded, but our data imply that
incidental Finno-Ugric link is rather expected in the European component if any.
Genetic relation of different Conqueror cemeteries
Archaeologist presume that the rich 10th century cemeteries of Karos and Kenézlő comprise the
Conqueror military elite, raising the question as to what extent can our findings be generalized to
the entire Conqueror population. Our fragmentary data from other cemeteries indicate the
presence of the same eastern and western genetic components (S1 Fig, Networks; 3, 4, 12, 36),
moreover [12] and [ 11] reported 91 other Conqueror HVR haplotypes from 24 cemeteries ( Fig
1), which show very similar major Hg distribution to our samples, with even larger proportion of
Asian major Hg components (Fig 6).
Download:
PPT
PowerPoint slide
PNG
larger image
TIFF
original image
Fig 6. Comparison of major Hg distributions from modern and ancient Hungarian
populations.
Asian main Hg-s are designated with brackets. Major Hg distribution of Conqueror
samples from this study are very similar to that of other 91 Conquerors taken from
previous studies [11,12]. Modern Hungarians have very small Asian components pointing
at small contribution from the Conquerors. Of the 289 modern Hungarian mitogenomes
272 are newly deposited [29].
https://doi.org/10.1371/journal.pone.0205920.g006
Thus our conclusions probably apply to the entire Conqueror population, but definitely to the 10 th
century immigrant military elite characterized with partial horse burials, though further
mitogenomic and genomic data are required for the accurate answer.
We have determined the maternal lineage of the majority of samples from the three neighboring
Karos cemeteries, and found likely maternal relatives with identical mtDNA genomes within
cemeteries allocated into the same circles on the phylogenetic trees in S1 Fig (summarized in
S1C Table), but surprisingly no identical haplotype was found between the three Karos
cemeteries. The only exceptions are the two chiefs; Karos2/52 and Karos3/11, who had identical
X2f maternal haplotypes and I2a1 Y chromosomal haplotypes (data not shown), so were
probably brothers. This indicates that these neighboring communities did not intermarry perhaps
because of different group-identity. Furthermore the east Eurasian haplogroup lineages from the
three Karos cemeteries indicate a discernible structuring (S1D Table); the Karos3 cemetery has
a definite south-east Chinese affinity, the Karos1 a North-East Siberian affinity, while the Karos2
lineages are widely distributed from East to Central Asia. In contrast, despite the low number of
samples analyzed from other Conqueror cemeteries we detected potential relatives with
identical mtDNA genomes between distant cemeteries (S1C Table). This suggests that
individual tribes might have been split and fragments of different tribes settled together upon the
conquest.
Relation of the Conquerors to modern Hungarians
Modern Hungarians are genetically very similar to their European neighbors [ 89] nevertheless
they contain some 3–5% east Eurasian components traceable with uniparental markers
[29,90,91], (Fig 6). Genome wide SNP data also detected the presence of 4% east Asian
component in modern Hungarians [92] with an approximate time of admixture dated to the first
millennium AD, corresponding to the invasions of Huns, Onogur-Bulgars, Avars and Hungarian
Conquerors from the Asian steppes, which are completely in line with our results.
Thus genetic heritage of the Conquerors definitely persists in modern Hungarians, but they
contributed to less than 10% of the recent Hungarian gene pool, as they were not alone to bring
in east Eurasian lineages. This dilution could have started with the conquer, as contemporary
local population size in the Carpathian Basin was estimated larger than that of the Conquerors
[93,94]. Anthropological data also have the same implication, as the Conquerors differed from
the subsequent Árpádian Age population, which was more similar to preconquest Avar Age
populations [95,96]. According to early anthropological studies people of the Avar and Conquest
age Carpathian Basin were very heterogeneous and immigrants arrived in several phases
between the 5 th and 9 th centuries [ 97], which in our view admixed with the autochthonous
population, of which genetic data are still barely available between the Bronze Age and Conquest
period.
The large genetic diversity of the Conquerors which seemingly assembled from multiple ethnic
sources and their relative low proportion, having no lasting effect on Hungarian ethnogenesis,
raises doubts about the Conqueror origin of the Hungarian language. Even if our samples
represent mainly the Conqueror elite, the “elite dominance” linguistic hypothesis seems
inconsistent when it presumes that the same Turkic elite was first readily assimilated linguistically
by Finno-Ugric groups, and then it assimilated locals of the Carpathian Basin. Turkic character of
the Conquerors is indicated by their “Turk” denomination in contemporary sources as well as
Turkic tribal names and person names of tribe leaders of the conquest-period [98]. Above data
infer that preconquest presence of the language in the Carpathian Basin, is an equally grounded
hypothesis, as had been proposed by several scientists (a summary in English is given in [99]).
Conclusions
The large diversity of Hg-s detected in the Conquerors reflects a quite complex genetic history,
which was summarized from our data on Fig 4. Their uniform archaeological findings and
predominantly Europid anthropological features (Fig 2, S1A Table) indicate a long lasting
admixture on the Pontic steppe, thus their final composition was likely formed there during the
last centuries prior to the conquest.
A significant fraction of their ancestors undoubtedly arrived from Asia, which probably originated
from Asian Scythians and Xiongnus. On the Pontic steppes Asian nomads assimilated with
descendants of the Srubnayas and this mixed population could have been the basis of many
medieval Pontic nomadic groups, including Conquerors. Their ancestors were certainly part of
the European Hun Empire, the succeeding Avar and Bulgar empires, and when they came into
power they very probably incorporated European Hun remains, as recognized previously [100].
Our genetic data seem to support the Hun-Conqueror connection which could have been the
basis of the historical-cultural Hungarian Hun tradition [3]. Direct genetic relation of the
Conquerors to medieval Onogur-Bulgars warrants further studies, as they are linked by
archaeological, anthropological and historical data as well as our population genetic indications.
Our conclusions are well supported by anthropological studies, which found analogies of the
lower class Conqueror individuals on the eastern European steppes, but parallels of the upper
warrior class were mainly found at the fringes of the Xiongnu empire, in South Siberia and
South-Central Asia [101]. Finally our data indicate that all potential ancestors of the Conquerors
were steppe nomadic people, which is in full agreement with their archaeological legacy.
Supporting information
S1 Appendix. Basic description of the studied cemeteries.
https://doi.org/10.1371/journal.pone.0205920.s001
(DOCX)
S1 Fig. Phylogenetic trees (1–58), made with Median-Joining Network, from mtDNA
sequences of the 102 Hungarian Conquerors.
Phylogenetic trees are arranged in alphabetic order according to haplogroups. The 67 subhaplogroups are depicted on 58 Networks. Samples falling into the same sub-haplogroup with
the studied sample are encircled. The smallest colored circles represent one individual; circle
sizes are proportional to the number of individuals with identical sequences. (When large number
of sequences with few phylogenetically informative SNP-s are aligned, the algorithm may force
the most similar but not identical sequences into the same large circle.) Green circles identify
Hungarian conqueror samples, red circles represent modern samples, and violet circles
correspond to ancient samples. A few ancient samples belonging to the shown haplogroup could
not be properly aligned due to incomplete sequences, and these were connected with dashed
line to the tree. Number of crosslines between neighboring circles denotes mutation distances.
Length of connecting lines is irrelevant, as they were modified in order to fit page. Genebank
accesion number and origin of samples closest to the studied conquerors are listed next to the
circles. Known Conqueror Y-chromosome haplogroups were added in blue color. We
summarized the probable origin of the samples’ Hg lineage in colored framed text. In some
cases comments are given next to the trees.
https://doi.org/10.1371/journal.pone.0205920.s002
(PDF)
S1 Table. Details of Hungarian Conqueror samples.
Table a. Description of samples including anthropological and archaeological details.
Haplogroups and closest matching sequences are also summarized from S1 Fig. Probable
origins of Hg lineages are color coded, codes and details of anthropological ages are given to the
right. Table b. Conqueror subpopulations considered in population genetic analysis. Sample
groupings are based on S1a Table. Table c. List of samples with Identical mtDNA sequences
indicating potential direct maternal relations. Table d. Distribution of the East Eurasian Hg-s in
the three Karos cemeteries.
https://doi.org/10.1371/journal.pone.0205920.s003
(XLSX)
S2 Table. Sequence data.
Table a. Details of NGS data for each samples. Samples highlighted with blue were published in
[15]. Lowest coverage sequences containing larger gaps were highlighted with pink.
Contamination was estimated with two methods; a) using the Schmutzi algorithm and b)
calculating the proportion of reads which did not correspond to the consensus sequence in
diagnostic positions as in [15]. All DNA extracts were partial UDG treated, except Karos2/52 for
which UDG treated and non treated libraries were merged to increase coverage and
misincorporation values of both libraries are provided (values of non treated labelled with *).
Table b. List of SNP-s provided against rCRS. Following the recommendations in [102], we
excluded common indels (hotspots) at nucleotide positions: 309.1C(C), 315.1C, 523-524del (or
522-523del), 3106del, 16182C, 16183C, 16193.1C(C), 16519C. Red numbers indicate SNP-s
missing due to lack of coverage. Haplogroup was determined both with HaploGrep (based on
SNP list) and HaploFind (based on Fasta sequences), haplogroups defined differently are
highlighted with yellow background. In the analysis we used consistently HapoFind defined Hg
assignments.
https://doi.org/10.1371/journal.pone.0205920.s004
(XLSX)
S3 Table. Population database.
Table a. Modern population database with NCBI GenBank accession numbers and haplogroups
given next to population abbreviations (provided in S3B Table). 314 newly deposited
mitogenomes from [29] including 272 Hungarian, 46 French, 18 Croatian, 12 Belgian, 11
German and 12 Romanian samples are provided at the end of the list. Table b. Summary of the
modern (top) and ancient (below) population database with abbreviations used in this study. In
case of low sample size related neighboring populations were merged whose list is provided
next to their group name. Table c. Ancient mtDNA genome database. Pink background highligts
samples, for which sequence was not available, so these were only included in SHD analysis.
Yellow background highlights haplogroups, which were classified differently by Haplofind than
published originally (Haplofind/original). Supplementary references are provided below the table.
Table d. Ancient samples considered only in the phylogenetic analysis but not used in population
genetic analysis. Supplementary references are provided below the table.
https://doi.org/10.1371/journal.pone.0205920.s005
(XLSX)
S4 Table. Population genetic data.
Table a. Pairwise Fst (top) and linearized Slatkin Fst (below) matrix of population distances
between all combinations of modern and ancient population. In the upper right part of the table
probability values are depicted, + correspond to significant P values (<0.05), while—means not
significant P values. Color code (provided below the table) highlights the best similarity levels.
Abbreviations of population names are given in S3B Table. Table b. Pair-wise Shared
Haplogroup Distance (SHD) values measured between all combinations of modern and ancient
populations. Color code (provided below the table) highlights the best similarity levels.
Abbreviations of population names are given in S3B Table. Table c. Comparison of population
genetic distance values measured with two different methods (Fst and SHD) between Hungarian
Conqueror subpopulations and all ancient (arch) and modern (rec) Eurasian populations. Color
code (provided below the table) highlights the best similarity levels. Abbreviations of population
names are given in S3B Table.
https://doi.org/10.1371/journal.pone.0205920.s006
(XLSX)
S5 Table. MITOMIX results.
Table a. Best MITOMIX results for the entire Conqueror population from available population Hg
frequency data (S3A and S3C Table ). Data were computed from all possible proportions of all
possible population combinations (top) or just from contemporary or older ancient populations
(below). The best 64/20 combinations giving the smallest SHD distances to the Conquerors are
listed. East Eurasian populations are highlighted with yellow, Volga Tatars are highlighted with
green. Finno-Ugric groups (Yug) do not appear among the potential sources. Table b. Best
MITOMIX results for modern Belarussians from available population Hg frequency data (S3A and
S3C Table). Data were computed from all possible proportions of all possible population
combinations. The best combinations giving the smallest SHD distances are listed. Table c. Best
MITOMIX results for modern Tuvans from available population Hg frequency data (S3A and S3C
Table). Data were computed from all possible proportions of all possible population
combinations. The best combinations giving the smallest SHD distances are listed. Table d. Best
MITOMIX results for the ancient Baltic Bronze Age population from available population Hg
frequency data (S3A and S3C Table ). Data were computed from all possible proportions of all
possible population combinations (top) or just from contemporary or older ancient populations
(below). The best combinations giving the smallest SHD distances are listed. Table e. Best
MITOMIX results for the ancient Srubnaya population from available population Hg frequency
data (S3A and S3C Table ). Data were computed from all possible proportions of all possible
population combinations (top) or just from contemporary or older ancient populations (below).
The best combinations giving the smallest SHD distances are listed. Table f. Best MITOMIX
results for Volga Tatars from available population Hg frequency data (S3A and S3C Table ). Data
were computed from all possible proportions of all possible population combinations. The best
106 combinations giving the smallest SHD distances to the Volga Tatars are listed. European
Scythian (Scy) and Tuvan (Tuv) admixture sources are highlighted with yellow. Finno-Ugric
groups (Yug) are among the major sources.
https://doi.org/10.1371/journal.pone.0205920.s007
(XLSX)
Acknowledgments
We would like to thank Alissa Mittnik, Mark Stoneking, Leyla Dzhansugurova and Eppie R. Jones
for providing unpublished sequences. We also thank László Révész, Kornél Bakay, Zoltán
Kristóf, Pál Sümegi and Balázs Tihanyi for their useful advices in topics of archaeology and
history. Erika Molnár, and András Bíró helped us with anthropological materials and background.
References
1. Szentpétery I. Scriptores rerum Hungaricarum tempore ducum regnumque stirpis
Arpadianae gestarum [Internet]. Budapest: Acad. Litter. Hungarica; 1937. Available:
https://books.google.hu/books?id=TcAVnQEACAAJ
2. Bakay K. Őstörténetünk régészeti forrásai. Miskolc: MBE; 1997.
3. Hóman B. A magyar hún-hagyomány és hún monda. Budepest: Studium; 1925.
4. Budenz J. A Comparative Dictionary of the Finno-Ugric Elements in the Hungarian
Vocabulary. University I, editor. Bloomington: Curzon Press Ltd.; 1997.
5. Róna-Tas A. Hungarians and Europe in the Early Middle Ages: An Introduction to Early
Hungarian History. Budapest: Central European University Press; 1999.
6. Lahermo P, Laitinen V, Sistonen P, Béres J, Karcagi V, Savontaus ML. MtDNA
polymorphism in the Hungarians: comparison to three other Finno-Ugric-speaking
populations. Hereditas. 2000;132: 35–42. pmid:10857257
View Article
PubMed/NCBI
Google Scholar
7. Semino O, Passarino G, Quintana-Murci L, Liu A, Béres J, Czeizel A, et al. MtDNA and
Y chromosome polymorphisms in Hungary: inferences from the palaeolithic, neolithic
and Uralic influences on the modern Hungarian gene pool. Eur J Hum Genet EJHG.
2000;8: 339–46. pmid:10854093
View Article
PubMed/NCBI
Google Scholar
8. Völgyi A, Zalán A, Szvetnik E, Pamjav H. Hungarian population data for 11 Y-STR and
49 Y-SNP markers. Forensic Sci Int Genet. 2009;3: e27–e28. pmid:19215861
View Article
PubMed/NCBI
Google Scholar
9. Fehér T, Németh E, Vándor A, Kornienko I V., Csáji LK, Pamjav H. Y-SNP L1034:
limited genetic link between Mansi and Hungarian-speaking populations. Mol Genet
Genomics. 2014;290: 377–386. pmid:25258186
View Article
PubMed/NCBI
Google Scholar
10. Pamjav H, Fóthi Á, Fehér T, Fóthi E. A study of the Bodrogköz population in northeastern Hungary by Y chromosomal haplotypes and haplogroups. Mol Genet
Genomics. 2017;292: 883–894. pmid:28409264
View Article
PubMed/NCBI
Google Scholar
11. Tömöry G, Csányi B, Bogácsi-Szabó E, Kalmár T, Czibula Á, Csősz A, et al.
Comparison of maternal lineage and biogeographic analyses of ancient and modern
Hungarian populations. Am J Phys Anthropol. 2007;134: 354–368. pmid:17632797
View Article
PubMed/NCBI
Google Scholar
12. Csősz A, Szécsényi-Nagy A, Csákyová V, Langó P, Bódis V, Köhler K, et al. Maternal
Genetic Ancestry and Legacy of 10th Century AD Hungarians. Sci Rep. 2016;6: 33446.
pmid:27633963
View Article
PubMed/NCBI
Google Scholar
13. Neparáczki E, Juhász Z, Pamjav H, Fehér T, Csányi B, Zink A, et al. Genetic structure
of the early Hungarian conquerors inferred from mtDNA haplotypes and Y-chromosome
haplogroups in a small cemetery. Molecular Genetics and Genomics. 2016: 1–14.
pmid:27803981
View Article
PubMed/NCBI
Google Scholar
14. Csányi B, Bogácsi-Szabo E, Tömöry G, Czibula Á, Priskin K, Csõsz A, et al. Ychromosome analysis of ancient Hungarian and two modern Hungarian-speaking
populations from the Carpathian Basin. Ann Hum Genet. 2008;72: 519–534.
pmid:18373723
View Article
PubMed/NCBI
Google Scholar
15. Neparáczki E, Kocsy K, Tóth GE, Maróti Z, Kalmár T, Bihari P, et al. Revising mtDNA
haplotypes of the ancient Hungarian conquerors with next generation sequencing.
Caramelli D, editor. PLoS One. Public Library of Science; 2017;12: e0174886.
pmid:28422985
View Article
PubMed/NCBI
Google Scholar
16. Stoljarova M, King JL, Takahashi M, Aaspõllu A, Budowle B. Whole mitochondrial
genome genetic diversity in an Estonian population sample. Int J Legal Med. 2016;130:
67–71. pmid:26289416
View Article
PubMed/NCBI
Google Scholar
17. Derenko M, Malyarchuk B, Grzybowski T, Denisova G, Rogalla U, Perkova M, et al.
Origin and post-glacial dispersal of mitochondrial DNA haplogroups C and D in Northern
Asia. PLoS One. 2010;5. pmid:21203537
View Article
PubMed/NCBI
Google Scholar
18. Derenko M, Malyarchuk B, Denisova G, Perkova M, Litvinov A, Grzybowski T, et al.
Western Eurasian ancestry in modern Siberians based on mitogenomic data. BMC Evol
Biol. 2014;14: 217. pmid:25301575
View Article
PubMed/NCBI
Google Scholar
19. Kivisild T. Maternal ancestry and population history from whole mitochondrial genomes.
Investig Genet. 2015;6: 3. pmid:25798216
View Article
PubMed/NCBI
Google Scholar
20. Achilli A, Perego U a, Lancioni H, Olivieri A, Gandini F, Hooshiar Kashani B, et al.
Reconciling migration models to the Americas with the variation of North American
native mitogenomes. Proc Natl Acad Sci U S A. 2013;110: 14308–13. pmid:23940335
View Article
PubMed/NCBI
Google Scholar
21. Pereira L, Richards M, Goios A, Alonso A, Albarrán C, Garcia O, et al. High-resolution
mtDNA evidence for the late-glacial resettlement of Europe from an Iberian refugium.
Genome Res. 2005;15: 19–24. pmid:15632086
View Article
PubMed/NCBI
Google Scholar
22. Posth C, Renaud G, Mittnik A, Drucker DG, Rougier H, Cupillard C, et al. Pleistocene
mitochondrial genomes suggest a single major dispersal of non-africans and a late
glacial population turnover in Europe. Curr Biol. 2016;26: 827–833. pmid:26853362
View Article
PubMed/NCBI
Google Scholar
23. Brotherton P, Haak W, Templeton J, Brandt G, Soubrier J, Jane Adler C, et al. Neolithic
mitochondrial haplogroup H genomes and the genetic origins of Europeans. Nat
Commun. 2013;4: 1764. pmid:23612305
View Article
PubMed/NCBI
Google Scholar
24. Becker RA, Wilks AR. maps: Draw Geographical Maps [Internet]. 2018. Available:
https://cran.r-project.org/web/packages/maps/index.html
View Article
Google Scholar
25. Révész L. A karosi honfoglaláskori temetők. Kovács L, Révész L, editors. Miskolc:
Herman Ottó Múzeum és Magyar Nemzeti Múzeum; 1996.
26. Llamas B, Valverde G, Fehren-Schmitz L, Weyrich LS, Cooper A, Haak W. From the
field to the laboratory: Controlling DNA contamination in human ancient DNA research
in the high-throughput sequencing era. STAR Sci Technol Archaeol Res. 2017;3: 1–14.
View Article
Google Scholar
27. Fortes GG, Paijmans JLA. Analysis of Whole Mitogenomes from Ancient Samples. In:
Kroneis T, editor. Whole Genome Amplification: Methods and Protocols. New York, NY:
Springer New York; 2015. pp. 179–195. https://doi.org/10.1007/978-1-4939-2990-0_13
28. Rohland N, Harney E, Mallick S, Nordenfelt S, Reich D. Partial uracil–DNA–glycosylase
treatment for screening of ancient DNA. Phil Trans R Soc B. 2015;370: 20130624.
pmid:25487342
View Article
PubMed/NCBI
Google Scholar
29. Maróti Z, Török T, Neparáczki E, Raskó I, Nagy I, Maróti M, et al. MITOMIX, an
Algorithm to Reconstruct Population Admixture Histories Indicates Ancient European
Ancestry of Modern Hungarians [Internet]. biRxiv. 2018. https://doi.org/10.1101/247395
View Article
Google Scholar
30. Vianello D, Sevini F, Castellani G, Lomartire L, Capri M, Franceschi C. HAPLOFIND: A
new method for high-throughput mtDNA haplogroup assignment. Hum Mutat. 2013;34:
1189–1194. pmid:23696374
View Article
PubMed/NCBI
Google Scholar
31. Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple
sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:
3059–3066. pmid:12136088
View Article
PubMed/NCBI
Google Scholar
32. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7:
Improvements in performance and usability. Mol Biol Evol. 2013;30: 772–780.
pmid:23329690
View Article
PubMed/NCBI
Google Scholar
33. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis
version 7.0 for bigger datasets. Mol Biol Evol. 2016; msw054. pmid:27004904
View Article
PubMed/NCBI
Google Scholar
34. Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific
phylogenies. Mol Biol Evol. 1999;16: 37–48. pmid:10331250
View Article
PubMed/NCBI
Google Scholar
35. Leigh JW, Bryant D. Popart: Full-Feature Software for Haplotype Network Construction.
Methods Ecol Evol. 2015;6: 1110–1116.
View Article
Google Scholar
36. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform
population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10.
pmid:21565059
View Article
PubMed/NCBI
Google Scholar
37. Slatkin M. A measure of population subdivision based on microsatellite allele
frequencies. Genetics. 1995;139: 457–462. pmid:7705646
View Article
PubMed/NCBI
Google Scholar
38. R Core Development Team. R: a language and environment for statistical computing,
3.2.1. Document freely available on the internet at: http://www.r-project.org. 2015.
https://doi.org/10.1017/CBO9781107415324.004
39. Kustár Á. A Karos-Eperjesszögi I.-II.-III. honfoglalás kori temetők embertani vizsgálata.
In: Kovács L, Révész L, editors. Magyarország honfoglalás kori és kora Árpád-kori
sírleletei 1. Miskolc: A Herman Ottó Múzeum és a Magyar Nemzeti Múzeum közös
kiadványa; 1996. pp. 395–456.
40. Kustár Á, Skultéty G. Facial reconstruction of a Hungarian conqueror male from the site
of Benepuszta. Savaria. 1996;22: 179–190.
View Article
Google Scholar
41. Gerasimov MM. The face finder. Hutchinson and Co, editor. London; 1971.
42. Prag J, Neave R. Making faces: Using forensic and archaeological evidence. London:
British Museum Press; 1997.
43. Taylor KT. Forensic Art and Illustration. London: CRC Press; 2001.
44. Sjovold T. Árpás anatomical method for face reconstruction. Ossa. 1981;7: 203–204.
View Article
Google Scholar
45. Röhrer-Ertl O, Helmer R. Zu Stand und Möglichkeiten der Erneut modifizierten
Kollmann–Methode. (Gesichtsrekonstruktion aufgrund des Schädels.). Gegenbaurs
morphologisches Jahrbuch 130. 1984. pp. 369–398.
View Article
Google Scholar
46. Palanichamy MG, Zhang C-L, Mitra B, Malyarchuk B, Derenko M, Chaudhuri TK, et al.
Mitochondrial haplogroup N1a phylogeography, with implication to the origin of
European farmers. BMC Evol Biol. 2010;10: 304. pmid:20939899
View Article
PubMed/NCBI
Google Scholar
47. Allentoft ME, Sikora M, Sjögren K-G, Rasmussen S, Rasmussen M, Stenderup J, et al.
Population genomics of Bronze Age Eurasia. Nature. 2015;522: 167–172.
pmid:26062507
View Article
PubMed/NCBI
Google Scholar
48. Ricaut FX, Keyser-Tracqui C, Bourgeois J, Crubezy E, Ludes B. Genetic analysis of a
Scytho-siberian skeleton and its implications for ancient Central Asian migrations. Hum
Biol. 2004;76: 109–125. pmid:15222683
View Article
PubMed/NCBI
Google Scholar
49. Unterländer M, Palstra F, Lazaridis I, Pilipenko A, Hofmanová Z, Groß M, et al. Ancestry
and demography and descendants of Iron Age nomads of the Eurasian Steppe. Nat
Commun. 2017;8: 14615. pmid:28256537
View Article
PubMed/NCBI
Google Scholar
50. Fedorova SA, Reidla M, Metspalu E, Metspalu M, Rootsi S, Tambets K, et al.
Autosomal and uniparental portraits of the native populations of Sakha (Yakutia):
implications for the peopling of Northeast Eurasia. BMC Evol Biol. 2013;13: 127.
pmid:23782551
View Article
PubMed/NCBI
Google Scholar
51. QGIS Development Team. QGIS Geographic Information System [Internet]. Open
Source Geospatial Foundation Project. 2015. http://www.qgis.org/
52. Pankratov V, Litvinov S, Kassian A, Shulhin D, Tchebotarev L, Yunusbayev B, et al.
East Eurasian ancestry in the middle of Europe: genetic footprints of Steppe nomads in
the genomes of Belarusian Lipka Tatars. Sci Rep. 2016;6: 30197. pmid:27453128
View Article
PubMed/NCBI
Google Scholar
53. Malyarchuk B, Derenko M, Denisova G, Kravtsova O. Mitogenomic diversity in Tatars
from the Volga-Ural region of Russia. Mol Biol Evol. 2010;27: 2220–2226.
pmid:20457583
View Article
PubMed/NCBI
Google Scholar
54. Türk A. Towards a classification of grave types and burial rites in the 10th– 11th century
Carpathian Basin–Some remarks and observations. In: Doncheva-Petkova L, Balogh C,
Türk A, editors. Avars, Bulgars and Magyars on the Middle and Lower Danube:
Proceedings of the Bulgarian-Hungarian Meeting, Sofia, May 27–28, 2009.
Archaeological Studies of PPCU Department of Archaeology; 2014. pp. 137–156.
55. Fóthi E. A Kárpát-medence X. századi antropológiai leletanyagának megoszlása. In:
Sziget A, editor. A Magyar Őstörténeti Munkaközösség Egyesület Jubileumi
Konferenciája. Budapest; 2006. pp. 49–59.
56. Bereczki Z, Molnár E, Marcsik A, Pálfi G. Rare Types of Trephination from Hungary
Shed New Light on Possible Cross-cultural Connections in the Carpathian Basin. Int J
Osteoarchaeol. 2015;25: 322–333.
View Article
Google Scholar
57. Jordanov J, Dimitrova B, Nikolov S. Symbolic trepanations of skulls from the middle
ages (IXth-Xth Century) in Bulgaria. Acta Neurochir (Wien). 1988;92: 15–18.
View Article
Google Scholar
58. Golden PB. Studies on the Peoples and Cultures of the Eurasian Steppes [Internet].
Cătălin H, editor. Bukarest: Editura Istros; 2011. Available:
https://www.academia.edu/9609971/Studies_on_the_Peoples_and_Cultures_of_the_Eurasian_Steppes
59. Sugar PF. A History of Hungary. Sugar PF, Hanák P, Frank T, editors. Indiana
University Press; 1990.
60. Runciman S. History of the First Bulgarian Empire, 482–1019 [Internet]. Bell; 1930.
Available: http://www.promacedonia.org/en/sr/index.html
View Article
Google Scholar
61. Tumen D. Anthropology of Archaeological populations from Inner Asia. Mong J
Anthropol Archaeol Ethnol. 2008;4: 162–183.
View Article
Google Scholar
62. González-Ruiz M, Santos C, Jordana X, Simón M, Lalueza-Fox C, Gigli E, et al. Tracing
the Origin of the East-West Population Admixture in the Altai Region (Central Asia).
PLoS One. 2012;7. pmid:23152818
View Article
PubMed/NCBI
Google Scholar
63. Molodin VI, Pilipenko AS, Romaschenko AG, Zhuravlev AA, Trapezov RO, Chikisheva
TA. Human migrations in the southern region of the West Siberian Plain during the
Bronze Age: Archaeological, palaeogenetic and anthropological data. Population
Dynamics in Prehistory and Early History: New Approaches Using Stable Isotopes and
Genetics. 2012. pp. 93–112.
View Article
Google Scholar
64. Hollard C, Keyser C, Giscard PH, Tsagaan T, Bayarkhuu N, Bemmann J, et al. Strong
genetic admixture in the Altai at the Middle Bronze Age revealed by uniparental and
ancestry informative markers. Forensic Sci Int Genet. 2014;12: 199–207.
pmid:25016250
View Article
PubMed/NCBI
Google Scholar
65. Lalueza-Fox C, Sampietro ML, Gilbert MTP, Castri L, Facchini F, Pettener D, et al.
Unravelling migrations in the steppe: mitochondrial DNA sequences from ancient
central Asians. Proc Biol Sci. 2004;271: 941–947. pmid:15255049
View Article
PubMed/NCBI
Google Scholar
66. Juras A, Krzewińska M, Nikitin AG, Ehler E, Chyleński M, Łukasik S, et al. Diverse
origin of mitochondrial lineages in Iron Age Black Sea Scythians. Sci Rep. 2017;7:
43950. pmid:28266657
View Article
PubMed/NCBI
Google Scholar
67. Ismagulov O. Ethnic anthropology of Kazakhstan. Institute of History, Archeaology and
Ethnography of the Academy of Science of Kazakhstan. 1982.
68. Pilipenko AS, Cherdantsev S V, Trapezov RO, Zhuravlev AA, Babenko VN,
Pozdnyakov D V, et al. Mitochondrial DNA diversity in a Transbaikalian Xiongnu
population. Archaeol Anthropol Sci. 2017;
View Article
Google Scholar
69. Rogers LL. Understanding ancient human population genetics of the eastern Eurasian
steppe through mitochondrial DNA analysis: Central Mongolian samples from the
Neolithic, Bronze Age, Iron Age and Mongol Empire periods [Internet]. Indiana
University. 2016. Available: http://pqdtopen.proquest.com/pubnum/10253175.html
70. Grousset R. The Empire of the Steppes: A History of Central Asia. Rutgers University
Press; 1970.
71. Botalov SG, Gutsalov S. Hunno-Sarmatians of the Ural-Kazakh steppes [Internet].
Chelyabinsk: Рифей; 2000. Available: https://books.google.hu/books?
id=m0MyAAAAMAAJ
72. Bóna I. A hunok és nagykirályaik. Corvina; 1993.
View Article
Google Scholar
73. Bakay K. Az Árpádok országa: őstörténetünk titkai [Internet]. Miskolc: MBE Nagy Lajos
Király Magánegyetem; 2002. Available: https://books.google.hu/books?
id=igfXAQAACAAJ
74. Anke B. Attila und die Hunnen. der Pfalz (Speyer), Historisches Museum; 2007.
75. Pohl W. Die Awaren: ein Steppenvolk im Mitteleuropa, 567–822 n. Chr [Internet]. Beck;
2002. Available: https://books.google.hu/books?id=qzlEdQRYB4IC
View Article
Google Scholar
76. Róna-Tas A, Berta Á, Károly L. West Old Turkic: Turkic Loanwords in Hungarian
[Internet]. Harrassowitz; 2011. Available: https://books.google.hu/books?
id=DtfVZwEACAAJ
View Article
Google Scholar
77. Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, et al.
Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015;528: 499–
503. pmid:26595274
View Article
PubMed/NCBI
Google Scholar
78. Lazaridis I, Nadel D, Rollefson G, Merrett DC, Rohland N, Mallick S, et al. Genomic
insights into the origin of farming in the ancient Near East. Nature. 2016;536.
pmid:27459054
View Article
PubMed/NCBI
Google Scholar
79. Mittnik A, Wang CC, Pfrengle S, Daubaras M, Zariņa G, Hallgren F, et al. The genetic
prehistory of the Baltic Sea region. Nat Commun. 2018;9. pmid:29382937
View Article
PubMed/NCBI
Google Scholar
80. Pimenoff VN, Comas D, Palo JU, Vershubsky G, Kozlov A, Sajantila A. Northwest
Siberian Khanty and Mansi in the junction of West and East Eurasian gene pools as
revealed by uniparental markers. Eur J Hum Genet. 2008;16: 1254–1264.
pmid:18506205
View Article
PubMed/NCBI
Google Scholar
81. Grosheva AN, Shneider Y V., Zhukova O V., Morozova IY, Rychkov SY. Features of the
Udmurt Mitochondrial Gene Pool in Relation to Tribal Structure. Russ J Genet. 2014;50:
975–986.
View Article
Google Scholar
82. Derbeneva O, Starikovskaya EB, Wallace DC, Sukernik RI. Traces of early Eurasians in
the Mansi of northwest Siberia revealed by mitochondrial DNA analysis. Am J Hum
Genet. 2002;70: 1009–1014. pmid:11845409
View Article
PubMed/NCBI
Google Scholar
83. Derenko M, Malyarchuk B, Denisova G, Perkova M, Rogalla U, Grzybowski T, et al.
Complete mitochondrial DNA analysis of eastern Eurasian haplogroups rarely found in
populations of northern Asia and eastern Europe. PLoS One. 2012;7. pmid:22363811
View Article
PubMed/NCBI
Google Scholar
84. Bermisheva MA, Tambets K, Villeins R, Khusnutdinova EK. Diversity of mitochondrial
DNA haplogroups in ethnic populations of the Volga-Ural region. Mol Biol. 2002;36:
802–812.
View Article
Google Scholar
85. Wong EHM, Khrunin A, Nichols L, Pushkarev D, Khokhrin D, Verbenko D, et al.
Reconstructing genetic history of Siberian and Northeastern European populations.
Genome Res. United States; 2017;27: 1–14. pmid:27965293
View Article
PubMed/NCBI
Google Scholar
86. Rootsi S, Zhivotovsky L a, Baldovic M, Kayser M, Kutuev I a, Khusainova R, et al. A
counter-clockwise northern route of the Y-chromosome haplogroup N from Southeast
Asia towards Europe. Eur J Hum Genet. 2007;15: 204–11. pmid:17149388
View Article
PubMed/NCBI
Google Scholar
87. Shi H, Qi X, Zhong H, Peng Y, Zhang X, Ma RZ, et al. Genetic Evidence of an East
Asian Origin and Paleolithic Northward Migration of Y-chromosome Haplogroup N.
PLoS One. 2013;8. pmid:23840409
View Article
PubMed/NCBI
Google Scholar
88. Ilumäe AM, Reidla M, Chukhryaeva M, Järve M, Post H, Karmin M, et al. Human Y
Chromosome Haplogroup N: A Non-trivial Time-Resolved Phylogeography that Cuts
across Language Families. Am J Hum Genet. 2016;99: 163–173. pmid:27392075
View Article
PubMed/NCBI
Google Scholar
89. Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Sudmant PH, et al. Ancient
human genomes suggest three ancestral populations for present-day Europeans.
Nature. 2014;513: 409–413. pmid:25230663
View Article
PubMed/NCBI
Google Scholar
90. Egyed B, Brandstatter A, Irwin JA, Pádár Z, Parsons TJ, Parson W. Mitochondrial
control region sequence variations in the Hungarian population: Analysis of population
samples from Hungary and from Transylvania (Romania). Forensic Sci Int Genet.
2007;1: 158–162. pmid:19083748
View Article
PubMed/NCBI
Google Scholar
91. Bíró A, Fehér T, Bárány G, Pamjav H. Testing Central and Inner Asian admixture
among contemporary Hungarians. Forensic Sci Int Genet. 2015;15: 121–126.
pmid:25468443
View Article
PubMed/NCBI
Google Scholar
92. Hellenthal G, Busby GBJ, Band G, Wilson JF, Capelli C, Falush D, et al. A genetic atlas
of human admixture history. Science. 2014;343: 747–51. pmid:24531965
View Article
PubMed/NCBI
Google Scholar
93. Kovacsics J. Magyarország történeti demográfiája (896–1995): millecentenáriumi
előadások [Internet]. Központi Statisztikai Hivatal; 1997. Available:
https://books.google.hu/books?id=dafsAAAAMAAJ
View Article
Google Scholar
94. Juhász P. A honfoglalók létszámának kérdéséhez (Population of a Hungarian tribes in
the 9th. century). Belvedere Meridionale. 2015;27: 66–76.
http://dx.doi.org/10.14232/belv.2015.3.7
View Article
Google Scholar
95. Lipták P. Avars and Ancient Hungarians [Internet]. Budapest: Akadémiai Kiadó; 1983.
Available: https://books.google.hu/books?id=dL6AAAAAMAAJ
96. Éry KK. A Kárpát-medence embertani képe a honfoglalás korában (Anthropological
visage of the Carpathian Basin population in the period of Hungarian conquest). In:
Györffy G, editor. Honfoglalás és régészet (Conquest and archeology). Budapest:
Balassi Kiadó; 1994. pp. 217–224, 291–302.
97. Lipták P. Physical Anthropology of the Finno-Ugric Peoples. In: Schwidetzky I, Chiarelli
B, Necrasov O, editors. Physical Anthropology of European Populations. Hague, Paris,
New York: Mouton Publishers; 1980. pp. 366–368.
98. Kristó G. Hungarian History in the Ninth Century [Internet]. Szegedi Középkorász
Műhely; 1996. Available: https://books.google.hu/books?id=N616AAAAMAAJ
View Article
Google Scholar
99. Dreisziger N. Ármin Vámbéry (1832–1913) as a Historian of Early Hungarian
Settlement in the Carpathian Basin. Hungarian Cult Stud. 2013;6: 18–39.
View Article
Google Scholar
100. Szász B. A húnok története: Attila nagykirály [Internet]. Budapest: Bartha Miklós
Társasaág Magyartoörténeti Szakosztály; 1943. Available:
http://mtdaportal.extra.hu/books/szasz_bela_a_hunok_tortenete.pdf
101. Fóthi E. A Kárpát-medence 6–11. századi történetének embertani vonatkozásai. In:
Sudár B, Szentpéteri J, Petkes Z, Lezsák G, Zsidai Z, editors. Magyar őstörténet
Tudomány és hagyományőrzés. Budapest: MTA Bölcsészettudományi Kutatóközpont;
2014. pp. 151–168.
102. Van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human
mitochondrial DNA variation. Hum Mutat. 2009;30. pmid:18853457
View Article
PubMed/NCBI
Google Scholar
Publications
PLOS Biology
PLOS Medicine
PLOS Computational Biology
PLOS Genetics
PLOS Pathogens
PLOS ONE
PLOS Neglected Tropical Diseases
Home
Blogs
Collections
Give feedback
LOCKSS
Privacy Policy
Terms of Use
Advertise
Media Inquiries
Contact
PLOS is a nonprofit 501(c)(3) corporation,
#C2354500, based in San Francisco,
California, US
Correction: Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya origin in the conquering Hungarians
(PLoS ONE (2018) 13: 10 (e0205920) DOI: 10.1371/journal.pone.0205920). Neparáczki, E., Maróti, Z., Kalmár, T., Kocsy, K., Maár, K.,
Bihari, P., Nagy, I., Fóthi, E., Pap, I., Kustár, Á., Pálfi, G., Raskó, I., Zink, A. & Török, T., Nov 1 2018, In : PloS one. 13, 11, e0208295.
Research output: Contribution to journal › Comment/debate. Licensure. Reproduction. Revising mtDNA haplotypes of the ancient
Hungarian conquerors with next generation sequencing. Neparáczki, E., Kocsy, K., Tóth, G. E., Maróti, Z., Kalmár, T., Bihari, P., Nagy, I.,
Pálfi, G., Molnár, E., Raskó, I. & Török, T., Apr 1 2017, In : PLoS One. 12, 4, e0174886. Mitogenomic data indicate admixture components
of Central-Inner Asian and Srubnaya origin in the conquering Hungarians. Endre Neparáczki, Zoltán Maróti, +11 authors Tibor Török. PloS
one. It has been widely accepted that the Finno-Ugric Hungarian language, originated from proto Uralic people, was brought into the
Carpathian Basin by the conquering Hungarians. From the middle of the… (More). Mitogenomic data indicate admixture components of
Central-Inner Asian and Srubnaya origin in the conquering Hungarians. Endre Neparáczki , Contributed equally to this work with: Endre
Neparáczki, Zoltán Maróti. Phylogenetic results reveal that more than one third of the Conqueror maternal lineages were derived from
Central-Inner Asia and their most probable ultimate sources were the Asian Scythians and Asian Huns, giving support to the Hungarian
Hun tradition. The rest of the lineages most likely originated from the Bronze Age Potapovka-Poltavka-Srubnaya cultures of the PonticCaspian steppe. Available data imply that the Conquerors did not have a major contribution to the gene pool of the Carpathian Basin.