Mutation analysis of 13 RM Y-STR loci in Han population from Beijing of China

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.

The sequence of the repetitive motif influences the frequency of multistep mutations in Short Tandem Repeats

Microsatellites, or Short Tandem Repeats (STRs), are subject to frequent length mutations that involve the loss or gain of an integer number of repeats. This work aimed to investigate the correlation between STRs' specific repetitive motif composition and mutational dynamics, specifically the occurrence of single- or multistep mutations. Allelic transmission data, comprising 323,818 allele transfers and 1,297 mutations, were gathered for 35 Y-chromosomal STRs with simple structure. Six structure groups were established: ATT, CTT, TCTA/GATA, GAAA/CTTT, CTTTT, and AGAGAT, according to the repetitive motif present in the DNA leading strand of the markers. Results show that the occurrence of multistep mutations varies significantly among groups of markers defined by the repetitive motif. The group of markers with the highest frequency of multistep mutations was the one with repetitive motif CTTTT (25% of the detected mutations) and the lowest frequency corresponding to the group with repetitive motifs TCTA/GATA (0.93%). Statistically significant differences (α = 0.05) were found between groups with repetitive motifs with different lengths, as is the case of TCTA/GATA and ATT (p = 0.0168), CTT (p < 0.0001) and CTTTT (p < 0.0001), as well as between GAAA/CTTT and CTTTT (p = 0.0102). The same occurred between the two tetrameric groups GAAA/CTTT and TCTA/GATA (p < 0.0001) - the first showing 5.7 times more multistep mutations than the second. When considering the number of repeats of the mutated paternal alleles, statistically significant differences were found for alleles with 10 or 12 repeats, between GATA and ATT structure groups. These results, which demonstrate the heterogeneity of mutational dynamics across repeat motifs, have implications in the fields of population genetics, epidemiology, or phylogeography, and whenever STR mutation models are used in evolutionary studies in general.

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.

Analysis of mutation rates and haplotypes of 23 Y-chromosomal STRs in Korean father-son pairs

Y-chromosomal short tandem repeats (Y-STRs) have been widely used in forensic genetics, and accurate knowledge of mutation rates at Y-STR loci is essential in kinship analysis. The main aim of this study was to estimate Y-STR mutation rates in Korean males. To obtain locus-specific mutations and haplotypes at 23 Y-STRs, we analyzed samples from 620 Korean father-son pairs. In addition, we also analyzed 476 unrelated individuals using the PowerPlex® Y23 System, with the aim of augmenting the available data for the Korean population. The PowerPlex® Y23 system facilitates analysis of the 23 Y-STR loci (DYS576, DYS570, DYS458, DYS635, DYS389 II, DYS549, DYS385, DYS481, DYS439, DYS456, DYS389 I, DYS19, DYS393, DYS391, DYS533, DYS437, DYS390, Y GATA H4, DYS448, DYS438, DYS392, and DYS643). Locus-specific mutation rate estimates varied from 0.00 to 8.06 × 10^-3 per generation, with an average mutation rate of 2.17 × 10^-3 (95% CI, 1.5-3.1 × 10^-3). To obtain comprehensive genetic values for the Korean population, we combined data obtained in this study with previously reported values, thereby enabling us to estimate the locus-specific mutation rates regarding 22,711 allele transmissions. By combining these data, we obtained an overall average mutation rate of 2.91 × 10^-3 (95% CI, 2.3-3.7 × 10^-3). In addition, among the 476 unrelated Korean males, we detected 467 different haplotypes, with an overall haplotype diversity value of 0.9999. By extracting haplotypes of Y-STRs described in previous literature on 23 Y-STR reported in Korea, we obtained gene diversities for 1133 Korean individuals. We believe that the values and characteristics of the 23 Y-STRs analyzed in this study will contribute to establishing criteria for forensic genetic interpretation, including kinship analysis.

Large-scale pedigree analysis highlights rapidly mutating Y-chromosomal short tandem repeats for differentiating patrilineal relatives and predicting their degrees of consanguinity

Rapidly mutating Y-chromosomal short tandem repeats (RM Y-STRs) were suggested for differentiating patrilineally related men as relevant in forensic genetics, anthropological genetics, and genetic genealogy. Empirical data are available for closely related males, while differentiation rates for more distant relatives are scarce. Available RM Y-STR mutation rate estimates are typically based on father-son pair data, while pedigree-based studies for efficient analysis requiring less samples are rare. Here, we present a large-scale pedigree analysis in 9379 pairs of men separated by 1-34 meioses on 30 Y-STRs with increased mutation rates including all known RM Y-STRs (RMplex). For comparison, part of the samples were genotyped at 25 standard Y-STRs mostly with moderate mutation rates (Yfiler Plus). For 43 of the 49 Y-STRs analyzed, pedigree-based mutation rates were similar to previous father-son based estimates, while for six markers significant differences were observed. Male relative differentiation rates from the 30 RMplex Y-STRs were 43%, 84%, 96%, 99%, and 100% for relatives separated by one, four, six, nine, and twelve meioses, respectively, which largely exceeded rates obtained by 25 standard Y-STRs. Machine learning based models for predicting the degree of patrilineal consanguinity yielded accurate and reasonably precise predictions when using RM Y-STRs. Fully matching haplotypes resulted in a 95% confidence interval of 1-6 meioses with RMplex compared to 1-25 with Yfiler Plus. Our comprehensive pedigree study demonstrates the value of RM Y-STRs for differentiating male relatives of various types, in many cases achieving individual identification, thereby overcoming the largest limitation of forensic Y-chromosome analysis.

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.

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.

Developmental Validation of a Rapidly Mutating Y-STR Panel Labeled by Six Fluoresceins for Forensic Research

The male-specific region of the human Y chromosome is a useful genetic marker for genealogical searching, male inheritance testing, and male DNA mixture deconvolution in forensic studies. However, the Y chromosomal short tandem repeats (Y-STRs) are difficult to distinguish among related males due to their low/medium mutation rate. In contrast, rapidly mutating (RM) Y-STRs exhibit unusually high mutation rates and possess great potential for differentiating male lineages. In this study, we developed a novel Y-STRs multiplex amplification assay of 32 RM Y-STRs by fragment analysis using six dye-labeled technologies (FAM, HEX, TAMRA, ROX, VIG, and SIZ). The development and the validation of the kit were carried out in accordance with the Scientific Working Group guidelines on DNA Analysis Methods. Identical allelic profiles of the 32 RM Y-STRs using a DNA 9948 sample as the positive control could be observed at different concentrations of PCR reagents. Further, the RM Y-STRs did not show cross-reactions with other common animal species, and the developed assay could tolerate interferences from common PCR inhibitors and mixed DNA samples. More importantly, the kit showed relatively high sensitivity and could detect trace DNA samples. Genetic distributions of 32 RM Y-STRs in the Guizhou Han population revealed that these RM Y-STRs showed relatively high genetic diversities. In conclusion, the RM Y-STR assay developed here showed good species specificity, high sensitivity, tolerance to inhibitors, and sample compatibility, which can be viewed as a highly efficient tool with high discrimination capacity for forensic male differentiation.

Improving the regional Y-STR haplotype resolution utilizing haplogroup-determining Y-SNPs and the application of machine learning in Y-SNP haplogroup prediction in a forensic Y-STR database: A pilot study on male Chinese Yunnan Zhaoyang Han population

Improving the resolution of the current widely used Y-chromosomal short tandem repeat (Y-STR) dataset is of great importance for forensic investigators, and the current approach is limited, except for the addition of more Y-STR loci. In this research, a regional Y-DNA database was investigated to improve the Y-STR haplotype resolution utilizing a Y-SNP Pedigree Tagging System that includes 24 Y-chromosomal single nucleotide polymorphism (Y-SNP) loci. This pilot study was conducted in the Chinese Yunnan Zhaoyang Han population, and 3473 unrelated male individuals were enrolled. Based on data on the male haplogroups under different panels, the matched or near-matching (NM) Y-STR haplotype pairs from different haplogroups indicated the critical roles of haplogroups in improving the regional Y-STR haplotype resolution. A classic median-joining network analysis was performed using Y-STR or Y-STR/Y-SNP data to reconstruct population substructures, which revealed the ability of Y-SNPs to correct misclassifications from Y-STRs. Additionally, population substructures were reconstructed using multiple unsupervised or supervised dimensionality reduction methods, which indicated the potential of Y-STR haplotypes in predicting Y-SNP haplogroups. Haplogroup prediction models were built based on nine publicly accessible machine-learning (ML) approaches. The results showed that the best prediction accuracy score could reach 99.71% for major haplogroups and 98.54% for detailed haplogroups. Potential influences on prediction accuracy were assessed by adjusting the Y-STR locus numbers, selecting Y-STR loci with various mutabilities, and performing data processing. ML-based predictors generally presented a better prediction accuracy than two available predictors (Nevgen and EA-YPredictor). Three tree models were developed based on the Yfiler Plus panel with unprocessed input data, which showed their strong generalization ability in classifying various Chinese Han subgroups (validation dataset). In conclusion, this study revealed the significance and application prospects of Y-SNP haplogroups in improving regional Y-STR databases. Y-SNP haplogroups can be used to discriminate NM Y-STR haplotype pairs, and it is important for forensic Y-STR databases to develop haplogroup prediction tools to improve the accuracy of biogeographic ancestry inferences.

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.

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.

Assessing 23 Y-STR loci mutation rates in Chinese Han father-son pairs from southwestern China

In this study, we have analyzed 23 Y-chromosomal short tandem repeats (Y-STRs) (DYS576, DYS389I, DYS389II, DYS448, DYS19, DYS391, DYS481, DYS549, DYS533, DYS438, DYS437, DYS570, DYS635, DYS390, DYS439, DYS392, DYS643, DYS393, DYS458, DYS460, DYS385ab, DYS456 and Y-GATA-H4) in 175 father-son sample pairs using a Microreader™ 24Y Direct ID system. Sixteen repeat mutations of father-son pairs at 10 loci, including three mutations at DYS570, 2 mutations at DYS549, DYS460, DYS458, and DYS576, and 1 mutation at other five loci, were revealed. Furthermore, all of the observed repeat mutations were single repeat changes with 5 (31.25%) repeat insertions and 11 (68.75%) repeat deletions. The deletion rate is more than two fold higher than of insertions (11:5 = 2.2-fold). Locus-specific mutation rates estimated varied between 5.71 × 10^-3 (CI from 0.1 × 10^-3 to 31.4 × 10^-3) and 1.71 × 10^-2 (CI from 3.6 × 10^-3 to 49.3 × 10^-3) for the 23 Y-STRs. An average mutation rate across all 23 Y-STR markers was estimated as 3.97 × 10^-3 (CI 2.3 × 10^-3 to 6.4 × 10^-3). Thus, locus-specific mutation rates in DYS460, DYS458, and DYS438, estimated are much higher than previously published comprehensive data, but an average mutation rate across all 23 Y-STR markers is similar to previous reports (3.97 × 10^-3 vs 4.34 × 10^-3). These results by characterizing Y-STR mutations will not only provided new information for Y-STR mutations but also might be important for paternal lineage identification, kinship analysis, and family relationship reconstruction in our forensic Y-STR analysis.

Analyzing genetic polymorphism and mutation of 44 Y-STRs in a Chinese Han population of Southern China

Short tandem repeat on the non-recombining part of chromosome Y with paternally inheritable capability is a valuable tool in the studies of forensic genetics, population genetics and anthropology. The mutation rate of Y-STR is an important parameter in the applications. A total of 629 haplotypes at 44 Y-STR markers were found in 629 unrelated males of our population. Mutation rates at 44 Y-STR loci ranged from 0 (CI: 0-5.70 × 10^-3) to 40.63 × 10^-3 (25.90 × 10^-3-57.2 × 10^-3) in our population. A higher mutation rate was noted at DYS612, DYS449, DYS547, DYS518, DYS576, DYS627, DYF403S1b, DYF387S1, DYS385a/b, DYS527a/b, DYF404S1, DYF403S1a and DYF399S1 in this population. The Y-STR set showed a higher discrimination capacity in forensic applications, and the present study provided reference data for the application of forensic and population genetics.