CSYseq: The first Y-chromosome sequencing tool typing a large number of Y-SNPs and Y-STRs to unravel worldwide human population genetics.

Male-specific Y-chromosome (chrY) polymorphisms are interesting components of the DNA for population genetics. While single nucleotide polymorphisms (Y-SNPs) indicate distant evolutionary ancestry, short tandem repeats (Y-STRs) are able to identify close familial kinships. Detailed chrY analysis provides thus both biogeographical background information as paternal lineage identification. The rapid advancement of high-throughput massive parallel sequencing (MPS) technology in the past decade has revolutionized genetic research. Using MPS, single-base information of both Y-SNPs as Y-STRs can be analyzed in a single assay typing multiple samples at once. In this study, we present the first extensive chrY-specific targeted resequencing panel, the 'CSYseq', which simultaneously identifies slow mutating Y-SNPs as evolution markers and rapid mutating Y-STRs as patrilineage markers. The panel was validated by paired-end sequencing of 130 males, distributed over 65 deep-rooted pedigrees covering 1,279 generations. The CSYseq successfully targets 15,611 Y-SNPs including 9,014 phylogenetic informative Y-SNPs to identify 1,443 human evolutionary Y-subhaplogroup lineages worldwide. In addition, the CSYseq properly targets 202 Y-STRs, including 81 slow, 68 moderate, 27 fast and 26 rapid mutating Y-STRs to individualize close paternal relatives. The targeted chrY markers cover a high average number of reads (Y-SNP = 717, Y-STR = 150), easy interpretation, powerful discrimination capacity and chrY specificity. The CSYseq is interesting for research on different time scales: to identify evolutionary ancestry, to find distant family and to discriminate closely related males. Therefore, this panel serves as a unique tool valuable for a wide range of genetic-genealogical applications in interdisciplinary research within evolutionary, population, molecular, medical and forensic genetics.

YMrCA: Improving Y-chromosomal ancestor time estimation for DNA kinship research.

The Y-chromosome is a valuable kinship indicator in family history and forensic research. To reconstruct genealogies, the time to the most recent common ancestor (tMRCA) between paternal relatives can be estimated through Y-STR analysis. Existing models are the stepwise mutation model (SMM, only one-step Y-STR changes) and the infinite allele model (IAM, new allele per Y-STR change). In this study, these mutation models and all existing tMRCA calculators were validated through a genetic-genealogy database containing 1,120 biologically related genealogical pairs confirmed by 46 Y-STRs with known tMRCA (18,109 generations). Consistent under- and overestimation and broad confidence intervals were observed, leading to dubious tMRCA estimates. This is because they do not include individual mutation rates or multi-step changes and ignore hidden multiple, back, or parallel modifications. To improve tMRCA estimation, we developed a user-friendly calculator, the "YMrCA", including all previously mentioned mutation characteristics. After extensive validation, we observed that the YMrCA calculator demonstrated a promising performance. The YMrCA yields a significantly higher tMRCA success rate (96%; +20%) and a lower tMRCA error (7; -3) compared to the mutation models and all online tMRCA calculators. Therefore, YMrCA offers the next step towards more objective tMRCA estimation for DNA kinship research.

The black legend on the Spanish presence in the low countries: Verifying shared beliefs on genetic ancestry.

OBJECTIVES: War atrocities committed by the Spanish army in the Low Countries during the 16th century are so ingrained in the collective memory of Belgian and Dutch societies that they generally assume a signature of this history to be present in their genetic ancestry. Historians claim this assumption is a consequence of the so-called "Black Legend" and negative propaganda portraying and remembering Spanish soldiers as extreme sexual aggressors. The impact of the presence of Spaniards during the Dutch Revolt on the genetic variation in the Low Countries has been verified in this study. MATERIALS AND METHODS: A recent population genetic analysis of Iberian-associated Y-chromosomal variation among Europe is enlarged with representative samples of Dutch (N = 250) and Flemish (N = 1,087) males. Frequencies of these variants are also compared between donors whose oldest reported paternal ancestors lived in-nowadays Flemish-cities affected by so-called Spanish Furies (N = 116) versus other patrilineages in current Flemish territory (N = 971). RESULTS: The frequencies of Y-chromosomal markers Z195 and SRY2627 decline steeply going north from Spain and the data for the Flemish and Dutch populations fits within this pattern. No trend of higher frequencies of these variants has been found within the well-ascertained samples associated with Spanish Fury cities. DISCUSSION: Although sexual aggression did occur in the 16th century, these activities did not leave a traceable "Spanish" genetic signature in the autochthonous genome of the Low Countries. Our results support the view that the 'Black Legend' and historical propaganda on sexual aggression have nurtured today's incorrect assumptions regarding genetic ancestry.

Bight of Benin: a Maternal Perspective of Four Beninese Populations and their Genetic Implications on the American Populations of African Ancestry.

The understanding of the first movements of the ancestral populations within the African continent is still unclear, particularly in West Africa, due to several factors that have shaped the African genetic pool across time. To improve the genetic representativeness of the Beninese population and to better understand the patterns of human settlement inside West Africa and the dynamics of peopling of the Democratic Republic of Benin, we analyzed the maternal genetic variation of 193 Beninese individuals belonging to Bariba, Berba, Dendi, and Fon populations. Results support the oral traditions indicating that the western neighbouring populations have been the ancestors of the first Beninese populations, and the extant genetic structure of the Beninese populations is most likely the result of admixture between populations from neighbouring countries and native people. The present findings highlight how the Beninese populations contributed to the gene pool of the extant populations of some American populations of African ancestry. This strengthens the hypothesis that the Bight of Benin was not only an assembly point for the slave trade during the Trans-Atlantic Slave Trade but also an important slave trapping area.

Genetic-genealogy approach reveals low rate of extrapair paternity in historical Dutch populations.

OBJECTIVE: Evolutionary theory has shown that seeking out extrapair paternity (EPP) can be a viable reproductive strategy for both sexes in pair-bonded species, also in humans. As yet, estimates of the contemporary or historical EPP rate in human population are still rare. In the present study, we estimated the historical EPP rate in the Dutch population over the last 400 years and compared the rate with those obtained for other human populations to determine the evolutionary, cultural, and socio-demographic factors that influence human cuckoldry behavior. METHODS: We estimated the historical EPP rate for the Dutch population via the "genealogical pair method", in which the EPP rate is derived from Y-chromosome mismatches between pairs of individuals that, based on genealogical evidence, share a common paternal ancestor. RESULTS: Based on the analysis of 68 representative genealogical pairs, separated by a total of 1013 fertilization events, we estimated that the historical EPP rate for the Dutch population over the last 400 years was 0.96% per generation (95% confidence interval 0.46%-1.76%). CONCLUSION: The Dutch EPP rate fits perfectly within the range reported for other contemporary and historical populations in Western Europe and was highly congruent with that estimated for neighboring Flanders, despite the socio-economic and religious differences between both populations. The estimated low EPP rate challenges the "dual mating strategy hypothesis" that states that women could obtain fitness benefits by securing investment from one man while cuckolding him to obtain good genes from an affair partner.

Pseudoautosomal region 1 length polymorphism in the human population.

The human sex chromosomes differ in sequence, except for the pseudoautosomal regions (PAR) at the terminus of the short and the long arms, denoted as PAR1 and PAR2. The boundary between PAR1 and the unique X and Y sequences was established during the divergence of the great apes. During a copy number variation screen, we noted a paternally inherited chromosome X duplication in 15 independent families. Subsequent genomic analysis demonstrated that an insertional translocation of X chromosomal sequence into the Y chromosome generates an extended PAR [corrected].The insertion is generated by non-allelic homologous recombination between a 548 bp LTR6B repeat within the Y chromosome PAR1 and a second LTR6B repeat located 105 kb from the PAR boundary on the X chromosome. The identification of the reciprocal deletion on the X chromosome in one family and the occurrence of the variant in different chromosome Y haplogroups demonstrate this is a recurrent genomic rearrangement in the human population. This finding represents a novel mechanism shaping sex chromosomal evolution.

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.

Toward male individualization with rapidly mutating y-chromosomal short tandem repeats.

Relevant for various areas of human genetics, Y-chromosomal short tandem repeats (Y-STRs) are commonly used for testing close paternal relationships among individuals and populations, and for male lineage identification. However, even the widely used 17-loci Yfiler set cannot resolve individuals and populations completely. Here, 52 centers generated quality-controlled data of 13 rapidly mutating (RM) Y-STRs in 14,644 related and unrelated males from 111 worldwide populations. Strikingly, >99% of the 12,272 unrelated males were completely individualized. Haplotype diversity was extremely high (global: 0.9999985, regional: 0.99836-0.9999988). Haplotype sharing between populations was almost absent except for six (0.05%) of the 12,156 haplotypes. Haplotype sharing within populations was generally rare (0.8% nonunique haplotypes), significantly lower in urban (0.9%) than rural (2.1%) and highest in endogamous groups (14.3%). Analysis of molecular variance revealed 99.98% of variation within populations, 0.018% among populations within groups, and 0.002% among groups. Of the 2,372 newly and 156 previously typed male relative pairs, 29% were differentiated including 27% of the 2,378 father-son pairs. Relative to Yfiler, haplotype diversity was increased in 86% of the populations tested and overall male relative differentiation was raised by 23.5%. Our study demonstrates the value of RM Y-STRs in identifying and separating unrelated and related males and provides a reference database.

Pig domestication and human-mediated dispersal in western Eurasia revealed through ancient DNA and geometric morphometrics.

Zooarcheological evidence suggests that pigs were domesticated in Southwest Asia ~8,500 BC. They then spread across the Middle and Near East and westward into Europe alongside early agriculturalists. European pigs were either domesticated independently or more likely appeared so as a result of admixture between introduced pigs and European wild boar. As a result, European wild boar mtDNA lineages replaced Near Eastern/Anatolian mtDNA signatures in Europe and subsequently replaced indigenous domestic pig lineages in Anatolia. The specific details of these processes, however, remain unknown. To address questions related to early pig domestication, dispersal, and turnover in the Near East, we analyzed ancient mitochondrial DNA and dental geometric morphometric variation in 393 ancient pig specimens representing 48 archeological sites (from the Pre-Pottery Neolithic to the Medieval period) from Armenia, Cyprus, Georgia, Iran, Syria, and Turkey. Our results reveal the first genetic signatures of early domestic pigs in the Near Eastern Neolithic core zone. We also demonstrate that these early pigs differed genetically from those in western Anatolia that were introduced to Europe during the Neolithic expansion. In addition, we present a significantly more refined chronology for the introduction of European domestic pigs into Asia Minor that took place during the Bronze Age, at least 900 years earlier than previously detected. By the 5th century AD, European signatures completely replaced the endemic lineages possibly coinciding with the widespread demographic and societal changes that occurred during the Anatolian Bronze and Iron Ages.

Recent radiation within Y-chromosomal haplogroup R-M269 resulted in high Y-STR haplotype resemblance.

Y-chromosomal short tandem repeats (Y-STRs) are often used in addition to Y-chromosomal single-nucleotide polymorphisms (Y-SNP) to detect subtle patterns in a population genetic structure. There are, however, indications for Y-STR haplotype resemblance across different subhaplogroups within haplogroup R1b1b2 (R-M269) which may lead to erosion in the observation of the population genetic pattern. Hence the question arises whether Y-STR haplotypes are still informative beyond high-resolution Y-SNP genotyping for population genetic studies. To address this question, we genotyped the Y chromosomes of more than 1000 males originating from the West-European regions of Flanders (Belgium), North-Brabant and Limburg (the Netherlands) at the highest resolution of the current Y-SNP tree together with 38 commonly used Y-STRs. We observed high resemblance of Y-STR haplotypes between males belonging to different subhaplogroups of haplogroup R-M269. Several subhaplogroups within R-M269 could not be distinguished from each other based on differences in Y-STR haplotype variation. The most likely hypothesis to explain this similarity of Y-STR haplotypes within the population of R-M269 members is a recent radiation where various subhaplogroups originated within a relatively short time period. We conclude that high-resolution Y-SNP typing rather than Y-STR typing might be more useful to study population genetic patterns in (Western) Europe.

In silico detection of phylogenetic informative Y-chromosomal single nucleotide polymorphisms from whole genome sequencing data.

A state-of-the-art phylogeny of the human Y-chromosome is an essential tool for forensic genetics. The explosion of whole genome sequencing (WGS) data due to the rapid progress of next-generation sequencing facilities is useful to optimize and to increase the resolution of the phylogenetic Y-chromosomal tree. The most interesting Y-chromosomal variants to increase the phylogeny are SNPs (Y-SNPs) especially since the software to call them in WGS data and to genotype them in forensic assays has been optimized over the past years. The PENNY software presented here detects potentially phylogenetic interesting Y-SNPs in silico based on SNP calling data files and classifies them into different types according to their position in the currently used Y-chromosomal tree. The software utilized 790 available male WGS samples of which 172 had a high SNP calling quality. In total, 1269 Y-SNPs potentially capable of increasing the resolution of the Y-chromosomal phylogenetic tree were detected based on a first run with PENNY. Based on a test panel of 57 high-quality and 618 low-quality WGS samples, we could prove that these newly added Y-SNPs indeed increased the resolution of the phylogenetic Y-chromosomal analysis substantially. Finally, we performed a second run with PENNY whereby all samples including those of the test panel are used and this resulted in 509 additional phylogenetic promising Y-SNPs. By including these additional Y-SNPs, a final update of the present phylogenetic Y-chromosomal tree which is useful for forensic applications was generated. In order to find more convincing forensic interesting Y-SNPs with this PENNY software, the number of samples and variety of the haplogroups to which these samples belong needs to increase. The PENNY software (inclusive the user manual) is freely available on the website http://bio.kuleuven.be/eeb/lbeg/software.

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.

Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications.

Nonrecombining Y-chromosomal microsatellites (Y-STRs) are widely used to infer population histories, discover genealogical relationships, and identify males for criminal justice purposes. Although a key requirement for their application is reliable mutability knowledge, empirical data are only available for a small number of Y-STRs thus far. To rectify this, we analyzed a large number of 186 Y-STR markers in nearly 2000 DNA-confirmed father-son pairs, covering an overall number of 352,999 meiotic transfers. Following confirmation by DNA sequence analysis, the retrieved mutation data were modeled via a Bayesian approach, resulting in mutation rates from 3.78 × 10(-4) (95% credible interval [CI], 1.38 × 10(-5) - 2.02 × 10(-3)) to 7.44 × 10(-2) (95% CI, 6.51 × 10(-2) - 9.09 × 10(-2)) per marker per generation. With the 924 mutations at 120 Y-STR markers, a nonsignificant excess of repeat losses versus gains (1.16:1), as well as a strong and significant excess of single-repeat versus multirepeat changes (25.23:1), was observed. Although the total repeat number influenced Y-STR locus mutability most strongly, repeat complexity, the length in base pairs of the repeated motif, and the father's age also contributed to Y-STR mutability. To exemplify how to practically utilize this knowledge, we analyzed the 13 most mutable Y-STRs in an independent sample set and empirically proved their suitability for distinguishing close and distantly related males. This finding is expected to revolutionize Y-chromosomal applications in forensic biology, from previous male lineage differentiation toward future male individual identification.

Deep into the roots of the Libyan Tuareg: a genetic survey of their paternal heritage.

Recent genetic studies of the Tuareg have begun to uncover the origin of this semi-nomadic northwest African people and their relationship with African populations. For centuries they were caravan traders plying the trade routes between the Mediterranean coast and south-Saharan Africa. Their origin most likely coincides with the fall of the Garamantes who inhabited the Fezzan (Libya) between the 1st millennium BC and the 5th century AD. In this study we report novel data on the Y-chromosome variation in the Libyan Tuareg from Al Awaynat and Tahala, two villages in Fezzan, whose maternal genetic pool was previously characterized. High-resolution investigation of 37 Y-chromosome STR loci and analysis of 35 bi-allelic markers in 47 individuals revealed a predominant northwest African component (E-M81, haplogroup E1b1b1b) which likely originated in the second half of the Holocene in the same ancestral population that contributed to the maternal pool of the Libyan Tuareg. A significant paternal contribution from south-Saharan Africa (E-U175, haplogroup E1b1a8) was also detected, which may likely be due to recent secondary introduction, possibly through slavery practices or fusion between different tribal groups. The difference in haplogroup composition between the villages of Al Awaynat and Tahala suggests that founder effects and drift played a significant role in shaping the genetic pool of the Libyan Tuareg.

At the Interface of Life and Death: Post-mortem and Other Applications of Vaginal, Skin, and Salivary Microbiome Analysis in Forensics.

Microbial forensics represents a promising tool to strengthen traditional forensic investigative methods and fill related knowledge gaps. Large-scale microbiome studies indicate that microbial fingerprinting can assist forensics in areas such as trace evidence, source tracking, geolocation, and circumstances of death. Nevertheless, the majority of forensic microbiome studies focus on soil and internal organ samples, whereas the microbiome of skin, mouth, and especially vaginal samples that are routinely collected in sexual assault and femicide cases remain underexplored. This review discusses the current and emerging insights into vaginal, skin, and salivary microbiome-modulating factors during life (e.g., lifestyle and health status) and after death (e.g., environmental influences and post-mortem interval) based on next-generation sequencing. We specifically highlight the key aspects of female reproductive tract, skin, and mouth microbiome samples relevant in forensics. To fill the current knowledge gaps, future research should focus on the degree to which the post-mortem succession rate and profiles of vaginal, skin, and saliva microbiota are sensitive to abiotic and biotic factors, presence or absence of oxygen and other gases, and the nutrient richness of the environment. Application of this microbiome-related knowledge could provide valuable complementary data to strengthen forensic cases, for example, to shed light on the circumstances surrounding death with (post-mortem) microbial fingerprinting. Overall, this review synthesizes the present knowledge and aims to provide a framework to adequately comprehend the hurdles and potential application of vaginal, skin, and salivary post-mortem microbiomes in forensic investigations.

Cell survival and DNA damage repair are promoted in the human blood thanatotranscriptome shortly after death.

RNA analysis of post-mortem tissues, or thanatotranscriptomics, has become a topic of interest in forensic science due to the essential information it can provide in forensic investigations. Several studies have previously investigated the effect of death on gene transcription, but it has never been conducted with samples of the same individual. For the first time, a longitudinal mRNA expression analysis study was performed with post-mortem human blood samples from individuals with a known time of death. The results reveal that, after death, two clearly differentiated groups of up- and down-regulated genes can be detected. Pathway analysis suggests active processes that promote cell survival and DNA damage repair, rather than passive degradation, are the source of early post-mortem changes of gene expression in blood. In addition, a generalized linear model with an elastic net restriction predicted post-mortem interval with a root mean square error of 4.75 h. In conclusion, we demonstrate that post-mortem gene expression data can be used as biomarkers to estimate the post-mortem interval though further validation using independent sample sets is required before use in forensic casework.

Ysurnames? The patrilineal Y-chromosome and surname correlation for DNA kinship research.

The Y-chromosome is a widely studied and useful small part of the genome providing different applications for interdisciplinary research. In many (Western) societies, the Y-chromosome and surnames are paternally co-inherited, suggesting a corresponding Y-haplotype for every namesake. While it has already been observed that this correlation may be disrupted by a false-paternity event, adoption, anonymous sperm donor or the co-founding of surnames, extensive information on the strength of the surname match frequency (SMF) with the Y-chromosome remains rather unknown. For the first time in Belgium and the Netherlands, we were able to study this correlation using 2,401 males genotyped for 46 Y-STRs and 183 Y-SNPs. The SMF was observed to be dependent on the number of Y-STRs analyzed, their mutation rates and the number of Y-STR differences allowed for a kinship. For a perfect match, the Yfiler® Plus and our in-house YForGen kit gave a similar high SMF of 98%, but for non-perfect matches, the latter could overall be identified as the best kit. The SMF generally increased due to less mismatches when encountering [1] deep Y-subhaplogroups, [2] less frequently occurring surnames, and [3] small geographical distances between relatives. This novel information enabled the design of a surname prediction model based on genetic and geographical distances of a kinship. The prediction model has an area under the curve (AUC) of 0.9 and is therefore useable for DNA kinship priority listing in estimation applications like forensic familial searching.

The Dutch Y-chromosomal landscape.

Previous studies indicated existing, albeit limited, genetic-geographic population substructure in the Dutch population based on genome-wide data and a lack of this for mitochondrial SNP based data. Despite the aforementioned studies, Y-chromosomal SNP data from the Netherlands remain scarce and do not cover the territory of the Netherlands well enough to allow a reliable investigation of genetic-geographic population substructure. Here we provide the first substantial dataset of detailed spatial Y-chromosomal haplogroup information in 2085 males collected across the Netherlands and supplemented with previously published data from northern Belgium. We found Y-chromosomal evidence for genetic-geographic population substructure, and several Y-haplogroups demonstrating significant clinal frequency distributions in different directions. By means of prediction surface maps we could visualize (complex) distribution patterns of individual Y-haplogroups in detail. These results highlight the value of a micro-geographic approach and are of great use for forensic and epidemiological investigations and our understanding of the Dutch population history. Moreover, the previously noted absence of genetic-geographic population substructure in the Netherlands based on mitochondrial DNA in contrast to our Y-chromosome results, hints at different population histories for women and men in the Netherlands.

Correction to: The Dutch Y-chromosomal landscape.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A game of hide and seq: Identification of parallel Y-STR evolution in deep-rooting pedigrees.

Short tandem repeats on the Y-chromosome (Y-STRs) are common DNA polymorphisms useful for genetic genealogy, population and evolutionary genetics, human genetics, pathology and forensic sciences. It is important to identify all Y-STR variants and to have knowledge of Y-STR mutation rates in order to correctly estimate the time to the most recent common ancestor (tMRCA) between paternally related individuals. When capillary electrophoresis (CE) is performed to analyze genealogical pairs, Y-STR sequence variations remain hidden when the number of repeats is identical. These hidden variations could be due to parallel Y-STR changes or modifications (PM) that occur independently in different lineages leading to alleles with identical number of repeats. In this study, we detect for the first time twelve PM by analyzing 133 males (960 meiosis) in extended deep-rooting family pedigrees on 42 Y-STRs. These PM were observed in nine Y-STR loci with mutation rates of at least 5.94 × 10-3 per generation. Sequencing analysis made it possible to distinguish insertions/deletions in different repeat regions revealing the presence of two unique changes in three PM on rapidly mutating and complex Y-STRs DYS724-ab and DYS518. Sequencing unraveled more information concerning the identity of alleles, and increased allelic discrimination possibilities which is of great importance in population genetics and forensic analysis. Limiting the analysis to CE could lead to wrong ancestral allele assumptions, to false negative interpretations and to tMRCA underestimations. These observations highlight the importance and added value of sequencing analysis and suggest a shift in genotyping methods from CE to next generation sequencing.