Differences in DYF387S1 copy number distribution among haplogroups caused by haplogroup-specific ancestral Y-chromosome mutations

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.

High incidence of AZF duplications in clan-structured Iranian populations detected through Y chromosome sequencing read depth analysis

The ampliconic region of the human Y chromosome consists of large duplicated sequences that can undergo non-allelic homologous recombination (NAHR), resulting in structural rearrangements that may cause infertility, especially when they occur in the azoospermia factor b/c (AZFb/c) region. Although AZF duplications have long been neglected due to the technical limitations of STS-based studies that focused mainly on deletions, recent next generation sequencing (NGS) technologies provided evidence for their importance in fertility. In this study, a NGS read depth approach was used to detect AZFb/c rearrangements in 87 Iranians from different ethnic groups. The duplication frequency in Iran proved to be twice as high as in the "1000 Genomes" dataset. Interestingly, most duplications were found in patrilineal ethnic groups, possibly as a consequence of their lower male effective population size which can counteract negative selection. Moreover, we found a large 8.0 Mb duplication, resulting in a fourfold increase in the copy number of AZFc genes, which to our knowledge is the largest duplication ever reported in this region. Overall, our results suggest that it is important to consider not only AZF deletions but also duplications to investigate the causes of male infertility, especially in patrilineal clan-based populations.

High-resolution genotyping of 58 STRs in 635 Northern Han Chinese with MiSeq FGx ® Forensic Genomics System

Sequence polymorphisms were characterized at 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), and 24 Y chromosomal STRs (Y-STRs) in 635 Northern Han Chinese with the ForenSeq DNA Signature Prep Kit on the MiSeq FGx Forensic Genomics System. Since repeat region (RR) and flanking region (FR) variation can be detected by massively parallel sequencing (MPS), the increase in the number of unique alleles and the average of gene diversity was 78.18% and 3.51% between sequence and length, respectively. A total of 74 novel RR variants were identified at 33 STRs compared with STRSeq and previous studies, and 13 FR variants (rs1770275883, rs2053373277, rs2082557941, rs1925525766, rs1926380862, rs1569322793, rs2051848492, rs2051848696, rs2016239814, rs2053269960, rs2044518192, rs2044536444, and rs2089968964) were first submitted to dbSNP. Also, 99.94% of alleles were concordant between the ForenSeq DNA Signature Prep Kit and commercial CE kits. Discordance resulted from the low performance at D22S1045 and occasionally at DYS392, flanking region deletions at D7S820 and DXS10074, and the strict alignment algorithm at DXS7132. Null alleles at DYS505 and DYS448 and multialleles at DYS387S1a/b, DYS385a/b, DYS448, DYS505, DXS7132, and HPRTB were validated with other MPS and CE kits. Thus, a high-resolution sequence-based (SB) and length-based (LB) allele frequencies dataset from Northern Han Chinese has been established already. As expected, forensic parameters increased significantly on combined power of discrimination (PD) and combined power of exclusion (PE) at A-STRs, mildly on combined PD and combined mean exclusion chance (MEC) at X-STRs, and barely on discrimination capacity (DC) at Y-STRs. Additionally, MiSeq FGx quality metrics and MPS performance were evaluated in this study, which presented the high-quality of the dataset at 20 consecutive runs, such as ≥ 60% bases with a quality score of 20 or higher (%≥ Q20), > 60% of effective reads, > 2000 × of depth of coverage (DoC), ≥ 60% of allele coverage ratio (ACR) or heterozygote balance, ≥ 70% of inter-locus balance, and ≤ 0.4 of the absolute value of observed minus expected heterozygosity (|H_exp - H_obs|). In conclusion, MiSeq FGx can help us generate a high-resolution and high-quality dataset for human identification and population genetic studies.

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.

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.

RMplex: An efficient method for analyzing 30 Y-STRs with high mutation rates

Y-chromosomal short tandem repeats (Y-STRs) with high mutation rates are recognized as valuable genetic markers for differentiating paternally related men, who typically cannot be separated with standard Y-STRs, and were shown to provide paternal lineage differentiation on a higher resolution level than standard Y-STRs. Both features make Y-STRs with high mutation rates relevant in criminal casework, particularly in sexual assault cases involving highly unbalanced male-female DNA mixtures that often fail autosomal forensic STR profiling for the male donor. Previously, the number of known Y-STRs with mutation rates higher than 10^-2 per locus per generation termed rapidly mutating Y-STRs (RM Y-STRs) was limited to 13, which has recently been overcome by the discovery and characterization of 12 additional RM Y-STRs. Here, we present the development and validation of RMplex, an efficient genotyping system for analyzing 30 Y-STRs with high mutation rates, including all currently known RM Y-STRs, using multiplex PCR with capillary electrophoresis (CE) or massively parallel sequencing (MPS), overall targeting a total of 44 male-specific loci. If previously unavailable, repeat number assignations were provided based on newly generated MPS data. Validation tests based on the CE method demonstrated that the results were both repeatable and reproducible, full profiles were achieved with minimal input DNA of 250 pg for RMplex 1 and 100 pg for RMplex 2, and in the presence of inhibitors, or with a surplus of female DNA, the assays performed reasonably well. Application of RMplex to differentiate between paternally related men was exemplified in 32 males belonging to five different paternal pedigrees. Given further successful forensic validation testing, we envision the future application of RMplex in criminal cases where it is suspected, or cannot be excluded, that the crime scene trace originated from a male relatives of the suspect who is highlighted with standard Y-STR matching. Other applications of RMplex are in criminal cases without known suspects to differentiate between male relatives highlighted in familial searching based on standard Y-STR matching.

Developmental validation of the Microreader™ RM-Y ID System: a new rapidly mutating Y-STR 17-plex system for forensic application

Y-chromosomal short tandem repeats (Y-STRs) are widely applied to evolutionary, genealogical, and kinship analyses of male linages in forensic studies, but these low to midrange mutated Y-STRs typically fail to separate related males from the same paternal lineage. Recently, rapidly mutating Y-STRs (RM Y-STRs) have been demonstrated to improve the differentiation of male relatives and individuals. The Microreader™ RM-Y ID System is a new RM Y-STR kit that is capable of simultaneously amplifying 17 RM Y-STRs. Herein, to verify the efficiency and accuracy of the Microreader™ RM-Y ID System, developmental validation was conducted, including PCR-based studies, sensitivity, stability, species specificity, mixture, stutter percentage, and precision studies. Full profiles could be obtained when the hematin concentration was 250 μM, humic acid concentration was 1500 ng/μl, and tannic acid concentration was 200 ng/μl. Full profiles of the mixture of males/males could be detected up to a ratio of 19:1, and full profiles of females/males could always be detected even at ratios up to 24,000:1. Moreover, the forensic characteristics of 250 DNA-confirmed father-son pairs were analysed. The results showed that these 17 RM Y-STRs had high power for forensic discrimination (HD = 1) in the Chinese Han population, and the mutation rates were in the range of 4 × 10^-3 (95% CI 1.00 × 10^-4 to 2.21 × 10^-2, DYS464) to 8.8 × 10^-2 (95% CI 5.60 × 10^-2 to 1.30 × 10^-1, DYF399S1), indicating that the kit was effective for RM Y-STR studies and absolute individualisation of interrelated male individuals.

Genetic insights into the paternal admixture history of Chinese Mongolians via high-resolution customized Y-SNP SNaPshot panels

The Mongolian people, one of the Mongolic-speaking populations, are native to the Mongolian Plateau in North China and southern Siberia. Many ancient DNA studies recently reported extensive population transformations during the Paleolithic to historic periods in this region, while little is known about the paternal genetic legacy of modern geographically different Mongolians. Here, we genotyped 215 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) and 37 Y-chromosomal short tandem repeats (Y-STRs) among 679 Mongolian individuals from Hohhot, Hulunbuir, and Ordos in North China using the AGCU Y37 kit and our developed eight Y-SNP SNaPshot panels (including two panels first reported herein). The C-M130 Y-SNP SNaPshot panel defines 28 subhaplogroups, and the N/O/Q complementary Y-SNP SNaPshot panel defines 30 subhaplogroups of N1b-F2930, N1a1a1a1a3-B197, Q-M242, and O2a2b1a1a1a4a-CTS4658, which improved the resolution our developed Y-SNP SNaPshot panel set and could be applied for dissecting the finer-scale paternal lineages of Mongolic speakers. We found a strong association between Mongolian-prevailing haplogroups and some observed microvariants among the newly generated Y-STR haplotype data, suggesting the possibility of haplogroup prediction based on the distribution of Y-STR haplotypes. We identified three main ancestral sources of the observed Mongolian-dominant haplogroups, including the local lineage of C2*-M217 and incoming lineages from other regions of southern East Asia (O2*-M122, O1b*-P31, and N1*-CTS3750) and western Eurasia (R1*-M173). We also observed DE-M145, D1*-M174, C1*-F3393, G*-M201, I-M170, J*-M304, L-M20, O1a*-M119, and Q*-M242 at relatively low frequencies (< 5.00%), suggesting a complex admixture history between Mongolians and other incoming Eurasians from surrounding regions. Genetic clustering analyses indicated that the studied Mongolians showed close genetic affinities with other Altaic-speaking populations and Sinitic-speaking Hui people. The Y-SNP haplotype/haplogroup-based genetic legacy not only revealed that the stratification among geographically/linguistically/ethnically different Chinese populations was highly consistent with the geographical division and language classification, but also demonstrated that patrilineal genetic materials could provide fine-scale genetic structures among geographically different Mongolian people, suggesting that our developed high-resolution Y-SNP SNaPshot panels have the potential for forensic pedigree searches and biogeographical ancestry inference.

Sequence Read Depth Analysis of a Monophyletic Cluster of Y Chromosomes Characterized by Structural Rearrangements in the AZFc Region Resulting in DYS448 Deletion and DYF387S1 Duplication

The azoospermia factor c region (AZFc), located in the long arm of the human Y chromosome, is frequently involved in chromosome rearrangements, mainly due to non-allelic homologous recombination events that occur between the nearly identical sequences (amplicon) that comprises it. These rearrangements may have major phenotypic effects like spermatogenic failure or other pathologies linked to male infertility. Moreover, they may also be relevant in forensic genetics, since some of the Y chromosome short tandem repeats (Y-STRs) commonly used in forensic analysis are located in amplicons or in inter-amplicon sequences of the AZFc. In a previous study, we identified four phylogenetically related samples with a null allele at DYS448 and a tetrallelic pattern at DYF387S1, two Y-STRs located in the AZFc. Through NGS read depth analysis, we found that the unusual Y-STR pattern may be due to a 1.6 Mb deletion arising concurrently or after a 3.5 Mb duplication event. The observed large genomic rearrangement results in copy number reduction for the RBMY gene family as well as duplication of other AZFc genes. Based on the diversity of 16 additional Y-STRs, we estimated that the duplication/deletion event occurred at least twenty generations ago, suggesting that it has not been affected by negative selection.

Developmental validation of the Yfiler Platinum PCR Amplification Kit for forensic genetic caseworks and databases

Y chromosome kits are successfully applied in cases where human biological material exists. With the development of genotyping ability, more Y chromosomal markers are needed for finer identification of male individuals and lineages. In this study, a developmental validation of a newly emerged Y chromosome kit that combines two different kinds of markers: 38 Y-STRs and 3 Y-indels are conducted. The results show that this kit has high sensitivity when there is a small amount of DNA (125 pg), more than one male (minor:major = 1:7), or a mixture of males and females (male:female = 125pg:1875pg), inhibited substances (800 μM hematin and more than 1600 ng/μL humic acid). The kit exhibits high precision level with a standard deviation of allele size no more than 0.14 nt. Locus DYS481 shows the largest stutter rate, with three stutters per true allele. Population samples are well identified (MP of 0.001106), and mutations can be observed in father-son pairs (46 mutations in 70 pairs, 10 in locus DYS627). Out of all the population samples, 13.2% belong to haplogroup M117-O2a2b1a1, with their ethnic group being Han Chinese. The results show that this kit can improve the performance of identifying male individuals, obtaining more unique haplotypes (increasing from 894 to 918 of 1000 male samples) and higher discrimination capacity (increasing from 0.942 to 0.955) in this study compared to previous widely used Yfiler Plus kit. Besides, it gives information about their paternal lineages in forensic genetic casework and genealogical database construction.