Reconstructing Roma history from genome-wide data

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.

Methods for Assessing Population Relationships and History Using Genomic Data

Genetic data contain a record of our evolutionary history. The availability of large-scale datasets of human populations from various geographic areas and timescales, coupled with advances in the computational methods to analyze these data, has transformed our ability to use genetic data to learn about our evolutionary past. Here, we review some of the widely used statistical methods to explore and characterize population relationships and history using genomic data. We describe the intuition behind commonly used approaches, their interpretation, and important limitations. For illustration, we apply some of these techniques to genome-wide autosomal data from 929 individuals representing 53 worldwide populations that are part of the Human Genome Diversity Project. Finally, we discuss the new frontiers in genomic methods to learn about population history. In sum, this review highlights the power (and limitations) of DNA to infer features of human evolutionary history, complementing the knowledge gleaned from other disciplines, such as archaeology, anthropology, and linguistics.

Founder lineages in the Iberian Roma mitogenomes recapitulate the Roma diaspora and show the effects of demographic bottlenecks

The Roma are the largest ethnic minority in Europe. With a Northwestern Indian origin around ~ 1.5 kya, they travelled throughout West Asia until their arrival in Europe around the eleventh century CE. Their diaspora through Europe is characterized by population bottlenecks and founder events which have contributed to their present day genetic and cultural diversity. In our study, we focus on the effects of founder effects in the mitochondrial DNA (mtDNA) pool of Iberian Roma by producing and analyzing 144 novel whole mtDNA sequences of Iberian Roma. Over 60% of their mtDNA pool is composed by founder lineages of South Asian origin or acquired by gene flow during their diaspora in the Middle East or locally in Europe in Europe. The TMRCA of these lineages predates the historical record of the Roma arrival in Spain. The abundance of founder lineages is in contrast with ~ 0.7% of autochthonous founder lineages present in the non-Roma Iberian population. Within those founder lineages, we found a substantial amount of South Asian M5a1b1a1 haplotypes and high frequencies of West Eurasian founder lineages (U3b1c, J2b1c, J1c1b, J1b3a, H88, among others), which we characterized phylogenetically and put in phylogeographical context. Besides, we found no evidence of genetic substructure of Roma within the Iberian Peninsula. These results show the magnitude of founder effects in the Iberian Roma and further explain the Roma history and genetic diversity from a matrilineal point of view.

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Reconstructing the history of founder events using genome-wide patterns of allele sharing across individuals

Founder events play a critical role in shaping genetic diversity, fitness and disease risk in a population. Yet our understanding of the prevalence and distribution of founder events in humans and other species remains incomplete, as most existing methods require large sample sizes or phased genomes. Thus, we developed ASCEND that measures the correlation in allele sharing between pairs of individuals across the genome to infer the age and strength of founder events. We show that ASCEND can reliably estimate the parameters of founder events under a range of demographic scenarios. We then apply ASCEND to two species with contrasting evolutionary histories: ~460 worldwide human populations and ~40 modern dog breeds. In humans, we find that over half of the analyzed populations have evidence for recent founder events, associated with geographic isolation, modes of sustenance, or cultural practices such as endogamy. Notably, island populations have lower population sizes than continental groups and most hunter-gatherer, nomadic and indigenous groups have evidence of recent founder events. Many present-day groups––including Native Americans, Oceanians and South Asians––have experienced more extreme founder events than Ashkenazi Jews who have high rates of recessive diseases due their known history of founder events. Using ancient genomes, we show that the strength of founder events differs markedly across geographic regions and time––with three major founder events related to the peopling of Americas and a trend in decreasing strength of founder events in Europe following the Neolithic transition and steppe migrations. In dogs, we estimate extreme founder events in most breeds that occurred in the last 25 generations, concordant with the establishment of many dog breeds during the Victorian times. Our analysis highlights a widespread history of founder events in humans and dogs and elucidates some of the demographic and cultural practices related to these events.

A founder event occurs when small numbers of ancestral individuals give rise to a large fraction of the population. Founder events reduce genetic variation and increase the risk of recessive diseases. Despite their importance in evolutionary and disease studies, we still only have a limited comprehension of their prevalence and properties in humans and other species, as most existing methods require large sample sizes or phased genomes. Here, we present a flexible method, ASCEND, to infer the timing and the strength of founder events that is suitable for sparse datasets with few samples or limited coverage. ASCEND provides reliable estimates across a wide range of demographic scenarios. By applying it to data from two species (humans and dogs), we document a widespread history of recent founder events in both species and provide insights about the demographic processes related to these events. Our analysis helps to identify groups with strong founder events that should be prioritized for future studies as they offer a unique opportunity for biological discovery and reducing disease burden through mapping of recessive disease-causing genes and pathways, as previously shown in studies of Ashkenazi Jews and Finns.

Paternal genetic structure of the Bosnian-Herzegovinian Roma: A Y-chromosomal STR study

OBJECTIVES

Studies indicate the complex nature of the genetic structure of the European Roma which has been shaped by different effects of their demographic history, while preserving their ancestral Indian origin. The primary aims of this study were to present for the first time the paternal profiles of the Roma from Bosnia and Herzegovina based on the data from Y-chromosome STR loci, identify the components of non-Roma paternal gene flow into the Roma, and evaluate the genetic relationships with other European Roma populations.

MATERIALS AND METHODS

In this study, 110 DNA samples of unrelated males from Roma populations residing in different regions of Bosnia and Herzegovina were genotyped using the 23 Y-STR loci included in the PowerPlex Y23 system.

RESULTS

The analysis of the genetic structure of the Bosnian-Herzegovinian Roma revealed intra-country population substructuring and indicated differing genetic affinities between the Bosnian-Herzegovinian Roma and other European Roma populations. The paternal genetic structure of the Bosnian-Herzegovinian Roma has two components: an ancestral component represented by haplogroup H1a1a-M82, and European component presented by haplogroups I1-M253, I2a1a2b-L621, J2a1a-L26, J2a1a1a2b2a3~Z7671, J2b2a-M241, G2a2b2a1a1b-L497, and E1b1b-M215.

CONCLUSION

Genetic relations between the Bosnian-Herzegovinian Roma and other European Roma are shaped by different influences on their demographic history. The data suggest that the paternal gene pool of the Roma from Bosnia and Herzegovina might be a consequence of an early separation of the proto-Roma population and the later gene flow as well as factors of the isolation that accompany the Roma populations in some Bosnian-Herzegovinian regions.

The genetic scenario of Mercheros: an under-represented group within the Iberian Peninsula

BACKGROUND

The general picture of human genetic variation has been vastly depicted in the last years, yet many populations remain broadly understudied. In this work, we analyze for the first time the Merchero population, a Spanish minority ethnic group that has been scarcely studied and historically persecuted. Mercheros have been roughly characterised by an itinerant history, common traditional occupations, and the usage of their own language.

RESULTS

Here, we examine the demographic history and genetic scenario of Mercheros, by using genome-wide array data, whole mitochondrial sequences, and Y chromosome STR markers from 25 individuals. These samples have been complemented with a wide-range of present-day populations from Western Eurasia and North Africa. Our results show that the genetic diversity of Mercheros is explained within the context of the Iberian Peninsula, evidencing a modest signal of Roma admixture. In addition, Mercheros present low genetic isolation and intrapopulation heterogeneity.

CONCLUSIONS

This study represents the first genetic characterisation of the Merchero population, depicting their fine-scale ancestry components and genetic scenario within the Iberian Peninsula. Since ethnicity is not only influenced by genetic ancestry but also cultural factors, other studies from multiple disciplines are needed to further explore the Merchero population. As with Mercheros, there is a considerable gap of underrepresented populations and ethnic groups in publicly available genetic data. Thus, we encourage the consideration of more ethnically diverse population panels in human genetic studies, as an attempt to improve the representation of human populations and better reconstruct their fine-scale history.

Admixture Has Shaped Romani Genetic Diversity in Clinically Relevant Variants

Genetic patterns of inter-population variation are a result of different demographic and adaptive histories, which gradually shape the frequency distribution of the variants. However, the study of clinically relevant mutations has a Eurocentric bias. The Romani, the largest transnational minority ethnic group in Europe, originated in South Asia and received extensive gene flow from West Eurasia. Most medical genetic studies have only explored founder mutations related to Mendelian disorders in this population. Here we analyze exome sequences and genome-wide array data of 89 healthy Spanish Roma individuals to study complex traits and disease. We apply a different framework and focus on variants with both increased and decreased allele frequencies, taking into account their local ancestry. We report several OMIM traits enriched for genes with deleterious variants showing increased frequencies in Roma or in non-Roma (e.g., obesity is enriched in Roma, with an associated variant linked to South Asian ancestry; while non-insulin dependent diabetes is enriched in non-Roma Europeans). In addition, previously reported pathogenic variants also show differences among populations, where some variants segregating at low frequency in non-Roma are virtually absent in the Roma. Lastly, we describe frequency changes in drug-response variation, where many of the variants increased in Roma are clinically associated with metabolic and cardiovascular-related drugs. These results suggest that clinically relevant variation in Roma cannot only be characterized in terms of founder mutations. Instead, we observe frequency differences compared to non-Roma: some variants are absent, while other have drifted to higher frequencies. As a result of the admixture events, these clinically damaging variants can be traced back to both European and South Asian-related ancestries. This can be attributed to a different prevalence of some genetic disorders or to the fact that genetic susceptibility variants are mostly studied in populations of European descent, and can differ in individuals with different ancestries.

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.

Clinical and Genomic Evaluation of 207 Genetic Myopathies in the Indian Subcontinent

Objective: Inherited myopathies comprise more than 200 different individually rare disease-subtypes, but when combined together they have a high prevalence of 1 in 6,000 individuals across the world. Our goal was to determine for the first time the clinical- and gene-variant spectrum of genetic myopathies in a substantial cohort study of the Indian subcontinent.

Methods: In this cohort study, we performed the first large clinical exome sequencing (ES) study with phenotype correlation on 207 clinically well-characterized inherited myopathy-suspected patients from the Indian subcontinent with diverse ethnicities.

Results: Clinical-correlation driven definitive molecular diagnosis was established in 49% (101 cases; 95% CI, 42–56%) of patients with the major contributing pathogenicity in either of three genes, GNE (28%; GNE-myopathy), DYSF (25%; Dysferlinopathy), and CAPN3 (19%; Calpainopathy). We identified 65 variant alleles comprising 37 unique variants in these three major genes. Seventy-eight percent of the DYSF patients were homozygous for the detected pathogenic variant, suggesting the need for carrier-testing for autosomal-recessive disorders like Dysferlinopathy that are common in India. We describe the observed clinical spectrum of myopathies including uncommon and rare subtypes in India: Sarcoglycanopathies (SGCA/B/D/G), Collagenopathy (COL6A1/2/3), Anoctaminopathy (ANO5), telethoninopathy (TCAP), Pompe-disease (GAA), Myoadenylate-deaminase-deficiency-myopathy (AMPD1), myotilinopathy (MYOT), laminopathy (LMNA), HSP40-proteinopathy (DNAJB6), Emery-Dreifuss-muscular-dystrophy (EMD), Filaminopathy (FLNC), TRIM32-proteinopathy (TRIM32), POMT1-proteinopathy (POMT1), and Merosin-deficiency-congenital-muscular-dystrophy-type-1 (LAMA2). Thirteen patients harbored pathogenic variants in >1 gene and had unusual clinical features suggesting a possible role of synergistic-heterozygosity/digenic-contribution to disease presentation and progression.

Conclusions: Application of clinically correlated ES to myopathy diagnosis has improved our understanding of the clinical and genetic spectrum of different subtypes and their overlaps in Indian patients. This, in turn, will enhance the global gene-variant-disease databases by including data from developing countries/continents for more efficient clinically driven molecular diagnostics.

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.

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.

Analysis of a DNA mixture involving Romani reference populations

Here we present an Italian criminal case that shows how statistical methods can be used to extract information from a series of mixed DNA profiles. The case involves several different individuals and a set of different DNA traces. The case possibly involves persons of interest of a small population of Romani origin. First, a brief description of the case is provided. Secondly, we introduce some heuristic tools that can be used to evaluate the data and we also briefly outline the statistical model used for analysing DNA mixtures. Finally, we illustrate some of the findings on the case and discuss further directions of research. The results show how the use of different population database allele frequencies for analysing the DNA mixtures can lead to very different results, some seemingly inculpatory and some seemingly exculpatory. We also illustrate the results obtained from combining the evidence from different samples.

Analysis of the mitochondrial CYTB gene sequence in human populations of northeastern Bosnia

This study offers the first report on variation sequence of the mitochondrial cytochrome b (MTCYTB) gene in populations from Bosnia (northeastern Bosnia). This study was designed on the analysis of the genetic diversity of two populations of different cultural-anthropological and genetic origin, Roma population and native/non-Roma population. The main aim of our study was to estimate the usefulness of the CYTB sequence in the analysis of genetic categorization of different populations and intergroup diversity, as well as to provide some additional information on haplogroup-associated polymorphisms within the CYTB region in defining haplogroup status. Estimation of the genetic diversity was done using intra and intergroup genetic indices. The population-specific polymorphisms have been found in both categories of the populations. The results of the analysis of genetic differentiation show significant pairwise Fst differences between the Romani and native populations. Also, registered significant genetic differentiation is illustrated on the level of genetic variation between subpopulations of the Roma and non-Roma origin. The important result in our study is the confirmation of the significance of the triad of polymorphisms T14783C-G15043A-G15301A, indicating the influence of Asian component of the maternal gene pool on the genetic structure of the studied population of the Roma. Our data show that the haplogroup polymorphisms exist in the CYTB region and can provide useful information on the haplogroups that were defined only by the control region of the mtDNA. The results of this study indicate the region of CYTB gene can be a benefit in providing some additional information in the analysis of genetic structure of human populations and can be additionally applied in population studies.

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.

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.

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.

A comprehensive portrait of Y-STR diversity of Indian populations and comparison with 129 worldwide populations

India, known for its rich cultural, linguistic and ethnic diversity, has attracted the attention of population geneticists to understand its genetic diversity employing autosomal, Y-chromosomal and mitochondrial DNA markers. Y-chromosomal short tandem repeats (Y-STRs) are useful in understanding population substructures and reveal the patrilineal affinities among populations. Previous studies on Indian populations based on Y-STR markers were either limited to restricted number of markers or focused on few selected populations. In this study we genotyped 407 unrelated male individuals from 12 states in India employing the suite of Y-STRs present in PowerPlex Y23 (Promega, Madison, WI, USA). These populations clustered genetically close to each other irrespective of their geographic co-ordinates and were characterized primarily by R1a, H and L haplogroups. Interestingly, comparison with 129 worldwide populations showed genetic affinity of the Indian populations with few populations from Europe and Levantine. This study presents the first pan-Indian landscape of 23 Y-STRs and serves as a useful resource for construction of an Indian Y-STR database.

Statistical methods for detecting admixture

The increasing availability of large-scale autosomal genetic variation data sampled from world-wide geographic areas, coupled with advances in the statistical methodology to analyse these data, is showcasing the power of DNA as a major tool to gain insights into the demographic history of humans and other organisms. Here we review statistical techniques that shed light on a specific aspect of demography: the detection and description of admixture events where two or more genetically distinct groups intermixed at one or more times in the past. In particular we give an overview of some of the widely used methods to identify and describe admixture events using autosomal DNA from unrelated individuals, with a particular focus on analysing biallelic Single-Nucleotide-Polymorphsim (SNP) markers.

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.

Target Organ Damage and Cardiovascular Risk in a Hypertensive Roma Sample Population in Romania

BACKGROUND

The largest European Roma community resides in Romania, but there is still little published data on cardiovascular (CV) risk factors and disease in this group. This study addresses the prevalence of arterial hypertension, associated CV disease risk, and target organ damage (TOD) in a Roma community from Bucharest, Romania.

METHODS

This is a cross-sectional community-based participatory research to assess for CV risk factors, TOD and CV disease, including 806 Roma ethnics (18-83 years) integrated in the local community, 36.16% males. Evaluation included physical examination with blood pressure, pulse wave velocity and anklebrachial measurements, laboratory tests, ECG, echocardiography and fundoscopy.

RESULTS

Prevalence of hypertension was 33.62%, awareness 76.38%, higher in females (p>0.01), and control rate 44.39%. Compared to age-matched normotensives, hypertensives had more left ventricle hypertrophy and more frequently increased pulse pressure. Differences in TOD were attenuated between newly and previously diagnosed, controlled and uncontrolled, hypertensives. Cardiovascular disease was almost absent in normotensives. Ten-year risk for fatal CV disease followed an increasing trend from normotension to long standing hypertension.

CONCLUSION

This is the first dedicated study to thoroughly assess TOD and risk for fatal CV disease in a Romanian Roma population. Hypertension was less prevalent than in the general population, with similar awareness, possibly as a consequence of integration in the surrounding community. Fatal CV disease risk followed the trend of increasing prevalence of risk factors, and hypertension played an important role in its modulation.

Haplotype Sharing Provides Insights into Fine-Scale Population History and Disease in Finland

Finland provides unique opportunities to investigate population and medical genomics because of its adoption of unified national electronic health records, detailed historical and birth records, and serial population bottlenecks. We assembled a comprehensive view of recent population history (≤100 generations), the timespan during which most rare-disease-causing alleles arose, by comparing pairwise haplotype sharing from 43,254 Finns to that of 16,060 Swedes, Estonians, Russians, and Hungarians from geographically and linguistically adjacent countries with different population histories. We find much more extensive sharing in Finns, with at least one ≥ 5 cM tract on average between pairs of unrelated individuals. By coupling haplotype sharing with fine-scale birth records from more than 25,000 individuals, we find that although haplotype sharing broadly decays with geographical distance, there are pockets of excess haplotype sharing; individuals from northeast Finland typically share several-fold more of their genome in identity-by-descent segments than individuals from southwest regions. We estimate recent effective population-size changes through time across regions of Finland, and we find that there was more continuous gene flow as Finns migrated from southwest to northeast between the early- and late-settlement regions than was dichotomously described previously. Lastly, we show that haplotype sharing is locally enriched by an order of magnitude among pairs of individuals sharing rare alleles and especially among pairs sharing rare disease-causing variants. Our work provides a general framework for using haplotype sharing to reconstruct an integrative view of recent population history and gain insight into the evolutionary origins of rare variants contributing to disease.

Tau haplotypes support the Asian ancestry of the Roma population settled in the Basque Country

We examined tau haplotype frequencies in two different ethnical groups from the Basque Country (BC): Roma people and residents of European ancestry (general population). In addition, we analyzed the spatial distribution of tau haplotypes in Eurasian populations to explore the genetic affinities of the Romani groups living in Europe in a broader scope. The 17q21.31 genomic region was characterized through the genotyping of two diagnostic single nucleotide polymorphisms, SNPs (rs10514879 and rs199451), which allow the identification of H1 and H2 haplotypes. A significant heterozygous deficit was detected in the Romani for rs10514879. The H2 haplotype frequency proved to be more than twice in the BC general population (0.283) than in the Roma people (0.127). In contrast, H2 frequency proved to be very similar between Basque and Hungarian Romani, and similar to the H2 frequencies found in northwestern India and Pakistan as well. Several statistical analyses unveiled genetic structuring for the MAPT diversity, mirrored in a significant association between geography and genetic distances, with an upward trend of H2 haplotype frequencies from Asia to Europe. Yet, Roma samples did not fit into this general spatial patterning because of their discrepancy between geographical position and H2 frequency. Despite the long spatial coexistence in the Basque region between the residents of European ancestry and the Roma, the latter have preserved their Asian genetic ancestry. Bearing in mind the lack of geographical barriers between both ethnical groups, these findings support the notion that sociocultural mores might promote assortative matings in human populations.

Refining the South Asian Origin of the Romani people

Background

 Recent genetic studies based on genome-wide Single Nucleotide Polymorphism (SNP) data further investigated the history of Roma and suggested that the source of South Asian ancestry in Roma originates most likely from the Northwest region of India.

Methods

In this study, based also on genome-wide SNP data, we attempted to refine these findings using significantly larger number of European Roma samples, an extended dataset of Indian groups and involving Pakistani groups into the analyses. Our Roma data contained 179 Roma samples. Our extended Indian data consisted of 51 distinct Indian ethnic groups, which provided us a higher resolution of the population living on the Indian subcontinent. We used in this study principal component analysis and other ancestry estimating methods for the study of population relationships, several formal tests of admixture and an improved algorithm for investigating shared IBD segments in order to investigate the main sources of Roma ancestry.

Results

According to our analyses, Roma showed significant IBD sharing of 0.132 Mb with Northwest Indian ethnic groups. The most significant IBD sharings included ethnic groups of Punjab, Rajasthan and Gujarat states. However, we found also significant IBD sharing of 0.087 Mb with ethnic groups living in Pakistan, such as Balochi, Brahui, Burusho, Kalash, Makrani, Pashtun and Sindhi.

Conclusion

Our results show that Northwest India could play an important role in the South Asian ancestry of Roma, however, the origin of Romani people might include the area of Pakistan as well.

Electronic supplementary material

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

Distinct Penetrance of Obesity-Associated Susceptibility Alleles in the Hungarian General and Roma Populations

AIMS

The aim of our study was to explore differences in genetic predisposition to obesity between the Hungarian general and Roma populations.

METHODS

A total of 1,152 samples from the Hungarian Roma population and 1,743 samples from the Hungarian general population were genotyped for 20 single nucleotide polymorphisms (SNPs) associated with the risk of obesity. Two types of multilocus genetic risk scores were constructed to estimate the combined effect of selected SNPs.

RESULTS

Risk allele frequencies differed significantly between the two populations for 11 SNPs, with no enrichment in any of the two study groups. Variants (rs1558902, rs1121980, rs9939609, and rs9941349) in the fat mass and obesity-associated (FTO) gene exhibited strong but ethnicity-independent association with obesity. Genetic risk scores showed stronger associations with obesity in the Roma population compared with the Hungarian general population; however, without significant gene-population interaction.

CONCLUSION

Differences in obesity prevalence between the Hungarian general and Hungarian Roma populations could not be explained by their distinct genetic susceptibility, rather by ethnicity-related environmental and behavioral factors. Nonetheless, particular gene-environment interactions might contribute to the distinct penetrance of the obesity-associated genetic factors in populations of different ethnic backgrounds.

Genetic diversity and divergence among Spanish beef cattle breeds assessed by a bovine high-density SNP chip

The availability of SNP chips for massive genotyping has proven to be useful to genetically characterize populations of domestic cattle and to assess their degree of divergence. In this study, the Illumina BovineHD BeadChip genotyping array was used to describe the genetic variability and divergence among 7 important autochthonous Spanish beef cattle breeds. The within-breed genetic diversity, measured as the marker expected heterozygosity, was around 0.30, similar to other European cattle breeds. The analysis of molecular variance revealed that 94.22% of the total variance was explained by differences within individuals whereas only 4.46% was the result of differences among populations. The degree of genetic differentiation was small to moderate as the pairwise fixation index of genetic differentiation among breeds (F) estimates ranged from 0.026 to 0.068 and the Nei's D genetic distances ranged from 0.009 to 0.016. A neighbor joining (N-J) phylogenetic tree showed 2 main groups of breeds: Pirenaica, Bruna dels Pirineus, and Rubia Gallega on the one hand and Avileña-Negra Ibérica, Morucha, and Retinta on the other. In turn, Asturiana de los Valles occupied an independent and intermediate position. A principal component analysis (PCA) applied to a distance matrix based on marker identity by state, in which the first 2 axes explained up to 17.3% of the variance, showed a grouping of animals that was similar to the one observed in the N-J tree. Finally, a cluster analysis for ancestries allowed assigning all the individuals to the breed they belong to, although it revealed some degree of admixture among breeds. Our results indicate large within-breed diversity and a low degree of divergence among the autochthonous Spanish beef cattle breeds studied. Both N-J and PCA groupings fit quite well to the ancestral trunks from which the Spanish beef cattle breeds were supposed to derive.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

The genetic history of Cochin Jews from India

Cochin Jews form a small and unique community on the Malabar coast in southwest India. While the arrival time of any putative Jewish ancestors of the community has been speculated to have taken place as far back as biblical times (King Solomon’s era), a Jewish community in the Malabar coast has been documented only since the 9th century CE. Here, we explore the genetic history of Cochin Jews by collecting and genotyping 21 community members and combining the data with that of 707 individuals from 72 other Indian, Jewish, and Pakistani populations, together with additional individuals from worldwide populations. We applied comprehensive genome-wide analyses based on principal component analysis, F_ST, ADMIXTURE, identity-by-descent sharing, admixture linkage disequilibrium decay, haplotype sharing, allele sharing autocorrelation decay and contrasting the X chromosome with the autosomes. We find that, as reported by several previous studies, the genetics of Cochin Jews resembles that of local Indian populations. However, we also identify considerable Jewish genetic ancestry that is not present in any other Indian or Pakistani populations (with the exception of the Jewish Bene Israel, which we characterized previously). Combined, Cochin Jews have both Jewish and Indian ancestry. Specifically, we detect a significant recent Jewish gene flow into this community 13–22 generations (~470–730 years) ago, with contributions from Yemenite, Sephardi, and Middle-Eastern Jews, in accordance with historical records. Genetic analyses also point to high endogamy and a recent population bottleneck in this population, which might explain the increased prevalence of some recessive diseases in Cochin Jews.

The Genetics of Bene Israel from India Reveals Both Substantial Jewish and Indian Ancestry

The Bene Israel Jewish community from West India is a unique population whose history before the 18^th century remains largely unknown. Bene Israel members consider themselves as descendants of Jews, yet the identity of Jewish ancestors and their arrival time to India are unknown, with speculations on arrival time varying between the 8th century BCE and the 6th century CE. Here, we characterize the genetic history of Bene Israel by collecting and genotyping 18 Bene Israel individuals. Combining with 486 individuals from 41 other Jewish, Indian and Pakistani populations, and additional individuals from worldwide populations, we conducted comprehensive genome-wide analyses based on F_ST, principal component analysis, ADMIXTURE, identity-by-descent sharing, admixture linkage disequilibrium decay, haplotype sharing and allele sharing autocorrelation decay, as well as contrasted patterns between the X chromosome and the autosomes. The genetics of Bene Israel individuals resemble local Indian populations, while at the same time constituting a clearly separated and unique population in India. They are unique among Indian and Pakistani populations we analyzed in sharing considerable genetic ancestry with other Jewish populations. Putting together the results from all analyses point to Bene Israel being an admixed population with both Jewish and Indian ancestry, with the genetic contribution of each of these ancestral populations being substantial. The admixture took place in the last millennium, about 19–33 generations ago. It involved Middle-Eastern Jews and was sex-biased, with more male Jewish and local female contribution. It was followed by a population bottleneck and high endogamy, which can lead to increased prevalence of recessive diseases in this population. This study provides an example of how genetic analysis advances our knowledge of human history in cases where other disciplines lack the relevant data to do so.

Origins, admixture and founder lineages in European Roma

The Roma, also known as 'Gypsies', represent the largest and the most widespread ethnic minority of Europe. There is increasing evidence, based on linguistic, anthropological and genetic data, to suggest that they originated from the Indian subcontinent, with subsequent bottlenecks and undetermined gene flow from/to hosting populations during their diaspora. Further support comes from the presence of Indian uniparentally inherited lineages, such as mitochondrial DNA M and Y-chromosome H haplogroups, in a significant number of Roma individuals. However, the limited resolution of most genetic studies so far, together with the restriction of the samples used, have prevented the detection of other non-Indian founder lineages that might have been present in the proto-Roma population. We performed a high-resolution study of the uniparental genomes of 753 Roma and 984 non-Roma hosting European individuals. Roma groups show lower genetic diversity and high heterogeneity compared with non-Roma samples as a result of lower effective population size and extensive drift, consistent with a series of bottlenecks during their diaspora. We found a set of founder lineages, present in the Roma and virtually absent in the non-Roma, for the maternal (H7, J1b3, J1c1, M18, M35b, M5a1, U3, and X2d) and paternal (I-P259, J-M92, and J-M67) genomes. This lineage classification allows us to identify extensive gene flow from non-Roma to Roma groups, whereas the opposite pattern, although not negligible, is substantially lower (up to 6.3%). Finally, the exact haplotype matching analysis of both uniparental lineages consistently points to a Northwestern origin of the proto-Roma population within the Indian subcontinent.

Interethnic variability of CYP4F2 (V433M) in admixed population of Roma and Hungarians

AIMS

Pharmacogenetic based dosing recommendations are provided in FDA-approved warfarin label for Caucasians. Evidence of notable difference in dosing algorithms of under-represented populations forced us to explore the genetic variability of CYP4F2 gene in Roma and Hungarian populations.

PATIENTS AND METHODS

484 Roma, 493 Hungarian untreated subjects were genotyped for the CYP4F2*3 (rs2108622) variant by PCR-RFLP assay.

RESULTS AND DISCUSSION

We firstly report, that frequencies of the CYP4F2 rs2108622 GG, GA, AA genotypes and A allele in the Roma population were 46.5%, 42.6%, 10.9% and 32.2%; in Hungarians 50.1%, 42.2%, 7.7% and 22.8%, respectively. Bearing of two minor alleles of CYP4F2 missense variant (AA genotype) modestly explains inter-ethnic differences of studied populations (p<0.08). CYP4F2*3 (V433M) risk allele frequency of Roma (0.32) was in higher range, and of Hungarians (0.23) in lower range, as compared with other world populations.

CONCLUSIONS

Roma have an elevated chance for higher mean warfarin dose, besides a decreased risk of major bleeding events in long-term warfarin use.

Functional variants of lipid level modifier MLXIPL, GCKR, GALNT2, CILP2, ANGPTL3 and TRIB1 genes in healthy Roma and Hungarian populations

The role of triglyceride metabolism in different diseases, such as cardiovascular or cerebrovascular diseases is still under extensive investigations. In genome-wide studies several polymorphisms have been reported, which are highly associated with plasma lipid level changes. Our goal was to examine eight variants: rs12130333 at the ANGPTL3, rs16996148 at the CILP2, rs17321515 at the TRIB1, rs17145738 and rs3812316 of the MLXIPL, rs4846914 at GALNT2, rs1260326 and rs780094 residing at the GCKR loci. A total of 399 Roma (Gypsy) and 404 Hungarian population samples were genotyped using PCR-RFLP method. Significant differences were found between Roma and Hungarian population samples in both MLXIPL variants (C allele frequency of rs17145738: 94.1% vs. 85.6%, C allele frequency of rs3812316: 94.2% vs. 86.8% in Romas vs. in Hungarians, p < 0.05), in ANGPTL3 (T allele frequency of rs1213033: 12.2% vs. 18.5% in Romas vs. Hungarians, p < 0.05) and GALNT2 (G allele frequency of rs4846914: 46.6% vs. 54.5% Romas vs. in Hungarians, p < 0.05), while no differences over SNPs could be verified and the known minor alleles showed no correlation with triglyceride levels in any population samples. The current study revealed fundamental differences of known triglyceride modifying SNPs in Roma population. Failure of finding evidence for affected triglyceride metabolism shows that these susceptibility genes are much less effective compared for example to the apolipoprotein A5 gene.

MARVELD2 (DFNB49) mutations in the hearing impaired Central European Roma population--prevalence, clinical impact and the common origin

Background

In the present study we aimed: 1) To establish the prevalence and clinical impact of DFNB49 mutations in deaf Roma from 2 Central European countries (Slovakia and Hungary), and 2) to analyze a possible common origin of the c.1331+2T>C mutation among Roma and Pakistani mutation carriers identified in the present and previous studies.

Methods

We sequenced 6 exons of the MARVELD2 gene in a group of 143 unrelated hearing impaired Slovak Roma patients. Simultaneously, we used RFLP to detect the c.1331+2T>C mutation in 85 Hungarian deaf Roma patients, control groups of 702 normal hearing Romanies from both countries and 375 hearing impaired Slovak Caucasians. We analyzed the haplotype using 21 SNPs spanning a 5.34Mb around the mutation c.1331+2T>C.

Results

One pathogenic mutation (c.1331+2T>C) was identified in 12 homozygous hearing impaired Roma patients. Allele frequency of this mutation was higher in Hungarian (10%) than in Slovak (3.85%) Roma patients. The identified common haplotype in Roma patients was defined by 18 SNP markers (3.89 Mb). Fourteen common SNPs were also shared among Pakistani and Roma homozygotes. Biallelic mutation carriers suffered from prelingual bilateral moderate to profound sensorineural hearing loss.

Conclusions

We demonstrate different frequencies of the c.1331+2T>C mutation in hearing impaired Romanies from 3 Central European countries. In addition, our results provide support for the hypothesis of a possible common ancestor of the Slovak, Hungarian and Czech Roma as well as Pakistani deaf patients. Testing for the c.1331+2T>C mutation may be recommended in GJB2 negative Roma cases with early-onset sensorineural hearing loss.

Genome-wide insights into the genetic history of human populations

Although mtDNA and the non-recombining Y chromosome (NRY) studies continue to provide valuable insights into the genetic history of human populations, recent technical, methodological and computational advances and the increasing availability of large-scale, genome-wide data from contemporary human populations around the world promise to reveal new aspects, resolve finer points, and provide a more detailed look at our past demographic history. Genome-wide data are particularly useful for inferring migrations, admixture, and fine structure, as well as for estimating population divergence and admixture times and fluctuations in effective population sizes. In this review, we highlight some of the stories that have emerged from the analyses of genome-wide SNP genotyping data concerning the human history of Southern Africa, India, Oceania, Island South East Asia, Europe and the Americas and comment on possible future study directions. We also discuss advantages and drawbacks of using SNP-arrays, with a particular focus on the ascertainment bias, and ways to circumvent it.

A renewal theory approach to IBD sharing

A long genomic segment inherited by a pair of individuals from a single, recent common ancestor is said to be identical-by-descent (IBD). Shared IBD segments have numerous applications in genetics, from demographic inference to phasing, imputation, pedigree reconstruction, and disease mapping. Here, we provide a theoretical analysis of IBD sharing under Markovian approximations of the coalescent with recombination. We describe a general framework for the IBD process along the chromosome under the Markovian models (SMC/SMC'), as well as introduce and justify a new model, which we term the renewal approximation, under which lengths of successive segments are independent. Then, considering the infinite-chromosome limit of the IBD process, we recover previous results (for SMC) and derive new results (for SMC') for the mean number of shared segments longer than a cutoff and the fraction of the chromosome found in such segments. We then use renewal theory to derive an expression (in Laplace space) for the distribution of the number of shared segments and demonstrate implications for demographic inference. We also compute (again, in Laplace space) the distribution of the fraction of the chromosome in shared segments, from which we obtain explicit expressions for the first two moments. Finally, we generalize all results to populations with a variable effective size.

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.

Genetic evidence for recent population mixture in India

Most Indian groups descend from a mixture of two genetically divergent populations: Ancestral North Indians (ANI) related to Central Asians, Middle Easterners, Caucasians, and Europeans; and Ancestral South Indians (ASI) not closely related to groups outside the subcontinent. The date of mixture is unknown but has implications for understanding Indian history. We report genome-wide data from 73 groups from the Indian subcontinent and analyze linkage disequilibrium to estimate ANI-ASI mixture dates ranging from about 1,900 to 4,200 years ago. In a subset of groups, 100% of the mixture is consistent with having occurred during this period. These results show that India experienced a demographic transformation several thousand years ago, from a region in which major population mixture was common to one in which mixture even between closely related groups became rare because of a shift to endogamy.

Inferring admixture histories of human populations using linkage disequilibrium

Long-range migrations and the resulting admixtures between populations have been important forces shaping human genetic diversity. Most existing methods for detecting and reconstructing historical admixture events are based on allele frequency divergences or patterns of ancestry segments in chromosomes of admixed individuals. An emerging new approach harnesses the exponential decay of admixture-induced linkage disequilibrium (LD) as a function of genetic distance. Here, we comprehensively develop LD-based inference into a versatile tool for investigating admixture. We present a new weighted LD statistic that can be used to infer mixture proportions as well as dates with fewer constraints on reference populations than previous methods. We define an LD-based three-population test for admixture and identify scenarios in which it can detect admixture events that previous formal tests cannot. We further show that we can uncover phylogenetic relationships among populations by comparing weighted LD curves obtained using a suite of references. Finally, we describe several improvements to the computation and fitting of weighted LD curves that greatly increase the robustness and speed of the calculations. We implement all of these advances in a software package, ALDER, which we validate in simulations and apply to test for admixture among all populations from the Human Genome Diversity Project (HGDP), highlighting insights into the admixture history of Central African Pygmies, Sardinians, and Japanese.