Molecular dissection of the basal clades in the human Y chromosome phylogenetic tree

An update and frequency distribution of Y chromosome haplogroups in modern Japanese males

Japanese males belong to the Y chromosome C1a1, C2, D1a2a, D1a2a-12f2b, O1b2, O1b2a1a1, O2a2b1, and O2a1b haplogroups. Notably, the regional frequency of each haplogroup is homogeneous. Owing to recent developments in genome sequencing technology, the phylogenetic tree of Y chromosome haplogroups is updated annually. Therefore, in this study, we aimed to provide an update on the Y chromosome haplogroups of modern Japanese males and examine their regional distributions. Using 1,640 samples of Japanese males from seven Japanese cities (Nagasaki, Fukuoka, Tokushima, Osaka, Kanazawa, Kawasaki, and Sapporo), haplogroups C1a1, C2, D1a2a, D1a2a-12f2b, O1b2, and O1b2a1a1 were updated based on the latest phylogenetic tree. Haplogroup C1a1 was mainly classified into C1a1a1a and C1a1a1b subgroups; C1a1a1b was more common in Tokushima and Osaka than in the other regions. Haplogroup C2 was mainly classified into C2a, C2b1a1a, C2b1a1b, C2b1a2, and C2b1b subgroups and exhibited frequency differences in Osaka. Haplogroup D1a2a was classified into D1a2a1c1 and D1a2a2 subgroups, and its frequency varied between Tokushima and Osaka. Haplogroup D1a2a-12f2b was classified into D1a2a1a2b1a1a and D1a2a1a3 subgroups; however, no significant frequency differences were observed. Haplogroup O1b2 was classified into O1b2a1a2a1a, O1b2a1a2a1b, and O1b2a1a3 subgroups, with frequency differences between Nagasaki and Kanazawa. Haplogroup O1b2a1a1 was mainly classified into O1b2a1a1a, O1b2a1a1b, and O1b2a1a1c subgroups; however, no significant frequency differences were observed. Our findings suggest that gene flow in the Kinki region is caused by human migration.

New insights into the evolution of human Y chromosome palindromes through mutation and gene conversion

About one-quarter of the euchromatic portion of the male-specific region of the human Y chromosome consists of large duplicated sequences that are organized in eight palindromes (termed P1-P8), which undergo arm-to arm gene conversion, a proposed mechanism for maintaining their sequence integrity. Although the relevance of gene conversion in the evolution of palindromic sequences has been profoundly recognized, the dynamic of this mechanism is still nuanced. To shed light into the evolution of these genomic elements, we performed a high-depth (50×) targeted next-generation sequencing of the palindrome P6 in 157 subjects belonging to the most divergent evolutionary lineages of the Y chromosome. We found 118 new paralogous sequence variants, which were placed into the context of a robust Y chromosome phylogeny based on 7240 SNPs of the X-degenerate region. We mapped along the phylogeny 80 gene conversion events that shaped the diversity of P6 arms during recent human history. In contrast to previous studies, we demonstrated that arm-to-arm gene conversion, which occurs at a rate of 6.01 × 10 -6 conversions/base/year, is not biased toward the retention of the ancestral state of sequences. We also found a significantly lower mutation rate of the arms (6.18 × 10-10 mutations/base/year) compared with the spacer (9.16 × 10-10 mutations/base/year), a finding that may explain the observed higher inter-species conservation of arms, without invoking any bias of conversion. Finally, by formally testing the mutation/conversion balance in P6, we found that the arms of this palindrome reached a steady-state equilibrium between mutation and gene conversion.

The unidirectional phylogeny of Homo sapiens anchors the origin of modern humans in Eurasia

BACKGROUND

The Out of Africa hypothesis, OOAH, was challenged recently in an extended mtDNA analysis, PPA (Progressive Phylogenetic Analysis), that identified the African human populations as paraphyletic, a finding that contradicted the common OOAH understanding that Hss had originated in Africa and invaded Eurasia from there. The results were consistent with the molecular Out of Eurasia hypothesis, OOEH, and Eurasian palaeontology, a subject that has been largely disregarded in the discussion of OOAH.

RESULTS

In the present study the mtDNA tree, a phylogeny based on maternal inheritance, was compared to the nuclear DNA tree of the paternally transmitted Y-chromosome haplotypes, Y-DNAs. The comparison showed full phylogenetic coherence between these two separate sets of data. The results were consistent with potentially four translocations of modern humans from Eurasia into Africa, the earliest taking place ≈ 250,000 years before present, YBP. The results were in accordance with the postulates behind OOEH at the same time as they lent no support to the OOAH.

CONCLUSIONS

The conformity between the mtDNA and Y-DNA phylogenies of Hss is consistent with the understanding that Eurasia was the donor and not the receiver in human evolution. The evolutionary problems related to OOAH became similarly exposed by the mtDNA introgression that took place from Hss into Neanderthals ≈ 500,000 YBP, a circumstance that demonstrated the early coexistence of the two lineages in Eurasia.

Searching for the roots of the first free African American community

San Basilio de Palenque is an Afro-descendant community near Cartagena, Colombia, founded in the sixteenth century. The recognition of the historical and cultural importance of Palenque has promoted several studies, namely concerning the African roots of its first inhabitants. To deepen the knowledge of the origin and diversity of the Palenque parental lineages, we analysed a sample of 81 individuals for the entire mtDNA Control Region as well as 92 individuals for 27 Y-STRs and 95 for 51 Y-SNPs. The results confirmed the strong isolation of the Palenque, with some degree of influx of Native American maternal lineages, and a European admixture exclusively mediated by men. Due to the high genetic drift observed, a pairwise F_ST analysis with available data on African populations proved to be inadequate for determining population affinities. In contrast, when a phylogenetic approach was used, it was possible to infer the phylogeographic origin of some lineages in Palenque. Contradicting previous studies indicating a single African origin, our results evidence parental genetic contributions from widely different African regions.

A Late Pleistocene human humerus from Rusinga Island, Lake Victoria, Kenya

In 2010, a hominin right humerus fragment (KNM-RU 58330) was surface collected in a small gully at Nyamita North in the Late Pleistocene Wasiriya Beds of Rusinga Island, Kenya. A combination of stratigraphic and geochronological evidence suggests the specimen is likely between ∼49 and 36 ka in age. The associated fauna is diverse and dominated by semiarid grassland taxa. The small sample of associated Middle Stone Age artifacts includes Levallois flakes, cores, and retouched points. The 139 mm humeral fragment preserves the shaft from distal to the lesser tubercle to 14 mm below the distal end of the weakly projecting deltoid tuberosity. Key morphological features include a narrow and weakly marked pectoralis major insertion and a distinctive medial bend in the diaphysis at the deltoid insertion. This bend is unusual among recent human humeri but occurs in a few Late Pleistocene humeri. The dimensions of the distal end of the fragment predict a length of 317.9 ± 16.4 mm based on recent samples of African ancestry. A novel method of predicting humeral length from the distance between the middle of the pectoralis major and the bottom of the deltoid insertion predicts a length of 317.3 mm ± 17.6 mm. Cross-sectional geometry at the midshaft shows a relatively high percentage of cortical bone and a moderate degree of flattening of the shaft. The Nyamita humerus is anatomically modern in its morphology and adds to the small sample of hominins from the Late Pleistocene associated with Middle Stone Age artifacts known from East Africa. It may sample a population closely related to the people of the out-of-Africa migration.

Y-Chromosome Variation in Southern African Khoe-San Populations Based on Whole-Genome Sequences

Although the human Y chromosome has effectively shown utility in uncovering facets of human evolution and population histories, the ascertainment bias present in early Y-chromosome variant data sets limited the accuracy of diversity and TMRCA estimates obtained from them. The advent of next-generation sequencing, however, has removed this bias and allowed for the discovery of thousands of new variants for use in improving the Y-chromosome phylogeny and computing estimates that are more accurate. Here, we describe the high-coverage sequencing of the whole Y chromosome in a data set of 19 male Khoe-San individuals in comparison with existing whole Y-chromosome sequence data. Due to the increased resolution, we potentially resolve the source of haplogroup B-P70 in the Khoe-San, and reconcile recently published haplogroup A-M51 data with the most recent version of the ISOGG Y-chromosome phylogeny. Our results also improve the positioning of tentatively placed new branches of the ISOGG Y-chromosome phylogeny. The distribution of major Y-chromosome haplogroups in the Khoe-San and other African groups coincide with the emerging picture of African demographic history; with E-M2 linked to the agriculturalist Bantu expansion, E-M35 linked to pastoralist eastern African migrations, B-M112 linked to earlier east-south gene flow, A-M14 linked to shared ancestry with central African rainforest hunter-gatherers, and A-M51 potentially unique to the Khoe-San.

The study of human Y chromosome variation through ancient DNA

High throughput sequencing methods have completely transformed the study of human Y chromosome variation by offering a genome-scale view on genetic variation retrieved from ancient human remains in context of a growing number of high coverage whole Y chromosome sequence data from living populations from across the world. The ancient Y chromosome sequences are providing us the first exciting glimpses into the past variation of male-specific compartment of the genome and the opportunity to evaluate models based on previously made inferences from patterns of genetic variation in living populations. Analyses of the ancient Y chromosome sequences are challenging not only because of issues generally related to ancient DNA work, such as DNA damage-induced mutations and low content of endogenous DNA in most human remains, but also because of specific properties of the Y chromosome, such as its highly repetitive nature and high homology with the X chromosome. Shotgun sequencing of uniquely mapping regions of the Y chromosomes to sufficiently high coverage is still challenging and costly in poorly preserved samples. To increase the coverage of specific target SNPs capture-based methods have been developed and used in recent years to generate Y chromosome sequence data from hundreds of prehistoric skeletal remains. Besides the prospects of testing directly as how much genetic change in a given time period has accompanied changes in material culture the sequencing of ancient Y chromosomes allows us also to better understand the rate at which mutations accumulate and get fixed over time. This review considers genome-scale evidence on ancient Y chromosome diversity that has recently started to accumulate in geographic areas favourable to DNA preservation. More specifically the review focuses on examples of regional continuity and change of the Y chromosome haplogroups in North Eurasia and in the New World.

Evidence of extensive non-allelic gene conversion among LTR elements in the human genome

Long Terminal Repeats (LTRs) are nearly identical DNA sequences found at either end of Human Endogenous Retroviruses (HERVs). The high sequence similarity that exists among different LTRs suggests they could be substrate of ectopic gene conversion events. To understand the extent to which gene conversion occurs and to gain new insights into the evolutionary history of these elements in humans, we performed an intra-species phylogenetic study of 52 LTRs on different unrelated Y chromosomes. From this analysis, we obtained direct evidence that demonstrates the occurrence of ectopic gene conversion in several LTRs, with donor sequences located on both sex chromosomes and autosomes. We also found that some of these elements are characterized by an extremely high density of polymorphisms, showing one of the highest nucleotide diversities in the human genome, as well as a complex patchwork of sequences derived from different LTRs. Finally, we highlighted the limits of current short-read NGS studies in the analysis of genetic diversity of the LTRs in the human genome. In conclusion, our comparative re-sequencing analysis revealed that ectopic gene conversion is a common event in the evolution of LTR elements, suggesting complex genetic links among LTRs from different chromosomes.

Refining the Y chromosome phylogeny with southern African sequences

The recent availability of large-scale sequence data for the human Y chromosome has revolutionized analyses of and insights gained from this non-recombining, paternally inherited chromosome. However, the studies to date focus on Eurasian variation, and hence the diversity of early-diverging branches found in Africa has not been adequately documented. Here, we analyze over 900 kb of Y chromosome sequence obtained from 547 individuals from southern African Khoisan- and Bantu-speaking populations, identifying 232 new sequences from basal haplogroups A and B. We identify new clades in the phylogeny, an older age for the root, and substantially older ages for some individual haplogroups. Furthermore, while haplogroup B2a is traditionally associated with the spread of Bantu speakers, we find that it probably also existed in Khoisan groups before the arrival of Bantu speakers. Finally, there is pronounced variation in branch length between major haplogroups; in particular, haplogroups associated with Bantu speakers have significantly longer branches. Technical artifacts cannot explain this branch length variation, which instead likely reflects aspects of the demographic history of Bantu speakers, such as recent population expansion and an older average paternal age. The influence of demographic factors on branch length variation has broader implications both for the human Y phylogeny and for similar analyses of other species.

Confirmation of Y haplogroup tree topologies with newly suggested Y-SNPs for the C2, O2b and O3a subhaplogroups

Y chromosome single nucleotide polymorphisms (Y-SNPs) are useful markers for reconstructing male lineages through hierarchically arranged allelic sets known as haplogroups, and are thereby widely used in the fields such as human evolution, anthropology and forensic genetics. The Y haplogroup tree was recently revised with newly suggested Y-SNP markers for designation of several subgroups of haplogroups C2, O2b and O3a, which are predominant in Koreans. Therefore, herein we analyzed these newly suggested Y-SNPs in 545 unrelated Korean males who belong to the haplogroups C2, O2b or O3a, and investigated the reconstructed topology of the Y haplogroup tree. We were able to confirm that markers L1373, Z1338/JST002613-27, Z1300, CTS2657, Z8440 and F845 define the C2 subhaplogroups, C2b, C2e, C2e1, C2e1a, C2e1b and C2e2, respectively, and that markers F3356, L682, F11, F238/F449 and F444 define the O subhaplogroups O2b1, O2b1b, O3a1c1, O3a1c2 and O3a2c1c, respectively. Among six C2 subhaplogroups (C2b, C2e, C2e1*, C2e1a, C2e1b and C2e2), the C2e haplogroup and its subhaplogroups were found to be predominant, and among the four O2b subhaplogroups (O2b*, O2b1*, O2b1a and O2b1b), O2b1b was most frequently observed. Among the O3a subhaplogroups, O3a2c1 was predominant and it was further divided into the subhaplogroups O3a2c1a and O3a2c1c with a newly suggested marker. However, the JST002613-27 marker, which had been known to define the haplogroup C2f, was found to be an ancestral marker of the C2e haplogroup, as is the Z1338 marker. Also, the M312 marker for the O2b1 haplogroup designation was replaced by F3356, because all of the O2b1 haplotypes showed a nucleotide change at F3356, but not at M312. In addition, the F238 marker was always observed to be phylogenetically equivalent to F449, while both of the markers were assigned to the O3a1c2 haplogroup. The confirmed phylogenetic tree of this study with the newly suggested Y-SNPs could be valuable for anthropological and forensic investigations of East Asians including Koreans.

Towards a consensus Y-chromosomal phylogeny and Y-SNP set in forensics in the next-generation sequencing era

Currently, several different Y-chromosomal phylogenies and haplogroup nomenclatures are presented in scientific literature and at conferences demonstrating the present diversity in Y-chromosomal phylogenetic trees and Y-SNP sets used within forensic and anthropological research. This situation can be ascribed to the exponential growth of the number of Y-SNPs discovered due to mostly next-generation sequencing (NGS) studies. As Y-SNPs and their respective phylogenetic positions are important in forensics, such as for male lineage characterization and paternal bio-geographic ancestry inference, there is a need for forensic geneticists to know how to deal with these newly identified Y-SNPs and phylogenies, especially since these phylogenies are often created with other aims than to carry out forensic genetic research. Therefore, we give here an overview of four categories of currently used Y-chromosomal phylogenies and the associated Y-SNP sets in scientific research in the current NGS era. We compare these categories based on the construction method, their advantages and disadvantages, the disciplines wherein the phylogenetic tree can be used, and their specific relevance for forensic geneticists. Based on this overview, it is clear that an up-to-date reduced tree with a consensus Y-SNP set and a stable nomenclature will be the most appropriate reference resource for forensic research. Initiatives to reach such an international consensus are therefore highly recommended.

New clustering methods for population comparison on paternal lineages

The goal of this study is to show two new clustering and visualising techniques developed to find the most typical clusters of 18-dimensional Y chromosomal haplogroup frequency distributions of 90 Western Eurasian populations. The first technique called "self-organizing cloud (SOC)" is a vector-based self-learning method derived from the Self Organising Map and non-metric Multidimensional Scaling algorithms. The second technique is a new probabilistic method called the "maximal relation probability" (MRP) algorithm, based on a probability function having its local maximal values just in the condensation centres of the input data. This function is calculated immediately from the distance matrix of the data and can be interpreted as the probability that a given element of the database has a real genetic relation with at least one of the remaining elements. We tested these two new methods by comparing their results to both each other and the k-medoids algorithm. By means of these new algorithms, we determined 10 clusters of populations based on the similarity of haplogroup composition. The results obtained represented a genetically, geographically and historically well-interpretable picture of 10 genetic clusters of populations mirroring the early spread of populations from the Fertile Crescent to the Caucasus, Central Asia, Arabia and Southeast Europe. The results show that a parallel clustering of populations using SOC and MRP methods can be an efficient tool for studying the demographic history of populations sharing common genetic footprints.

Inter- and intraspecies phylogenetic analyses reveal extensive X-Y gene conversion in the evolution of gametologous sequences of human sex chromosomes

It has long been believed that the male-specific region of the human Y chromosome (MSY) is genetically independent from the X chromosome. This idea has been recently dismissed due to the discovery that X–Y gametologous gene conversion may occur. However, the pervasiveness of this molecular process in the evolution of sex chromosomes has yet to be exhaustively analyzed. In this study, we explored how pervasive X–Y gene conversion has been during the evolution of the youngest stratum of the human sex chromosomes. By comparing about 0.5 Mb of human–chimpanzee gametologous sequences, we identified 19 regions in which extensive gene conversion has occurred. From our analysis, two major features of these emerged: 1) Several of them are evolutionarily conserved between the two species and 2) almost all of the 19 hotspots overlap with regions where X–Y crossing-over has been previously reported to be involved in sex reversal. Furthermore, in order to explore the dynamics of X–Y gametologous conversion in recent human evolution, we resequenced these 19 hotspots in 68 widely divergent Y haplogroups and used publicly available single nucleotide polymorphism data for the X chromosome. We found that at least ten hotspots are still active in humans. Hence, the results of the interspecific analysis are consistent with the hypothesis of widespread reticulate evolution within gametologous sequences in the differentiation of hominini sex chromosomes. In turn, intraspecific analysis demonstrates that X–Y gene conversion may modulate human sex-chromosome-sequence evolution to a greater extent than previously thought.

Y-chromosome E haplogroups: their distribution and implication to the origin of Afro-Asiatic languages and pastoralism

Archeological and paleontological evidences point to East Africa as the likely area of early evolution of modern humans. Genetic studies also indicate that populations from the region often contain, but not exclusively, representatives of the more basal clades of mitochondrial and Y-chromosome phylogenies. Most Y-chromosome haplogroup diversity in Africa, however, is present within macrohaplogroup E that seem to have appeared 21 000-32 000 YBP somewhere between the Red Sea and Lake Chad. The combined analysis of 17 bi-allelic markers in 1214 Y chromosomes together with cultural background of 49 populations displayed in various metrics: network, multidimensional scaling, principal component analysis and neighbor-joining plots, indicate a major contribution of East African populations to the foundation of the macrohaplogroup, suggesting a diversification that predates the appearance of some cultural traits and the subsequent expansion that is more associated with the cultural and linguistic diversity witnessed today. The proto-Afro-Asiatic group carrying the E-P2 mutation may have appeared at this point in time and subsequently gave rise to the different major population groups including current speakers of the Afro-Asiatic languages and pastoralist populations.

An unbiased resource of novel SNP markers provides a new chronology for the human Y chromosome and reveals a deep phylogenetic structure in Africa

Sequence diversity and the ages of the deepest nodes of the MSY phylogeny remain largely unexplored due to the severely biased collection of SNPs available for study. We characterized 68 worldwide Y chromosomes by high-coverage next-generation sequencing, including 18 deep-rooting ones, and identified 2386 SNPs, 80% of which were novel. Many aspects of this pool of variants resembled the pattern observed among genome-wide de novo events, suggesting that in the MSY, a large proportion of newly arisen alleles has survived in the phylogeny. Some degree of purifying selection emerged in the form of an excess of private missense variants. Our tree recapitulated the previously known topology, but the relative lengths of major branches were drastically modified and the associated node ages were remarkably older. We found significantly different branch lengths when comparing the rare deep-rooted A1b African lineage with the rest of the tree. Our dating results and phylogeography led to the following main conclusions: (1) Patrilineal lineages with ages approaching those of early AMH fossils survive today only in central-western Africa; (2) only a few evolutionarily successful MSY lineages survived between 160 and 115 kya; and (3) an early exit out of Africa (before 70 kya), which fits recent western Asian archaeological evidence, should be considered. Our experimental design produced an unbiased resource of new MSY markers informative for the initial formation of the anatomically modern human gene pool, i.e., a period of our evolution that had been previously considered to be poorly accessible with paternally inherited markers.

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.

Multiplex genotyping assays for fine-resolution subtyping of the major human Y-chromosome haplogroups E, G, I, J, and R in anthropological, genealogical, and forensic investigations

Inherited DNA polymorphisms located within the nonrecombing portion of the human Y chromosome provide a powerful means of tracking the patrilineal ancestry of male individuals. Recently, we introduced an efficient genotyping method for the detection of the basal Y-chromosome haplogroups A to T, as well as an additional method for the dissection of haplogroup O into its sublineages. To further extend the use of the Y chromosome as an evolutionary marker, we here introduce a set of genotyping assays for fine-resolution subtyping of haplogroups E, G, I, J, and R, which make up the bulk of Western Eurasian and African Y chromosomes. The marker selection includes a total of 107 carefully selected bi-allelic polymorphisms that were divided into eight hierarchically organized multiplex assays (two for haplogroup E, one for I, one for J, one for G, and three for R) based on the single-base primer extension (SNaPshot) technology. Not only does our method allow for enhanced Y-chromosome lineage discrimination, the more restricted geographic distribution of the subhaplogroups covered also enables more fine-scaled estimations of patrilineal bio-geographic origin. Supplementing our previous method for basal Y-haplogroup detection, the currently introduced assays are thus expected to be of major relevance for future DNA studies targeting male-specific ancestry for forensic, anthropological, and genealogical purposes.

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.