A novel forensic panel of 186-plex SNPs and 123-plex STR loci based on massively parallel sequencing

Investigating SNP typing using alternative reference materials with the FORCE panel and QIAseq® chemistry

Advancements in forensic science have introduced single nucleotide polymorphism (SNP) markers as crucial tools in kinship, ancestry, and identity testing, and in predicting phenotypic traits. The emergence of Forensic Genetic Genealogy (FGG) and massively parallel sequencing (MPS) technologies have further enhanced the utility of SNP markers, generating interest in their application within the forensic community. The FORensic Capture Enrichment (FORCE) panel, targeting 5497 SNPs including those associated with ancestry, phenotype, identity, and kinship, is specifically designed for direct kinship comparisons and the recovery of autosomal markers for direct identification in cases where traditional short tandem repeat (STR) typing is unsuitable. This study aimed to evaluate alternative reference materials-hair roots, hair shafts, and fingernail clippings-using the FORCE panel and QIAseq® chemistry for direct identification. We assessed SNP recovery and concordance between these materials and buccal swabs, and we compared predictions for phenotype, Y-haplogroup, and biogeographic ancestry. Additionally, the study examined the concordance and performance on two sequencing platforms: MiSeq FGx and NextSeq 550. Our results demonstrated high (99.62-100 %) SNP concordance rates between alternative reference materials and buccal swabs, with fingernail samples showing the highest SNP recovery and concordance of the alternate reference materials. Phenotypic, ancestry, and Y-haplogroup ancestry predictions from alternative materials were 100 % consistent with those from buccal samples. However, some discrepancies in phenotype predictions were noted when comparing predictions to self-reported data. Both sequencing platforms provided comparable results, with NextSeq 550 mid-output kit offering 3.5X higher coverage and potential cost-saving advantages through greater sample multiplexing. In conclusion, the FORCE panel, combined with QIAseq® chemistry, effectively profiles SNPs from various alternative reference materials, offering reliable concordance and comprehensive genotype recovery even from low-input and degraded samples. This highlights its potential application in forensic investigations where traditional reference samples are unavailable.

Improved second-degree kinship analysis using the FGID forensic four-in-one DNA typing kit

In general, using additional genetic markers can greatly improve the system power and accuracy of complex kinship identification. In this study, we employed the FGID Forensic Four-in-one DNA typing kit on the DNBSEQ-G99RS platform to gather sequence information for 66 autosomal STRs (A-STRs) and 132 autosomal SNPs (A-SNPs) of 24 pairs of second-degree relatives (13 grandparent-grandchild pairs and 11 uncle/aunt-nephew/niece pairs) from 9 families, along with 100 unrelated Han Chinese individuals. The likelihood ratio (LR) and identical-by-state (IBS) methods were applied to the kinship analysis, respectively. As anticipated, compared to the results using the ForenSeq™ DNA Signature Prep Kit, the LR values for second-degree relatives improved with the addition of more STR and SNP markers in this study, while LR values for unrelated individuals decreased. Sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) all exceeded 0.9600 at the cut-off values of t1 = -3 and t2 = 3, with 96.80% of simulated relationship pairs accurately classified as either second-degree relationships or unrelated individuals. Furthermore, IBS values based on sequence-based STRs were significantly lower than those based on length-based STRs. Overall, these results indicated that the FGID kit, which incorporates an increasing number of genetic markers, could effectively enhance the system performance in second-degree kinship analysis.

Massively parallel sequencing of a forensic combined panel of 107-plex STR loci and 292-plex SNP loci in the Han Chinese population

Massively parallel sequencing (MPS), a well-established strategy for forensic DNA profiling, enables the simultaneous sequencing of multiple targeted loci of multiple samples at a single-base resolution with high coverage. In this study, we developed a novel typing system by combining solution-based hybrid capture methods with MPS to target as many as 107 short tandem repeats (STRs) and 292 single nucleotide polymorphisms (SNPs) in the Han Chinese population. Completely accurate and concordant STR genotypes were obtained when compared to typing results generated from conventional capillary electrophoresis analysis, with six loci exhibiting inferior performance due to allele dropout or even locus dropout. The locus detection success reached 85.2 % for STRs at a DNA input of 10 ng and 95.61 % for SNPs at a DNA input of 5 ng. Mixture studies substantiated the considerable potential of our system in identifying minor contributor alleles at both STR and SNP loci. Additionally, the system demonstrated full inferential abilities in distinguishing first-degree kinship from unrelated individual pairs and achieved significant effectiveness of 99.78 % and 80.2 % for the identification of second- and third-degree kinship, respectively. These findings indicated that our novel typing system is highly discriminative and informative when used in the Han Chinese population and would be highly efficient for use in paternity testing and complex kinship analysis.

Forensic efficiency evaluation of a mtDNA whole genome sequencing system constructed with long fragment amplification strategy on DNA nanoball sequencing platform

Mitochondrial DNA (mtDNA) is an important genetic marker for degraded biological sample identification, maternal pedigree tracing, and population genetic structure study owing to its characteristics of high copy number, anti-degradable ring structure, and maternal inheritance. Whole mtDNA genome sequencing is an optimal method for the analysis of mtDNA polymorphism and heterogeneity because it allows for the comprehensive use of maternal genetic information. However, because of lacking quantitative evaluations for sequencing data, the scientific interpretation standards for mtDNA sequencing results of the previously used sequencing systems are often different, and false positive or false negative results are prone to occur when faced with the interference of nuclear genomic DNA, or the heterogeneities of mtDNA sequence and structure. In this study, we evaluated a novel mtDNA whole genome sequencing system using long fragment amplification strategy on the DNA nanoball (DNB) sequencing platform. This system demonstrated high sequencing quality and specific mtDNA sequencing efficiencies on positive control DNA and FTA bloodstain samples, as the average Q20 and Q30 values of the corresponding samples were 97.17 % and 91.93 %; 97.37 % and 92.48 %, respectively. The mean mapping percentages for the reference sequences of whole genome DNA (wgDNA), mtDNA, and nuclear genomic DNA (ngDNA) in the corresponding samples were 99.98 %, 99.97 %, 0.03 %, and 99.91 %, 99.40 %, 0.60 %; respectively. The average error calling rates for the bases A, C, G, and T of the whole mtDNA genome were 0.2519 %, 0.2550 %, 0.2906 %; and 0.2392 %, respectively. The efficacy of heteroplasmy identification was assessed using a set of theoretical sites with predetermined rates. These sites were created by combining the samples with known mtDNA haplotypes in certain proportions. The absolute errors between observed and theoretical heteroplasmy values were 89.59 %, 74.68 %, 50.20 %, 12.65 %, 8.31 %, and 4.85 %, while the theoretical heteroplasmy values were 5 %, 10 %, 20 %, 80 %, 90 %, and 95 %, respectively. The absolute error exhibited relative stability when the mtDNA sequencing depth exceeded 500×. Furthermore, the system sequencing efficiency was also confirmed among different kinds of samples, and these samples included natural samples (e.g., peripheral blood samples preserved on FTA cards for 2 and 11 years, and on filter paper for 6 and 9 years), degraded samples, sensitivity samples, samples derived from various bodily fluids, and maternal pedigree samples. In summary, the whole mtDNA genome sequencing system used for forensic identification demonstrated high performance in analyzing mtDNA sequence information, and showed significant prospects for forensic application and maternal genetic research.

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.

Comprehensive landscape of non-CODIS STRs in global populations provides new insights into challenging DNA profiles

The worldwide implementation of short tandem repeats (STR) profiles in forensic genetics necessitated establishing and expanding the CODIS core loci set to facilitated efficient data management and exchange. Currently, the mainstay CODIS STRs are adopted in most general-purpose forensic kits. However, relying solely on these loci failed to yield satisfactory results for challenging tasks, such as bio-geographical ancestry inference, complex DNA mixture profile interpretation, and distant kinship analysis. In this context, non-CODIS STRs are potent supplements to enhance the systematic discriminating power, particularly when combined with the high-throughput next-generation sequencing (NGS) technique. Nevertheless, comprehensive evaluation on non-CODIS STRs in diverse populations was scarce, hindering their further application in routine caseworks. To address this gap, we investigated genetic variations of 178 historically available non-CODIS STRs from ethnolinguistically different worldwide populations and studied their characteristics and forensic potentials via high-coverage whole genome sequencing (WGS) data. Initially, we delineated the genomic properties of these non-CODIS markers through sequence searching, repeat structure scanning, and manual inspection. Subsequent population genetics analysis suggested that these non-CODIS STRs had comparable polymorphism levels and forensic utility to CODIS STRs. Furthermore, we constructed a theoretical next-generation sequencing (NGS) panel comprising 108 STRs (20 CODIS STRs and 88 non-CODIS STRs), and evaluated its performance in inferring bio-geographical ancestry origins, deconvoluting complex DNA mixtures, and differentiating distant kinships using real and simulated datasets. Our findings demonstrated that incorporating supplementary non-CODIS STRs enabled the extrapolation of multidimensional information from a single STR profile, thereby facilitating the analysis of challenging forensic tasks. In conclusion, this study presents an extensive genomic landscape of forensic non-CODIS STRs among global populations, and emphasized the imperative inclusion of additional polymorphic non-CODIS STRs in future NGS-based forensic systems.

A novel 193-plex MPS panel integrating STRs and SNPs highlights the application value of forensic genetics in individual identification and paternity testing

Massively parallel sequencing (MPS) has emerged as a promising technology for targeting multiple genetic loci simultaneously in forensic genetics. Here, a novel 193-plex panel was designed to target 28 A-STRs, 41 Y-STRs, 21 X-STRs, 3 sex-identified loci, and 100 A-SNPs by employing a single-end 400 bp sequencing strategy on the MGISEQ-2000™ platform. In the present study, a series of validations and sequencing of 1642 population samples were performed to evaluate the overall performance of the MPS-based panel and its practicality in forensic application according to the SWGDAM guidelines. In general, the 193-plex markers in our panel showed good performance in terms of species specificity, stability, and repeatability. Compared to commercial kits, this panel achieved 100% concordance for standard gDNA and 99.87% concordance for 14,560 population genotypes. Moreover, this panel detected 100% of the loci from 0.5 ng of DNA template and all unique alleles at a 1:4 DNA mixture ratio (0.2 ng minor contributor), and the applicability of the proposed approach for tracing and degrading DNA was further supported by case samples. In addition, several forensic parameters of STRs and SNPs were calculated in a population study. High CPE and CPD values greater than 0.9999999 were clearly demonstrated and these results could be useful references for the application of this panel in individual identification and paternity testing. Overall, this 193-plex MPS panel has been shown to be a reliable, repeatable, robust, inexpensive, and powerful tool sufficient for forensic practice.

A pipeline for sample tagging of whole genome bisulfite sequencing data using genotypes of whole genome sequencing

BACKGROUND

In large-scale high-throughput sequencing projects and biobank construction, sample tagging is essential to prevent sample mix-ups. Despite the availability of fingerprint panels for DNA data, little research has been conducted on sample tagging of whole genome bisulfite sequencing (WGBS) data. This study aims to construct a pipeline and identify applicable fingerprint panels to address this problem.

RESULTS

Using autosome-wide A/T polymorphic single nucleotide variants (SNVs) obtained from whole genome sequencing (WGS) and WGBS of individuals from the Third China National Stroke Registry, we designed a fingerprint panel and constructed an optimized pipeline for tagging WGBS data. This pipeline used Bis-SNP to call genotypes from the WGBS data, and optimized genotype comparison by eliminating wildtype homozygous and missing genotypes, and retaining variants with identical genomic coordinates and reference/alternative alleles. WGS-based and WGBS-based genotypes called from identical or different samples were extensively compared using hap.py. In the first batch of 94 samples, the genotype consistency rates were between 71.01%-84.23% and 51.43%-60.50% for the matched and mismatched WGS and WGBS data using the autosome-wide A/T polymorphic SNV panel. This capability to tag WGBS data was validated among the second batch of 240 samples, with genotype consistency rates ranging from 70.61%-84.65% to 49.58%-61.42% for the matched and mismatched data, respectively. We also determined that the number of genetic variants required to correctly tag WGBS data was on the order of thousands through testing six fingerprint panels with different orders for the number of variants. Additionally, we affirmed this result with two self-designed panels of 1351 and 1278 SNVs, respectively. Furthermore, this study confirmed that using the number of genetic variants with identical coordinates and ref/alt alleles, or identical genotypes could not correctly tag WGBS data.

CONCLUSION

This study proposed an optimized pipeline, applicable fingerprint panels, and a lower boundary for the number of fingerprint genetic variants needed for correct sample tagging of WGBS data, which are valuable for tagging WGBS data and integrating multi-omics data for biobanks.

Recent advances in forensic biology and forensic DNA typing: INTERPOL review 2019-2022

This review paper covers the forensic-relevant literature in biological sciences from 2019 to 2022 as a part of the 20th INTERPOL International Forensic Science Managers Symposium. Topics reviewed include rapid DNA testing, using law enforcement DNA databases plus investigative genetic genealogy DNA databases along with privacy/ethical issues, forensic biology and body fluid identification, DNA extraction and typing methods, mixture interpretation involving probabilistic genotyping software (PGS), DNA transfer and activity-level evaluations, next-generation sequencing (NGS), DNA phenotyping, lineage markers (Y-chromosome, mitochondrial DNA, X-chromosome), new markers and approaches (microhaplotypes, proteomics, and microbial DNA), kinship analysis and human identification with disaster victim identification (DVI), and non-human DNA testing including wildlife forensics. Available books and review articles are summarized as well as 70 guidance documents to assist in quality control that were published in the past three years by various groups within the United States and around the world.

Systematic selections and forensic application evaluations of 111 individual identification SNPs in the Chinese Inner Mongolia Manchu group

Single nucleotide polymorphism (SNP) possesses a promising application in forensic individual identification due to its wide distribution in the human genome and the ability to carry out the genotyping of degraded biological samples by designing short amplicons. Some commonly used individual identification SNPs are less polymorphic in East Asian populations. In order to improve the individual identification efficiencies in East Asian populations, SNP genetic markers with relatively higher polymorphisms were selected from the 1,000 Genome Project phase III database in East Asian populations. A total of 111 individual identification SNPs (II-SNPs) with the observed heterozygosity values greater than 0.4 were screened in East Asian populations, and then, the forensic efficiencies of these selected SNPs were also evaluated in Chinese Inner Mongolia Manchu group. The observed heterozygosity and power of discrimination values at 111 II-SNPs in the Inner Mongolia Manchu group ranged from 0.4011 to 0.7005, and 0.5620 to 0.8025, respectively, and the average value of polymorphism information content was greater than 0.3978. The cumulative match probability and combined probability of exclusion values at II-SNPs were 7.447E^-51 and 1-4.17E^-12 in the Inner Mongolia Manchu group, respectively. The accumulative efficiency results indicated that the set of II-SNPs could be used as a potential tool for forensic individual identification and parentage testing in the Manchu group. The sequencing depths ranged from 781× to 12374×. And the mean allele count ratio and noise level were 0.8672 and 0.0041, respectively. The sequencing results indicated that the SNP genetic marker detection based on the massively parallel sequencing technology for SNP genetic markers had high sequencing performance and could meet the sequencing requirements of II-SNPs in the studied group.

Development a multiplex panel of AISNPs, multi-allelic InDels, microhaplotypes and Y-SNP/InDel loci for multiple forensic purposes via the NGS

Recently, next generation sequencing showed the promising application value in forensic research. In this study, we constructed a multiplex system of different molecular genetic markers based on the previous selected AISNPs, multi-allelic InDels, microhaplotypes and Y-SNP/InDel loci and evaluated forensic efficiencies of the system in Chinese Shaanxi Han, Hui and Mongolian groups via the NGS platform. Ancestry informative analyses of Shaanxi Han, Hui and Mongolian groups revealed that most Mongolian individuals could be differentiated from Shaanxi Hans and Huis based on the selected AISNPs. Multi-allelic InDels and microhaplotypes showed the multiple allele variations and possessed relatively high genetic polymorphisms in these three groups, indicating these loci could also provide higher forensic efficiencies for individual identification and paternity testing. Based on Y-SNPs, different haplogroup distributions were observed among Shaanxi Han, Hui and Mongolian groups. In conclusion, the self-developed system could be used to simultaneously carry out the individual identification, paternity analysis, mixture deconvolution, forensic ancestry information analysis and Y chromosomal haplogroup inference, which could provide more investigative clues in forensic practices. This article is protected by copyright. All rights reserved.

Development and validation of a novel 133-plex forensic STR panel (52 STRs and 81 Y-STRs) using single-end 400 bp massive parallel sequencing

Short tandem repeats (STRs) are the preferred genetic markers in forensic DNA analysis, routinely measured by capillary electrophoresis (CE) method based on the fragment length features. While, the massive parallel sequencing (MPS) technology could simultaneously target a large number of intriguing forensic STRs, bypassing the intrinsic limitations of amplicon size separation and accessible fluorophores in CE, which is efficient and promising for enabling the identification of forensic biological evidence. Here, we developed a novel MPS-based Forensic Analysis System Multiplecues SetB Kit of 133-plex forensic STR markers (52 STRs and 81 Y-STRs) and one Y-InDel (M175) based on multiplex PCR and single-end 400 bp sequencing strategy. This panel was subjected to developmental validation studies according to the SWGDAM Validation Guidelines. Approximately 2185 MPS-based reactions using 6 human DNA standards and 8 male donors were conducted for substrate studies (filter paper, gauze, cotton swab, four different types of FTA cards, peripheral venous blood, saliva, and exfoliated cells), sensitivity studies (from 2 ng down to 0.0625 ng), mixture studies (two-person DNA mixtures), PCR inhibitor studies (seven commonly encountered PCR inhibitors), species specificity studies (11 non-human species), and repeatability studies. Results of concordance studies (413 Han males and 6 human DNA standards) generated by STRait Razor and in-house Python scripts indicated 99.98% concordance rate in STR calling relative to CE for STRs between 41,900 genotypes at 100 STR markers. Moreover, the limitations of present studies, the nomenclature rules and forensic MPS applications were also described. In conclusion, the validation studies based on ~ 2200 MPS-based and ~ 2500 CE-based DNA profiles demonstrated that the novel MPS-based panel meets forensic DNA quality assurance guidelines with robust, reliable, and reproducible performance on samples of various quantities and qualities, and the STR nomenclature rules should be further regulated to integrate the inconformity between MPS-based and CE-based methods.

BGISEQ-500RS sequencing of a 448-plex SNP panel for forensic individual identification and kinship analysis

Next generation sequencing (NGS)-based single nucleotide polymorphism (SNP) genotyping is widely used in the field of forensics. SNP genotyping data from several NGS platforms have been published, but forensic application trials of DNA nanoball sequencing platforms have been very limited. In this work, we developed a 448-plex SNP panel on the BGISEQ-500RS platform. The sequencing metrics of a total of 261 samples that were sequenced with this panel are reported in detail. The average sequencing depth was 8373 × and the average heterozygosity of the 448-plex assay was 0.85. Sensitivity analysis showed that 325 SNPs were successfully genotyped with as little as 50 pg of genomic DNA, with the mean quality score of the sequencing data above Q30. Forensic parameters were calculated based on the data of 142 unrelated Chinese Han individuals and the combined matching probability was as low as 5.21 × 10^-101. Kinship analyses based on experiments and computer simulations showed that the 448-panel was as effective as the ForenSeq™ DNA Signature Prep Kit for second-degree kinship identification, and when the two panels were merged, the related pairs were almost completely distinguished from unrelated pairs. The 448-plex SNP panel on the BGISEQ-500RS platform provides a powerful tool for forensic individual identification and kinship analysis.