New clues to the evolutionary history of the main European paternal lineage M269: dissection of the Y-SNP S116 in Atlantic Europe and Iberia

de Pancorbo M, Calafell F. Y-chromosome target enrichment reveals rapid expansion of haplogroup R1b-DF27 in Iberia during the Bronze Age transition

The Y chromosome can yield a unique perspective into the study of human demographic history. However, due to the repetitive nature of part of its sequence, only a small set of regions are suitable for variant calling and discovery from short-read sequencing data. These regions combined represent 8.9 Mbp or 0.14% of a diploid human genome. Consequently, investing in whole-genome sequencing to resolve Y-chromosome questions is poorly efficient. Here we use, as an alternative, target enrichment technology to greatly increase sequencing effectiveness, validating and applying the technique to 181 males, for 162 of whom we obtained a positive result. Additionally, 75 samples sequenced for the whole genome were also included, for a total sample size of 237. These samples were chosen for their Y chromosome haplogroup: R1b-DF27. In the context of European populations, and particularly in Iberia, this haplogroup stands out for its high frequency and its demographic history. Current evidence indicates that the diffusion of this haplogroup is related to the population movements that mark the cultural Bronze Age transition, making it remarkably interesting for population geneticists. The results of this study show the effects of the rapid radiation of the haplogroup in Spain, as even with the higher discriminating power of whole sequences, most haplotypes still fall within the R1b-DF27* paragroup rather than in the main derived branches. However, we were able to refine the ISOGG 2019-2020 phylogeny, and its two main subbranches, namely L176.2 and Z272, which present geographical differentiation between the Atlantic and Mediterranean coasts of Iberia.

Post-Austronesian migrational wave of West Polynesians to Micronesia

This study examines Y-chromosome and mtDNA markers in the population of the island of Kiritimati in the context of geographically targeted reference populations from the Pacific. Kiritimati derives its population from the atoll islands of the Gilbert Archipelago and representsa geographicaltransitional region between Micronesia, Polynesia and Melanesia that likely played a critical role during theAustronesian expansion. The large presence(84.1%)of individuals withO-M175, O2a-M324 and O2a2b-P164 sub-haplogroups, 69.9% being O2a2b-P164, the Y-STR homogeneity within O2a2b-P164 and the very recent age of the sub-haplogroup(363-548 years ago)inKiritimati suggestthe arrival ofa genetically homogenous population to the Gilberteses followed by a population expassion.The close Y-STR haplotype affinities with profiles from the Samoa and Tonga Archipelagos point to an unprecedented massive post-Austronesian expansionexodus from West Polynesia.Contrasting the abundance of AustronesianO2a2b-P164 sub-haplogroup, the most abundantMelanesian/Papuansub-haplogroup,C-M130is present at a frequency of 13.5%. Thenetwork topology suggests that C-M130 arrived to theKiribati Archipelago from West Polynesia, specifically from West Samoa, Tonga and/or Tutuila subsequent to the Austronesian expansion about 832-1408 years ago. The haplotype affinities withinO2a2b-P164 argue for anoriginal source in Taiwan and its dispersal to West Polynesia and then to Southeast Micronesia. The present investigation provides an understanding of the genetic composition and complex migration history of an understudied region of the Pacific and provides evidence for recent dispersals towards Micronesia from West Polynesia subsequent to the initial Austronesian expansion.

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas

The Isthmus of Panama was a crossroads between North and South America during the continent’s first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

Phylogeographic review of Y chromosome haplogroups in Europe

The Y chromosome has been widely explored for the study of human migrations. Due to its paternal inheritance, the Y chromosome polymorphisms are helpful tools for understanding the geographical distribution of populations all over the world and for inferring their origin, which is really useful in forensics. The remarkable historical context of Europe, with numerous migrations and invasions, has turned this continent into a melting pot. For this reason, it is interesting to study the Y chromosome variability and how it has contributed to improving our knowledge of the distribution and development of European male genetic pool as it is today. The analysis of Y lineages in Europe shows the predominance of four haplogroups, R1b-M269, I1-M253, I2-M438 and R1a-M420. However, other haplogroups have been identified which, although less frequent, provide significant evidence about the paternal origin of the populations. In addition, the study of the Y chromosome in Europe is a valuable tool for revealing the genetic trace of the different European colonizations, mainly in several American countries, where the European ancestry is mostly detected by the presence of the R1b-M269 haplogroup. Therefore, the objective of this review is to compile the studies of the Y chromosome haplogroups in current European populations, in order to provide an outline of these haplogroups which facilitate their use in forensic studies.

The Y chromosome of autochthonous Basque populations and the Bronze Age replacement

Here we report on the Y haplogroup and Y-STR diversity of the three autochthonous Basque populations of Alava (n = 54), Guipuzcoa (n = 30) and Vizcaya (n = 61). The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STR loci (DYS19, DYS385a/b, DYS398I/II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, Y-GATA H4) using the AmpFlSTR Yfiler system. Six major haplogroups (R, I, E, J, G, and DE) were detected, being R-S116 (P312) haplogroup the most abundant at 75.0% in Alava, 86.7% in Guipuzcoa and 87.3% in Vizcaya. Age estimates for the R-S116 mutation in the Basque Country are 3975 ± 303, 3680 ± 345 and 4553 ± 285 years for Alava, Guipuzcoa and Vizcaya, respectively. Pairwise Rst genetic distances demonstrated close Y-chromosome affinities among the three autochthonous Basque populations and between them and the male population of Ireland and Gascony. In a MDS plot, the population of Ireland segregates within the Basque cluster and closest to the population of Guipuzcoa, which plots closer to Ireland than to any of the other Basque populations. Overall, the results support the notion that during the Bronze Age a dispersal of individuals carrying the R-S116 mutation reached the Basque Country replacing the Paleolithic/Neolithic Y chromosome of the region.

A Y-chromosomal survey of Ecuador's multi-ethnic population reveals new insights into the tri-partite population structure and supports an early Holocene age of the rare Native American founder lineage C3-MPB373

Ecuador is a multiethnic and pluricultural country with a complex history defined by migration and admixture processes. The present study aims to increase our knowledge on the Ecuadorian Native Amerindian groups and the unique South American Y-chromosome haplogroup C3-MPB373 through the analysis of up to 23 Y-chromosome STRs (Y-STRs) and several Y-SNPs in a sample of 527 Ecuadorians from 7 distinct populations and geographic areas, including Kichwa and non-Kichwa Native Amerindians, Mestizos and Afro-Ecuadorians. Our results reveal the presence of C3-MPB373 both in the Amazonian lowland Kichwa with frequencies up to 28 % and, for the first time, in notable proportions in Kichwa populations from the Ecuadorian highlands. The substantially higher frequencies of C3-MPB373 in the Amazonian lowlands found in Kichwa and Waorani individuals suggest a founder effect in that area. Notably, estimates for the time to the most recent common ancestor (TMRCA) in the range of 7.2-9.0 kya point to an ancient origin of the haplogroup and suggest an early Holocene expansion of C3-MPB373 into South America. Finally, the pairwise genetic distances (R_ST) separate the Kichwa Salasaka from all the other Native Amerindian and Ecuadorian groups, indicating a so far hidden diversity among the Kichwa-speaking populations and suggesting a more southern origin of this population. In sum, our study provides a more in-depth knowledge of the male genetic structure of the multiethnic Ecuadorian population, as well as a valuable reference dataset for forensic use.

The Y-chromosomal haplotype and haplogroup distribution of modern Switzerland still reflects the alpine divide as a geographical barrier for human migration

A sample of 606 Swiss individuals has been characterized for 27 Y-STR and 34 Y-SNPs, defining major European haplogroups. For the first time, a subsample from the southernmost part of Switzerland, the Italian speaking canton Ticino, has been included. The data reveals significant intra-national differences in the distribution of haplogroups R1b-U106, R1b-U152, I1 and J2a north and south of the alpine divide, with R1b-U152 being the most frequent haplogroup among all Swiss subpopulations, reaching 26 % in average and 53 % in the Ticino sample. In addition, a high percentage of haplogroup E1b1b-M35 in Eastern Switzerland corresponds well with data reported from Western Austria. In general, we detected a low level of differentiation between the subgroups north of the alpine divide. The dataset also revealed a variety of microvariants. Some of them were previously known to be associated with particular haplogroups. However, we discovered one microvariant in DYS533 that seems to be closely associated with haplogroup I2-P215 (xM223). This association had not yet been reported to date. The concordance study with two STR-kits suggests that the DYS533 microvariant is due to an InDel in the flanking regions of the marker. One individual carried a large deletion, frequently detected in people of East Asian ancestry, encompassing the amelogenin locus. To our knowledge, this is the first time that such a deletion has been observed within European haplogroup R1b-U152. This is the first comprehensive Y chromosomal dataset for Switzerland, demonstrating significant population substructure due to an intra-national geographical barrier.

A Statistical Method to Enhance the Analysis of the Differences Among High-Resolution Melting (HRM) Curves of PCR-Amplified DNA Fragments

Consistent differences among melting curves of PCR-amplified DNA fragments are treated by normalizing the relative fluorescence units (RFU) and performing a clustering analysis, but statistically significant differences among curves are not usually determined. In the present study, an analysis based on functional data analysis (FDA) was implemented to evaluate the existence of statistically significant differences between normalized RFU curves obtained from PCR-HRM (high-resolution melting) analysis by using ANOVA for functional data. The effectiveness of the FDA method was analyzed with data from a set of samples of eight animal species of interest in food analysis, as well as mixtures of DNA from these species, analyzed by PCR-HRM to differentiate them. The statistical method described in this study has been demonstrated to be a robust and precise tool to discriminate among melting curves derived from HRM analysis. This method has advantages over the current comparison methods. PRACTICAL APPLICATION: As long as food fraud and mislabeling exist, new techniques for species identification are needed. High-resolution melting (HRM) has been shown to be a rapid, reliable and inexpensive species identification method. In the present study, functional data analysis (FDA) was applied to HRM curves of DNA from eight animal species used for food, as well as to mixtures of these species in different proportions. FDA has advantages over the usual methods, providing a deeper statistical analysis and facilitating the data interpretation as shown by the HRM analysis for a clearer comparison among individual species and mixtures of species.

The genomic history of the Iberian Peninsula over the past 8000 years

We assembled genome-wide data from 271 ancient Iberians, of whom 176 are from the largely unsampled period after 2000 BCE, thereby providing a high-resolution time transect of the Iberian Peninsula. We document high genetic substructure between northwestern and southeastern hunter-gatherers before the spread of farming. We reveal sporadic contacts between Iberia and North Africa by ~2500 BCE and, by ~2000 BCE, the replacement of 40% of Iberia's ancestry and nearly 100% of its Y-chromosomes by people with Steppe ancestry. We show that, in the Iron Age, Steppe ancestry had spread not only into Indo-European-speaking regions but also into non-Indo-European-speaking ones, and we reveal that present-day Basques are best described as a typical Iron Age population without the admixture events that later affected the rest of Iberia. Additionally, we document how, beginning at least in the Roman period, the ancestry of the peninsula was transformed by gene flow from North Africa and the eastern Mediterranean.

Paternal lineages in southern Iberia provide time frames for gene flow from mainland Europe and the Mediterranean world

BACKGROUND

The geography of southern Iberia and an abundant archaeological record of human occupation are ideal conditions for a full understanding of scenarios of genetic history in the area. Recent advances in the phylogeography of Y-chromosome lineages offer the opportunity to set upper bounds for the appearance of different genetic components.

AIM

To provide a global knowledge on the Y haplogroups observed in Andalusia with their Y microsatellite variation. Preferential attention is given to the vehement debate about the age, origin and expansion of R1b-M269 clade and sub-lineages.

SUBJECT AND METHODS

Four hundred and fourteen male DNA samples from western and eastern autochthonous Andalusians were genotyped for a set of Y-SNPs and Y-STRs. Gene diversity, potential population genetic structures and coalescent times were assessed.

RESULTS

Most of the analysed samples belong to the European haplogroup R1b1a1a2-M269, whereas haplogroups E, J, I, G and T show lower frequencies. A phylogenetic dissection of the R1b-M269 was performed and younger time frames than those previously reported in the literature were obtained for its sub-lineages.

CONCLUSION

The particular Andalusian R1b-M269 assemblage confirms the shallow topology of the clade. Moreover, the sharing of lineages with the rest of Europe indicates the impact in Iberia of an amount of pre-existing diversity, with the possible exception of R1b-DF27. Lineages such as J2-M172 and G-M201 highlight the importance of maritime travels of early farmers who reached the Iberian Peninsula.

The Marquesans at the fringes of the Austronesian expansion

In the present study, 87 unrelated individuals from the Marquesas Archipelago in French Polynesia were typed using mtDNA, Y-chromosome and autosomal (STRs) markers and compared to key target populations from Island South East Asia (ISEA), Taiwan, and West and East Polynesia to investigate their genetic relationships. The Marquesas, located at the eastern-most fringes of the Austronesian expansion, offer a unique opportunity to examine the effects of a protracted population expansion wave on population structure. We explore the contribution of Melanesian, Asian and European heritage to the Marquesan islands of Nuku-Hiva, Hiva-Oa and Tahuata. Overall, the Marquesas Islands are genetically homogeneous. In the Marquesan Archipelago all of the mtDNA haplogroups are of Austronesian origin belonging to the B4a1 subhaplogroup as the region marks the end of a west to east decreasing cline of Melanesian mtDNA starting with the West Polynesian population of Tonga. Genetic discrepancies are less pronounced between the Marquesan and Society islands, and among the Marquesan islands. Interestingly, a number of Melanesian, Polynesian and European Y-chromosome haplogroups exhibit very different distribution between the Marquesan islands of Nuku Hiva and Hiva Oa, likely resulting from drift, differential migration involving various source populations and/or unique trading routes.

Paternal portrait of populations of the middle Magdalena River region (Tolima and Huila, Colombia): New insights on the peopling of Central America and northernmost South America

The valley of the Magdalena River is one of the main population pathways in Colombia. The gene pool and spatial configuration of human groups in this territory have been outlined throughout three historical stages: the Native pre-Hispanic world, Spanish colonization, and XIX century migrations. This research was designed with the goal of characterizing the diversity and distribution pattern of Y-chromosome lineages that are currently present in the Tolima and Huila departments (middle Magdalena River region). Historic cartography was used to identify the main geographic sites where the paternal lineages belonging to this area have gathered. Twelve municipalities were chosen, and a survey that included genealogical information was administered. Samples collected from 83 male volunteers were analyzed for 48 Y-SNPs and 17 Y-STRs. The results showed a highly diverse region characterized by the presence of 16 sublineages within the major clades R, Q, J, G, T and E and revealed that 93% (n = 77) of haplotypes were different. Among these haplogroups, European-specific R1b-M269 lineages were the most representative (57.83%), with six different subhaplogroups and 43 unique haplotypes. Native American paternal ancestry was also detected based on the presence of the Q1a2-M3*(xM19, M194, M199) and Q1a2-M346*(xM3) lineages. Interestingly, all Q1a2-M346*(xM3) samples (n = 7, with five different haplotypes) carried allele six at the DYS391 locus. This allele has a worldwide frequency of 0.169% and was recently associated with a new Native subhaplogroup. An in-depth phylogenetic analysis of these samples suggests the Tolima and Huila region to be the principal area in all Central and South America where this particular Native lineage is found. This lineage has been present in the region for at least 1,809 (+/- 0,5345) years.

Effective resolution of the Y chromosome sublineages of the Iberian haplogroup R1b-DF27 with forensic purposes

Single-nucleotide polymorphisms (SNPs) found within the non-recombining region of the Y chromosome (NRY) represent a powerful tool in forensic genetics for inferring the paternal ancestry of a vestige and complement the determination of biogeographical origin in combination with other markers like AIMs. In the present study, we introduce a panel of 15 Y-SNPs for a fine-resolution subtyping of the haplogroup R1b-DF27, in a single minisequencing reaction. This is the first minisequencing panel that allows a fine subtyping of R1b-DF27, which displays high frequencies in Iberian and Iberian-influenced populations. This panel includes subhaplogroups of DF27 that display moderate geographical differentiation, of interest to link a sample with a specific location of the Iberian Peninsula or with Iberian ancestry. Conversely, part of the intricacy of a new minisequencing panel is to have all the included variants available to test the effectiveness of the analysis method. We have overcome the absence of the least common variants through site-directed mutagenesis. Overall, the results show that our panel is a robust and effective method for subtyping R1b-DF27 lineages from a minimal amount of DNA, and its high resolution enables to improve male lineage discrimination in Iberian and Southwest European descent individuals. The small length of the amplicons and its reproducibility makes this assay suitable for forensic and population genetics purposes.

Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective

BACKGROUND

Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then.

AIM

This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation.

SUBJECTS AND METHODS

High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals.

RESULTS

Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones.

CONCLUSIONS

A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

The Beaker phenomenon and the genomic transformation of northwest Europe

From around 2750 to 2500 bc, Bell Beaker pottery became widespread across western and central Europe, before it disappeared between 2200 and 1800 bc. The forces that propelled its expansion are a matter of long-standing debate, and there is support for both cultural diffusion and migration having a role in this process. Here we present genome-wide data from 400 Neolithic, Copper Age and Bronze Age Europeans, including 226 individuals associated with Beaker-complex artefacts. We detected limited genetic affinity between Beaker-complex-associated individuals from Iberia and central Europe, and thus exclude migration as an important mechanism of spread between these two regions. However, migration had a key role in the further dissemination of the Beaker complex. We document this phenomenon most clearly in Britain, where the spread of the Beaker complex introduced high levels of steppe-related ancestry and was associated with the replacement of approximately 90% of Britain's gene pool within a few hundred years, continuing the east-to-west expansion that had brought steppe-related ancestry into central and northern Europe over the previous centuries.

Analysis of the R1b-DF27 haplogroup shows that a large fraction of Iberian Y-chromosome lineages originated recently in situ

Haplogroup R1b-M269 comprises most Western European Y chromosomes; of its main branches, R1b-DF27 is by far the least known, and it appears to be highly prevalent only in Iberia. We have genotyped 1072 R1b-DF27 chromosomes for six additional SNPs and 17 Y-STRs in population samples from Spain, Portugal and France in order to further characterize this lineage and, in particular, to ascertain the time and place where it originated, as well as its subsequent dynamics. We found that R1b-DF27 is present in frequencies ~40% in Iberian populations and up to 70% in Basques, but it drops quickly to 6–20% in France. Overall, the age of R1b-DF27 is estimated at ~4,200 years ago, at the transition between the Neolithic and the Bronze Age, when the Y chromosome landscape of W Europe was thoroughly remodeled. In spite of its high frequency in Basques, Y-STR internal diversity of R1b-DF27 is lower there, and results in more recent age estimates; NE Iberia is the most likely place of origin of DF27. Subhaplogroup frequencies within R1b-DF27 are geographically structured, and show domains that are reminiscent of the pre-Roman Celtic/Iberian division, or of the medieval Christian kingdoms.

Male Lineages in Brazil: Intercontinental Admixture and Stratification of the European Background

The non-recombining nature of the Y chromosome and the well-established phylogeny of Y-specific Single Nucleotide Polymorphisms (Y-SNPs) make them useful for defining haplogroups with high geographical specificity; therefore, they are more apt than the Y-STRs to detect population stratification in admixed populations from diverse continental origins. Different Y-SNP typing strategies have been described to address issues of population history and movements within geographic territories of interest. In this study, we investigated a set of 41 Y-SNPs in 1217 unrelated males from the five Brazilian geopolitical regions, aiming to disclose the genetic structure of male lineages in the country. A population comparison based on pairwise FST genetic distances did not reveal statistically significant differences in haplogroup frequency distributions among populations from the different regions. The genetic differences observed among regions were, however, consistent with the colonization history of the country. The sample from the Northern region presented the highest Native American ancestry (8.4%), whereas the more pronounced African contribution could be observed in the Northeastern population (15.1%). The Central-Western and Southern samples showed the higher European contributions (95.7% and 93.6%, respectively). The Southeastern region presented significant European (86.1%) and African (12.0%) contributions. The subtyping of the most frequent European lineage in Brazil (R1b1a-M269) allowed differences in the genetic European background of the five Brazilian regions to be investigated for the first time.

The Paternal Landscape along the Bight of Benin - Testing Regional Representativeness of West-African Population Samples Using Y-Chromosomal Markers

Patterns of genetic variation in human populations across the African continent are still not well studied in comparison with Eurasia and America, despite the high genetic and cultural diversity among African populations. In population and forensic genetic studies a single sample is often used to represent a complete African region. In such a scenario, inappropriate sampling strategies and/or the use of local, isolated populations may bias interpretations and pose questions of representativeness at a macrogeographic-scale. The non-recombining region of the Y-chromosome (NRY) has great potential to reveal the regional representation of a sample due to its powerful phylogeographic information content. An area poorly characterized for Y-chromosomal data is the West-African region along the Bight of Benin, despite its important history in the trans-Atlantic slave trade and its large number of ethnic groups, languages and lifestyles. In this study, Y-chromosomal haplotypes from four Beninese populations were determined and a global meta-analysis with available Y-SNP and Y-STR data from populations along the Bight of Benin and surrounding areas was performed. A thorough methodology was developed allowing comparison of population samples using Y-chromosomal lineage data based on different Y-SNP panels and phylogenies. Geographic proximity turned out to be the best predictor of genetic affinity between populations along the Bight of Benin. Nevertheless, based on Y-chromosomal data from the literature two population samples differed strongly from others from the same or neighbouring areas and are not regionally representative within large-scale studies. Furthermore, the analysis of the HapMap sample YRI of a Yoruban population from South-western Nigeria based on Y-SNPs and Y-STR data showed for the first time its regional representativeness, a result which is important for standard population and forensic genetic applications using the YRI sample. Therefore, the uniquely and powerful geographical information carried by the Y-chromosome makes it an important locus to test the representativeness of a certain sample even in the genomic era, especially in poorly investigated areas like Africa.