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.

Representing vulnerable populations in genetic studies: The case of the Roma

Moreau () has raised concerns about the use of DNA data obtained from vulnerable populations, such as the Uighurs in China. We discuss another case, situated in Europe and with a research history dating back 100 years: genetic investigations of Roma. In our article, we focus on problems surrounding representativity in these studies. We claim that many of the circa 440 publications in our sample neglect the methodological and conceptual challenges of representativity. Moreover, authors do not account for problematic misrepresentations of Roma resulting from the conceptual frameworks and sampling schemes they use. We question the representation of Roma as a "genetic isolate" and the underlying rationales, with a strong focus on sampling strategies. We discuss our results against the optimistic prognosis that the "new genetics" could help to overcome essentialist understandings of 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.

Mitochondrial DNA diversity of the Roma from northeastern Bosnia, Bosnia and Herzegovina

This study is the first report on the mtDNA diversity in the Roma population from Bosnia. The main aim of this study was to analyse the mtDNA diversity in the studied population, evaluate the genetic relations with other European Roma populations, and analyse the influences of the Roma gene flow on the mitochondrial profile of the Roma from northeastern Bosnia. MtDNA variability in the analysed population has been studied by means of hypervariable segment I and II (HVSI/II) of the control region sequencing and analysis of restriction fragment-length polymorphisms of the coding region. Our results show that genetic structure of the Roma from northeastern Bosnia has a combination of lineages of three main layers: specific founder of Indian origin (M5a1 and M35b) and founder non-M lineages of Indian/European origin (H7a, X2b and X2d) and lineages of European/Middle East origin (H, H1, H11a, V, T2b, K1b and W). The distribution of the haplogroups in the gene pool of the comparative European Roma populations indicates the separate origin of the Bosnian Roma and the Bulgarian Vlax and Croatian Vlax Roma. The data suggest that mitochondrial gene pool of the Roma population from northeastern Bosnia might be a consequence of early parting and the later different migration routes that are part of their demographic history. Our data confirm the high genetic heterogeneity of the Roma populations that can be shaped by effects of genetic drift, isolation and low effective population size, and this correlates with the migratory history of the Roma.

Characterization of ADME genes variation in Roma and 20 populations worldwide

The products of the polymorphic ADME genes are involved in Absorption, Distribution, Metabolism, and Excretion of drugs. The pharmacogenetic data have been studied extensively due to their clinical importance in the appropriate drug prescription, but such data from the isolated populations are rather scarce. We analyzed the distribution of 95 polymorphisms in 31 core ADME genes in 20 populations worldwide and in newly genotyped samples from the Roma (Gypsy) population living in Croatia. Global distribution of ADME core gene loci differentiated three major clusters; (1) African, (2) East Asian, and (3) joint European, South Asian and South American cluster. The SLCO1B3 (rs4149117) and CYP3A4 (rs2242480) genes differentiated at the highest level the African group of populations, while NAT2 gene loci (rs1208, rs1801280, and rs1799929) and VKORC1 (rs9923231) differentiated East Asian populations. The VKORC1 rs9923231 was among the investigated loci the one with the largest global minor allele frequency (MAF) range; its MAF ranged from 0.027 in Nigeria to 0.924 in Han Chinese. The distribution of the investigated gene loci positions Roma population within the joined European and South Asian clusters, suggesting that their ADME gene pool is a combination of ancestral (Indian) and more recent (European) surrounding, as it was already implied by other genetic markers. However, when compared to the populations worldwide, the Croatian Roma have extreme MAF values in 10 out of the 95 investigated ADME core gene loci. Among loci which have extraordinary MAFs in Roma population two have strong proof of clinical importance: rs1799853 (CYP2C9) for warfarin dosage, and rs12248560 (CYP2C19) for clopidogrel dosage, efficacy and toxicity. This finding confirms the importance of taking the Roma as well as the other isolated populations`genetic profiles into account in pharmaco-therapeutic practice.

Distinctiveness of the Roma population within CYP2B6 worldwide variation

AIM

To determine variation of CYP2B6 gene within the genetically specific Croatian Roma (Gypsy) population originating from India and to examine it in the worldwide perspective.

MATERIALS & METHODS

Seven SNP loci (rs12721655, rs2279343, rs28399499, rs34097093, rs3745274, rs7260329 and rs8192709) were genotyped in 439 subjects using Kompetitive Allele Specific PCR (KASP) method.

RESULTS

The Croatian Roma took an outlying position in CYP2B6 variation from the worldwide perspective mainly due to their exceptionally high minor allele frequency (MAF) for rs8192709 (12.8%), and lower for rs2279343 (21.1%) compared with south Asian populations.

CONCLUSION

This study provides the first data of several CYP2B6 polymorphisms in Roma population and indicates the need for systematic investigation of the most important pharmacogenes' variants in this large, transnationally isolated population worldwide.

Maternal genetic heritage of Southeastern Europe reveals a new Croatian isolate and a novel, local sub-branching in the x2 haplogroup

High mtDNA variation in Southeastern Europe (SEE) is a reflection of the turbulent and complex demographic history of this area, influenced by gene flow from various parts of Eurasia and a long history of intermixing. Our results of 1035 samples (488 from Croatia, 239 from Bosnia and 130 from Herzegovina, reported earlier, and 97 Slovenians and 81 individuals from Žumberak, reported here for the first time) show that the SEE maternal genetic diversity fits within a broader European maternal genetic landscape. The study also shows that the population of Žumberak, located in the continental part of Croatia, developed some unique mtDNA haplotypes and elevated haplogroup frequencies due to distinctive demographic history and can be considered a moderate genetic isolate. We also report seven samples from the Bosnian population and one Herzegovinian sample designated as X2* individuals that could not be assigned to any of its sublineages (X2a'o) according to the existing X2 phylogeny. In an attempt to clarify the phylogeny of our X2 samples, their mitochondrial DNA has been completely sequenced. We suppose that these lineages are signs of local microdifferentiation processes that occurred in the recent demographic past in this area and could possibly be marked as SEE-specific X2 sublineages.

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.

Ancestral modal Y-STR haplotype shared among Romani and South Indian populations

One of the primary unanswered questions regarding the dispersal of Romani populations concerns the geographical region and/or the Indian caste/tribe that gave rise to the proto-Romani group. To shed light on this matter, 161 Y-chromosomes from Roma, residing in two different provinces of Serbia, were analyzed. Our results indicate that the paternal gene pool of both groups is shaped by several strata, the most prominent of which, H1-M52, comprises almost half of each collection's patrilineages. The high frequency of M52 chromosomes in the two Roma populations examined may suggest that they descend from a single founder that has its origins in the Indian subcontinent. Moreover, when the Y-STR profiles of haplogroup H derived individuals in our Roma populations were compared to those typed in the South Indian emigrants from Malaysia and groups from Madras, Karnataka (Lingayat and Vokkaliga castes) and tribal Soligas, sharing of the two most common haplotypes was observed. These similarities suggest that South India may have been one of the contributors to the proto-Romanis. European genetic signatures (i.e., haplogroups E1b1b1a1b-V13, G2a-P15, I-M258, J2-M172 and R1-M173), on the other hand, were also detected in both groups, but at varying frequencies. The divergent European genetic signals in each collection are likely the result of differential gene flow and/or admixture with the European host populations but may also be attributed to dissimilar endogamous practices following the initial founder effect. Our data also support the notion that a number of haplogroups including G2a-P15, J2a3b-M67(xM92), I-M258 and E1b1b1-M35 were incorporated into the proto-Romani paternal lineages as migrants moved from northern India through Southwestern Asia, the Middle East and/or Anatolia into the Balkans.