Mitochondrial control region sequence variations in the Hungarian population: Analysis of population samples from Hungary and from Transylvania (Romania)

High Coverage Mitogenomes and Y-Chromosomal Typing Reveal Ancient Lineages in the Modern-Day Székely Population in Romania

Here we present 115 whole mitogenomes and 92 Y-chromosomal Short Tandem Repeat (STR) and Single Nucleotide Polymorphism (SNP) profiles from a Hungarian ethnic group, the Székelys (in Romanian: Secuii, in German: Sekler), living in southeast Transylvania (Romania). The Székelys can be traced back to the 12th century in the region, and numerous scientific theories exist as to their origin. We carefully selected sample providers that had local ancestors inhabiting small villages in the area of Odorheiu Secuiesc/Székelyudvarhely in Romania. The results of our research and the reported data signify a qualitative leap compared to previous studies since it presents the first complete mitochondrial DNA sequences and Y-chromosomal profiles of 23 STRs from the region. We evaluated the results with population genetic and phylogenetic methods in the context of the modern and ancient populations that are either geographically or historically related to the Székelys. Our results demonstrate a predominantly local uniparental make-up of the population that also indicates limited admixture with neighboring populations. Phylogenetic analyses confirmed the presumed eastern origin of certain maternal (A, C, D) and paternal (Q, R1a) lineages, and, in some cases, they could also be linked to ancient DNA data from the Migration Period (5th-9th centuries AD) and Hungarian Conquest Period (10th century AD) populations.

Y-chromosomal connection between Hungarians and geographically distant populations of the Ural Mountain region and West Siberia

Hungarians who live in Central Europe today are one of the westernmost Uralic speakers. Despite of the proposed Volga-Ural/West Siberian roots of the Hungarian language, the present-day Hungarian gene pool is highly similar to that of the surrounding Indo-European speaking populations. However, a limited portion of specific Y-chromosomal lineages from haplogroup N, sometimes associated with the spread of Uralic languages, link modern Hungarians with populations living close to the Ural Mountain range on the border of Europe and Asia. Here we investigate the paternal genetic connection between these spatially separated populations. We reconstruct the phylogeny of N3a4-Z1936 clade by using 33 high-coverage Y-chromosomal sequences and estimate the coalescent times of its sub-clades. We genotype close to 5000 samples from 46 Eurasian populations to show the presence of N3a4-B539 lineages among Hungarians and in the populations from Ural Mountain region, including Ob-Ugric-speakers from West Siberia who are geographically distant but linguistically closest to Hungarians. This sub-clade splits from its sister-branch N3a4-B535, frequent today among Northeast European Uralic speakers, 4000-5000 ya, which is in the time-frame of the proposed divergence of Ugric languages.

Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya origin in the conquering Hungarians

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.

Maternal genetic features of the Iron Age Tagar population from Southern Siberia (1st millennium BC)

Early nomads in the Eurasian steppes since the beginning of the 1st millennium BC played a key role in the formation of the cultural and genetic landscape of populations of a significant part of Eurasia, from Eastern Europe to Eastern Central Asia. Numerous archaeological cultures associated with early nomads have been discovered throughout the Eurasian steppe belt. The Tagar archaeological culture existed in the Minusinsk basin (Sayan Mountains, Southern Siberia, Russia) in the northeastern periphery of the Eurasian steppe belt from the 8th to 1st century BC during the pre-Scythian, Scythian, and Early Xiongnu-Sarmatian periods. In this study, we evaluated mtDNA diversity in the Tagar population based on representative series (N = 79) belonging to all chronological stages of the culture. The Tagar population had a mixed mtDNA pool dominated by Western Eurasian haplogroups and subgroups (H, HV6, HV*, I, K, T, U2e, U4, U5a, and U*) and, to a lesser degree, Eastern Eurasian haplogroups (A*, A8, C*, C5, D, G2a, and F1b). The Tagar population showed a similar mtDNA pool structure to those of other Iron Age populations representing the "Scythian World." We observed particularly high similarity between the Tagar and Classic Scythians from the North Pontic region. Our results support the assumption that genetic components introduced by Bronze Age migrants from Western Eurasia contributed to the formation of the genetic composition of Scythian period populations in Southern Siberia. Another important component of the Tagar mtDNA pool was autochthonous East Eurasian lineages, some of which (A8 and C4a2a) are potential markers of the westward genetic influence of the eastern populations of the Scythian period. Our results suggest a genetic continuity (at least partial) between the Early, Middle, and Late Tagar populations.

Whole mitochondrial genome diversity in two Hungarian populations

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

Forensic genetic analyses in isolated populations with examples of central European Valachs and Roma

Isolated populations present a constant threat to the correctness of forensic genetic casework. In this review article we present several examples of how analyzing samples from isolated populations can bias the results of the forensic statistics and analyses. We select our examples from isolated populations from central and southeastern Europe, namely the Valachs and the European Roma. We also provide the reader with general strategies and principles to improve the laboratory practice (best practice) and reporting of samples from supposedly isolated populations. These include reporting the precise population data used for computing the forensic statistics, using the appropriate θ correction factor for calculating allele frequencies, typing ancestry informative markers in samples of unknown or uncertain ethnicity and establishing ethnic-specific forensic databases.

Genetic affinities among the historical provinces of Romania and Central Europe as revealed by an mtDNA analysis

Background

As a major crossroads between Asia and Europe, Romania has experienced continuous migration and invasion episodes. The precise routes may have been shaped by the topology of the territory and had diverse impacts on the genetic structure of mitochondrial DNA (mtDNA) in historical Romanian provinces. We studied 714 Romanians from all historical provinces, Wallachia, Dobrudja, Moldavia, and Transylvania, by analyzing the mtDNA control region and coding markers to encompass the complete landscape of mtDNA haplogroups.

Results

We observed a homogenous distribution of the majority of haplogroups among the Romanian provinces and a clear association with the European populations. A principal component analysis and multidimensional scaling analysis supported the genetic similarity of the Wallachia, Moldavia, and Dobrudja groups with the Balkans, while the Transylvania population was closely related to Central European groups. These findings could be explained by the topology of the Romanian territory, where the Carpathian Arch played an important role in migration patterns. Signals of Asian maternal lineages were observed in all Romanian historical provinces, indicating gene flow along the migration routes through East Asia and Europe.

Conclusions

Our current findings based on the mtDNA analysis of populations in historical provinces of Romania suggest similarity between populations in Transylvania and Central Europe, supported both by the observed clines in haplogroup frequencies for several European and Asian maternal lineages and MDS analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0487-5) contains supplementary material, which is available to authorized users.

Maternal Genetic Ancestry and Legacy of 10(th) Century AD Hungarians

The ancient Hungarians originated from the Ural region in today’s central Russia and migrated across the Eastern European steppe, according to historical sources. The Hungarians conquered the Carpathian Basin 895–907 AD, and admixed with the indigenous communities. Here we present mitochondrial DNA results from three datasets: one from the Avar period (7^th–9^th centuries) of the Carpathian Basin (n = 31); one from the Hungarian conquest-period (n = 76); and a completion of the published 10^th–12^th century Hungarian-Slavic contact zone dataset by four samples. We compare these mitochondrial DNA hypervariable segment sequences and haplogroup results with published ancient and modern Eurasian data. Whereas the analyzed Avars represents a certain group of the Avar society that shows East and South European genetic characteristics, the Hungarian conquerors’ maternal gene pool is a mixture of West Eurasian and Central and North Eurasian elements. Comprehensively analyzing the results, both the linguistically recorded Finno-Ugric roots and historically documented Turkic and Central Asian influxes had possible genetic imprints in the conquerors’ genetic composition. Our data allows a complex series of historic and population genetic events before the formation of the medieval population of the Carpathian Basin, and the maternal genetic continuity between 10^th–12^th century and modern Hungarians.

Inter-laboratory evaluation of the EUROFORGEN Global ancestry-informative SNP panel by massively parallel sequencing using the Ion PGM™

The EUROFORGEN Global ancestry-informative SNP (AIM-SNPs) panel is a forensic multiplex of 128 markers designed to differentiate an individual's ancestry from amongst the five continental population groups of Africa, Europe, East Asia, Native America, and Oceania. A custom multiplex of AmpliSeq™ PCR primers was designed for the Global AIM-SNPs to perform massively parallel sequencing using the Ion PGM™ system. This study assessed individual SNP genotyping precision using the Ion PGM™, the forensic sensitivity of the multiplex using dilution series, degraded DNA plus simple mixtures, and the ancestry differentiation power of the final panel design, which required substitution of three original ancestry-informative SNPs with alternatives. Fourteen populations that had not been previously analyzed were genotyped using the custom multiplex and these studies allowed assessment of genotyping performance by comparison of data across five laboratories. Results indicate a low level of genotyping error can still occur from sequence misalignment caused by homopolymeric tracts close to the target SNP, despite careful scrutiny of candidate SNPs at the design stage. Such sequence misalignment required the exclusion of component SNP rs2080161 from the Global AIM-SNPs panel. However, the overall genotyping precision and sensitivity of this custom multiplex indicates the Ion PGM™ assay for the Global AIM-SNPs is highly suitable for forensic ancestry analysis with massively parallel sequencing.

Indian signatures in the westernmost edge of the European Romani diaspora: new insight from mitogenomes

In agreement with historical documentation, several genetic studies have revealed ancestral links between the European Romani and India. The entire mitochondrial DNA (mtDNA) of 27 Spanish Romani was sequenced in order to shed further light on the origins of this population. The data were analyzed together with a large published dataset (mainly hypervariable region I [HVS-I] haplotypes) of Romani (N=1,353) and non-Romani worldwide populations (N>150,000). Analysis of mitogenomes allowed the characterization of various Romani-specific clades. M5a1b1a1 is the most distinctive European Romani haplogroup; it is present in all Romani groups at variable frequencies (with only sporadic findings in non-Romani) and represents 18% of their mtDNA pool. Its phylogeographic features indicate that M5a1b1a1 originated 1.5 thousand years ago (kya; 95% CI: 1.3-1.8) in a proto-Romani population living in Northwest India. U3 represents the most characteristic Romani haplogroup of European/Near Eastern origin (12.4%); it appears at dissimilar frequencies across the continent (Iberia: ≈ 31%; Eastern/Central Europe: ≈ 13%). All U3 mitogenomes of our Iberian Romani sample fall within a new sub-clade, U3b1c, which can be dated to 0.5 kya (95% CI: 0.3-0.7); therefore, signaling a lower bound for the founder event that followed admixture in Europe/Near East. Other minor European/Near Eastern haplogroups (e.g. H24, H88a) were also assimilated into the Romani by introgression with neighboring populations during their diaspora into Europe; yet some show a differentiation from the phylogenetically closest non-Romani counterpart. The phylogeny of Romani mitogenomes shows clear signatures of low effective population sizes and founder effects. Overall, these results are in good agreement with historical documentation, suggesting that cultural identity and relative isolation have allowed the Romani to preserve a distinctive mtDNA heritage, with some features linking them unequivocally to their ancestral Indian homeland.

Introducing the Algerian mitochondrial DNA and Y-chromosome profiles into the North African landscape

North Africa is considered a distinct geographic and ethnic entity within Africa. Although modern humans originated in this Continent, studies of mitochondrial DNA (mtDNA) and Y-chromosome genealogical markers provide evidence that the North African gene pool has been shaped by the back-migration of several Eurasian lineages in Paleolithic and Neolithic times. More recent influences from sub-Saharan Africa and Mediterranean Europe are also evident. The presence of East-West and North-South haplogroup frequency gradients strongly reinforces the genetic complexity of this region. However, this genetic scenario is beset with a notable gap, which is the lack of consistent information for Algeria, the largest country in the Maghreb. To fill this gap, we analyzed a sample of 240 unrelated subjects from a northwest Algeria cosmopolitan population using mtDNA sequences and Y-chromosome biallelic polymorphisms, focusing on the fine dissection of haplogroups E and R, which are the most prevalent in North Africa and Europe respectively. The Eurasian component in Algeria reached 80% for mtDNA and 90% for Y-chromosome. However, within them, the North African genetic component for mtDNA (U6 and M1; 20%) is significantly smaller than the paternal (E-M81 and E-V65; 70%). The unexpected presence of the European-derived Y-chromosome lineages R-M412, R-S116, R-U152 and R-M529 in Algeria and the rest of the Maghreb could be the counterparts of the mtDNA H1, H3 and V subgroups, pointing to direct maritime contacts between the European and North African sides of the western Mediterranean. Female influx of sub-Saharan Africans into Algeria (20%) is also significantly greater than the male (10%). In spite of these sexual asymmetries, the Algerian uniparental profiles faithfully correlate between each other and with the geography.

Genetic structure of the paternal lineage of the Roma people

According to written sources, Roma (Romanies, Gypsies) arrived in the Balkans around 1,000 years ago from India and have subsequently spread through several parts of Europe. Genetic data, particularly from the Y chromosome, have supported this model, and can potentially refine it. We now provide an analysis of Y-chromosomal markers from five Roma and two non-Roma populations (N = 787) in order to investigate the genetic relatedness of the Roma population groups to one another, and to gain further understanding of their likely Indian origins, the genetic contribution of non-Roma males to the Roma populations, and the early history of their splits and migrations in Europe. The two main sources of the Roma paternal gene pool were identified as South Asian and European. The reduced diversity and expansion of H1a-M82 lineages in all Roma groups imply shared descent from a single paternal ancestor in the Indian subcontinent. The Roma paternal gene pool also contains a specific subset of E1b1b1a-M78 and J2a2-M67 lineages, implying admixture during early settlement in the Balkans and the subsequent influx into the Carpathian Basin. Additional admixture, evident in the low and moderate frequencies of typical European haplogroups I1-M253, I2a-P37.2, I2b-M223, R1b1-P25, and R1a1-M198, has occurred in a more population-specific manner.

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

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

Paternal genetic history of the Vlax Roma

Romanies constitute the largest minority group belonging to different subgroups in Hungary. Vlax Romanies are one of these Romani subgroups. The Gypsies came to Hungary from the Balkans in two large migrations. The Carpathian Romanies arrived in the 15th century and the Vlax Romanies came in the 19th century. The Carpathian Gypsies speak Hungarian and the Vlax Romanies speak Hungarian and Romani languages. Only a limited number of genetic studies of Y-chromosomal haplotypes/haplogroups have been done before, moreover most studies did not contain information regarding the investigated Roma populations which subgroups belong to. In the present study, we analyzed a wide set of Y-chromosomal markers to do comparable studies of the Vlax Roma in eastern Hungarian regions. The results can be compared in the context of previously published data on other Romani groups, Indian and Hungarian reference populations. Haplogroups H1a-M82 and J2a2-M67 were most common in the investigated population groups. A median-joining network of haplogroup H1a-M82 has demonstrated the sharing of identical Indian specific Y-chromosomal lineages between all Romani populations including Malaysian Indians as well as the Vlax Romanies. This common lineage of haplogroup H1a-M82 represents a common descent from a single ancestor provides a strong genetic link to the ancestral geographical origin of the proto-Gypsies. The detected haplogroups in the Vlax Romani population groups can be classified into two different Y-chromosomal lineages based on their putative origin. These lineages include ancestral Indian (H1a-M82), present-day Eurasian (J2a2-M67, J2*-M172, E1b1b1a-M78, I1-M253, R1a1-M198 and R1b1-P25) Y-chromosome lineages. Presence of these lineages in the paternal gene pool of the Roma people is illustrative of the Gypsy migration route from India through the Balkan to the Carpathian Basin.

Mitochondrial DNA phylogeny in Eastern and Western Slavs

To resolve the phylogeny of certain mitochondrial DNA (mtDNA) haplogroups in eastern Europe and estimate their evolutionary age, a total of 73 samples representing mitochondrial haplogroups U4, HV*, and R1 were selected for complete mitochondrial genome sequencing from a collection of about 2,000 control region sequences sampled in eastern (Russians, Belorussians, and Ukrainians) and western (Poles, Czechs, and Slovaks) Slavs. On the basis of whole-genome resolution, we fully characterized a number of haplogroups (HV3, HV4, U4a1, U4a2, U4a3, U4b, U4c, U4d, and R1a) that were previously described only partially. Our findings demonstrate that haplogroups HV3, HV4, and U4a1 could be traced back to the pre-Neolithic times ( approximately 12,000-19,000 years before present [YBP]) in eastern Europe. In addition, an ancient connection between the Caucasus/Europe and India has been revealed by analysis of haplogroup R1 diversity, with a split between the Indian and Caucasus/European R1a lineages occurring about 16,500 years ago. Meanwhile, some mtDNA subgroups detected in Slavs (such as U4a2a, U4a2*, HV3a, and R1a1) are definitely younger being dated between 6,400 and 8,200 YBP. However, robust age estimations appear to be problematic due to the high ratios of nonsynonymous to synonymous substitutions found in young mtDNA subclusters.

Mitochondrial DNA variability in Slovaks, with application to the Roma origin

To gain insight into the mitochondrial gene pool diversity of European populations, we studied mitochondrial DNA (mtDNA) variability in 207 subjects from western and eastern areas of Slovakia. Sequencing of two hypervariable segments, HVS I and HVS II, in combination with screening of coding region haplogroup-specific RFLP-markers, revealed that the majority of Slovak mtDNAs belong to the common West Eurasian mitochondrial haplogroups (HV, J, T, U, N1, W, and X). However, a few sub-Saharan African (L2a) mtDNAs were detected in a population from eastern part of Slovakia. In addition, about 3% of mtDNAs from eastern Slovakia encompass Roma-specific lineages. By means of complete mtDNA sequencing we demonstrate here that the Roma-specific M-lineages observed in gene pools of different Slavonic populations (Slovaks, Poles and Russians), belong to Indian-specific haplogroups M5a1 and M35. Moreover, we show that haplogroup J lineages found in gene pools of the Roma and some Slavonic populations (Czechs and Slovaks) belong to new subhaplogroup J1a, which is defined by coding region mutation at position 8460.