Increased Resolution of Y Chromosome Haplogroup T Defines Relationships among Populations of the Near East, Europe, and Africa

High Coverage Mitogenomes and Y-Chromosomal Typing Reveal Ancient Lineages in the Modern-Day Székely Population in Romania

Here we present 115 whole mitogenomes and 92 Y-chromosomal Short Tandem Repeat (STR) and Single Nucleotide Polymorphism (SNP) profiles from a Hungarian ethnic group, the Székelys (in Romanian: Secuii, in German: Sekler), living in southeast Transylvania (Romania). The Székelys can be traced back to the 12th century in the region, and numerous scientific theories exist as to their origin. We carefully selected sample providers that had local ancestors inhabiting small villages in the area of Odorheiu Secuiesc/Székelyudvarhely in Romania. The results of our research and the reported data signify a qualitative leap compared to previous studies since it presents the first complete mitochondrial DNA sequences and Y-chromosomal profiles of 23 STRs from the region. We evaluated the results with population genetic and phylogenetic methods in the context of the modern and ancient populations that are either geographically or historically related to the Székelys. Our results demonstrate a predominantly local uniparental make-up of the population that also indicates limited admixture with neighboring populations. Phylogenetic analyses confirmed the presumed eastern origin of certain maternal (A, C, D) and paternal (Q, R1a) lineages, and, in some cases, they could also be linked to ancient DNA data from the Migration Period (5th-9th centuries AD) and Hungarian Conquest Period (10th century AD) populations.

Distribution of Y-chromosome haplogroups in Serbian population groups originating from historically and geographically significant distinct parts of the Balkan Peninsula

Our study enrolled 1200 Serbian males originating from three geographical regions in the Balkan Peninsula inhabited by Serbs: present-day Serbia, regions of Old Herzegovina and Kosovo and Metohija. These samples were genotyped using the combination of 23 Y-chromosomal short tandem repeats (Y-STRs) loci and 17 Ychromosomal single nucleotide polymorphisms (Y-SNPs) loci for the haplotype and haplogroup analysis in order to characterize in detail Y chromosome flow in the recent history. Serbia's borders have changed through history, forcing Serbs constantly to migrate to different regions of Balkan Peninsula. The most significant migration waves in the recent history towards present-day Serbia occurred from the regions of Old- Herzegovina and Kosovo and Metohija that lie in the south-west/south. High haplotype diversity and discrimination capacity were observed in all three datasets, with the highest number of unique haplotypes (381) and discrimination capacity (0.97) detected in the samples originating from the present-day Serbia. Haplogroup composition didn't differ significantly among datasets, with three dominant haplogroups (I-M170, E-P170 and R-M198), and haplogroup I-M170 being the most frequent in all three datasets. Haplogroup E-P170 was the second most dominant in the dataset originating from geographical region of Kosovo and Metohija, whereas haplogroup R-M198 was the second most prevalent in the dataset from historical region of Old Herzegovina. Based on the phylogenetic three for haplogroup I constructed within this study, haplogroup I2a1-P37.2 was the most dominant within all three datasets, especially in the dataset from historical region of Old Herzegovina, where 182 out of 400 samples were derived for SNP P37.2. Genetic distances between three groups of samples, evaluated by the Fst and Rst statistical values, and further visualized through multidimensional scaling plot, showed great genetic similarity between datasets from Old Herzegovina and present-day Serbia. Genetic difference in the haplogroup distribution and frequency between datasets from historical region of Old Herzegovina and from geographical region of Kosovo and Metohija was confirmed with highest Fst and Rst vaules. In this study we have distinguished genetic structure, diversity and haplogroup frequencies within 1200 Serbian males from three datasets, relationships among them as well as with other Balkan and European populations, which is useful for studying recent demographic history.

Croatian genetic heritage: an updated Y-chromosome story

AIM

To analyze an additional set of ˝Y-chromosome genetic markers to acquire a more detailed insight into the diversity of the Croatian population.

METHODS

A total of 518 Yfiler Plus profiles were genotyped. Allele frequencies, haplotype frequencies, and haplotype diversity were calculated by using the STRAF software v. 2.0.4. Genetic distances were quantified by Rst with AMOVA online tool from the YHRD. The evolutionary history was inferred with the neighbor-joining method of phylogenetic tree construction in the MEGAX software. Whit Athey's Haplogroup Predictor v. 5 was used for additional comparison with regional and other European populations.

RESULTS

A total of 507 haplotypes were used for genetic STR analysis. An interpopulation study on 17 Y-STR markers showed the lowest genetic diversity between the Croatian and Bosnian-Herzegovinian populations and the highest between the Croatian and Irish populations. Additional interpopulation comparison with the original 27 Y-STR markers (for the population with available data) was also performed. A total of 518 haplotypes were used in the determination of haplogroup diversity. Haplogroup I with its sublineage I2a expressed the highest prevalence. The second most prevalent haplogroup was R, with its major sublineage R1a, except for the subpopulation of Hvar, where E1b1b was the second most prevalent haplogroup. Rare haplogroups also confirmed in this study were L, T, and Q. G1 was detected for the first time in the Croatian population.

CONCLUSION

We obtained a new insight into the differences between examined subpopulations of Croatia and their possible (dis)similarities with neighboring and distant populations.

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.

Ancient genomes from a rural site in Imperial Rome (1st-3rd cent. CE): a genetic junction in the Roman Empire

BACKGROUND

Rome became the prosperous Capital of the Roman Empire through the political and military conquests of neighbouring areas. People were able to move Romeward modifying the Rome area's demographic structure. However, the genomic evidence for the population of one of the broadest Empires in antiquity has been sparse until recently.

AIM

The genomic analysis of people buried in Quarto Cappello del Prete (QCP) necropolis was carried out to help elucidate the genomic structure of Imperial Rome inhabitants.

SUBJECTS AND METHODS

We recruited twenty-five individuals from QCP for ancient DNA analysis through whole-genome sequencing. Multiple investigations were carried out to unveil the genetic components featuring in the studied samples and the community's putative demographic structure.

RESULTS

We generated reliable whole-genome data for 7 samples surviving quality controls. The distribution of Imperial Romans from QCP partly overlaps with present-day Southern Mediterranean and Southern-Near Eastern populations.

CONCLUSION

The genomic legacy with the south-eastern shores of the Mediterranean Sea and the Central and Western Northern-African coast funerary influence pave the way for considering people buried in QCP as resembling a Punic-derived human group.

Insights into the Middle Eastern paternal genetic pool in Tunisia: high prevalence of T-M70 haplogroup in an Arab population

To obtain refreshed insights into the paternal lineages of Tunisian populations, Y-chromosome diversity was assessed in two populations belonging to an Arab genealogical lineage, Kairouan and Wesletia, as well as in four Tunisian Andalusian populations, Testour, Slouguia, Qalaat-El-Andalous and El Alia. The Arabs from Kairouan revealed 73.47% of E-M81 and close affinities with Berber groups, indicating they are likely arabized Berbers, clearly differentiated from the Arabs from Wesletia, who harbored the highest frequency (71.8%) of the Middle Eastern component ever observed in North Africa. In the Tunisian Andalusians, the North African component largely prevailed, followed by the Middle Eastern contribution. Global comparative analysis highlighted the heterogeneity of Tunisian populations, among which, as a whole, dominated a set of lineages ascribed to be of autochthonous Berber origin (71.67%), beside a component of essentially Middle Eastern extraction (18.35%), and signatures of Sub-Saharan (5.2%), European (3.45%) and Asiatic (1.33%) contributions. The remarkable frequency of T-M70 in Wesletia (17.4%) prompted to refine its phylogeographic analysis, allowing to confirm its Middle Eastern origin, though signs of local evolution in Northern Africa were also detected. Evidence was clear on the ancient introduction of T lineages into the region, probably since Neolithic times associated to spread of agriculture.

Haplogroup Prediction Using Y-Chromosomal Short Tandem Repeats in the General Population of Bosnia and Herzegovina

Human Y-chromosomal haplogroups are an important tool used in population genetics and forensic genetics. A conventional method used for Y haplogroup assignment is based on a set of Y-single nucleotide polymorphism (SNP) markers deployed, which exploits the low mutation rate nature of these markers. Y chromosome haplogroups can be successfully predicted from Y-short tandem repeat (STR) markers using different software packages, and this method gained much attention recently due to its labor-, time-, and cost-effectiveness. The present study was based on the analysis of a total of 480 adult male buccal swab samples collected from different regions of Bosnia and Herzegovina. Y haplogroup prediction was performed using Whit Athey’s Haplogroup Predictor, based on haplotype data on 23 Y-STR markers contained within the PowerPlex® Y23 kit. The results revealed the existence of 14 different haplogroups, with I2a, R1a, and E1b1b being the most prevalent with frequencies of 43.13, 14.79, and 14.58%, respectively. Compared to the previously published studies on Bosnian-Herzegovinian population based on Y-SNP and Y-STR data, this study represents an upgrade of molecular genetic data with a significantly larger number of samples, thus offering more accurate results and higher probability of detecting rare haplogroups.

The assembly of caprine Y chromosome sequence reveals a unique paternal phylogenetic pattern and improves our understanding of the origin of domestic goat

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y-linked scaffolds using existing data. A total of 11.5 Mb Y-linked sequences were clustered into 33 scaffolds, and 187 protein-coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole-genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post-Neolithic period.

Y-chromosome variability and genetic history of Commons from Northern Italy

Objectives

Genetic drift and admixture are driving forces in human evolution, but their concerted impact to population evolution in historical times and at a micro‐geographic scale is poorly assessed. In this study we test a demographic model encompassing both admixture and drift to the case of social‐cultural isolates such as the so‐called “Commons.”

Materials and methods

Commons are peculiar institutions of medieval origins whose key feature is the tight relationship between population and territory, mediated by the collective property of shared resources. Here, we analyze the Y‐chromosomal genetic structure of four Commons (for a total of 366 samples) from the Central and Eastern Padana plain in Northern Italy.

Results

Our results reveal that all these groups exhibit patterns of significant diversity reduction, peripheral/outlier position within the Italian/European genetic space and high frequency of Common‐specific haplogroups. By explicitly testing different drift‐admixture models, we show that a drift‐only model is more probable for Central Padana Commons, while additional admixture (~20%) from external population around the same time of their foundation cannot be excluded for the Eastern ones.

Discussion

Building on these results, we suggest central Middle Ages as the most probable age of foundation for three of the considered Commons, the remaining one pointing to late antiquity. We conclude that an admixture‐drift model is particularly useful for interpreting the genetic structure and recent demographic history of small‐scale populations in which social‐cultural features play a significant role.

Y Haplogroup Diversity of the Dominican Republic: Reconstructing the Effect of the European Colonization and the Trans-Atlantic Slave Trades

The Dominican Republic is one of the two countries on the Hispaniola island, which is part of the Antilles. Hispaniola was affected by the European colonization and massive deportation of African slaves since the XVI century and these events heavily shaped the genetic composition of the present-day population. To shed light about the effect of the European rules, we analyzed 92 single nucleotide polymorphisms on the Y chromosome in 182 Dominican individuals from three different locations. The Dominican Y haplogroup composition was characterized by an excess of northern African/European lineages (59%), followed by the African clades (38%), whereas the Native-American lineages were rare (3%). The comparison with the mitochondrial DNA variability, dominated by African clades, revealed a sex-biased admixture pattern, in line with the colonial society dominated by European men. When other Caribbean and non-Caribbean former colonies were also considered, we noted a difference between territories under a Spanish rule (like the Dominican Republic) and British/French rule, with the former characterized by an excess of European Y lineages reflecting the more permissive Iberian legislation about mixed people and slavery. Finally, we analyzed the distribution in Africa of the Dominican lineages with a putative African origin, mainly focusing on central and western Africa, which were the main sources of African slaves. We found that most (83%) of the African lineages observed in Santo Domingo have a central African ancestry, suggesting that most of the slaves were deported from regions.

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.

Ancient DNA from Chalcolithic Israel reveals the role of population mixture in cultural transformation

The material culture of the Late Chalcolithic period in the southern Levant (4500–3900/3800 BCE) is qualitatively distinct from previous and subsequent periods. Here, to test the hypothesis that the advent and decline of this culture was influenced by movements of people, we generated genome-wide ancient DNA from 22 individuals from Peqi’in Cave, Israel. These individuals were part of a homogeneous population that can be modeled as deriving ~57% of its ancestry from groups related to those of the local Levant Neolithic, ~17% from groups related to those of the Iran Chalcolithic, and ~26% from groups related to those of the Anatolian Neolithic. The Peqi’in population also appears to have contributed differently to later Bronze Age groups, one of which we show cannot plausibly have descended from the same population as that of Peqi’in Cave. These results provide an example of how population movements propelled cultural changes in the deep past.

The Late Chalcolithic material culture in the southern Levant has unique attributes that suggest spread of people or culture. Here, the authors use genome-wide ancient DNA data from 22 individuals from a Chalcolithic site and show evidence of complex population movements and turnovers.

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 genetic variation in the R1a clade among the Ashkenazi Levites' Y chromosome

Approximately 300,000 men around the globe self-identify as Ashkenazi Levites, of whom two thirds were previously shown to descend from a single male. The paucity of whole Y-chromosome sequences precluded conclusive identification of this ancestor's age, geographic origin and migration patterns. Here, we report the variation of 486 Y-chromosomes within the Ashkenazi and non-Ashkenazi Levite R1a clade, other Ashkenazi Jewish paternal lineages, as well as non-Levite Jewish and non-Jewish R1a samples. Cumulatively, the emerging profile is of a Middle Eastern ancestor, self-affiliating as Levite, and carrying the highly resolved R1a-Y2619 lineage, which was likely a minor haplogroup among the Hebrews. A star-like phylogeny, coalescing similarly to other Ashkenazi paternal lineages, ~1,743 ybp, suggests it to be one of the Ashkenazi paternal founders; to have expanded as part of the overall Ashkenazi demographic expansion, without special relation to the Levite affiliation; and to have subsequently spread to non-Ashkenazi Levites.

Mapping human dispersals into the Horn of Africa from Arabian Ice Age refugia using mitogenomes

Rare mitochondrial lineages with relict distributions can sometimes be disproportionately informative about deep events in human prehistory. We have studied one such lineage, haplogroup R0a, which uniquely is most frequent in Arabia and the Horn of Africa, but is distributed much more widely, from Europe to India. We conclude that: (1) the lineage ancestral to R0a is more ancient than previously thought, with a relict distribution across the Mediterranean/Southwest Asia; (2) R0a has a much deeper presence in Arabia than previously thought, highlighting the role of at least one Pleistocene glacial refugium, perhaps on the Red Sea plains; (3) the main episode of dispersal into Eastern Africa, at least concerning maternal lineages, was at the end of the Late Glacial, due to major expansions from one or more refugia in Arabia; (4) there was likely a minor Late Glacial/early postglacial dispersal from Arabia through the Levant and into Europe, possibly alongside other lineages from a Levantine refugium; and (5) the presence of R0a in Southwest Arabia in the Holocene at the nexus of a trading network that developed after ~3 ka between Africa and the Indian Ocean led to some gene flow even further afield, into Iran, Pakistan and India.

Y chromosome diversity in a linguistic isolate (Mirandese, NE Portugal)

OBJECTIVES

The purpose of this study was to genetically characterize the male lineages of people who speak Mirandese, an interesting case of a linguistic relict that can still be found in the municipality of Miranda do Douro, NE Portugal. This region lies within the area of the Leonese dialects, which are remnants of the Romance dialects spoken in the Kingdom of Leon currently grouped in the Astur-Leonese linguistic continuum. We intended to disclose affinities with surrounding populations, namely from Spain where the Astur-Leonese is also spoken.

METHODS

Eighty-eight unrelated males (58 from Miranda and 30 from Bragança, the broad Portuguese region where Miranda is located) were genotyped with the combined use of 17 Y chromosome short tandem repeats (Y-STRs) and a high resolution Y chromosome single nucleotide polymorphism (Y-SNPs) strategy. Moreover, 236 males from Miranda and neighboring regions, previously classified as R-M269, were also genotyped.

RESULTS

R-P312 was the most frequent haplogroup in the Mirandese, followed by J-12f2.1 and T-M70. The male lineages J-12f2.1 and T-M70 were also well represented, and both were shared with descendants of Sephardic Jews. No signs of diversity reduction were detected.

CONCLUSIONS

Mirandese speakers display a Y chromosome gene pool that shows a subtle differentiation from neighboring populations, mainly attributable to the assimilation of lineages ascribed to be of Jewish ancestry. Although not revealing signs of geographic/linguistic isolation, no clear affinities with other Astur-Leonese populations were detected. The results suggest that in Miranda language sharing is not accompanied by significant gene flow between populations from both sides of the political border. Am. J. Hum. Biol. 28:671-680, 2016. © 2016 Wiley Periodicals, Inc.

A time- and cost-effective strategy to sequence mammalian Y Chromosomes: an application to the de novo assembly of gorilla Y

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

Exploring the Y Chromosomal Ancestry of Modern Panamanians

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama’s population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.

Genetic population study of Y-chromosome markers in Benin and Ivory Coast ethnic groups

Ninety-six single nucleotide polymorphisms (SNPs) and seventeen short tandem repeat (STRs) were investigated on the Y-chromosome of 288 unrelated healthy individuals from populations in Benin (Bariba, Yoruba, and Fon) and the Ivory Coast (Ahizi and Yacouba). We performed a multidimensional scaling analysis based on FST and RST genetic distances using a large extensive database of sub-Saharan African populations. There is more genetic homogeneity in Ivory Coast populations compared with populations from Benin. Notably, the Beninese Yoruba are significantly differentiated from neighbouring groups, but also from the Yoruba from Nigeria (FST>0.05; P<0.01). The Y-chromosome dataset presented here provides new valuable data to understand the complex genetic diversity and human male demographic events in West Africa.

Carriers of Mitochondrial DNA Macrohaplogroup N Lineages Reached Australia around 50,000 Years Ago following a Northern Asian Route

Background

The modern human colonization of Eurasia and Australia is mostly explained by a single-out-of-Africa exit following a southern coastal route throughout Arabia and India. However, dispersal across the Levant would better explain the introgression with Neanderthals, and more than one exit would fit better with the different ancient genomic components discovered in indigenous Australians and in ancient Europeans. The existence of an additional Northern route used by modern humans to reach Australia was previously deduced from the phylogeography of mtDNA macrohaplogroup N. Here, we present new mtDNA data and new multidisciplinary information that add more support to this northern route.

Methods

MtDNA hypervariable segments and haplogroup diagnostic coding positions were analyzed in 2,278 Saudi Arabs, from which 1,725 are new samples. Besides, we used 623 published mtDNA genomes belonging to macrohaplogroup N, but not R, to build updated phylogenetic trees to calculate their coalescence ages, and more than 70,000 partial mtDNA sequences were screened to establish their respective geographic ranges.

Results

The Saudi mtDNA profile confirms the absence of autochthonous mtDNA lineages in Arabia with coalescence ages deep enough to support population continuity in the region since the out-of-Africa episode. In contrast to Australia, where N(xR) haplogroups are found in high frequency and with deep coalescence ages, there are not autochthonous N(xR) lineages in India nor N(xR) branches with coalescence ages as deep as those found in Australia. These patterns are at odds with the supposition that Australian colonizers harboring N(xR) lineages used a route involving India as a stage. The most ancient N(xR) lineages in Eurasia are found in China, and inconsistently with the coastal route, N(xR) haplogroups with the southernmost geographical range have all more recent radiations than the Australians.

Conclusions

Apart from a single migration event via a southern route, phylogeny and phylogeography of N(xR) lineages support that people carrying mtDNA N lineages could have reach Australia following a northern route through Asia. Data from other disciplines also support this scenario.

Mitochondrial DNA and the Y chromosome suggest the settlement of Madagascar by Indonesian sea nomad populations

Background

Linguistic, cultural and genetic characteristics of the Malagasy suggest that both Africans and Island Southeast Asians were involved in the colonization of Madagascar. Populations from the Indonesian archipelago played an especially important role because linguistic evidence suggests that the Malagasy language branches from the Southeast Barito language family of southern Borneo, Indonesia, with the closest language spoken today by the Ma’anyan. To test for a genetic link between Malagasy and these linguistically related Indonesian populations, we studied the Ma’anyan and other Indonesian ethnic groups (including the sea nomad Bajo) that, from their historical and linguistic contexts, may be modern descendants of the populations that helped enact the settlement of Madagascar.

Result

A combination of phylogeographic analysis of genetic distances, haplotype comparisons and inference of parental populations by linear optimization, using both maternal and paternal DNA lineages, suggests that Malagasy derive from multiple regional sources in Indonesia, with a focus on eastern Borneo, southern Sulawesi and the Lesser Sunda islands.

Conclusion

Settlement may have been mediated by ancient sea nomad movements because the linguistically closest population, Ma’anyan, has only subtle genetic connections to Malagasy, whereas genetic links with other sea nomads are more strongly supported. Our data hint at a more complex scenario for the Indonesian settlement of Madagascar than has previously been recognized.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1394-7) contains supplementary material, which is available to authorized users.

Simultaneous analysis of hundreds of Y-chromosomal SNPs for high-resolution paternal lineage classification using targeted semiconductor sequencing

SNPs from the non-recombining part of the human Y chromosome (Y-SNPs) are informative to classify paternal lineages in forensic, genealogical, anthropological, and evolutionary studies. Although thousands of Y-SNPs were identified thus far, previous Y-SNP multiplex tools target only dozens of markers simultaneously, thereby restricting the provided Y-haplogroup resolution and limiting their applications. Here, we overcome this shortcoming by introducing a high-resolution multiplex tool for parallel genotyping-by-sequencing of 530 Y-SNPs using the Ion Torrent PGM platform, which allows classification of 432 worldwide Y haplogroups. Contrary to previous Y-SNP multiplex tools, our approach covers branches of the entire Y tree, thereby maximizing the paternal lineage classification obtainable. We used a default DNA input amount of 10 ng per reaction but preliminary sensitivity testing revealed positive results from as little as 100 pg input DNA. Furthermore, we demonstrate that sample pooling using barcodes is feasible, allowing increased throughput for lower per-sample costs. In addition to the wetlab protocol, we provide a software tool for automated data quality control and haplogroup classification. The unique combination of ultra-high marker density and high sensitivity achievable from low amounts of potentially degraded DNA makes this new multiplex tool suitable for a wide range of Y-chromosome applications.

Improved phylogenetic resolution and rapid diversification of Y-chromosome haplogroup K-M526 in Southeast Asia

The highly structured distribution of Y-chromosome haplogroups suggests that current patterns of variation may be informative of past population processes. However, limited phylogenetic resolution, particularly of subclades within haplogroup K, has obscured the relationships of lineages that are common across Eurasia. Here we genotype 13 new highly informative single-nucleotide polymorphisms in a worldwide sample of 4413 males that carry the derived allele at M526, and reconstruct an NRY haplogroup tree with significantly higher resolution for the major clade within haplogroup K, K-M526. Although K-M526 was previously characterized by a single polytomy of eight major branches, the phylogenetic structure of haplogroup K-M526 is now resolved into four major subclades (K2a-d). The largest of these subclades, K2b, is divided into two clusters: K2b1 and K2b2. K2b1 combines the previously known haplogroups M, S, K-P60 and K-P79, whereas K2b2 comprises haplogroups P and its subhaplogroups Q and R. Interestingly, the monophyletic group formed by haplogroups R and Q, which make up the majority of paternal lineages in Europe, Central Asia and the Americas, represents the only subclade with K2b that is not geographically restricted to Southeast Asia and Oceania. Estimates of the interval times for the branching events between M9 and P295 point to an initial rapid diversification process of K-M526 that likely occurred in Southeast Asia, with subsequent westward expansions of the ancestors of haplogroups R and Q.

Early back-to-Africa migration into the Horn of Africa

Genetic studies have identified substantial non-African admixture in the Horn of Africa (HOA). In the most recent genomic studies, this non-African ancestry has been attributed to admixture with Middle Eastern populations during the last few thousand years. However, mitochondrial and Y chromosome data are suggestive of earlier episodes of admixture. To investigate this further, we generated new genome-wide SNP data for a Yemeni population sample and merged these new data with published genome-wide genetic data from the HOA and a broad selection of surrounding populations. We used multidimensional scaling and ADMIXTURE methods in an exploratory data analysis to develop hypotheses on admixture and population structure in HOA populations. These analyses suggested that there might be distinct, differentiated African and non-African ancestries in the HOA. After partitioning the SNP data into African and non-African origin chromosome segments, we found support for a distinct African (Ethiopic) ancestry and a distinct non-African (Ethio-Somali) ancestry in HOA populations. The African Ethiopic ancestry is tightly restricted to HOA populations and likely represents an autochthonous HOA population. The non-African ancestry in the HOA, which is primarily attributed to a novel Ethio-Somali inferred ancestry component, is significantly differentiated from all neighboring non-African ancestries in North Africa, the Levant, and Arabia. The Ethio-Somali ancestry is found in all admixed HOA ethnic groups, shows little inter-individual variance within these ethnic groups, is estimated to have diverged from all other non-African ancestries by at least 23 ka, and does not carry the unique Arabian lactase persistence allele that arose about 4 ka. Taking into account published mitochondrial, Y chromosome, paleoclimate, and archaeological data, we find that the time of the Ethio-Somali back-to-Africa migration is most likely pre-agricultural.

The history of the North African mitochondrial DNA haplogroup U6 gene flow into the African, Eurasian and American continents

Background

Complete mitochondrial DNA (mtDNA) genome analyses have greatly improved the phylogeny and phylogeography of human mtDNA. Human mitochondrial DNA haplogroup U6 has been considered as a molecular signal of a Paleolithic return to North Africa of modern humans from southwestern Asia.

Results

Using 230 complete sequences we have refined the U6 phylogeny, and improved the phylogeographic information by the analysis of 761 partial sequences. This approach provides chronological limits for its arrival to Africa, followed by its spreads there according to climatic fluctuations, and its secondary prehistoric and historic migrations out of Africa colonizing Europe, the Canary Islands and the American Continent.

Conclusions

The U6 expansions and contractions inside Africa faithfully reflect the climatic fluctuations that occurred in this Continent affecting also the Canary Islands. Mediterranean contacts drove these lineages to Europe, at least since the Neolithic. In turn, the European colonization brought different U6 lineages throughout the American Continent leaving the specific sign of the colonizers origin.

Tracing Arab-Islamic inheritance in Madagascar: study of the Y-chromosome and mitochondrial DNA in the Antemoro

Madagascar is located at the crossroads of the Asian and African worlds and is therefore of particular interest for studies on human population migration. Within the large human diversity of the Great Island, we focused our study on a particular ethnic group, the Antemoro. Their culture presents an important Arab-Islamic influence, but the question of an Arab biological inheritance remains unresolved. We analyzed paternal (n=129) and maternal (n=135) lineages of this ethnic group. Although the majority of Antemoro genetic ancestry comes from sub-Saharan African and Southeast Asian gene pools, we observed in their paternal lineages two specific haplogroups (J1 and T1) linked to Middle Eastern origins. This inheritance was restricted to some Antemoro sub-groups. Statistical analyses tended to confirm significant Middle Eastern genetic contribution. This study gives a new perspective to the large human genetic diversity in Madagascar.

Seeing the wood for the trees: a minimal reference phylogeny for the human Y chromosome

During the last few decades, a wealth of studies dedicated to the human Y chromosome and its DNA variation, in particular Y-chromosome single-nucleotide polymorphisms (Y-SNPs), has led to the construction of a well-established Y-chromosome phylogeny. Since the recent advent of new sequencing technologies, the discovery of additional Y-SNPs is exploding and their continuous incorporation in the phylogenetic tree is leading to an ever higher resolution. However, the large and increasing amount of information included in the "complete" Y-chromosome phylogeny, which now already includes many thousands of identified Y-SNPs, can be overwhelming and complicates its understanding as well as the task of selecting suitable markers for genotyping purposes in evolutionary, demographic, anthropological, genealogical, medical, and forensic studies. As a solution, we introduce a concise reference phylogeny whereby we do not aim to provide an exhaustive tree that includes all known Y-SNPs but, rather, a quite stable reference tree aiming for optimal global discrimination capacity based on a strongly reduced set that includes only the most resolving Y-SNPs. Furthermore, with this reference tree, we wish to propose a common standard for Y-marker as well as Y-haplogroup nomenclature. The current version of our tree is based on a core set of 417 branch-defining Y-SNPs and is available online at http://www.phylotree.org/Y.

Increasing phylogenetic resolution still informative for Y chromosomal studies on West-European populations

Many Y-chromosomal lineages which are defined in the latest phylogenetic tree of the human Y chromosome by the Y Chromosome Consortium (YCC) in 2008 are distributed in (Western) Europe due to the fact that a large number of phylogeographic studies focus on this area. Therefore, the question arises whether newly discovered polymorphisms on the Y chromosome will still be interesting to study Western Europeans on a population genetic level. To address this question, the West-European region of Flanders (Belgium) was selected as study area since more than 1000 Y chromosomes from this area have previously been genotyped at the highest resolution of the 2008 YCC-tree and coupled to in-depth genealogical data. Based on these data the temporal changes of the population genetic pattern over the last centuries within Flanders were studied and the effects of several past gene flow events were identified. In the present study a set of recently reported novel Y-SNPs were genotyped to further characterize all those Flemish Y chromosomes that belong to haplogroups G, R-M269 and T. Based on this extended Y-SNP set the discrimination power increased drastically as previous large (sub-)haplogroups are now subdivided in several non-marginal groups. Next, the previously observed population structure within Flanders appeared to be the result of different gradients of independent sub-haplogroups. Moreover, for the first time within Flanders a significant East-West gradient was observed in the frequency of two R-M269 lineages, and this gradient is still present when considering the current residence of the DNA donors. Our results thus suggest that an update of the Y-chromosomal tree based on new polymorphisms is still useful to increase the discrimination power based on Y-SNPs and to study population genetic patterns in more detail, even in an already well-studied region such as Western Europe.

Updating the Y-chromosomal phylogenetic tree for forensic applications based on whole genome SNPs

The Y-chromosomal phylogenetic tree has a wide variety of important forensic applications and therefore it needs to be state-of-the-art. Nevertheless, since the last 'official' published tree many publications reported additional Y-chromosomal lineages and other phylogenetic topologies. Therefore, it is difficult for forensic scientists to interpret those reports and use an up-to-date tree and corresponding nomenclature in their daily work. Whole genome sequencing (WGS) data is useful to verify and optimise the current phylogenetic tree for haploid markers. The AMY-tree software is the first open access program which analyses WGS data for Y-chromosomal phylogenetic applications. Here, all published information is collected in a phylogenetic tree and the correctness of this tree is checked based on the first large analysis of 747 WGS samples with AMY-tree. The obtained result is one phylogenetic tree with all peer-reviewed reported Y-SNPs without the observed recurrent and ambiguous mutations. Nevertheless, the results showed that currently only the genomes of a limited set of Y-chromosomal (sub-)haplogroups is available and that many newly reported Y-SNPs based on WGS projects are false positives, even with high sequencing coverage methods. This study demonstrates the usefulness of AMY-tree in the process of checking the quality of the present Y-chromosomal tree and it accentuates the difficulties to enlarge this tree based on only WGS methods.

AMY-tree: an algorithm to use whole genome SNP calling for Y chromosomal phylogenetic applications

Background

Due to the rapid progress of next-generation sequencing (NGS) facilities, an explosion of human whole genome data will become available in the coming years. These data can be used to optimize and to increase the resolution of the phylogenetic Y chromosomal tree. Moreover, the exponential growth of known Y chromosomal lineages will require an automatic determination of the phylogenetic position of an individual based on whole genome SNP calling data and an up to date Y chromosomal tree.

Results

We present an automated approach, ‘AMY-tree’, which is able to determine the phylogenetic position of a Y chromosome using a whole genome SNP profile, independently from the NGS platform and SNP calling program, whereby mistakes in the SNP calling or phylogenetic Y chromosomal tree are taken into account. Moreover, AMY-tree indicates ambiguities within the present phylogenetic tree and points out new Y-SNPs which may be phylogenetically relevant. The AMY-tree software package was validated successfully on 118 whole genome SNP profiles of 109 males with different origins. Moreover, support was found for an unknown recurrent mutation, wrong reported mutation conversions and a large amount of new interesting Y-SNPs.

Conclusions

Therefore, AMY-tree is a useful tool to determine the Y lineage of a sample based on SNP calling, to identify Y-SNPs with yet unknown phylogenetic position and to optimize the Y chromosomal phylogenetic tree in the future. AMY-tree will not add lineages to the existing phylogenetic tree of the Y-chromosome but it is the first step to analyse whole genome SNP profiles in a phylogenetic framework.

A calibrated human Y-chromosomal phylogeny based on resequencing

We have identified variants present in high-coverage complete sequences of 36 diverse human Y chromosomes from Africa, Europe, South Asia, East Asia, and the Americas, representing eight major haplogroups. After restricting our analysis to 8.97 Mb of the unique male-specific Y sequence, we identified 6662 high-confidence variants, including single-nucleotide polymorphisms (SNPs), multi-nucleotide polymorphisms (MNPs), and indels. We constructed phylogenetic trees using these variants, or subsets of them, and recapitulated the known structure of the tree. Assuming a male mutation rate of 1 × 10⁻⁹ per base pair per year, the time depth of the tree (haplogroups A3-R) was ∼101,000–115,000 yr, and the lineages found outside Africa dated to 57,000–74,000 yr, both as expected. In addition, we dated a striking Paleolithic male lineage expansion to 41,000–52,000 yr ago and the node representing the major European Y lineage, R1b, to 4000–13,000 yr ago, supporting a Neolithic origin for these modern European Y chromosomes. In all, we provide a nearly 10-fold increase in the number of Y markers with phylogenetic information, and novel historical insights derived from placing them on a calibrated phylogenetic tree.

A haplotype at STAT2 Introgressed from neanderthals and serves as a candidate of positive selection in Papua New Guinea

Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (~80 thousand years ago) and is found throughout Eurasia at an average frequency of ~5%. Interestingly, N is found in Melanesian populations at ~10-fold higher frequency (~54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2-all of which are part of the same 250 kb introgressive haplotype-as good candidates.

Global genetic variation at OAS1 provides evidence of archaic admixture in Melanesian populations

Recent analysis of DNA extracted from two Eurasian forms of archaic human shows that more genetic variants are shared with humans currently living in Eurasia than with anatomically modern humans in sub-Saharan Africa. Although these genome-wide average measures of genetic similarity are consistent with the hypothesis of archaic admixture in Eurasia, analyses of individual loci exhibiting the signal of archaic introgression are needed to test alternative hypotheses and investigate the admixture process. Here, we provide a detailed sequence analysis of the innate immune gene OAS1, a locus with a divergent Melanesian haplotype that is very similar to the Denisova sequence from the Altai region of Siberia. We resequenced a 7-kb region encompassing the OAS1 gene in 88 individuals from six Old World populations (San, Biaka, Mandenka, French Basque, Han Chinese, and Papua New Guineans) and discovered previously unknown and ancient genetic variation. The 5' region of this gene has unusual patterns of diversity, including 1) higher levels of nucleotide diversity in Papuans than in sub-Saharan Africans, 2) very deep ancestry with an estimated time to the most recent common ancestor of >3 myr, and 3) a basal branching pattern with Papuan individuals on either side of the rooted network. A global geographic survey of >1,500 individuals showed that the divergent Papuan haplotype is nearly restricted to populations from eastern Indonesia and Melanesia. Polymorphic sites within this haplotype are shared with the draft Denisova genome over a span of ∼90 kb and are associated with an extended block of linkage disequilibrium, supporting the hypothesis that this haplotype introgressed from an archaic source that likely lived in Eurasia.

Neolithic patrilineal signals indicate that the Armenian plateau was repopulated by agriculturalists

Armenia, situated between the Black and Caspian Seas, lies at the junction of Turkey, Iran, Georgia, Azerbaijan and former Mesopotamia. This geographic position made it a potential contact zone between Eastern and Western civilizations. In this investigation, we assess Y-chromosomal diversity in four geographically distinct populations that represent the extent of historical Armenia. We find a striking prominence of haplogroups previously implicated with the Agricultural Revolution in the Near East, including the J2a-M410-, R1b1b1(*)-L23-, G2a-P15- and J1-M267-derived lineages. Given that the Last Glacial Maximum event in the Armenian plateau occured a few millennia before the Neolithic era, we envision a scenario in which its repopulation was achieved mainly by the arrival of farmers from the Fertile Crescent temporally coincident with the initial inception of farming in Greece. However, we detect very restricted genetic affinities with Europe that suggest any later cultural diffusions from Armenia to Europe were not associated with substantial amounts of paternal gene flow, despite the presence of closely related Indo-European languages in both Armenia and Southeast Europe.

An efficient multiplex genotyping approach for detecting the major worldwide human Y-chromosome haplogroups

The Y chromosome is paternally inherited and therefore serves as an evolutionary marker of patrilineal descent. Worldwide DNA variation within the non-recombining portion of the Y chromosome can be represented as a monophyletic phylogenetic tree in which the branches (haplogroups) are defined by at least one SNP. Previous human population genetics research has produced a wealth of knowledge about the worldwide distribution of Y-SNP haplogroups. Here, we apply previous and very recent knowledge on the Y-SNP phylogeny and Y-haplogroup distribution by introducing two multiplex genotyping assays that allow for the hierarchical detection of 28 Y-SNPs defining the major worldwide Y haplogroups. PCR amplicons were kept small to make the method sensitive and thereby applicable to DNA of limited amount and/or quality such as in forensic settings. These Y-SNP assays thus form a valuable tool for researchers in the fields of forensic genetics and genetic anthropology to infer a man's patrilineal bio-geographic ancestry from DNA.