Microhaplotype and Y-SNP/STR (MY): A novel MPS-based system for genotype pattern recognition in two-person DNA mixtures

Nanopore sequencing of MiniHap biomarkers for forensic DNA mixture deconvolution: A proof-of-principle study

Mixture deconvolution remains one of the major challenges in the field of forensic science. Currently, genetic markers are used and studied in the field of forensic genetics, including short tandem repeat (STR), insertion/deletion polymorphism (InDel), single nucleotide polymorphism (SNP), InDel closely linked to STR (DIP-STR), SNP closely linked to STR (SNP-STR), InDel closely linked to SNP (DIP-SNP) and microhaplotype (MH), all of which have been studied for DNA mixture analysis and have their own advantages and disadvantages. Mini-haplotype (MiniHap), as a novel haplotype genetic marker, contains 5 or more SNPs. A previous study has substantiated its significant high polymorphic characteristics, and it is expected to have potential applications in individual identification, paternity testing, ancestry inference, and mixture deconvolution. In this study, we first screened 22 MiniHaps with high polymorphism potential and constructed a panel based on the QNome nanopore sequencing device. Subsequently, we tested 100 unrelated Chinese Han individuals to evaluate the sequencing performance, allele (haplotype) frequencies, effective number of alleles (Ae) and forensic parameters of the 22 MiniHaps markers included in this novel assay. Next, a series of mixture simulations (two- or three-person mixtures with mixing ratios of 1:1-1:99 or 1:1:1-1:8:1) based on three standard materials (9947 A, 9948 and 2800 M) were detected by this MiniHap panel to explore its potential for DNA mixture deconvolution. The average Ae value was 10.9574, and 52.38 % of MiniHap loci had Ae values greater than 12.0000. The mean values of genetic diversity (GD) and power of discrimination (PD) were 0.8717 and 0.9457, respectively. Notably, most MiniHaps (85.71 %) had PD values exceeding 0.9000. The combined match probability (CMP) and combined power of exclusion (CPE) of this MiniHap panel were 4.4505 × 10-31 and 0.999999999999999996653, respectively. Moreover, the results of mixture analysis demonstrated that this MiniHap panel allowed detecting the components of minor contributor (s) even in imbalanced mixture samples, with detection limits of 1:39 and 1:8:1 for two- and three-person mixtures, respectively. In summary, MiniHap markers have remarkable application potential in mixture deconvolution, and it is necessary to conduct in-depth research on MiniHap markers for mixture analysis in the future.

MhGeneS: An Analytical Pipeline to Allow for Robust Microhaplotype Genotyping

Microhaplotypes are small linked genomic regions comprising two or more single-nucleotide polymorphisms (SNPs) that are being applied in forensics and are emerging in wildlife monitoring studies and genomic epidemiology. Typically, targeted in non-coding regions, microhaplotypes in exonic regions can be designed with larger amplicons to capture functional non-synonymous sites and minimise insertion/deletion (indel) polymorphisms. Quality control is an important first step for high-confidence genotyping to counteract such false-positive variants. As genetic markers with higher polymorphism compared to biallelic SNPs, it is critical to ensure sequencing errors across the microhaplotype amplicon are filtered out to avoid introducing false-haplotypes. We developed the MhGeneS pipeline which works in tandem with Seq2Sat to help validate microhaplotype genotyping of the coding region of genes, with broader applicability to any microhaplotype profiling. We genotyped microhaplotype regions of the Zfx (≅ 160 bp) and Zfy (≅ 140 bp) genes, as well as an exon of the prion protein (Prnp) gene (≅ 370 bp) in caribou (Rangifer tarandus) using paired-end Illumina technology. As important quality metrics affecting microhaplotype calling, we identified the sequencing error rate profile related to the overlap or non-overlap of paired-end reads as well as the read depth as significant. In the case of Prnp, we achieved confident microhaplotype calling through MhGeneS by removing small sections of the 5' and 3' amplicons and using a minimum read depth of 20. Read depth and sequence trimming may be locus-specific, and validation of these parameters is recommended before the high-throughput profiling of samples.

Discovery unbalanced DNA mixtures and evaluation mixing ratio via a droplet digital PCR method

BACKGROUND

Small amounts of DNA from a perpetrator collected during crime-scene investigations can be masked by large amounts of DNA from the victim. These samples can provide important information for the perpetrator's conviction. Short tandem repeat (STR) detection system is not sensitive enough to detect trace amounts of minor components in unbalanced mixed DNA. We developed a system using droplet digital polymerase chain reaction (ddPCR) capable of discovering trace components and accurately determining the ratio of mixed DNA in extremely unbalanced mixtures.

METHODS

The non-recombining regions of the X chromosome and Y chromosome were quantified in the DNA of male and female mixtures using duplex ddPCR. Absolute quantification of low-abundance portions of trace samples and unbalanced mixtures was done using different mixing ratios.

RESULTS

The ddPCR system could be used to detect low-abundance samples with < 5 copies of DNA components in an extremely unbalanced mixture at a mixing ratio of 10000:1. The high sensitivity and specificity of the system could identify the mixing ratio of mixed DNA accurately.

CONCLUSIONS

A ddPCR system was developed for evaluation of mixed samples of male DNA and female DNA. Our system could detect DNA quantities as low as 5 copies in extremely unbalanced mixed samples with good specificity and applicability. This method could assist forensic investigators in avoiding the omission of important physical evidence, and evaluating the ratio of mixed male/female trace samples.

Developmental and validation of a novel small and high-efficient panel of microhaplotypes for forensic genetics by the next generation sequencing

BACKGROUND

In the domain of forensic science, the application of kinship identification and mixture deconvolution techniques are of critical importance, providing robust scientific evidence for the resolution of complex cases. Microhaplotypes, as the emerging class of genetic markers, have been widely studied in forensics due to their high polymorphisms and excellent stability.

RESULTS AND DISCUSSION

In this research, a novel and high-efficient panel integrating 33 microhaplotype loci along with a sex-determining locus was developed by the next generation sequencing technology. In addition, we also assessed its forensic utility and delved into its capacity for kinship analysis and mixture deconvolution. The average effective number of alleles (Ae) of the 33 microhaplotype loci in the Guizhou Han population was 6.06, and the Ae values of 30 loci were greater than 5. The cumulative power of discrimination and cumulative power of exclusion values of the novel panel in the Guizhou Han population were 1-5.6 × 10- 43 and 1-1.6 × 10- 15, respectively. In the simulated kinship analysis, the panel could effectively distinguish between parent-child, full-sibling, half-sibling, grandfather-grandson, aunt-nephew and unrelated individuals, but uncertainty rates clearly increased when distinguishing between first cousins and unrelated individuals. For the mixtures, the novel panel had demonstrated excellent performance in estimating the number of contributors of mixtures with 1 to 5 contributors in combination with the machine learning methods.

CONCLUSIONS

In summary, we have developed a small and high-efficient panel for forensic genetics, which could provide novel insights into forensic complex kinships testing and mixture deconvolution.

Application of a newly constructed NGS panel with 45 X-linked microhaplotypes demonstrates the unique value of X-MH for kinship testing and mixture analysis

X-linked microhaplotypes (X-MHs) have the potential to be a valuable supplementary tool in complex kinship identification or the resolution of DNA mixtures, because they bring together the distinctive genetic pattern of X chromosomal markers and the benefits of microhaplotypes (MHs). In this study, we used the 1000 Genome database to screen and select 63 X-MHs; 18 MHs were filtered out though a batch sequencing assessment of the DNA samples collected from 112 unrelated Chinese Han individuals. The resulting 45-plex panel performed well in comprehensive assessments including repeatability, sensitivity, species specificity, resistance to PCR inhibitors or degradation, mutation rate, and accuracy in detecting DNA mixture samples. The minimum amount of DNA template that can be tested with this panel is 0.5 ng. Additionally, the alleles of the minor contributor can be accurately detected when the mixture rate is larger than 1:9 in female-male mixture or 1:19 in male-male mixture. Then, we calculated population parameters on each MH based on the allele frequency data obtained from the sequence results of the aforementioned 112 unrelated samples. Combining these parameters on each MH, it can be calculated that TDPm, TDPf, CPET, CPEDFM, CPEDFF and CNCEP3 of the 45-plex system were 1-8.99×10-13, 1-1.62×10-19, 0.9999999995, 0.9999981, 0.9955, 0.9999971 and 0.99940, respectively, indicating that the panel is capable in personal identification and parentage testing. To reveal the unique advantage of X-MHs in the analyses of complex kinship and male DNA mixture, further assessments were made. For complex kinship identification, 22 types of individual pairs with different second-degree kinship were simulated and different types of likelihood ratios (LR) were calculated for each. The results revealed that the panel can achieve accuracy of approximately 70 %∼80 % when dividing each of the three types of second-degree kinships into three or four groups. Theoretically, such sub-division cannot be done by using independent autosomal markers. For male DNA mixture analysis without suspects, the maximum likelihood ratio strategy was derived and employed in the estimation of the number of male contributors (NOMC). Simulations were conducted to verify the efficacy of the 45-plex panel in the field and to compare it with autosomal markers by assuming the 45 MHs as autosomal ones. The results showed that X-MHs can achieve higher accuracy in the estimation of NOMC than autosomal ones when the mixed males were unrelated. The results highlighted the unique value of X-linked MHs in complex kinship and male mixture analyses.

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.

Developmental validation of the STRSeqTyper122 kit for massively parallel sequencing of forensic STRs

Massively parallel sequencing allows for integrated genotyping of different types of forensic markers, which reduces DNA consumption, simplifies experimental processes, and provides additional sequence-based genetic information. The STRseqTyper122 kit genotypes 63 autosomal STRs, 16 X-STRs, 42 Y-STRs, and the Amelogenin locus. Amplicon sizes of 117 loci were below 300 bp. In this study, MiSeq FGx sequencing metrics for STRseqTyper122 were presented. The genotyping accuracy of this kit was examined by comparing to certified genotypes of NIST standard reference materials and results from five capillary electrophoresis-based kits. The sensitivity of STRseqTyper122 reached 125 pg, and > 80% of the loci were correctly called with 62.5 pg and 31.25 pg input genomic DNA. Repeatability, species specificity, and tolerance for DNA degradation and PCR inhibitors of this kit were also evaluated. STRseqTyper122 demonstrated reliable performance with routine case-work samples and provided a powerful tool for forensic applications.

Deconvoluting multi-person biological mixtures and accurate characterization and identification of separated contributors using non-targeted single-cell DNA sequencing

The genetic characterization and identification of individuals who contributed to biological mixtures are complex and mostly unresolved tasks. These tasks are relevant in various fields, particularly in forensic investigations, which frequently encounters crime scene stains generated by more than one person. Currently, forensic mixture deconvolution is mostly performed subsequent to forensic DNA profiling at the level of the mixed DNA profiles, which comes with several limitations. Some previous studies attempted at separating single cells prior to forensic DNA profiling. However, these approaches are biased at selection of the cells and, due to their targeted DNA analysis on low template DNA, provide incomplete and unreliable forensic DNA profiles. We recently demonstrated the feasibility of performing mixture deconvolution prior to forensic DNA profiling through the utilization of a non-targeted single-cell transcriptome sequencing (scRNA-seq). In addition to individual-specific mixture deconvolution, this approach also allowed accurate characterisation of biological sex, biogeographic ancestry and individual identification of the separated mixture contributors. However, RNA has the forensic disadvantage of being prone to degradation, and sequencing RNA - focussing on coding regions - limits the number of single nucleotide polymorphisms (SNPs) utilized for genetic mixture deconvolution, characterization, and identification. These limitations can be overcome by performing single-cell sequencing on the level of DNA instead of RNA. Here, for the first time, we applied non-targeted single-cell DNA sequencing (scDNA-seq) by applying the scATAC-seq (Assay for Transposase-Accessible Chromatin with sequencing) technique to address the challenges of mixture deconvolution in the forensic context. We demonstrated that scATAC-seq, together with our recently developed De-goulash data analysis pipeline, is capable of deconvoluting complex blood mixtures of five individuals from both sexes with varying biogeographic ancestries. We further showed that our approach achieved correct genetic characterization of the biological sex and the biogeographic ancestry of each of the separated mixture contributors and established their identity. Furthermore, by analysing in-silico generated scATAC-seq data mixtures, we demonstrated successful individual-specific mixture deconvolution of i) highly complex mixtures of 11 individuals, ii) balanced mixtures containing as few as 20 cells (10 per each individual), and iii) imbalanced mixtures with a ratio as low as 1:80. Overall, our proof-of-principle study demonstrates the general feasibility of scDNA-seq in general, and scATAC-seq in particular, for mixture deconvolution, genetic characterization and individual identification of the separated mixture contributors. Furthermore, it shows that compared to scRNA-seq, scDNA-seq detects more SNPs from fewer cells, providing higher sensitivity, that is valuable in forensic genetics.

Development and evaluation of a novel panel containing 188 microhaplotypes for 2nd-degree kinship testing in the Hebei Han population

Distant kinship identification is one of the critical problems in forensic genetics. As a new type of genetic marker defined and discussed in the last decade, the microhaplotype (MH) has drawn much attention in such identification owing to its specific advantages to traditional short tandem repeat (STR) or single nucleotide polymorphism (SNP) markers. In this study, MH markers were screened step by step from the 1000 Genomes Project database, and a novel multiplex panel containing 188 MHs (in which 181 are reported the first time, while 1 was reported in a previous study and the other 6 have partial overlaps with known markers) was constructed for application in 2nd- and 3rd-degree kinship identification. Along with the construction, a novel MH nomenclature was proposed, in which the SNP position information they contained was taken into account to eliminate the possibility that the same locus was named differently interlaboratory. After a series of evaluations, the panel was shown to have good sequencing accuracy, high sensitivity, species specificity, and resistance to anti-PCR inhibitors or degradation. Population data of the 188 MHs were calculated based on the genetic information of 221 unrelated Hebei Han individuals, and the effective number of alleles (Ae) ranged from 2.0925 to 8.2634 (with an average of 2.9267). For the whole system, the cumulative matching probability (CMP), the cumulative power of exclusion in paternity testing of duos (CPEduo) and that of trios (CPEtrio) reached 2.8422 × 10^-137, 1-1.3109 × 10^-21, and 1-2.8975 × 10^-39, respectively, indicating that this panel was satisfactory for individual identification and paternity testing. Then, the efficiency of the 188 MHs in 2nd- and 3rd-degree kinship testing was studied based on 30 extended families consisting of 179 2nd-degree and 121 3rd-degree relatives, as well as simulations of 0.5 million pairs of those two kinships. The results showed that clear opinions would be given in 83.36% of 2nd-degree identifications with a false rate less than 10^-5, when the confirming and excluding thresholds of cumulative likelihood ratio (CLR) were set as 10⁴ and 10^-4, respectively. This panel is still not sufficient to solve the problem of 3rd-degree kinship identification alone, and approximately 300 or 870 MH loci would be needed in 2nd- or 3rd-degree kinship identification, respectively, to achieve a system efficiency not less than 0.99 with such a threshold set; such necessary numbers would be used only as a reference in further research.

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

State of the Art for Microhaplotypes

In recent years, the number of publications on microhaplotypes has averaged more than a dozen papers annually. Many have contributed to a significant increase in the number of highly polymorphic microhaplotype loci. This increase allows microhaplotypes to be very informative in four main areas of forensic uses of DNA: individualization, ancestry inference, kinship analysis, and mixture deconvolution. The random match Probability (RMP) can be as small as 10−100 for a large panel of microhaplotypes. It is possible to measure the heterozygosity of an MH as the effective number of alleles (Ae). Ae > 7.5 exists for African populations and >4.5 exists for Native American populations for a smaller panel of two dozen selected microhaplotypes. Using STRUCTURE, at least 10 different ancestral clusters can be defined by microhaplotypes. The Ae for a locus is also identical to the Paternity Index (PI), the measure of how informative a locus will be in parentage testing. High Ae loci can also be useful in missing persons cases. Finally, high Ae microhaplotypes allow the near certainty of seeing multiple additional alleles in a mixture of two or more individuals in a DNA sample. In summary, a panel of higher Ae microhaplotypes can outperform the standard CODIS markers.