Understanding the genomic heterogeneity of North African Imazighen: from broad to microgeographical perspectives

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare.

The Impact of Recent Demography on Functional Genetic Variation in North African Human Groups

The strategic location of North Africa has made the region the core of a wide range of human demographic events, including migrations, bottlenecks, and admixture processes. This has led to a complex and heterogeneous genetic and cultural landscape, which remains poorly studied compared to other world regions. Whole-exome sequencing is particularly relevant to determine the effects of these demographic events on current-day North Africans' genomes, since it allows to focus on those parts of the genome that are more likely to have direct biomedical consequences. Whole-exome sequencing can also be used to assess the effect of recent demography in functional genetic variation and the efficacy of natural selection, a long-lasting debate. In the present work, we use newly generated whole-exome sequencing and genome-wide array genotypes to investigate the effect of demography in functional variation in 7 North African populations, considering both cultural and demographic differences and with a special focus on Amazigh (plur. Imazighen) groups. We detect genetic differences among populations related to their degree of isolation and the presence of bottlenecks in their recent history. We find differences in the functional part of the genome that suggest a relaxation of purifying selection in the more isolated groups, allowing for an increase of putatively damaging variation. Our results also show a shift in mutational load coinciding with major demographic events in the region and reveal differences within and between cultural and geographic groups.

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.

Whole mitogenomes reveal that NW Africa has acted both as a source and a destination for multiple human movements

Despite being enclosed between the Mediterranean Sea and the Sahara Desert, North Africa has been the scenario of multiple human migrations that have shaped the genetic structure of its present-day populations. Despite its richness, North Africa remains underrepresented in genomic studies. To overcome this, we have sequenced and analyzed 264 mitogenomes from the Algerian Chaoui-speaking Imazighen (a.k.a. Berbers) living in the Aurès region. The maternal genetic composition of the Aurès is similar to Arab populations in the region, dominated by West Eurasian lineages with a moderate presence of M1/U6 North African and L sub-Saharan lineages. When focusing on the time and geographic origin of the North African specific clades within the non-autochthonous haplogroups, different geographical neighboring regions contributed to the North African maternal gene pool during time periods that could be attributed to previously suggested admixture events in the region, since Paleolithic times to recent historical movements such as the Arabization. We have also observed the role of North Africa as a source of geneflow mainly in Southern European regions since Neolithic times. Finally, the present work constitutes an effort to increase the representation of North African populations in genetic databases, which is key to understand their history.

Population analysis of complete mitogenomes for 334 samples from El Salvador

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980-1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war.

Identifying signatures of positive selection in human populations from North Africa

Because of its location, North Africa (NA) has witnessed continuous demographic movements with an impact on the genomes of present-day human populations. Genomic data describe a complex scenario with varying proportions of at least four main ancestry components: Maghrebi, Middle Eastern-, European-, and West-and-East-African-like. However, the footprint of positive selection in NA has not been studied. Here, we compile genome-wide genotyping data from 190 North Africans and individuals from surrounding populations, investigate for signatures of positive selection using allele frequencies and linkage disequilibrium-based methods and infer ancestry proportions to discern adaptive admixture from post-admixture selection events. Our results show private candidate genes for selection in NA involved in insulin processing (KIF5A), immune function (KIF5A, IL1RN, TLR3), and haemoglobin phenotypes (BCL11A). We also detect signatures of positive selection related to skin pigmentation (SLC24A5, KITLG), and immunity function (IL1R1, CD44, JAK1) shared with European populations and candidate genes associated with haemoglobin phenotypes (HPSE2, HBE1, HBG2), other immune-related (DOCK2) traits, and insulin processing (GLIS3) traits shared with West and East African populations. Finally, the SLC8A1 gene, which codifies for a sodium-calcium exchanger, was the only candidate identified under post-admixture selection in Western NA.

Indigenous people from Amazon show genetic signatures of pathogen-driven selection

Ecological conditions in the Amazon rainforests are historically favorable for the transmission of numerous tropical diseases, especially vector-borne diseases. The high diversity of pathogens likely contributes to the strong selective pressures for human survival and reproduction in this region. However, the genetic basis of human adaptation to this complex ecosystem remains unclear. This study investigates the possible footprints of genetic adaptation to the Amazon rainforest environment by analyzing the genomic data of 19 native populations. The results based on genomic and functional analysis showed an intense signal of natural selection in a set of genes related to Trypanosoma cruzi infection, which is the pathogen responsible for Chagas disease, a neglected tropical parasitic disease native to the Americas that is currently spreading worldwide.

Genomic Insights into the Population History of the Resande or Swedish Travelers

The Resande are a minority ethnic group in Sweden, who were characterized by an itinerant way of life, and they have been suggested to originate from the mixture between Swedish and Romani populations. Because the population history of the Resande has been scarcely studied, we analyzed genome-wide genotype array data from unrelated Resande individuals in order to shed light on their origins and demographic history for the first time from a genetic perspective. Our results confirm the Romani-related ancestry of this population and suggest an admixture event between a Romani-like population and a general Swedish-like population that occurred approximately between the mid-18th and mid-19th centuries, two centuries after the arrival of the first historically reported Romani families in Sweden. This inferred date suggests that the Romani group involved in the admixture is related to the pre-18th-century arrivals of Romani in Scandinavia. In addition, a reduction in the population size is detected previous to the admixture event, suggesting a subtle signal of isolation. The present work constitutes a step forward toward a better representation of ethnic minorities and underrepresented groups in population genetic analyses. In order to know in more detail the complete history of human populations, it is time to focus on studying populations that have not been previously considered for a general scenario and that can provide valuable information to fill in the gaps that still remain uncovered.

The genomic analysis of current-day North African populations reveals the existence of trans-Saharan migrations with different origins and dates

The Sahara Desert has acted as a barrier to human gene-flow between the northern and central parts of Africa since its aridification. Nonetheless, some contacts between both sides of the desert have occurred throughout history, mainly driven by commercial activity. Part of this was the infamous trans-Saharan slave trade, which forcedly brought peoples from south of the Sahara to North Africa from Roman times until the nineteenth century. Although historical records exist, the genetic aspects of these trans-Saharan migrations have not been deeply studied. In the present study, we assess the genetic influence of trans-Saharan migrations in current-day North Africa and characterize its amount, geographical origin, and dates. We confirm the heterogeneous and generally low-frequency presence of genomic segments of sub-Saharan origin in present-day North Africans acquired in recent historical times, and we show evidence of at least two admixture events: one dated around the thirteenth-fourteenth centuries CE between North Africans and a Western-sub-Saharan-like source similar to current-day Senegambian populations, and another one dated around the seventeenth century CE involving Tunisians and an Eastern-sub-Saharan-like source related to current-day south-Sudan and Kenyan populations. Time and location coincide with the peak of trans-Saharan slave-trade activity between Western African empires and North African powers, and are also concordant with the possibility of continuous recent south-to-north gene-flow. These findings confirm the trans-Saharan human genetic contacts, providing new and precise evidence about its possible dates and geographical origins, which are pivotal to understanding the genomic composition of an underrepresented region such as North Africa.

Evidence of selection in the ectodysplasin pathway among endangered aquatic mammals

Synopsis

The ectodysplasin pathway has been a target of evolution repeatedly. Genetic variation in the key genes of this pathway (EDA, EDAR, and EDARADD) results in a rich source of pleiotropic effects across ectodermally-derived structures, including teeth, hair, sweat glands, and mammary glands. In addition, a non-canonical Wnt pathway has a very similar functional role, making variation in the WNT10A gene also of evolutionary significance. The adaptation of mammals to aquatic environments has occurred independently in at least 4 orders, whose species occupy a wide geographic range (from equatorial to polar regions) and exhibit great phenotypic variation in ectodermally-derived structures, including the presence or absence of fur and extreme lactational strategies. The role of the ectodysplasin pathway in the adaptation to aquatic environments has been never explored in mammalian species. In the present study, we analyze the genetic variation in orthologous coding sequences from EDA, EDAR, EDARADD, and WNT10A genes together with ectodermally-derived phenotypic variation from 34 aquatic and non-aquatic mammalian species to assess signals of positive selection, gene-trait coevolution, and genetic convergence. Our study reveals strong evidence of positive selection in a proportion of coding sites in EDA and EDAR genes in 3 endangered aquatic mammals (the Hawaiian monk seal, the Yangtze finless porpoise, and the sea otter). We hypothesize functional implications potentially related to the adaptation to the low-latitude aquatic environment in the Hawaiian monk seal and the freshwater in the Yangtze finless porpoise. The signal in the sea otter is likely the result of an increased genetic drift after an intense bottleneck and reduction of genetic diversity. Besides positive selection, we have not detected robust signals of gene-trait coevolution or convergent amino acid shifts in the ectodysplasin pathway associated with shared phenotypic traits among aquatic mammals. This study provides new evidence of the evolutionary role of the ectodysplasin pathway and encourages further investigation, including functional studies, to fully resolve its relationship with mammalian aquatic adaptation.

Spanish

La vía de la ectodisplasina ha sido objeto de la evolución repetidamente. La variación genética en los principales genes de esta vía (EDA, EDAR y EDARADD) da como resultado una gran diversidad de efectos pleiotrópicos en las estructuras derivadas del ectodermo, incluidos los dientes, el cabello, las glándulas sudoríparas y las glándulas mamarias. Además, una vía wnt no canónica tiene un papel funcional muy similar, por lo que la variación en el gen WNT10A también tiene importancia evolutiva. La adaptación de los mamíferos a los entornes acuáticos se ha producido de forma independiente en al menos cuatro órdenes, cuyas especies ocupan un amplio rango geográfico (desde regiones ecuatoriales a polares) y presentan una gran variación fenotípica en las estructuras derivadas del ectodermo, incluyendo la presencia o ausencia de pelaje y estrategias de lactancia muy diferentes. El papel de la vía de la ectodisplasina en la adaptación a entornos acuáticos no se ha explorado nunca en especies de mamíferos. En este estudio, analizamos la variación genética en las secuencias codificantes ortólogas de los genes EDA, EDAR, EDARADD y WNT10A junto con la variación fenotípica derivada del ectodermo de 34 especies de mamíferos acuáticos y no acuáticos para evaluar señales de selección positiva, coevolución gen-rasgo y convergencia genética. Nuestro estudio revela señales de selección positiva en regiones de las secuencias codificantes de los genes EDA y EDAR en tres mamíferos acuáticos en peligro de extinción (la foca monje de Hawái, la marsopa lisa y la nutria marina). Estas señales podrían tener implicaciones funcionales potencialmente relacionadas con la adaptación al entorno acuático de baja latitud en la foca monje de Hawái y el agua dulce en la marsopa lisa. La señal en la nutria marina es probablemente el resultado de una mayor deriva genética tras un intenso un cuello de botella y una reducción de la diversidad genética. A parte de selección positiva, no hemos detectado señales sólidas de coevolución gen-rasgo o cambios convergentes de aminoácidos en la vía de la ectodisplasina asociados a rasgos fenotípicos compartidos entre mamíferos acuáticos. Este estudio proporciona nuevas evidencias del papel evolutivo de la vía de la ectodisplasina y quiere promover futuras investigaciones con estudios funcionales para acabar de resolver la relación de esta vía con la adaptación acuática de los mamíferos.

Pharmacological characterization of GLPG1972/S201086, a potent and selective small-molecule inhibitor of ADAMTS5

OBJECTIVE

A disintegrin and metalloproteinase with thrombospondin motifs 5 (ADAMTS5) is a key enzyme in degradation of cartilage in osteoarthritis (OA). We report the pharmacological characterization of GLPG1972/S201086, a new, potent and selective small-molecule ADAMTS5 inhibitor.

METHODS

Potency and selectivity of GLPG1972/S201086 for ADAMTS5 were determined using fluorescently labeled peptide substrates. Inhibitory effects of GLPG1972/S201086 on interleukin-1α-stimulated glycosaminoglycan release in mouse femoral head cartilage explants and on interleukin-1β-stimulated release of an ADAMTS5-derived aggrecan neoepitope (quantified with ELISA) in human articular cartilage explants were determined. In the destabilization of the medial meniscus (DMM) mouse and menisectomized (MNX) rat models, effects of oral GLPG1972/S201086 on relevant OA histological and histomorphometric parameters were evaluated.

RESULTS

GLPG1972/S201086 inhibited human and rat ADAMTS5 (IC50 ± SD: 19 ± 2 nM and <23 ± 1 nM, respectively), with 8-fold selectivity over ADAMTS4, and 60->5,000-fold selectivity over other related proteases in humans. GLPG1972/S201086 dose-dependently inhibited cytokine-stimulated aggrenolysis in mouse and human cartilage explants (100% at 20 μM and 10 μM, respectively). In DMM mice, GLPG1972/S201086 (30-120 mg/kg b.i.d) vs vehicle reduced femorotibial cartilage proteoglycan loss (23-37%), cartilage structural damage (23-39%) and subchondral bone sclerosis (21-36%). In MNX rats, GLPG1972/S201086 (10-50 mg/kg b.i.d) vs vehicle reduced cartilage damage (OARSI score reduction, 6-23%), and decreased proteoglycan loss (∼27%) and subchondral bone sclerosis (77-110%).

CONCLUSIONS

GLPG1972/S201086 is a potent, selective and orally available ADAMTS5 inhibitor, demonstrating significant protective efficacy on both cartilage and subchondral bone in two relevant in vivo preclinical OA models.

Population Histories and Genomic Diversity of South American Natives

South America is home to one of the most culturally diverse present-day native populations. However, the dispersion pattern, genetic substructure, and demographic complexity within South America are still poorly understood. Based on genome-wide data of 58 native populations, we provide a comprehensive scenario of South American indigenous groups considering the genomic, environmental, and linguistic data. Clear patterns of genetic structure were inferred among the South American natives, presenting at least four primary genetic clusters in the Amazonian and savanna regions and three clusters in the Andes and Pacific coast. We detected a cline of genetic variation along a west-east axis, contradicting a hard Andes-Amazon divide. This longitudinal genetic variation seemed to have been shaped by both serial population bottlenecks and isolation by distance. Results indicated that present-day South American substructures recapitulate ancient macroregional ancestries and western Amazonia groups show genetic evidence of cultural exchanges that led to language replacement in precontact times. Finally, demographic inferences pointed to a higher resilience of the western South American groups regarding population collapses caused by the European invasion and indicated precontact population reductions and demic expansions in South America.

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.

Ancient and modern mitogenomes from Central Argentina: new insights into population continuity, temporal depth and migration in South America

The inverted triangle shape of South America places Argentina territory as a geographical crossroads between the two principal peopling streams that followed either the Pacific or the Atlantic coasts, which could have then merged in Central Argentina (CA). Although the genetic diversity from this region is therefore crucial to decipher past population movements in South America, its characterization has been overlooked so far. We report 92 modern and 22 ancient mitogenomes spanning a temporal range of 5000 years, which were compared with a large set of previously reported data. Leveraging this dataset representative of the mitochondrial diversity of the subcontinent, we investigate the maternal history of CA populations within a wider geographical context. We describe a large number of novel clades within the mitochondrial DNA tree, thus providing new phylogenetic interpretations for South America. We also identify several local clades of great temporal depth with continuity until the present time, which stem directly from the founder haplotypes, suggesting that they originated in the region and expanded from there. Moreover, the presence of lineages characteristic of other South American regions reveals the existence of gene flow to CA. Finally, we report some lineages with discontinuous distribution across the Americas, which suggest the persistence of relic lineages likely linked to the first population arrivals. The present study represents to date the most exhaustive attempt to elaborate a Native American genetic map from modern and ancient complete mitochondrial genomes in Argentina and provides relevant information about the general process of settlement in South America.

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.

Population history of North Africa based on modern and ancient genomes

Compared with the rest of the African continent, North Africa has provided limited genomic data. Nonetheless, the genetic data available show a complex demographic scenario characterized by extensive admixture and drift. Despite the continuous gene flow from the Middle East, Europe and sub-Saharan Africa, an autochthonous genetic component that dates back to pre-Holocene times is still present in North African groups. The comparison of ancient and modern genomes has evidenced a genetic continuity in the region since Epipaleolithic times. Later population movements, especially the gene flow from the Middle East associated with the Neolithic, have diluted the genetic autochthonous component, creating an east to west gradient. Recent historical movements, such as the Arabization, have also contributed to the genetic landscape observed currently in North Africa and have culturally transformed the region. Genome analyses have not shown evidence of a clear correlation between cultural and genetic diversity in North Africa, as there is no genetic pattern of differentiation between Tamazight (i.e. Berber) and Arab speakers as a whole. Besides the gene flow received from neighboring areas, the analysis of North African genomes has shown that the region has also acted as a source of gene flow since ancient times. As a result of the genetic uniqueness of North African groups and the lack of available data, there is an urgent need for the study of genetic variation in the region and its implications in health and disease.

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

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

Sex-biased patterns shaped the genetic history of Roma

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

Dissecting human North African gene-flow into its western coastal surroundings

North African history and populations have exerted a pivotal influence on surrounding geographical regions, although scant genetic studies have addressed this issue. Our aim is to understand human historical migrations in the coastal surroundings of North Africa. We built a refined genome-wide dataset of North African populations to unearth the fine-scale genetic structure of the region, using haplotype information. The results suggest that the gene-flow from North Africa into the European Mediterranean coast (Tuscany and the Iberian Peninsula) arrived mainly from the Mediterranean coast of North Africa. In Tuscany, this North African admixture date estimate suggests the movement of peoples during the fall of the Roman Empire around the fourth century. In the Iberian Peninsula, the North African component probably reflects the impact of the Arab expansion since the seventh century and the subsequent expansion of the Christian Kingdoms. By contrast, the North African component in the Canary Islands has a source genetically related to present-day people from the Atlantic North African coast. We also find sub-Saharan gene-flow from the Senegambia region in the Canary Islands. Specifically, we detect a complex signal of admixture involving Atlantic, Senegambian and European sources intermixing around the fifteenth century, soon after the Castilian conquest. Our results highlight the differential genetic influence of North Africa into the surrounding coast and show that specific historical events have not only had a socio-cultural impact but additionally modified the gene pool of the populations.

Reconstructed Lost Native American Populations from Eastern Brazil Are Shaped by Differential Jê/Tupi Ancestry

After the colonization of the Americas by Europeans and the consequent Trans-Atlantic Slave Trade, most Native American populations in eastern Brazil disappeared or went through an admixture process that configured a population composed of three main genetic components: the European, the sub-Saharan African, and the Native American. The study of the Native American genetic history is challenged by the lack of availability of genome-wide samples from Native American populations, the technical difficulties to develop ancient DNA studies, and the low proportions of the Native American component in the admixed Brazilian populations (on average 7%). We analyzed genome-wide data of 5,825 individuals from three locations of eastern Brazil: Salvador (North-East), Bambui (South-East), and Pelotas (South) and we reconstructed populations that emulate the Native American groups that were living in the 16th century around the sampling locations. This genetic reconstruction was performed after local ancestry analysis of the admixed Brazilian populations, through the rearrangement of the Native American haplotypes into reconstructed individuals with full Native American ancestry (51 reconstructed individuals in Salvador, 45 in Bambui, and 197 in Pelotas). We compared the reconstructed populations with nonadmixed Native American populations from other regions of Brazil through haplotype-based methods. Our results reveal a population structure shaped by the dichotomy of Tupi-/Jê-speaking ancestry related groups. We also show evidence of a decrease of the diversity of nonadmixed Native American groups after the European contact, in contrast with the reconstructed populations, suggesting a reservoir of the Native American genetic diversity within the admixed Brazilian population.

Distribution of local ancestry and evidence of adaptation in admixed populations

Admixed American populations have different global proportions of European, Sub-Saharan African, and Native-American ancestry. However, individuals who display the same global ancestry could exhibit remarkable differences in the distribution of local ancestry blocks. We studied for the first time the distribution of local ancestry across the genome of 264 Brazilian admixed individuals, ascertained within the scope of the Brazilian Initiative on Precision Medicine. We found a decreased proportion of European ancestry together with an excess of Native-American ancestry on chromosome 8p23.1 and showed that this is due to haplotypes created by chromosomal inversion events. Furthermore, Brazilian non-inverted haplotypes were more similar to Native-American haplotypes than to European haplotypes, in contrast to what was found in other American admixed populations. We also identified signals of recent positive selection on chromosome 8p23.1, and one gene within this locus, PPP1R3B, is related to glycogenesis and has been associated with an increased risk of type 2 diabetes and obesity. These findings point to a selection event after admixture, which is still not entirely understood in recent admixture events.

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.

Patterns of genetic structure and adaptive positive selection in the Lithuanian population from high-density SNP data

The analysis of geographically specific regions and the characterization of fine-scale patterns of genetic diversity may facilitate a much better understanding of the microevolutionary processes affecting local human populations. Here we generated genome-wide high-density SNP genotype data in 425 individuals from six geographical regions in Lithuania and combined our dataset with available ancient and modern data to explore genetic population structure, ancestry components and signatures of natural positive selection in the Lithuanian population. Our results show that Lithuanians are a homogenous population, genetically differentiated from neighbouring populations but within the general expected European context. Moreover, we not only confirm that Lithuanians preserve one of the highest proportions of western, Scandinavian and eastern hunter-gather ancestry components found in European populations but also that of an steppe Early to Middle Bronze Age pastoralists, which together configure the genetic distinctiveness of the Lithuanian population. Finally, among the top signatures of positive selection detected in Lithuanians, we identified several candidate genes related with diet (PNLIP, PPARD), pigmentation (SLC24A5, TYRP1 and PPARD) and the immune response (BRD2, HLA-DOA, IL26 and IL22).

Spatially explicit analysis reveals complex human genetic gradients in the Iberian Peninsula

The Iberian Peninsula is a well-delimited geographic region with a rich and complex human history. However, the causes of its genetic structure and past migratory dynamics are not yet fully understood. In order to shed light on them, here we evaluated the gene flow and genetic structure throughout the Iberian Peninsula with spatially explicit modelling applied to a georeferenced genetic dataset composed of genome-wide SNPs from 746 individuals belonging to 17 different regions of the Peninsula. We found contrasting patterns of genetic structure throughout Iberia. In particular, we identified strong patterns of genetic differentiation caused by relevant barriers to gene flow in northern regions and, on the other hand, a large genetic similarity in central and southern regions. In addition, our results showed a preferential north to south migratory dynamics and suggest a sex-biased dispersal in Mediterranean and southern regions. The estimated genetic patterns did not fit with the geographical relief of the Iberian landscape and they rather seem to follow political and linguistic territorial boundaries.

Whole-genome sequence analysis of a Pan African set of samples reveals archaic gene flow from an extinct basal population of modern humans into sub-Saharan populations

Background

Population demography and gene flow among African groups, as well as the putative archaic introgression of ancient hominins, have been poorly explored at the genome level.

Results

Here, we examine 15 African populations covering all major continental linguistic groups, ecosystems, and lifestyles within Africa through analysis of whole-genome sequence data of 21 individuals sequenced at deep coverage. We observe a remarkable correlation among genetic diversity and geographic distance, with the hunter-gatherer groups being more genetically differentiated and having larger effective population sizes throughout most modern-human history. Admixture signals are found between neighbor populations from both hunter-gatherer and agriculturalists groups, whereas North African individuals are closely related to Eurasian populations. Regarding archaic gene flow, we test six complex demographic models that consider recent admixture as well as archaic introgression. We identify the fingerprint of an archaic introgression event in the sub-Saharan populations included in the models (~ 4.0% in Khoisan, ~ 4.3% in Mbuti Pygmies, and ~ 5.8% in Mandenka) from an early divergent and currently extinct ghost modern human lineage.

Conclusion

The present study represents an in-depth genomic analysis of a Pan African set of individuals, which emphasizes their complex relationships and demographic history at population level.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1684-5) contains supplementary material, which is available to authorized users.

People from Ibiza: an unexpected isolate in the Western Mediterranean

In this study, we seek to understand and to correlate the genetic patterns observed in the population of the island of Ibiza in the Western Mediterranean basin with past events. Genome-wide genotypes of 189 samples representing 13 of 17 regions in Spain have been analyzed, in addition to 105 samples from the Levant, 157 samples from North Africa, and one ancient sample from the Phoenician Cas Molí site in Ibiza. Before the Catalans conquered the island in 1235 CE, Ibiza (Eivissa) had already been influenced by several cultures, starting with the Phoenicians, then the Carthaginians, followed by the Umayyads. The impact of these various cultures on the genetic structure of the islanders is still unexplored. Our results show a clear distinction between Ibiza and the rest of Spain. To investigate whether this was due to the Phoenician colonization or to more recent events, we compared the genomes of current Ibizans to that of an ancient Phoenician sample from Ibiza and to both modern Levantine and North African genomes. We did not identify any trace of Phoenician ancestry in the current Ibizans. Interestingly, the analysis of runs of homozygosity and changes in the effective population size through time support the idea that drift has shaped the genetic structure of current Ibizans. In addition to the small carrying capacity of the island, Ibiza experienced a series of dramatic demographic changes due to several instances of famine, war, malaria and plague that could have significantly contributed to its current genetic differentiation.

Sequence diversity of the Rh blood group system in Basques

Basques show specific cultural, demographic, and genetic characteristics that have placed them as an isolated and unique population within Europe, such as their non-Indo-European language, Euskara. They have historically lived along the Western Pyrenees, between Spain and France, in one of the most important European glacial refugia during the Last Glacial Maximum. The most striking genetic characteristic is their highest frequency of the RhD blood group negative allele, a variant related to the hemolytic disease of the newborn. Both demographic and adaptive processes have been suggested as possible causes of the high frequency of RhD negative in Basques, but neither hypothesis has been clearly demonstrated. While previous studies on the Rh system in Basques have been mostly focused on serological and genotyping diversity, in this work we analyze genotyping and next generation sequencing data in order to provide a general framework of the genetic scenario of the system in Basques. In particular, we genotyped the most relevant variants of the system (D/d, E/e, and C/c), and sequenced three ~6 kb flanking regions surrounding the Rh genes in Basques and also in other populations for comparison. Our results are in agreement with previous studies, with Basques presenting the highest frequency of the RHD deletion (47.2%). Haplotype analyses of D/d, E/e, and C/c variants confirmed an association between the RhC allele, previously suggested to be under positive selection, and the RhD positive variant in non-sub-Saharan populations, including Basques. We also found extreme differentiation for the C/c variant when comparing sub-Saharan to non-sub-Saharan populations.

A Roma founder BIN1 mutation causes a novel phenotype of centronuclear myopathy with rigid spine

OBJECTIVE

To describe a large series of BIN1 patients, in which a novel founder mutation in the Roma population of southern Spain has been identified.

METHODS

Patients diagnosed with centronuclear myopathy (CNM) at 5 major reference centers for neuromuscular disease in Spain (n = 53) were screened for BIN1 mutations. Clinical, histologic, radiologic, and genetic features were analyzed.

RESULTS

Eighteen patients from 13 families carried the p.Arg234Cys variant; 16 of them were homozygous for it and 2 had compound heterozygous p.Arg234Cys/p.Arg145Cys mutations. Both BIN1 variants have only been identified in Roma, causing 100% of CNM in this ethnic group in our cohort. The haplotype analysis confirmed all families are related. In addition to clinical features typical of CNM, such as proximal limb weakness and ophthalmoplegia, most patients in our cohort presented with prominent axial weakness, often associated with rigid spine. Severe fat replacement of paravertebral muscles was demonstrated by muscle imaging. This phenotype seems to be specific to the p.Arg234Cys mutation, not reported in other BIN1 mutations. Extreme clinical variability was observed in the 2 compound heterozygous patients for the p.Arg234Cys/p.Arg145Cys mutations, from a congenital onset with catastrophic outcome to a late-onset disease. Screening of European Roma controls (n = 758) for the p.Arg234Cys variant identified a carrier frequency of 3.5% among the Spanish Roma.

CONCLUSION

We have identified a BIN1 founder Roma mutation associated with a highly specific phenotype, which is, from the present cohort, the main cause of CNM in Spain.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Recent Historical Migrations Have Shaped the Gene Pool of Arabs and Berbers in North Africa

North Africa is characterized by its diverse cultural and linguistic groups and its genetic heterogeneity. Genomic data has shown an amalgam of components mixed since pre-Holocean times. Though no differences have been found in uniparental and classical markers between Berbers and Arabs, the two main ethnic groups in the region, the scanty genomic data available have highlighted the singularity of Berbers. We characterize the genetic heterogeneity of North African groups, focusing on the putative differences of Berbers and Arabs, and estimate migration dates. We analyze genome-wide autosomal data in five Berber and six Arab groups, and compare them to Middle Easterns, sub-Saharans, and Europeans. Haplotype-based methods show a lack of correlation between geographical and genetic populations, and a high degree of genetic heterogeneity, without strong differences between Berbers and Arabs. Berbers enclose genetically diverse groups, from isolated endogamous groups with high autochthonous component frequencies, large homozygosity runs and low effective population sizes, to admixed groups with high frequencies of sub-Saharan and Middle Eastern components. Admixture time estimates show a complex pattern of recent historical migrations, with a peak around the 7th century C.E. coincident with the Arabization of the region; sub-Saharan migrations since the 1st century B.C. in agreement with Roman slave trade; and a strong migration in the 17th century C.E., coincident with a huge impact of the trans-Atlantic and trans-Saharan trade of sub-Saharan slaves in the Modern Era. The genetic complexity found should be taken into account when selecting reference groups in population genetics and biomedical studies.

A tale of agriculturalists and hunter-gatherers: Exploring the thrifty genotype hypothesis in native South Americans

OBJECTIVES

To determine genetic differences between agriculturalist and hunter-gatherer southern Native American populations for selected metabolism-related markers and to test whether Neel's thrifty genotype hypothesis (TGH) could explain the genetic patterns observed in these populations.

MATERIALS AND METHODS

375 Native South American individuals from 17 populations were genotyped using six markers (APOE rs429358 and rs7412; APOA2 rs5082; CD36 rs3211883; TCF7L2 rs11196205; and IGF2BP2 rs11705701). Additionally, APOE genotypes from 39 individuals were obtained from the literature. AMOVA, main effects, and gene-gene interaction tests were performed.

RESULTS

We observed differences in allele distribution patterns between agriculturalists and hunter-gatherers for some markers. For instance, between-groups component of genetic variance (F_CT ) for APOE rs429358 showed strong differences in allelic distributions between hunter-gatherers and agriculturalists (p = 0.00196). Gene-gene interaction analysis indicated that the APOE E4/CD36 TT and APOE E4/IGF2BP2 A carrier combinations occur at a higher frequency in hunter-gatherers, but this combination is not replicated in archaic (Neanderthal and Denisovan) and ancient (Anzick, Saqqaq, Ust-Ishim, Mal'ta) hunter-gatherer individuals.

DISCUSSION

A complex scenario explains the observed frequencies of the tested markers in hunter-gatherers. Different factors, such as pleotropic alleles, rainforest selective pressures, and population dynamics, may be collectively shaping the observed genetic patterns. We conclude that although TGH seems a plausible hypothesis to explain part of the data, other factors may be important in our tested populations.

Complement factor 5 (C5) p.A252T mutation is prevalent in, but not restricted to, sub-Saharan Africa: implications for the susceptibility to meningococcal disease

Complement C5 deficiency (C5D) is a rare primary immunodeficiency associated with recurrent infections, particularly meningitis, by Neisseria species. To date, studies to elucidate the molecular basis of hereditary C5D have included fewer than 40 families, and most C5 mutations (13 of 17) have been found in single families. However, the recently described C5 p.A252T mutation is reported to be associated with approximately 7% of meningococcal disease cases in South Africa. This finding raises the question of whether the mutation may be prevalent in other parts of Africa or other continental regions. The aim of this study was to investigate the prevalence of C5 p.A252T in Africa and other regions and discuss the implications for prophylaxis against meningococcal disease. In total, 2710 samples from healthy donors within various populations worldwide were analysed by quantitative polymerase chain reaction (qPCR) assay to detect the C5 p.A252T mutation. Eleven samples were found to be heterozygous for p.A252T, and nine of these samples were from sub-Saharan African populations (allele frequency 0·94%). Interestingly, two other heterozygous samples were from individuals in populations outside Africa (Israel and Pakistan). These findings, together with data from genomic variation databases, indicate a 0·5-2% prevalence of the C5 p.A252T mutation in heterozygosity in sub-Saharan Africa. Therefore, this mutation may have a relevant role in meningococcal disease susceptibility in this geographical area.

The Simons Genome Diversity Project: 300 genomes from 142 diverse populations

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans.

Genomic insights into the origin of farming in the ancient Near East

We report genome-wide ancient DNA from 44 ancient Near Easterners ranging in time between ~12,000 and 1,400 bc, from Natufian hunter-gatherers to Bronze Age farmers. We show that the earliest populations of the Near East derived around half their ancestry from a 'Basal Eurasian' lineage that had little if any Neanderthal admixture and that separated from other non-African lineages before their separation from each other. The first farmers of the southern Levant (Israel and Jordan) and Zagros Mountains (Iran) were strongly genetically differentiated, and each descended from local hunter-gatherers. By the time of the Bronze Age, these two populations and Anatolian-related farmers had mixed with each other and with the hunter-gatherers of Europe to greatly reduce genetic differentiation. The impact of the Near Eastern farmers extended beyond the Near East: farmers related to those of Anatolia spread westward into Europe; farmers related to those of the Levant spread southward into East Africa; farmers related to those of Iran spread northward into the Eurasian steppe; and people related to both the early farmers of Iran and to the pastoralists of the Eurasian steppe spread eastward into South Asia.