Refining the South Asian Origin of the Romani people

Population structure of wild soybean (Glycine soja) based on SLAF-seq have implications for its conservation

Background

Glycine soja Sieb. & Zucc. is the wild ancestor from which the important crop plant soybean was bred. G. soja provides important germplasm resources for the breeding and improvement of cultivated soybean crops, however the species is threatened by habitat loss and fragmentation, and is experiencing population declines across its natural range. Understanding the patterns of genetic diversity in G. soja populations can help to inform conservation practices.

Methods

In this study, we analyzed the genetic diversity and differentiation of G. soja at different sites and investigated the gene flow within the species. We obtained 147 G. soja accessions collected from 16 locations across the natural range of the species from China, Korea and Japan. Samples were analyzed using SLAF-seq (Specific-Locus Amplified Fragment Sequencing).

Results

We obtained a total of 56,489 highly consistent SNPs. Our results suggested that G. soja harbors relatively high diversity and that populations of this species are highly differentiated. The populations harboring high genetic diversity, especially KR, should be considered first when devising conservation plans for the protection of G. soja, and in situ protection should be adopted in KR. G. soja populations from the Yangtze River, the Korean peninsula and northeastern China have a close relationship, although these areas are geographically disconnected. Other populations from north China clustered together. Analysis of gene flow suggested that historical migrations of G. soja may have occurred from the south northwards across the East-Asia land-bridge, but not across north China. All G. soja populations could be divided into one of two lineages, and these two lineages should be treated separately when formulating protection policies.

Population history modulates the fitness effects of Copy Number Variation in the Roma

We provide the first whole genome Copy Number Variant (CNV) study addressing Roma, along with reference populations from South Asia, the Middle East and Europe. Using CNV calling software for short-read sequence data, we identified 3171 deletions and 489 duplications. Taking into account the known population history of the Roma, as inferred from whole genome nucleotide variation, we could discern how this history has shaped CNV variation. As expected, patterns of deletion variation, but not duplication, in the Roma followed those obtained from single nucleotide polymorphisms (SNPs). Reduced effective population size resulting in slightly relaxed natural selection may explain our observation of an increase in intronic (but not exonic) deletions within Loss of Function (LoF)-intolerant genes. Over-representation analysis for LoF-intolerant gene sets hosting intronic deletions highlights a substantial accumulation of shared biological processes in Roma, intriguingly related to signaling, nervous system and development features, which may be related to the known profile of private disease in the population. Finally, we show the link between deletions and known trait-related SNPs reported in the genome-wide association study (GWAS) catalog, which exhibited even frequency distributions among the studied populations. This suggests that, in general human populations, the strong association between deletions and SNPs associated to biomedical conditions and traits could be widespread across continental populations, reflecting a common background of potentially disease/trait-related CNVs.

Genetic Susceptibility to Insulin Resistance and Its Association with Estimated Longevity in the Hungarian General and Roma Populations

Diabetes mellitus is a major public health problem with a wide range of prevalence among different ethnic groups. Early recognition of pre-diabetes is important to prevent the development of the disease, its complications, co-morbidities, and consequently early death. Insulin resistance (IR) is considered a condition that precedes type 2 diabetes; thus, understanding its underlying causes (genetic and non-genetic factors) will bring us closer to preventing it. The present study aimed to investigate the genetic susceptibility to IR and its impact on estimated longevity in populations with different ethnic origins using randomly selected samples of 372 Hungarian general (HG, as a reference with Caucasian origin) and 334 Roma participants (largest ethnic minority in Europe, with a northern India origin). In the present study, we used the Homeostasis Model Assessment—Insulin Resistance (HOMA—IR) to identify people with IR (>3.63) at the population level. To investigate the genetic predisposition to IR, 29 single nucleotide polymorphisms (SNPs) identified in a systematic literature search were selected and genotyped in sample populations. In the analyses, the adjusted p < 0.0033 was considered significant. Of these 29 SNPs, the commutative effects of 15 SNPs showing the strongest association with HOMA—IR were used to calculate an optimized genetic risk score (oGRS). The oGRS was found nominally significantly (p = 0.019) higher in the Roma population compared to HG one, and it was more strongly correlated with HOMA—IR. Therefore, it can be considered as a stronger predictor of the presence of IR among the Roma (AUCRoma = 0.673 vs. AUCHG = 0.528). Furthermore, oGRS also showed a significant correlation with reduced estimated longevity in the Roma population (β = −0.724, 95% CI: −1.230−−0.218; p = 0.005), but not in the HG one (β = 0.065, 95% CI: −0.388−0.518; p = 0.779). Overall, IR shows a strong correlation with a genetic predisposition among Roma, but not in the HG population. Furthermore, the increased genetic risk of Roma is associated with shorter estimated longevity, whereas this association is not observed in the HG one. Increased genetic susceptibility of Roma to IR should be considered in preventive programs targeting the development of type 2 diabetes, which may also reduce the risk of preventable premature death among them.

The Counteracting Effects of Demography on Functional Genomic Variation: The Roma Paradigm

Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations.

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.

The shaping of immunological responses through natural selection after the Roma Diaspora

The Roma people are the largest transnational ethnic minority in Europe and can be considered the last human migration of South Asian origin into the continent. They left Northwest India approximately 1,000 years ago, reaching the Balkan Peninsula around the twelfth century and Romania in the fourteenth century. Here, we analyze whole-genome sequencing data of 40 Roma and 40 non-Roma individuals from Romania. We performed a genome-wide scan of selection comparing Roma, their local host population, and a Northwestern Indian population, to identify the selective pressures faced by the Roma mainly after they settled in Europe. We identify under recent selection several pathways implicated in immune responses, among them cellular metabolism pathways known to be rewired after immune stimulation. We validated the interaction between PIK3-mTOR-HIF-1α and cytokine response influenced by bacterial and fungal infections. Our results point to a significant role of these pathways for host defense against the most prevalent pathogens in Europe during the last millennium.

Prevalence of Insulin Resistance in the Hungarian General and Roma Populations as Defined by Using Data Generated in a Complex Health (Interview and Examination) Survey

Data mainly from one-off surveys clearly show that the health of Roma, the largest ethnic minority of Europe, is much worse than that of the general population. However, results from comprehensive exploratory studies are missing. The aim of our study was to create a complex database for comparative and association studies to better understand the background of the very unfavourable health of Roma, especially the high burden of cardiometabolic diseases. A three-pillar (questionnaire-based, physical and laboratory examinations) health survey was carried out on randomly selected samples of the Hungarian general (HG, n = 417) and Roma (HR, n = 415) populations, and a database consisting of more than half a million datapoints was created. Using selected data, the prevalence rates of metabolic syndrome (MetS) and of its components were determined, and to estimate the risk of insulin resistance (IR), surrogate measures (the homeostasis model assessment of insulin resistance index, quantitative insulin sensitivity check index, McAuley and TyG indices and the TG/HDL-C ratio) were calculated. Receiver operating characteristic curve analysis and Youden’s method were used to define the optimal cut-off values of each IR index. The prevalence of MetS was very high in both study populations (HG: 39.8%, HR: 44.0%) with no statistically significant difference between the two groups in females or males. The prevalence of MetS showed a very marked increase in the HR 35–49 years age group. Among surrogate measures, the TyG index showed the greatest power for predicting IR/MetS at a cut-off value of 4.69 (77% sensitivity, 84% specificity) and indicated a 42.3% (HG) and 40.5% (HR) prevalence of IR. The prevalence of MetS and IR is almost equally very unfavourable in both groups; thus, the factors underlying the high premature mortality burden of Roma should be further clarified by investigating the full spectrum of risk factors available in the database, with a special focus on the access of Roma people to preventive and curative health services.

Applicability of Obesity-Related SNPs and their Effect Size Measures Defined on Populations with European Ancestry for Genetic Risk Estimation among Roma

Investigations on the impact of genetic factors on the development of obesity have been limited regarding the Roma population—the largest and most vulnerable ethnic minority in Europe of Asian origin. Genetic variants identified from genetic association studies are primarily from European populations. With that in mind, we investigated the applicability of data on selected obesity-related single nucleotide polymorphisms (SNPs), obtained from the Hungarian general (HG) population of European origin, on the Hungarian Roma (HR) population. Twenty preselected SNPs in susceptible alleles, known to be significantly associated with obesity-related phenotypes, were used to estimate the effect of these SNPs on body mass index (BMI) and waist circumference (WC) in HG (N = 1783) and HR (N = 1225) populations. Single SNP associations were tested using linear and logistic regression models, adjusted for known covariates. Out of 20 SNPs, four located in FTO (rs1121980, rs1558902, rs9939609, and rs9941349) showed strong association with BMI and WC as continuous variables in both samples. Computations based on Adult Treatment Panel III (ATPIII) and the International Diabetes Federation’s (IDF) European and Asian criteria showed rs9941349 in FTO to be associated only with WC among both populations, and two SNPs (rs2867125, rs6548238) in TMEM18 associated with WC only in HG population. A substantial difference (both in direction and effect size) was observed only in the case of rs1801282 in PPARγ on WC as a continuous outcome. Findings suggest that genetic risk scores based on counting SNPs with relatively high effect sizes, defined based on populations with European ancestry, can sufficiently allow estimation of genetic susceptibility for Roma. Further studies are needed to clarify the role of SNP(s) with protective effect(s).

Revealing the Genetic Impact of the Ottoman Occupation on Ethnic Groups of East-Central Europe and on the Roma Population of the Area

History of East-Central Europe has been intertwined with the history of Turks in the past. A significant part of this region of Europe has been fallen under Ottoman control during the 150 years of Ottoman occupation in the 16–17th centuries. The presence of the Ottoman Empire affected this area not only culturally but also demographically. The Romani people, the largest ethnic minority of the East-Central European area, share an even more eventful past with Turkish people from the time of their migration throughout Eurasia and they were a notable ethnic group in East-Central Europe in the Ottoman era already. The relationship of Turks with East-Central European ethnic groups and with regional Roma ethnicity was investigated based on genome-wide autosomal single nucleotide polymorphism data. Population structure analysis, ancestry estimation, various formal tests of admixture and DNA segment analyses were carried out in order to shed light to the conclusion of these events on a genome-wide basis. Analyses show that the Ottoman occupation of Europe left detectable impact in the affected East-Central European area and shaped the ancestry of the Romani people as well. We estimate that the investigated European populations have an average identity-by-descent share of 0.61 with Turks, which is notable, compared to other European populations living in West and North Europe far from the affected area, and compared to the share of Sardinians, living isolated from these events. Admixture of Roma and Turks during the Ottoman rule show also high extent.

Hearing impairment locus heterogeneity and identification of PLS1 as a new autosomal dominant gene in Hungarian Roma

Roma are a socially and culturally distinct isolated population with genetically divergent subisolates, residing mainly across Central, Southern, and Eastern Europe. We evaluated the genetic etiology of hearing impairment (HI) in 15 Hungarian Roma families through exome sequencing. A family with autosomal dominant non-syndromic HI segregating a rare variant in the Calponin-homology 2 domain of PLS1, or Plastin 1 [p.(Leu363Phe)] was identified. Young adult Pls1 knockout mice have progressive HI and show morphological defects to their inner hair cells. There is evidence that PLS1 is important in the preservation of adult stereocilia and normal hearing. Four families segregated the European ancestral variant c.35delG [p.(Gly12fs)] in GJB2, and one family was homozygous for p.(Trp24*), an Indian subcontinent ancestral variant which is common amongst Roma from Slovakia, Czech Republic, and Spain. We also observed variants in known HI genes USH1G, USH2A, MYH9, MYO7A, and a splice site variant in MANBA (c.2158-2A>G) in a family with HI, intellectual disability, behavioral problems, and respiratory inflammation, which was previously reported in a Czech Roma family with similar features. Lastly, using multidimensional scaling and ADMIXTURE analyses, we delineate the degree of Asian/European admixture in the HI families understudy, and show that Roma individuals carrying the GJB2 p.(Trp24*) and MANBA c.2158-2A>G variants have a more pronounced South Asian background, whereas the other hearing-impaired Roma display an ancestral background similar to Europeans. We demonstrate a diverse genetic HI etiology in the Hungarian Roma and identify a new gene PLS1, for autosomal dominant human non-syndromic HI.

Genetic profiling revealed an increased risk of venous thrombosis in the Hungarian Roma population

BACKGROUND

Besides modifiable risk factors, genetic susceptibility may also explain the high cardiovascular disease burden of the Roma population.

OBJECTIVES

Aim of this study was to define the genetic susceptibility of Hungarian Roma to venous thrombosis (VT) and comparing it to that of the general population.

METHODS

Fifty-two SNPs associated with VT (in F2, F5, F9, F11, F15, FGA, FGB, FGG, CYP4V2, KLKB1 and vWF) were selected and analyzed in the group of Roma (N = 962) and general (N = 1492) subjects collected by cross-sectional studies. Allele frequencies and genetic risk scores (GRS, unweighted and weighted) were computed for the study groups and compared to estimate the joint effects of SNPs.

RESULTS

The majority of the susceptible alleles were more prevalent in the Roma population, and both GRS and wGRS were found to be significantly higher in Roma than in the general population (GRS: 41.83 ± 5.78 vs. 41.04 ± 6.04; wGRS: 7.78 ± 1.28 vs. 7.46 ± 1.33, p = .001). Only 2.39% of subjects in the Roma population were in the bottom fifth of the wGRS (wGRS≤0.19) compared with 3.62% of those in the general population (p = .080); 2.88% of the general subjects were in the top fifth of the wGRS (wGRS≥10.02), while 4.26% of the Roma population were (p = .066).

CONCLUSION

In conclusion, the Roma population seems to have increased genetic susceptibility to VT. This might have important implications in the future in identifying possible new opportunities for targeted prevention and treatment for those subgroups in the populations who are at greater risk for VT development.

The Origin of a Coastal Indigenous Horse Breed in China Revealed by Genome-Wide SNP Data

The Jinjiang horse is a unique Chinese indigenous horse breed distributed in the southern coastal areas, but the ancestry of Jinjiang horses is not well understood. Here, we used Equine SNP70 Bead Array technology to genotype 301 horses representing 10 Chinese indigenous horse breeds, and we integrated the published genotyped data of 352 individuals from 14 foreign horse breeds to study the relationships between Jinjiang horses and horse breeds from around the world. Principal component analysis (PCA), linkage disequilibrium (LD), runs of homozygosity (ROH) analysis, and ancestry estimating methods were conducted to study the population relationships and the ancestral sources and genetic structure of Jinjiang horses. The results showed that there is no close relationship between foreign horse breeds and Jinjiang horses, and Jinjiang horses shared a similar genetic background with Baise horses. TreeMix analysis revealed that there was gene flow from Chakouyi horses to Jinjiang horses. The ancestry analysis showed that Baise horses and Chakouyi horses are the most closely related ancestors of Jinjiang horses. In conclusion, our results showed that Jinjiang horses have a native origin and that Baise horses and Chakouyi horses were key ancestral sources of Jinjiang horses. The study also suggested that ancient trade activities and the migration of human beings had important effects on indigenous horse breeds in China.

Generalizability and applicability of results obtained from populations of European descent regarding the effect direction and size of HDL-C level-associated genetic variants to the Hungarian general and Roma populations

OBJECTIVES

Large-scale association studies that mainly involve European populations identified many genetic loci related to high-density lipoprotein cholesterol (HDL-C) levels, one of the most important indicators of the risk for cardiovascular diseases. The question with intense speculation of whether the effect estimates obtained from European populations for different HDL-C level-related SNPs are applicable to the Roma ethnicity, the largest minority group in Europe with a South Asian origin, was addressed in the present study.

DESIGN

The associations between 21 SNPs (in the genes LIPC(G), CETP, GALNT2, HMGCP, ABCA1, KCTD10 and WWOX) and HDL-C levels were examined separately in adults of the Hungarian general (N = 1542) and Roma (N = 646) populations by linear regression. Individual effects (direction and size) of single SNPs on HDL-C levels were computed and compared between the study groups and with data published in the literature.

RESULTS

Significant associations between SNPs and HDL-C levels were more frequently found in general subjects than in Roma subjects (11 SNPs in general vs. 4 SNPs in Roma). The CETP gene variants rs1532624, rs708272 and rs7499892 consistently showed significant associations with HDL-C levels across the study groups (p ˂ 0.05), indicating a possible causal variant(s) in this region. Although nominally significant differences in effect size were found for three SNPs (rs693 in gene APOB, rs9989419 in gene CETP, and rs2548861 in gene WWOX) by comparing the general and Roma populations, most of these SNPs did not have a significant effect on HDL-C levels. The β coefficients for SNPs in the Roma population were found to be identical both in direction and magnitude to the effect obtained previously in large-scale studies on European populations.

CONCLUSIONS

The effect of the vast majority of the SNPs on HDL-C levels could be replicated in the Hungarian general and Roma populations, which indicates that the effect size measurements obtained from the literature can be used for risk estimation for both populations.

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.

Revealing the impact of the Caucasus region on the genetic legacy of Romani people from genome-wide data

Romani people are a significant minority in Europe counting about 10 million individuals scattered throughout the continent. They are a migratory group originating from Northwestern India. Their exodus from India occurred approximately 1000–1500 years ago. The migration route of the Romani people was reconstructed with the help of cultural anthropology, linguistics and historical records. Their migration made them through Central Asia, Middle East and the Caucasus region, prior to the arriving into Europe. Yet the significance of these regions, especially of the Caucasus, in Roma ancestry was a rather neglected topic. Contribution of the Caucasus and further affected regions to the ancestry of Roma was investigated based on genome-wide autosomal marker data. 158 European Roma samples and 41 populations from the Caucasus region, from Middle East, Central Asia and from South Asia were considered in our tests. Population structure and ancestry analysis algorithms were applied to investigate the relationship of Roma with these populations. Identical by descent DNA segment analyses and admixture linkage disequilibrium based tests were also applied. Our results suggest that the Caucasus region plays also a significant role in the genetic legacy of Romani people besides the main sources, Europe and South Asia, previously investigated by other population genetic studies. The Middle East and Central Asia seems slightly less important but far from negligible in connection with the sources of Roma ancestry. Our results point out that the Caucasus region and altogether the area of the Caspian and Black Seas had a significant role in the migration of Romani people towards Europe and contributed significantly to the genetic legacy of Roma rival to the European and Indian main sources.