Vlax Roma history: what do coalescent-based methods tell us?

From Croatian Roma to 1000 Genomes: The Story of the CYP2D6 Gene Promoter and Enhancer SNPs

The CYP2D6 gene encodes an enzyme responsible for the metabolism of ~20% of clinically prescribed drugs. In this study, 18 SNPs from the enhancer and promoter regions of CYP2D6 in 323 Roma from Croatia were genotyped, to find out whether the demographic history of Roma affected the distribution of the studied SNPs and their linkage disequilibrium (LD) values, with the major SNPs defining the CYP2D6 star alleles. No differences were found between the three Roma groups in allele and genotype frequencies. The distribution of LD values of Roma was compared with LD values of European and Asian populations. Regulatory CYP2D6 SNPs (rs5758550, rs28624811, rs1080985 and rs1080983) showed similar distribution and the highest LDs with rs16947 from the gene-coding region in all populations. In the promoter region, a complete LD between rs1080989 and rs28588594, and between rs1080983 and rs28624811, was found in Croatian Roma and investigated populations from 1000 genomes. A high LD was also found between rs1080985 from the promoter and rs5758550 from the enhancer region. SNP rs28735595 from the gene promoter region had the highest LD, with two gene region SNPs, rs1058164 and rs1135840. To conclude, the Croatian Roma population shows an LD pattern of the CYP2D6 gene region similar to the 1000 Genomes European and Asian populations.

Untangling SNP Variations within CYP2D6 Gene in Croatian Roma

CYP2D6 is a highly polymorphic gene whose variations affect its enzyme activity. To assess whether the specific population history of Roma, characterized by constant migrations and endogamy, influenced the distribution of alleles and thus phenotypes, the CYP2D6 gene was sequenced using NGS (Next Generation Sequencing) method-targeted sequencing in three groups of Croatian Roma (N = 323) and results were compared to European and Asian populations. Identified single nucleotide polymorphisms (SNPs) were used to reconstruct haplotypes, which were translated into the star-allele nomenclature and later into phenotypes. A total of 43 polymorphic SNPs were identified. The three Roma groups differed significantly in the frequency of alleles of polymorphisms 6769 A > G, 6089 G > A, and 5264 A > G (p < 0.01), as well as in the prevalence of the five most represented star alleles: *1, *2, *4, *10, and *41 (p < 0.0001). Croatian Roma differ from the European and Asian populations in the accumulation of globally rare SNPs (6089 G > A, 4589 C > T, 4622 G > C, 7490 T > C). Our results also show that demographic history influences SNP variations in the Roma population. The three socio-culturally different Roma groups studied differ significantly in the distribution of star alleles, which confirms the importance of a separate study of different Roma groups.

Recent Common Origin, Reduced Population Size, and Marked Admixture Have Shaped European Roma Genomes

The Roma Diaspora-traditionally known as Gypsies-remains among the least explored population migratory events in historical times. It involved the migration of Roma ancestors out-of-India through the plateaus of Western Asia ultimately reaching Europe. The demographic effects of the Diaspora-bottlenecks, endogamy, and gene flow-might have left marked molecular traces in the Roma genomes. Here, we analyze the whole-genome sequence of 46 Roma individuals pertaining to four migrant groups in six European countries. Our analyses revealed a strong, early founder effect followed by a drastic reduction of ∼44% in effective population size. The Roma common ancestors split from the Punjabi population, from Northwest India, some generations before the Diaspora started, <2,000 years ago. The initial bottleneck and subsequent endogamy are revealed by the occurrence of extensive runs of homozygosity and identity-by-descent segments in all Roma populations. Furthermore, we provide evidence of gene flow from Armenian and Anatolian groups in present-day Roma, although the primary contribution to Roma gene pool comes from non-Roma Europeans, which accounts for >50% of their genomes. The linguistic and historical differentiation of Roma in migrant groups is confirmed by the differential proportion, but not a differential source, of European admixture in the Roma groups, which shows a westward cline. In the present study, we found that despite the strong admixture Roma had in their diaspora, the signature of the initial bottleneck and the subsequent endogamy is still present in Roma genomes.

Sex-biased patterns shaped the genetic history of Roma

The Roma population is a European ethnic minority characterized by recent and multiple dispersals and founder effects. After their origin in South Asia around 1,500 years ago, they migrated West. In Europe, they diverged into ethnolinguistically distinct migrant groups that spread across the continent. Previous genetic studies based on genome-wide data and uniparental markers detected Roma founder events and West-Eurasian gene flow. However, to the best of our knowledge, it has not been assessed whether these demographic processes have equally affected both sexes in the population. The present study uses the largest and most comprehensive dataset of complete mitochondrial and Y chromosome Roma sequences to unravel the sex-biased patterns that have shaped their genetic history. The results show that the Roma maternal genetic pool carries a higher lineage diversity from South Asia, as opposed to a single paternal South Asian lineage. Nonetheless, the European gene flow events mainly occurred through the maternal lineages; however, a signal of this gene flow is also traceable in the paternal lineages. We also detect a higher female migration rate among European Roma groups. Altogether, these results suggest that sociocultural factors influenced the emergence of sex-biased genetic patterns at global and local scales in the Roma population through time.

European Roma groups show complex West Eurasian admixture footprints and a common South Asian genetic origin

The Roma population is the largest transnational ethnic minority in Europe, characterized by a linguistic, cultural and historical heterogeneity. Comparative linguistics and genetic studies have placed the origin of European Roma in the Northwest of India. After their migration across Persia, they entered into the Balkan Peninsula, from where they spread into Europe, arriving in the Iberian Peninsula in the 15th century. Their particular demographic history has genetic implications linked to rare and common diseases. However, the South Asian source of the proto-Roma remains still untargeted and the West Eurasian Roma component has not been yet deeply characterized. Here, in order to describe both the South Asian and West Eurasian ancestries, we analyze previously published genome-wide data of 152 European Roma and 34 new Iberian Roma samples at a fine-scale and haplotype-based level, with special focus on the Iberian Roma genetic substructure. Our results suggest that the putative origin of the proto-Roma involves a Punjabi group with low levels of West Eurasian ancestry. In addition, we have identified a complex West Eurasian component (around 65%) in the Roma, as a result of the admixture events occurred with non-proto-Roma populations between 1270–1580. Particularly, we have detected the Balkan genetic footprint in all European Roma, and the Baltic and Iberian components in the Northern and Western Roma groups, respectively. Finally, our results show genetic substructure within the Iberian Roma, with different levels of West Eurasian admixture, as a result of the complex historical events occurred in the Peninsula.

Human demographic processes and admixture events leave traceable footprints in the genomes of the populations and they can modulate the genetic architecture of complex diseases. Here, we aim to study the Roma people, an admixed population with a particular demographic history recognized as the largest ethnic minority in Europe. Previous studies suggest that they originated in South Asia 1,500 years ago and followed a diaspora towards Europe with extensive admixture with non-Roma West Eurasian groups. However, the genetic components of the Roma have not been deeply characterized. Our study reveals a common South Asian origin of all European Roma, closely related to a Punjabi group from Northwestern India. Through fine-scale haplotype-based methods, we describe a complex West Eurasian genetic component in the Roma groups, identifying a common Balkan ancestry and country-specific admixture footprints consistent with the dispersion through Europe. Our findings provide new insights into the demographic history and recent admixture events that have shaped the genetic composition of European Roma groups and could enable a better genetic characterization of complex disease in this population.

Carrier rates of four single-gene disorders in Croatian Bayash Roma

To assess how specific population history, different migration routes, isolation, and endogamy practices contributed to the distribution of several rare diseases found in specific Roma groups, we conducted a population-based research study of rare disease mutations in Croatian Vlax Roma. We tested a total of 427 subjects from Baranja and Međimurje for the presence of four mutations causing hereditary motor and sensory neuropathy type Lom (HMSNL), GM1 gangliosidosis (GM1), congenital cataracts, facial dysmorphism and neuropathy (CCFDN), and limb girdle muscle dystrophy type 2C (LGMD2C), using the RFLP-PCR method to estimate carrier frequencies. We identified a total of four individuals heterozygous for the mutation causing HMSNL in the Baranja population, with a carrier rate amounting to 1.5%. Carriers for other three mutations causing GM1, CCFDN, and LGMD2C were not found in our sample. The carrier rate for the HMSNL mutation in Baranja is lower than in other Vlax Roma groups. In addition, distinct differences in carrier rates between the Croatian Vlax groups point to different genetic history, despite their belonging to the same Roma migration category and subgroup. The difference in carrier rates is either the result of admixture or the reflection of a greater extent of genetic drift since recent founding, maintained by a high degree of endogamy.

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.

Different patterns of Toll-like receptor 2 polymorphisms in populations of various ethnic and geographic origins

Upon the invasion of the host by microorganisms, innate immunity is triggered through pathogen recognition by pattern recognition receptors (PRRs). Toll-like receptors (TLRs) are the best-studied class of PRRs, and they recognize specific pathogen-associated molecular patterns (PAMPs) from various microorganisms. A large number of studies have shown that genetic variation in TLRs may influence susceptibility to infections. We assessed the genetic variation of TLR2, which encodes one of the most important TLRs, in various populations around the globe and correlated it with changes in the function of the molecule. The three best-known nonsynonymous TLR2 polymorphisms (1892C>A, 2029C>T, and 2258G>A) were assessed in different populations from the main continental masses: Romanians, Vlax-Roma, Dutch (European populations), Han Chinese (East Asia), Dogon, Fulani (Africa), and Trio Indians (America). The 2029C>T polymorphism was absent in both European and non-European populations, with the exception of the Vlax-Roma, suggesting that this polymorphism most likely arose in Indo-Aryan people after migration into South Asia. The 1892C>A polymorphism that was found exclusively in European populations, but not in Asian, African, or American volunteers, probably occurred in proto-Indo-Europeans. Interestingly, 2258G>A was present only in Europeans, including Vlax-Roma, but at a very low frequency. The differential pattern of the TLR2 polymorphisms in various populations may explain some of the differences in susceptibility to infections between these populations.

Unique frequencies of HFE gene variants in Roma/Gypsies

The aim of this study was to assess the frequencies of three hemochromatosis gene (HFE) mutations in ethnic Roma/Gypsies in Slovakia. A cohort of 367 individuals representing general population and not preselected for health status was genotyped by TaqMan real-time PCR assay for C282Y, H63D and S65C mutations in HFE gene. A unique genetic profile was revealed: C282Y is found in the highest frequency of all Central European countries (4.90%), while the frequency of H63D mutation (4.09%) is lower than any reported in Europe so far. S65C mutation was not present in the cohort. These mutation frequencies can be explained rather by gene influx and genetic isolation than by genetic inheritance from a former Roma/Gypsy homeland.

E2 allele of the apolipoprotein E gene polymorphism is predictive for obesity status in Roma minority population of Croatia

Background and Aims

The Roma (Gypsies) are a transnational minority, founder population characterized by unique genetic background modeled by culturally determined endogamy. The present study explores whether the widely found cardiovascular diseases (CVD) risk effects of ACE I/D, APOE (ε2, ε3, ε4), eNOS-VNTR and LEP G2548A polymorphisms can be replicated in this specific population.

Methods and Results

The community-based study was carried on 208 adult Bayash Roma living in rural settlements of eastern and northern Croatia. Risk effect of four CVD candidate polymorphisms are related to the most prominent classical CVD risk phenotypes: obesity indicators (body mass index and waist circumference), hypertension and hyperlipidemia (triglycerides, HDL and LDL cholesterol). For all of them the standard risk cut-offs were applied. The extent to which the phenotypic status is related to genotype was assessed by logistic regression analysis. The strongest associations were found for ε2 allele of the APOE as a predictor of waist circumference (OR 3.301; 95%CI 1.254-8.688; p = 0.016) as well as for BMI (OR 3.547; 95%CI 1.471-8.557; p = 0.005). It is notable that ε3 allele of APOE gene turned out to be a protective genetic factor determining low lipid levels.

Conclusion

The strength of the relation and the similarity of the results obtained for both tested indicators of obesity provide firm evidence that APOE plays an important role in obesity development in the Roma population.

Y chromosome variation in Europe: Continental and local processes in the formation of the extant gene pool

The polymorphism of the male-specific portion of the Y chromosome has been increasingly used to describe the composition of the European gene pool and to reconstruct its formation. Here the theoretical grounds and the limitations of this approach are presented, together with the different views on debated issues. The emerging picture for the composition of the male gene pool of the continent is illustrated, but local peculiarities that represent departures from the main trends are also highlighted, in order to illustrate the main unifying feature, i.e. the overlay of recent patterns onto more ancient ones. A synopsis of the main findings and conclusions obtained in regional studies has also been compiled.

Inter- and extra-Indian admixture and genetic diversity in reunion island revealed by analysis of mitochondrial DNA

Reunion Island is a French territory located in the western Indian Ocean. The genetic pattern of the Reunionese population has been shaped by contributions from highly contrasting regions of the world. Over the last 350 years, several migration waves and cultural and socio-economic factors have led to the emergence of six main ethnic groups in Reunion. India is one of the principal regions that contributed to the setting up of the Reunionese population. Diversity, demographic and admixture analyses were performed on mtDNA variation of the Reunionese of Indian ancestry, including the Malbar and Zarab ethnic groups, in order to question their history. Using a phylogeographical approach, we generated and analysed quantitative data on the contribution of the Indian parental populations. Furthermore, we showed that the settlement of Reunion Island by Indians did not involve a founder effect, except in the very beginning of the Reunionese settlement (at the end of the 17(th) century). The accuracy of our results was optimised by a re-evaluation of the classification of the Southern Asian mtDNA haplogroups. Finally, by comparing our results to a previous study dealing with the Reunionese population, we highlighted how ethno-historical data are critical for reconstructing the complex history of multiethnic populations.

Investigations of the Y chromosome, male founder structure and YSTR mutation rates in the Old Order Amish

OBJECTIVES

Using Y chromosome short tandem repeat (YSTR) genotypes, (1) evaluate the accuracy and completeness of the Lancaster County Old Order Amish (OOA) genealogical records and (2) estimate YSTR mutation rates.

METHODS

Nine YSTR markers were genotyped in 739 Old Order Amish males who participated in several ongoing genetic studies of complex traits and could be connected into one of 28 all-male lineage pedigrees constructed using the Anabaptist Genealogy Database and the query software Ped-Hunter. A putative founder YSTR haplotype was constructed for each pedigree, and observed and inferred father-son transmissions were used to estimate YSTR mutation rates.

RESULTS

We inferred 27 distinct founder Y chromosome haplotypes in the 28 male lineages, which encompassed 27 surnames accounting for 98% of Lancaster OOA households. Nearly all deviations from founder haplotypes were consistent with mutation events rather than errors. The estimated marker-specific mutation rates ranged from 0 to 1.09% (average 0.33% using up to 283 observed meioses only and 0.28% using up to 1,232 observed and inferred meioses combined).

CONCLUSIONS

These data confirm the accuracy and completeness of the male lineage portion of the Anabaptist Genealogy Database and contribute mutation rate estimates for several commonly used Y chromosome STR markers.

Hungarian mtDNA population databases from Budapest and the Baranya county Roma

To facilitate forensic mtDNA testing in Hungary, we have generated control region databases for two Hungarian populations: 211 individuals were sampled from the urban Budapest population and 208 individuals were sampled from a Romani ("gypsy") population in Baranya county. Sequences were generated using a highly redundant approach to minimize potential database errors. The Budapest population had high sequence diversity with 180 lineages, 183 polymorphic positions, and a random match probability of 1%. In contrast, the Romani population exhibited low sequence diversity, with only 56 lineages, 109 segregating sites, and a random match probability of 8.8%. The mtDNA haplogroup compositions of the two populations were also distinct, with the large proportion of haplogroup M samples (35%) in the Roma the most obvious difference between the two populations. These factors highlight the importance of considering population structure when generating reference databases for forensic testing purposes. Comparisons between our Romani population sample and other published data indicate the need for heightened caution when sampling and using mtDNA databases of small endogamous populations. The Romani populations that we compared showed significant departures from genetic uniformity.

Mitochondrial DNA diversity in the Polish Roma

Mitochondrial DNA variability in the Polish Roma population has been studied by means of hypervariable segment I and II (HVS I and II) sequencing and restriction fragment-length polymorphism analysis of the mtDNA coding region. The mtDNA haplotypes detected in the Polish Roma fall into the common Eurasian mitochondrial haplogroups (H, U3, K, J1, X, I, W, and M*). The results of complete mtDNA sequencing clearly indicate that the Romani M*-lineage belongs to the Indian-specific haplogroup M5, which is characterized by three transitions in the coding region, at sites 12477, 3921 and 709. Molecular variance analysis inferred from mtDNA data reveals that genetic distances between the Roma groups are considerably larger than those between the surrounding European populations. Also, there are significant differences between the Bulgarian Roma (Balkan and Vlax groups) and West European Roma (Polish, Lithuanian and Spanish groups). Comparative analysis of mtDNA haplotypes in the Roma populations shows that different haplotypes appear to demonstrate impressive founder effects: M5 and H (16261-16304) in all Romani groups; U3, I and J1 in some Romani groups. Interestingly, haplogroup K (with HVS I motif 16224-16234-16311) found in the Polish Roma sample seems to be specific for Ashkenazi Jewish populations.

A newly discovered founder population: the Roma/Gypsies

The Gypsies (a misnomer, derived from an early legend about Egyptian origins) defy the conventional definition of a population: they have no nation-state, speak different languages, belong to many religions and comprise a mosaic of socially and culturally divergent groups separated by strict rules of endogamy. Referred to as "the invisible minority", the Gypsies have for centuries been ignored by Western medicine, and their genetic heritage has only recently attracted attention. Common origins from a small group of ancestors characterise the 8-10 million European Gypsies as an unusual trans-national founder population, whose exodus from India played the role of a profound demographic bottleneck. Social and economic pressures within Europe led to gradual fragmentation, generating multiple genetically differentiated subisolates. The string of population bottlenecks and founder effects have shaped a unique genetic profile, whose potential for genetic research can be met only by study designs that acknowledge cultural tradition and self-identity.

Mutation history of the roma/gypsies

The 8-10 million European Roma/Gypsies are a founder population of common origins that has subsequently split into multiple socially divergent and geographically dispersed Gypsy groups. Unlike other founder populations, whose genealogy has been extensively documented, the demographic history of the Gypsies is not fully understood and, given the lack of written records, has to be inferred from current genetic data. In this study, we have used five disease loci harboring private Gypsy mutations to examine some missing historical parameters and current structure. We analyzed the frequency distribution of the five mutations in 832-1,363 unrelated controls, representing 14 Gypsy populations, and the diversification of chromosomal haplotypes in 501 members of affected families. Sharing of mutations and high carrier rates supported a strong founder effect, and the identity of the congenital myasthenia 1267delG mutation in Gypsy and Indian/Pakistani chromosomes provided the best evidence yet of the Indian origins of the Gypsies. However, dramatic differences in mutation frequencies and haplotype divergence and very limited haplotype sharing pointed to strong internal differentiation and characterized the Gypsies as a founder population comprising multiple subisolates. Using disease haplotype coalescence times at the different loci, we estimated that the entire Gypsy population was founded approximately 32-40 generations ago, with secondary and tertiary founder events occurring approximately 16-25 generations ago. The existence of multiple subisolates, with endogamy maintained to the present day, suggests a general approach to complex disorders in which initial gene mapping could be performed in large families from a single Gypsy group, whereas fine mapping would rely on the informed sampling of the divergent subisolates and searching for the shared genomic region that displays the strongest linkage disequilibrium with the disease.