Determining Y-STR mutation rates in deep-routing genealogies: Identification of haplogroup differences

Evaluation of RMplex system for differentiating father-son pairs using Y-STRs in a Korean population

Y-chromosomal short tandem repeats (Y-STRs) at rapidly mutating (RM) loci have been suggested as tools for differentiating paternally related males. RMplex is a recently developed system that incorporates 26 RM loci and four fast-mutating (FM) loci, targeting 44 male-specific loci. Here, we evaluated the RMplex by estimating Y-STR mutation rates and the overall differentiation rates for 542 Korean father-son pairs, as well as the genetic population values for 409 unrelated males. RMplex performed well, distinguishing 50.7 % of the father-son pairs by at least one mutation, a value 10 times higher than the previously reported differentiation rate achieved using the PowerPlex® Y23 System. Of the 369 mutations, 361 (97.8 %) were single-step mutations, with locus-specific mutation rates varying from 1.8 × 10-3 to 1.1 × 10-1 mutations per generation, and an average mutation rate of 2.3 × 10-2. Gene diversity values ranged from 0.5696 for DYS442 to 0.9970 for DYF1000, and the haplotype discrimination capacity of unrelated males was 100 %. Among the loci studied, DYS712 exhibited the highest mutation rate in this study of the Korean population. Similarly, the mutation rate of this locus is reported to be substantially higher for the Japanese and Chinese populations than for European populations. These findings suggest that DYS712 mutations are relatively frequent in East Asian populations. Although we did not detect significant relationships among the Y-chromosome single nucleotide polymorphism-based haplogroups, allele length was strongly correlated with the mutation rate at DYS712, which is consistent with previous studies. Although the incorporation of multi-copy loci into RMplex contributed significantly to the high mutation rates detected and to its discrimination capacity, this requires careful interpretation, owing to the potential for duplications. Nonetheless, these findings provide evidence regarding the suitability of the RMplex for distinguishing paternally related males in the Korean population.

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Opt-in or out? Public perspectives on forensic DNA kinship investigations within the Dutch-speaking community

Forensic DNA kinship investigation involves analyzing genetic relationships between individuals to offer new leads for solving (cold) cases. Familial DNA matching has become a valuable asset in criminal case investigations, especially when traditional DNA methods hit dead ends. However, concerns surrounding ethical and privacy implications raised questions about its implementation and acceptance among the general public. The present study investigated the public perspectives regarding forensic DNA kinship investigations among 1710 Dutch-speaking Belgians using an online cross-sectional survey. The questionnaire consisted of three categories, including personal information, DNA knowledge, and their opinion on several familial DNA searching and investigative genetic genealogy related questions. The participants' average DNA knowledge score was 71 %, indicating a relatively high level of understanding of DNA-related concepts. Remarkably, the study revealed that 92 % of the participants expressed willingness to cooperate as a volunteer in a forensic DNA kinship investigation, irrespective of their scientific background or educational level. Key factors influencing participation included assurance of painless sampling and robust privacy safeguards. Participants lacking familiarity with DNA hesitated more towards participating in forensic DNA analysis, referring to "the fear of the unknown". Despite ethical and privacy concerns, the highly positive attitude towards forensic DNA analysis reflects a level of empathy and willingness to contribute to the pursuit of justice. Nearly all participants (95 %) agreed to use online DNA databases for resolving violent crimes with forensic genetic genealogy, but half emphasized the need for prior informed consent, referring to the current "opt-in" system. The results underscore the need for stringent regulations and ethical oversight to ensure the responsible use of genetic data while striking a balance between public safety and the protection of individuals' privacy rights. These findings add to the growing body of evidence regarding the potential benefits of forensic DNA kinship matching as a tool in criminal investigations, suggesting its potential future utilization and legalization.

Y chromosome haplogroup N in a Japanese population is classified into three subclades, and two DYS385 loci, a duplicated Y-STR, are duplicated again in subclade N-M1819

DYS385 is one of the major Y chromosome short tandem repeats (Y-STRs) in forensic genetics and exists as 2 copies in the human Y chromosome palindrome P4 region. In this study, we found that some samples were estimated to have ≥ 4 copies of DYS385 in Y chromosome haplogroup N in a Japanese population. Y chromosome haplogroup N is distributed widely in eastern/central Asia, Siberia, and eastern/northern Europe, and is also observed in Japan; however, little is known about haplogroup N subclades in the Japanese population. To reveal the link between increased DYS385 copy number and haplogroup N subclades, we sequenced single nucleotide polymorphisms to classify the subclades. As a result, the Japanese Y chromosomes of haplogroup N were classified into three subclades, and an increased DYS385 copy number was specific to subclade N-M1819* (N1b2*). These results are of use in forensic DNA analysis because Y-STR copy number is important for the interpretation of Y-STR typing results of male DNA mixtures and kinship analysis. We also found that DYS458.1 microvariants (DYS458 intermediate alleles with single-nucleotide insertion) were observed only in subclade N-CTS962 (N1b1b∼) samples. Given that previous studies reported that DYS458.1 microvariants are observed in Y chromosomes of haplogroup N in other populations, DYS458.1 might be used to infer haplogroup N subclades without limitation to the Japanese population. The results of this study will be beneficial not only to forensic genetics but also to anthropological studies.

Microsatellites' mutation modeling through the analysis of the Y-chromosomal transmission: Results of a GHEP-ISFG collaborative study

The Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) organized a collaborative study on mutations of Y-chromosomal short tandem repeats (Y-STRs). New data from 2225 father-son duos and data from 44 previously published reports, corresponding to 25,729 duos, were collected and analyzed. Marker-specific mutation rates were estimated for 33 Y-STRs. Although highly dependent on the analyzed marker, mutations compatible with the gain or loss of a single repeat were 23.2 times more likely than those involving a greater number of repeats. Longer alleles (relatively to the modal one) showed to be nearly twice more mutable than the shorter ones. Within the subset of longer alleles, the loss of repeats showed to be nearly twice more likely than the gain. Conversely, shorter alleles showed a symmetrical trend, with repeat gains being twofold more frequent than reductions. A positive correlation between the paternal age and the mutation rate was observed, strengthening previous findings. The results of a machine learning approach, via logistic regression analyses, allowed the establishment of algebraic formulas for estimating the probability of mutation depending on paternal age and allele length for DYS389I, DYS393 and DYS627. Algebraic formulas could also be established considering only the allele length as predictor for DYS19, DYS389I, DYS389II-I, DYS390, DYS391, DYS393, DYS437, DYS439, DYS449, DYS456, DYS458, DYS460, DYS481, DYS518, DYS533, DYS576, DYS626 and DYS627 loci. For the remaining Y-STRs, a lack of statistical significance was observed, probably as a consequence of the small effective size of the subsets available, a common difficulty in the modeling of rare events as is the case of mutations. The amount of data used in the different analyses varied widely, depending on how the data were reported in the publications analyzed. This shows a regrettable waste of produced data, due to inadequate communication of the results, supporting an urgent need of publication guidelines for mutation studies.

Considerations on the application of a mutation model for Y-STR interpretation

If Y-STR profiling is to be more effective in criminal casework, the methods used to evaluate evidential weight require improvement. Many forensic scientists assign an evidential weight by estimating the number of times a Y-STR profile obtained from a questioned sample has been observed in YHRD datasets. More sophisticated models have been suggested but not yet implemented into routine casework, e.g. Andersen & Balding [1]. Mutation is inherent to STR meiosis (or inheritance) and is encountered in practice. We evaluated a mutation model that can be incorporated into a method for assigning evidential weight to Y-STR profiles, an essential part of bringing any method into practice. Since an important part of implementation to casework is communication, the article is written in an accessible format for practitioners as well as statisticians. The mutation component within the MUTEA model by Willems et al. [2] incorporates the potential for multistep mutations and a tendency for alleles to revert towards a central length, reflecting observed mutation data, e.g. [3]. We have estimated the parameters in this model and in a simplified symmetric version of this model, using sequence data from father/son pairs [4] and deep-rooted pedigrees [5]. Both datasets contain multistep mutations, which may have an effect on models based on simulations [1]. We introduce Beta-Binomial and Beta-Geometric conjugate analyses for estimating rate and step parameters for the mutation models presented here, which require only summations and multiplications. We proved mathematically that the parameters can be estimated independently. We show the importance of reporting the variability of the parameters and not only a point estimate. The parameters can be easily incorporated into statistical models, and updated sequentially as more data becomes available. We recommend fuller publication of data to enable the development and evaluation of a wider range of mutation models.

Improving the strategy to identify historical military remains: a literature review and Y-STR meta-analysis

The identification of historical military remains by Unrecovered War Casualties-Army (UWC-A) currently relies on Y-chromosome Short Tandem Repeat (Y-STR) testing when maternal relatives are not available, or when a mitochondrial DNA match does not provide sufficient certainty of identification. However, common Y-STR profiles (using Yfiler™) between sets of remains or families often prevent identification. To resolve these cases, an investigation of additional Y-DNA markers is needed for their potential inclusion into the DNA identification strategy. The number of genetic transmissions between missing soldiers and their living relatives needs to be considered to avoid false exclusions between paternal relatives. Analysis of 236 World War I/II (WWI/II) era pairs of relatives identified up to seven genetic transmissions between WWII soldiers and their living relatives, and nine for WWI. Previous Y-STR meta-analyses were published approximately 10 years ago when rapidly mutating markers were relatively new. This paper reports a contemporary literature review and meta-analysis of 35 studies (which includes 23 studies not previously used in meta-analysis) and 23 commonly used Y-STR's mutation rates to inform the inclusion of additional loci to UWC-A's DNA identification strategy. Meta-analysis found mutation data for a given Y-STR locus could be pooled between studies and that the mutation rates were significantly different between some loci (at P < 0.05). Based on this meta-analysis, we have identified two additional markers from PowerPlex® Y23 for potential inclusion in UWC-A's identification strategy. Further avenues for potential experimental exploration are discussed.

Key points

From 236 UWC-A pairs of relatives, we observed up to nine genetic transmissions between WWI soldiers and their living relatives, and seven for WWII.MedCalc® software for meta-analysis utilizing the Freeman-Tukey transformation was run, which analysed 35 published studies and 23 commonly used loci. Previous Y-STR mutation rate meta-analyses are now 10 years old; this paper includes 23 studies that were not included in previous meta-analyses.Through meta-analysis, we identify two markers from PowerPlex® Y23 for potential inclusion in UWC-A's historical remains identification strategy (alongside Yfiler™). We discuss potential next steps for experimental exploration of additional Y-DNA markers.

Genetic Analysis of 27 Y-STR Haplotypes in 11 Iranian Ethnic Groups

BACKGROUND

The study of Y-chromosomal variations provides valuable insights into male susceptibility in certain diseases like cardiovascular disease (CVD). In this study, we analyzed paternal lineage in different Iranian ethnic groups, not only to identify developing medical etiology, but also to pave the way for gender-specific targeted strategies and personalized medicine in medical genetic research studies.

METHODS

The diversity of eleven Iranian ethnic groups was studied using 27 Y-chromosomal short tandem repeat (Y-STR) haplotypes from Y-filer® Plus kit. Analysis of molecular variance (AMOVA) based on pair-wise RST along with multidimensional scaling (MDS) calculation and Network phylogenic analysis was employed to quantify the differences between 503 unrelated individuals from each ethnicity.

RESULTS

Results from AMOVA calculation confirmed that Gilaks and Azeris showed the largest genetic distance (RST=0.35434); however, Sistanis and Lurs had the smallest considerable genetic distance (RST=0.00483) compared to other ethnicities. Although Azeris had a considerable distance from other ethnicities, they were still close to Turkmens. MDS analysis of ethnic groups gave the indication of lack of similarity between different ethnicities. Besides, network phylogenic analysis demonstrated insignificant clustering between samples.

CONCLUSION

The AMOVA analysis results explain that the close distance of Azeris and Turkmens may be the effect of male-dominant expansions across Central Asia that contributed to historical and demographics of populations in the region. Insignificant differences in network analysis could be the consequence of high mutation events that happened in the Y-STR regions over the years. Considering the ethnic group affiliations in medical research, our results provided an understanding and characterization of Iranian male population for future medical and population genetics studies.

Sequence-based mutation patterns at 41 Y chromosomal STRs in 2 548 father-son pairs

A total of 2 548 unrelated healthy father-son pairs from a Northern Han Chinese population were genotyped at 41 Y chromosomal short tandem repeat (Y-STRs) including DYS19, DYS388, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS444, DYS447, DYS448, DYS449, DYS456, DYS458, DYS460, DYS481, DYS518, DYS522, DYS549, DYS533, DYS557, DYS570, DYS576, DYS593, DYS596, DYS627, DYS635, DYS643, DYS645, Y-GATA-H4, DYF387S1a/b, DYF404S1a/b, DYS385a/b, and DYS527a/b. In 2 548 father samples, 2 387 unique haplotypes were detected with the haplotype diversity and discrimination capacity values of 0.999 956 608 and 0.96 741 007. The average gene diversity (GD) value was 0.6934 with a range from 0.1051 at DYS645 to 0.9657 at DYS385a/b. When comparing alleles at 24 overlapped Y-STRs between the ForenSeq™ deoxyribonucleic acid (DNA) Signature Prep Kit on the MiSeq FGx® Forensic Genomics System and the Goldeneye® DNA ID Y Plus Kit on the Applied Biosystems™ 3730 DNA Analyzer from 308 father samples in mutational pairs, 258 alleles were detected by massively parallel sequencing (MPS) typing including 156 length-based alleles that could be obtained by capillary electrophoresis (CE) typing, 95 repeat region (RR) variant alleles and seven flanking region variant alleles. Hereof, we found 16 novel RR variant alleles and firstly identified two SNPs (rs2016239814 at DYS19 and rs2089968964 at DYS448) and one 4-bp deletion (rs2053269960 at DYS439) that had been validated by the Database of Short Genetic Variation. Sanger sequencing or MPS was employed to confirm 356 mutations from 104 468 allele transfers generated from CE, where 96.63% resulted in one-step mutations, 2.25% in two-step, and 1.12% in multi-step, and the overall ratio of repeat gains versus losses was balanced (173 gains vs. 183 losses). In 308 father-son pairs, 268 pairs occurred mutations at a single locus, 33 pairs at two loci, six pairs at three loci, and one pair at four loci. The average Y-STR mutation rate at 41 Y-STRs was ⁓3.4 × 10-3 (95% confidence intervals: 3.1 × 10-3-3.8 × 10-3). The mutation rates at DYS576 and DYS627 were higher than 1 × 10-2 in Northern Han Chinese, whilst the mutation rates at DYF387S1a/b, DYF404S1a/b, DYS449, DYS518, and DYS570 were lower than initially defined. In this study, the classical molecular factors (the longer STR region, the more complex motif and the order father) were confirmed to drive Y-STR mutation rates increased, but the length of repeat unit did not conform to the convention. Lastly, the interactive graphical and installable StatsY was developed to facilitate forensic scientists to automatically calculate allele and haplotype frequencies, forensic parameters, and mutation rates at Y-STRs.

Key points

308 of 2 548 father-son pairs from Northern Han Chinese occurred at least one mutation(s) across 41 Y-STRs.Sanger sequencing or MPS was employed to confirm those mutations generated from CE.The longer STR region, the more complex motif and the order father drove Y-STR mutation rates increased.StatsY was developed to calculate allele and haplotype frequencies, forensic parameters and mutation rates at Y-STRs.

Mutation analysis for newly suggested 30 Y-STR loci with high mutation rates in Chinese father-son pairs

Rapidly mutating Y-STRs (RM Y-STRs) harbor great potential to distinguish male relatives and achieve male identification. However, forensic applications were greatly limited by the small number of the initially identified 14 RM Y-STRs. Recently, with the emergence of 12 novel RM Y-STRs, an integrated panel named RMplex was introduced, which contains all 26 RM Y-STRs and four fast mutating Y-STRs (FM Y-STRs). To obtain the first data on the mutation rates and father-son differentiation rates of the 30 newly proposed Y-STRs in Chinese populations, we performed an empirical mutation study on 307 DNA-confirmed Chinese paternal pairs. Previously reported mutation rates for 14 RM Y-STRs in Chinese and European populations were pooled and merged with our data. The highest meiosis number for the two groups reached 4771 and 2687, respectively. Five loci showed significant differences between the populations (DYS570, DYS399S1, DYS547, DYS612, and DYF403S1b). For the new panel covering 30 Y-STR loci, our results show extensive differences in the mutation rates between the two populations, as well. 10 RM Y-STR loci showed relatively low mutation rates (10^-3-10^-2 per meiosis) and 2 FM Y-STR loci had rapid mutation rates (> 10^-2 per meiosis) in the Chinese population. Several-fold differences in mutation rates were found in nine Y-STR loci between the Chinese and reference populations, with two loci having significantly higher mutation rates and one locus with a significantly lower mutation rate in the Chinese population (P < 0.05). Eighteen RM Y-STRs (> 10^-2 per meiosis), 8 FM Y-STR loci (5×10^-3-10^-2 per meiosis), 3 moderately mutating Y-STRs (MM Y-STRs, 10^-3-5×10^-3 per meiosis), and one locus with no observed mutation events were identified in the Chinese population. 40.06% of the Chinese paternity pairs were discriminated with RMplex while only 20.84% with the initial 14 RM Y-STRs, indicating that RMplex is beneficial for distinguishing paternally related males. Future studies on populations of different genetic backgrounds are necessary to obtain comprehensive estimates of mutation rates at these new loci.

Assessing the factors influencing the performance of machine learning for classifying haplogroups from Y-STR haplotypes

Two distinct genetic markers, single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs), exist simultaneously in the non-recombining portion of the Y chromosome. Because of their different rates of mutation, Y-STRs and Y-SNPs play distinct roles in forensic and evolutionary genetics. Current approaches to infer haplogroup status rely on genotyping lots of Y-SNP loci. Given the relationship between haplotype and haplogroup of a Y chromosome, a cost-effective strategy of Y-STRs typing had an advantage in haplogroup prediction. Many machine learning algorithms have sprung up for assigning a Y-STR haplotype to a haplogroup. However, a series of issues must be solved before the using of machine learning method in practice. Thus, the k-nearest neighbor (kNN) classifier was built respectively based on different situations in this study. We assessed different factors which may influence the performance of the kNN prediction model for classifying haplogroups. The training set was based on a diverse ground-truth data set comprising Y-STR haplotypes and corresponding Y-SNP haplogroups. Our results showed that combining different levels of haplogroups into the observations or transracial prediction was impractical. Moreover, using more slow mutation Y-STR loci in the category is good for promoting classification accuracy. The preconditions for an effective and accurate haplogroup assignment by the kNN classifier were revealed.

RMplex reveals population differences in RM Y-STR mutation rates and provides improved father-son differentiation in Japanese

Rapidly mutating Y chromosomal short tandem repeat markers (RM Y-STRs) -characterized by at least one mutation per 100 generations- are suitable for differentiating both related and unrelated males. The recently introduced multiplex method RMplex allows for the efficient analysis of 30 Y-STRs with increased mutation rates, including all 26 currently known RM Y-STRs. While currently available RM Y-STR mutation rates were established mostly from European individuals, here we applied RMplex to DNA samples of 178 genetically confirmed father-son pairs from East Asia. For several Y-STRs, we found significantly higher mutation rates in Japanese compared to previous estimates. The consequent father-son differentiation rate based on RMplex was significantly higher (52%) in Japanese than previously reported for Europeans (42%), and much higher than with Yfiler Plus in both sample sets (14% and 13%, respectively). Further analysis suggests that the higher mutation and relative differentiation rates in Japanese can in part be explained by on average longer Y-STR alleles relative to Europeans. Moreover, we show that the most striking difference, which was found in DYS712, could be linked to a Y-SNP haplogroup (O1b2-P49) that is common in Japanese and rare in other populations. We encourage the forensic Y-STR community to generate more RMplex data from more population samples of sufficiently large sample size in combination with Y-SNP data to further investigate population effects on mutation and relative differentiation rates. Until more RMplex data from more populations become available, caution shall be placed when applying RM Y-STR mutation rate estimates established in one population, such as Europeans, to forensic casework involving male suspects of paternal origin from other populations, such as non-Europeans.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

Empirical Evidence on Enhanced Mutation Rates of 19 RM-YSTRs for Differentiating Paternal Lineages

Rapidly mutating Y-chromosomal short tandem repeats (RM Y STRs) with mutation rates ≥ 10⁻² per locus per generation are valuable for differentiating amongst male paternal relatives where standard Y STRs with mutation rates of ≤10⁻³ per locus per generation may not. Although the 13 RM Y STRs commonly found in commercial assays provide higher levels of paternal lineage differentiation than conventional Y STRs, there are many male paternal relatives that still cannot be differentiated. This can be improved by increasing the number of Y STRs or choosing those with high mutation rates. We present a RM Y STR multiplex comprising 19 loci with high mutation rates and its developmental validation (repeatability, sensitivity and male specificity). The multiplex was found to be robust, reproducible, specific and sensitive enough to generate DNA profiles from samples with inhibitors. It was also able to detect all contributor alleles of mixtures in ratios up to 9:1. We provide preliminary evidence for the ability of the multiplex to discriminate between male paternal relatives by analyzing large numbers of male relative pairs (536) separated by one to seven meioses. A total of 96 mutations were observed in 162 meioses of father–son pairs, and other closely related male pairs were able to be differentiated after 1, 2, 3, 4, 5, 6 and 7 meiosis in 44%, 69%, 68%, 85%, 0%, 100% and 100% of cases, respectively. The multiplex offers a noticeable enhancement in the ability to differentiate paternally related males compared with the 13 RM Y STR set. We envision the future application of our 19 RM Yplex in criminal cases for the exclusion of male relatives possessing matching standard Y STR profiles and in familial searching with unknown suspects. It represents a step towards the complete individualization of closely related males.

On the Forensic Use of Y-Chromosome Polymorphisms

Nowadays, the use of Y-chromosome polymorphisms forms an essential part of many forensic DNA investigations. However, this was not always the case. Only since 1992 have we seen that some forensic scientists started to have an interest in this chromosome. In this review, I will sketch a brief history focusing on the forensic use of Y-chromosome polymorphisms. Before describing the various applications of short-tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) on the Y-chromosome, I will discuss a few often ignored aspects influencing proper use and interpretation of Y-chromosome information: (i) genotyping Y-SNPs and Y-STRs, (ii) Y-STR haplotypes shared identical by state (IBS) or identical by descent (IBD), and (iii) Y-haplotype database frequencies.

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.

Forensic characteristics and population construction of two major minorities from southwest China revealed by a novel 37 Y-STR loci system

Y-chromosome short tandem repeats (Y-STRs) have become important supplementary evidence in forensic science. Nowadays, the Y-chromosome STR haplotype reference database (YHRD) contains abundant Y-STR haplotype data from all over the world, while haplotype data of Guizhou Miao and Tujia are scarce. Hence, genetic polymorphisms of 37 Y-STRs were investigated in 446 unrelated males (206 Miao males and 246 Tujia males) residing in Guizhou Province. A total of 206 and 242 unique haplotypes with the highest diversity value of 0.9665 and 0.9470 were obtained. The heatmap, multidimensional scaling (MDS), the unweighted pair-group method with arithmetic means (UPGMA) tree and principal component analysis (PCA) based on the genetic distance (Rst) value within our studied populations and other 26 populations indicated that population structures follow the boundary of the continent. Guizhou Miao and Guizhou Tujia populations have intimate relationships with East Asian populations, especially the geographically close, similar history and the same language family populations.

Improved Models of Coalescence Ages of Y-DNA Haplogroups

Databases of commercial DNA-testing companies now contain more customers with sequenced DNA than any completed academic study, leading to growing interest from academic and forensic entities. An important result for both these entities and the test takers themselves is how closely two individuals are related in time, as calculated through one or more molecular clocks. For Y-DNA, existing interpretations of these clocks are insufficiently accurate to usefully measure relatedness in historic times. In this article, I update the methods used to calculate coalescence ages (times to most-recent common ancestor, or TMRCAs) using a new, probabilistic statistical model that includes Y-SNP, Y-STR and ancilliary historical data, and provide examples of its use.

The Y-STR landscape of coastal southeastern Han: Forensic characteristics, haplotype analyses, mutation rates, and population genetics

The Y-STR landscape of Coastal Southeastern Han (CSEH) living in Chinese southeast areas (including Guangdong, Fujian, and Zhejiang provinces) is still unclear. We investigated 62 Y-STR markers in a reasonably large number of 1021 unrelated males and 1027 DNA-confirmed father-son pairs to broaden the genetic backgrounds of CSEH. In total, 85 null alleles, 121 off-ladder alleles, and 95 copy number variants were observed, and 1012 distinct haplotypes were determined with the overall HD and DC values of 0.999974 and 0.9912. We observed 369 mutations in 76 099 meiotic transfers, and the average estimated Y-STR mutation rate was 4.85 × 10^-3 (95% CI, 4.4 × 10^-3 -5.4 × 10^-3 ). The Spearman correlation analyses indicated that GD values (R² = 0.6548) and average allele sizes (R² = 0.5989) have positive correlations with Y-STR mutation rates. Our RM Y-STR set including 8 candidate RM Y-STRs, of which DYS534, DYS630, and DYS713 are new candidates in CSEH, distinguished 18.52% of father-son pairs. This study also clarified the population structures of CSEH which isolated in population-mixed South China relatively. The strategy, SM Y-STRs for familial searching and RM Y-STRs for individual identification regionally, could be applicable based on enough knowledge of the Y-STR mutability of different populations.

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.

Mutation analysis for 25 Y-STR markers in Japanese population

In this study, we analyzed DNA samples from 213 Japanese father son pairs with 25 Y-chromosome short tandem repeat (Y-STR) (DYS576, DYS389I, DYS635, DYS389II, DYS627, DYS460, DYS458, DYS19, YGATAH4, DYS448, DYS391, DYS456, DYS390, DYS438, DYS392, DYS518, DYS570, DYS437, DYS385, DYS449, DYS393, DYS439, DYS481, DYF387S1, and DYS533) markers using the Yfiler™ Plus PCR amplification kit. We calculated Y-STR mutation rates for each locus to evaluate the efficacy of the 25 Y-STR markers for paternity testing and forensic identification using samples from male relatives. Six rapidly mutating Y-STR markers (DYS576, DYS627, DYS518, DYS570, DYS449 and DYF387S1), previously reported to have high mutation rates (>1.0 × 10^-2), are included in the 25 Y-STR markers, but our findings revealed that the mutation rates for all Y-STR markers except for DYS576 and DYS458 were lower than 1.0 × 10^-2. Therefore, the use of these 25 Y-STR markers may be useful for forensic identification in the Japanese population.

Genetic polymorphisms and phylogenetic analyses of the Ü-Tsang Tibetan from Lhasa based on 30 slowly and moderately mutated Y-STR loci

As a result of the expansion of old Tibet on the Qinghai-Tibet Plateau, Tibetans diverged into three main branches, Ü-Tsang, Amdo, and Kham Tibetan. Ü-Tsang Tibetans are geographically distributed across the wide central and western portions of the Qinghai-Tibet Plateau while Lhasa is the central gathering place for Tibetan culture. The AGCU Y30, a 6-dye fluorescence kit including 30 slowly and moderately mutated Y-STR loci, has been validated for its stability and sensitivity in different biomaterials and diverse Chinese populations (Han and other minorities), and widely used in the practical work of forensic science. However, the 30 Y-STR profiling of Tibetan, especially for Ü-Tsang Tibetan, were insufficient. We utilized the AGCU Y30 to genotype 577 Ü-Tsang Tibetan unrelated males from Lhasa in the Tibet Autonomous Region of China to fill up the full and accurate Y-STR profiles. A total of 552 haplotypes were observed, 536 (97.10%) of which were unique. One hundred and ninety-four alleles were observed at 26 single copy loci and the allelic frequencies ranged from 0.0017 to 0.8180. For the two multi-copy loci DYS385a/b and DYS527a/b, 64 and 36 allelic combinations were observed, respectively. The gene diversity (GD) values ranged from 0.3079 at DYS391 to 0.9142 at DYS385a/b and the overall haplotype diversity (HD) was 0.9998, and its discrimination capacity (DC) was 0.9567. The population genetic analyses demonstrated that Lhasa Ü-Tsang Tibetan had close relationships with other Tibetan populations from Tibet and Qinghai, especially with Ü-Tsang Tibetan. From the perspective of Y haplogroups, the admixture of the southward Qiang people with dominant haplogroup O-M122 and the northward migrations of the initial settlers of East Asia with haplogroup D-M175 hinted the Sino-Tibetan homologous, thus, we could not ignore the gene flows with other Sino-Tibetan populations, especially for Han Chinese, to characterize the forensic genetic landscape of Tibetan.

Mutational data and population profiling of 23 Y-STRs in three Brazilian populations

Y-chromosomal STRs are important markers in forensic genetics, due to some peculiar characteristics. The absence of recombination makes them a useful tool to infer kinship in complex cases involving distant paternal relatives, or to infer paternal bio-geographic ancestry. The presence of a single copy, being transmitted from father to son, allow tracing mutational events in Y-STRs without ambiguity. For the statistical interpretation of forensic evidences based on Y-STR profiles, it is necessary to have estimates on both mutation rates and haplotype frequencies. In this work, 407 father-son duos from São Paulo and Rio de Janeiro states and 204 unrelated individuals from Manaus were analyzed. Haplotype frequencies and mutation rates for the Y-STRs from the PowerPlex Y23 commercial kit were estimated. Thirty-six mutations were observed in 15 of the 22 Y-STRs analyzed, for an average mutation rate of 3.84 × 10^-3 (95 % CI 2.69 × 10^-3 to 5.32 × 10^-3). All mutations in GAAA repeats occurred in alleles with 13 or more uninterrupted units. Mutations in GATA repeats were observed in alleles with 9-17 uninterrupted units. An analysis carried out in different father's age groups showed an increase of 2.48 times the mutation rate in the age group of 40-50 years, when compared to the 20-30 age group, in agreement with the described for autosomal STRs. A high haplotype diversity was found in the three Brazilian populations. Pairwise genetic distance analysis (F_ST) showed no significant differences between the three populations in this study, which were also close to populations with strong European influence. The highest distances among the Brazilian populations were with São Gabriel da Cachoeira, which has a high Native American ancestry.

The construction and application of a new 17-plex Y-STR system using universal fluorescent PCR

Y-chromosomal short tandem repeat (Y-STR) polymorphisms are useful in forensic identification, population genetics, and human structures. However, the current Y-STR systems are limited in discriminating distant relatives in a family with a low discrimination power. Increasing the capacity of detecting Y chromosomal polymorphisms will drastically narrow down the matching number of genealogy populations or pedigrees. In this study, we developed a system containing 17 Y-STRs that are complementary to the current commercially available Y-STR kits. This system was constructed by multiplex PCR with expected sizes of 126-400 bp labeled by different fluorescence molecules (DYS715, DYS709, DYS716, DYS713, and DYS607 labeled by FAM; DYS718, DYS723, DYS708, and DYS714 labeled by JOE; DYS712, DYS717, DYS721, and DYS605 labeled by TAMRA; and DYS719, DYS726, DYS598, and DYS722 labeled by ROX). The system was extensively tested for sensitivity, male specificity, species specificity, mixture, population genetics, and mutation rates following the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines. The genetic data were obtained from eight populations with a total of 1260 individuals. Our results showed that all the 17 Y-STRs are human- and male-specific and include only one copy of the Y-chromosome. The 17 Y-STR system detects 143 alleles and has a high discrimination power (0.996031746). Mutation rates were different among the 17 Y-STRs, ranging from 0.30 to 3.03%. In conclusion, our study provides a robust, sensitive, and cost-effective genotyping method for human identification, which will be beneficial for narrowing the search scope when applied to genealogy searching with the Y-STR DNA databank.

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.

Estimating Y-Str Mutation Rates and Tmrca Through Deep-Rooting Italian Pedigrees

In the population genomics era, the study of Y-chromosome variability is still of the greatest interest for several fields ranging from molecular anthropology to forensics and genetic genealogy. In particular, mutation rates of Y-chromosomal Short Tandem Repeats markers (Y-STRs) are key parameters for different interdisciplinary applications. Among them, testing the patrilineal relatedness between individuals and calculating their Time of Most Recent Common Ancestors (TMRCAs) are of the utmost importance. To provide new valuable estimates and to address these issues, we typed 47 Y-STRs (comprising Yfiler, PowerPlex23 and YfilerPlus loci, the recently defined Rapidly Mutating [RM] panel and 11 additional markers often used in genetic genealogical applications) in 135 individuals belonging to 66 deep-rooting paternal genealogies from Northern Italy. Our results confirmed that the genealogy approach is an effective way to obtain reliable Y-STR mutation rate estimates even with a limited number of samples. Moreover, they showed that the impact of multi-step mutations and backmutations is negligible within the temporal scale usually adopted by forensic and genetic genealogy analyses. We then detected a significant association between the number of mutations within genealogies and observed TMRCAs. Therefore, we compared observed and expected TMRCAs by implementing a Bayesian procedure originally designed by Walsh (2001) and showed that the method yields a good performance (up to 96.72%), especially when using the Infinite Alleles Model (IAM).

Relating Clans Ao and Aisin Gioro from northeast China by whole Y-chromosome sequencing

The Y-chromosome haplogroup C2b1a3a2-F8951 is the paternal lineage of the Aisin Gioro clan, the most important brother branch of the famous Mongolic-speaking population characteristic haplogroup C2*-Star Cluster (C2b1a3a1-F3796). However, investigations on its internal phylogeny are still limited. In this study, we used whole Y-chromosome sequencing to update its phylogenetic tree. In the revised tree, C2b1a3a2-F8951 and C2*-Star Cluster differentiated 3852 years ago (95% CI = 3295-4497). Approximately 3558 years ago (95% CI = 3013-4144), C2b1a3a2-F8951 was divided into two main subclades, C2b1a3a2a-F14753 and C2b1a3a2b-F5483. Currently, samples of C2b1a3a2-F8951 were mainly from the House of Aisin Gioro clan, the Ao family from Daur and some individuals mainly from northeast China. Although other haplogroups are also found in the Ao family, including C2b1a2-M48, C2b1a3a1-F3796, C2a1b-F845, and N1c-M178, the haplogroup C2b1a3a2-F8951 is still the most distinct genetic component. For haplogroup C2b1a3a2-F8951, the time of the most recent common ancestor of the House of Aisin Gioro clan and the Ao family were both very late, just a few hundred years ago. Some family-specific Y-SNPs of the House of Aisin Gioro and the Ao family were also discovered. This revision evidently improved the resolving power of Y-chromosome phylogeny in northeast Asia, deepening our understanding of the origin of these two families, even the Mongolic-speaking population.

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.