Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 27 September 2024
Sec. Statistical Genetics and Methodology
This article is part of the Research Topic Forensic Investigative Genetic Genealogy and Fine-Scale Structure of Human Populations, Volume II View all articles

Human complex mixture analysis by “FD Multi-SNP Mixture Kit”

Anqi ChenAnqi Chen1Lun Li,Lun Li2,3Junfei Zhou,Junfei Zhou2,3Tiantian Li,Tiantian Li2,3Chunyan YuanChunyan Yuan4Hai Peng,
Hai Peng2,3*Chengtao Li,
Chengtao Li1,4*Suhua Zhang,
Suhua Zhang1,4*
  • 1Institute of Forensic Science, Fudan University, Shanghai, China
  • 2Institute for Systems Biology, Jianghan University, Wuhan, Hubei, China
  • 3School of Life Sciences, Jianghan University, Wuhan, Hubei, China
  • 4Shanghai Key Laboratory of Forensic Medicine, Shanghai Forensic Service Platform, Ministry of Justice, Academy of Forensic Science, Shanghai, China

Introduction: Multiple linked single nucleotide polymorphisms (SNPs) have shown potential in personal identification and mixture detection. However, the limited number of marker and sequencing errors have obstructed accurate DNA typing.

Methods: To develop more candidate loci, the diversity value (D-value) was introduced as a new parameter for screening the novel polymorphic multiple linked-SNP markers, referred to as multi-SNP. In this study, a “FD Multi-SNP Mixture Kit” comprising 567 multi-SNPs was developed for mixture detection. Additionally, a new computational error correction method was applied as a quality control approach for sequencing data.

Results: The results demonstrated higher typing success rates than the conventional CE typing method. For single-source DNA, approximately 70–80 loci were detected with a DNA input of 0.009765625 ng. More than 65% of the minor alleles were distinguishable at 1 ng DNA with a frequency of 0.5% in 2- to 4-person mixtures.

Conclusion: This study offers a polymorphic and high-resolution detection method for DNA genotyping and complex mixture detection, providing an alternative strategy for addressing challenging mixed DNA traces.

1 Introduction

Mixture identification is a demanding task in the field of DNA analysis. The detection and analysis of complex mixtures have been among the most challenging problems faced by forensic DNA laboratories (Butler, 2023). Over the past decade, researchers have conducted studies on the collection, detection, and analysis of mixed samples, achieving some success in identifying two-person mixed samples. However, the examination and analysis of complex mixtures using existing techniques remain unsatisfactory in terms of their efficiency.

The conventional capillary electrophoresis-based short tandem repeat (CE-STR) typing method is the most commonly used approach for DNA identification (Chen et al., 2021a). However, STR profiles of mixed samples often show the coexistence and superposition of alleles from each contributor (Chen et al., 2021b; Chen et al., 2020). The generation of stutter, imbalance of alleles, random amplification effects, and interpretation of background signals all affect the possible typing combinations, making it more challenging for forensic scientists (Wu et al., 2019; Chen et al., 2022). Studies have shown that the STR profile might be incomplete if the minor contributor’s DNA proportion is less than 5%–20% (Chen et al., 2022; Zheng et al., 2019; Chen et al., 2019), which has limited its application in the identification of low template DNA. The International Society of Forensic Genetics (ISFG) also states that detection using autosomal STR locus is not completely reliable when the proportion of the minor contributor in a mixed two-person sample is less than 5% (Crespillo et al., 2014).

Compared to STR, single nucleotide polymorphism (SNP) is an attractive genetic marker because of its simplistic and easily readable typing results, and the lack of stutter peaks in the typing process (Babol-Pokora and Berent, 2008; Freire-Aradas et al., 2012; Brown et al., 2017). Although the SNP loci have limitations in their discriminatory power, several studies revealed that SNPs have some usefulness in the analysis of mixed samples (Isaacson et al., 2015; Kaur et al., 2014). Kidd et al. (Kidd et al., 2017; Kidd and Speed, 2015) first introduced the concept of microhaplotype loci (microhaps, MHs), which contain two or more closely linked SNPs. MHs provide significantly more information than single SNP markers (Aly and Sabri, 2015). They have demonstrated their utility in various fields, such as identifying missing persons, determining paternity, and detecting mixed samples. The innovations in haplotype phasing methods are essential to the effectiveness of MH analysis. Statistical phasing allows for the inference of haplotypes from population data (Kidd et al., 2022; Wu et al., 2021), family-based phasing leverages known familial relationships to improve accuracy (Wei et al., 2024), and long-read physical phasing directly sequences longer DNA fragments, preserving the linkage between variants (Kureshi et al., 2020). These advances collectively enhance the resolution and reliability of MH analysis, making it a powerful tool for complex genetic investigations. To date, an increasing number of MHs have been identified, but distinguishing minor alleles with frequencies below 1.5% remains challenging (Tao et al., 2022). Furthermore, the error rates of Illumina sequencing machines are approximately 0.1% per bases sequenced (Stoler and Nekrutenko, 2021), which suggests the low-abundant signals are likely to be masked. To further address the challenges mentioned above, it is necessary to develop methods for marker screening as well as accurate allele calling.

In this study, we developed a new solution named “FD Multi-SNP Mixture Kit”, which screened a specific category of SNP-based genetic markers and applied a novel strategy to call SNP alleles. Unlike the classic MHs selection approach (Kidd et al., 2017; Kidd and Speed, 2015), the new markers (referred to as multi-SNP) include combinations of two or more SNP loci that possess the D (diversity) value within 75 base pairs (bp). The good polymorphism and short amplification fragments might make them suitable for forensic samples. In addition to marker screening, new strategies for analyzing multi-SNP marker data were employed to decode the identity information of minor contributors. The results demonstrated that the approach was highly sensitive and effective for analyzing mixed samples.

2 Materials and methods

2.1 Genome-wide screening of multi-SNPs

The multi-SNPs markers were selected in accordance with our previous research with a few minor modifications (Fang. et al., 2021). The SNPs were obtained from the publicly accessible whole-genome sequence data of 91 Chinese Han individuals from the 1,000 Genomes Project (BioProject accession: PRJEB11005). A window size of 75 bp was chosen to accommodate complex DNA characteristics, such as trace amounts or high degradation. To select highly discriminatory multi-SNPs, the polymorphism within the window was evaluated using the D-value, which was calculated by t⁄((N¦2), where (N¦2) represented the total number of pairwise combinations that could be made within a population of N individuals, and t denoted the number of distinguishable pairings based on at least two single nucleotide variations found within the window. The windows susceptible to amplification, sequencing, and computational analysis errors, such as sequences comprising PolyN or tandem repeats, were filtered out. Only windows with a D-value of ≥ 0.6 were submitted to https://ampliseq.com for multiplex primer design. The amplicon length was set to be below 140 bp so that a single sequence read could cover the entire amplicon.

2.2 Library construction and sequencing

For each sample, the sequencing library was constructed using 5 μL of DNA following the protocol of the MGIEasy Universal DNA Library Prep Set (MGI-Tech, China), with the number of PCR cycles set to 28. Ten additional PCR cycles were then performed to anneal eight-nucleotide barcode sequences specific to each sample for multiplexing during sequencing. Finally, the pooled libraries were sequenced on the Illumina NovaSeq X platform to obtain 150 bp paired-end reads.

2.3 DNA samples

A total of 409 unrelated individuals from the Chinese Han population were recruited for the study. All samples were collected upon the approval of the Ethics Committee at the Academy of Forensic Science, Ministry of Justice, China. Each participant provided written informed consent to participate in the study. Genomic DNA (gDNA) samples were extracted from 200 μL of peripheral blood by using a QIAamp DNA Blood Kit (Qiagen, Netherlands). The gDNA was extracted according to the manufacturer’s instructions and quantified using a Qubit fluorometer (Life Technologies, United States). All extracted DNA was stored at −80 °C until use.

Two single-source DNA samples (one male and one female) with known genotypes were employed for the sensitivity study of the “FD Multi-SNP Mixture Kit”. The gDNA samples were diluted with molecular biology grade water (Phenix, United States). Ten different DNA inputs were prepared as follows: 5 ng, 2.5 ng; 1.25 ng, 0.625 ng, 0.3125 ng, 0.15625 ng, 0.078125 ng, 0.0390625 ng, 0.01953125 ng and 0.00976563 ng. In total, 60 samples (2 single-source DNA samples × 10 input amounts × 3 replicates) were sequenced in the sensitivity study.

In this study, two batches of artificial DNA mixtures were prepared. The first batch was prepared to compare the conventional CE-STR method and the next-generation sequencing (NGS) method of the “FD Multi-SNP Mixture Kit”. This batch included thirteen 2-person mixtures, eight 3-person mixtures, six 4-person mixtures and ten 5-person mixtures (Supplementary Table S3). The second batch was prepared to evaluate the mixture detection performance of the “FD Multi-SNP Mixture Kit”. It included seventy-five contrived mixture samples in total: twenty-five 2-person mixtures, twenty-five 3-person mixtures, twenty-three 4-person mixtures, one 5-person mixture and one 10-person mixture (Table 1). All participants were initially tested for all multi-SNPs using the “FD Multi-SNP Mixture Kit” and the allele typing results were obtained.

Table 1
www.frontiersin.org

Table 1. Summary of the DNA mixtures prepared for the mixture study (n = 75).

2.4 Bioinformatics

For each sample, the raw reads were mapped to the human reference genome using bowtie2 (Langmead and Salzberg, 2012) (version 2.1.0), and the unmapped or partially mapped reads were discarded in the further studies. For the fully mapped read, the nucleotide sequence spanning a multi-SNP locus was taken as its allele.

To lessen technical errors, several measures were employed as follows:

(1) Any variants located within or two base pairs from repeat sequences, consecutive mismatches, or indels were ignored to reduce the errors due to improper alignment. Since sequencing errors happen independently between reads, it’s unlikely that both reads of a pair have incorrect but identical alleles. Therefore, only the paired-end reads with identical alleles were retained to minimize false alleles brought about by sequencing errors.

(2) An iterative procedure was performed to identify true alleles from PCR artifacts: 1) All detected alleles in one multi-SNP locus were gathered, and the allele with the greatest read number was designated as the major allele. 2) For the remaining alleles (denoted as candidate alleles), they were temporary presumed as products of the major allele (e.g., PCR errors). The reads number of each candidate allele was modeled as a binomial distribution of total read counts of the major and candidate alleles and the error rate (e). If the counts number (c) of a candidate allele was excessive, resulting a significantly low P (X ≥ c), the null hypothesis that the allele was merely an artifact should be rejected. When more than one candidate alleles were tested, the multiple test correction was conducted, and the candidate alleles with false discovery rate (FDR) < 0.5% were deemed as the true minor alleles. 3) In the subsequent iterations, the authenticity of the remaining candidate alleles was assessed similarly to the first round, with the exception that for each candidate allele, the true allele that was the most similar was taken as the template of the erroneous amplification. The procedure terminated when no more minor alleles were found.

The PCR error rates may be affected by the nucleotide composition of locus, thus locus-specific error rates should be employed in the statistical models. Besides, since the errors happened independently across reads, the odds of multiple errors occurring within one read sequence was lower than that of a single one, making it necessary to use an error rate specific to on the number of single nucleotide variants (n) between the candidate allele and true allele. All alleles except the two alleles with most supporting reads at the multi-SNP loci of individual samples were produced by technical errors, and the variants on these alleles were regarded as false. In this study, for a multi-SNP locus l, an error profile was estimated from 409 Chinese Han individual samples, and the parameter el(n) was the average ratio of reads that allele bearing n false SNPs.

For each individual sample, loci with major allele read depths greater than 50 were further genotyped. The homozygous locus was recognized if the minor alleles with α values (the ratio of supporting read count of the minor allele to that of the major allele) less than 0.05, and the major allele was taken as the genotype. In contrary, if a locus had only one allele with an α value between 0.05 and 0.2, it was determined as heterozygous with the major allele and the minor allele as its genotype. Otherwise, the locus failed genotyping and was discarded from further studies.

2.5 Statistical analysis

The ACR value at heterozygous genotypes was calculated as the smaller number of depths at one allele compared to the other. Allele frequencies were calculated using the counting method. The other forensic statistical parameters, including expected heterozygosity (He), observed heterozygosity (Ho), match probability (MP), discrimination power (DP) and probability of exclusion (PE), were calculated using Power Stats v12 software (Promega, United States). Hardy–Weinberg Equilibrium (HWE) was estimated using Arlequin 3.533 software. Detection rate was applied to assess the concordance for both the single- and multi-sourced DNA samples, calculated as the ratio of the number of alleles matched to the total number of alleles detected.

2.6 Conventional CE-STR typing

The STR profiles of the first batch of artificial DNA mixtures were determined using a PowerPlex®21 kit (Promega, United States), which includes 20 autosomal STRs and one sex-linked polymorphism Amelogenin. Fluorescence multiplex PCR was performed according to the manufacturer’s instructions and the amplified products were separated on a 3100 ABI Prism Genetic Analyzer (Applied Biosystems, United States). Genotyping was per-formed using GeneMapper software (Applied Biosystems, United States), and all results were checked by an experienced technician.

3 Results and discussion

The MHs, which exhibit a high degree of polymorphism, have been shown to have significant potential in mixture deconvolutions (Chen et al., 2019). To identify more polymorphic multi-allelic markers, we proposed the D-value as a new parameter. Using this approach, a total of 519 multi-SNP markers were identified, including 157 two-linked SNPs, 246 three-linked SNPs, 91 four-linked SNPs and 25 five-linked SNPs. To further enhance polymorphic diversity, an additional 48 single SNP markers were also incorporated (Table 2). As a result, a toolkit named the “FD Multi-SNP Mixture Kit” was developed for high-throughput DNA profiling and mixture detection.

Table 2
www.frontiersin.org

Table 2. Loci constitution of the “FD multi-SNP Mixture Kit”.

3.1 Sequencing performance of the “FD Multi-SNP mixture kit”

To evaluate the performance of the panel, we sequenced the DNA samples of 409 unrelated individuals from the Chinese Han population. For single-source sample, the panel achieved an allele coverage ratio of 98.21%, indicating that there were sufficient loci for forensic identification. The average sequencing depth was 1150.87 ± 716.64×, ranging from 99.19 ± 35.53× to 5841.08 ± 929.49×. The depth of CHM527 was the highest, followed by CHM277 (3917.78 ± 2041.55×), CHM511 (3856.50 ± 1857.29×), CHM528 (3712.10 ± 1723.12×) and CHM101 (3629.66 ± 1502.32×). The depths of CHM578, CHM044, CHM303, CHM416 and CHM085 were the lowest, with depths below 200×. Although sequencing depths varied across loci, the lowest depth remained at approximately 100× (Figure 1A). Jiang et al. (2019) argued that deeper sequencing depths lead to higher typing accuracy. In this study, the sequencing depths of the markers were comparable to those in published SNP-based sequencing studies (Guo et al., 2017; Tao et al., 2021), ensuring good data quality for further analysis.

Figure 1
www.frontiersin.org

Figure 1. Sequencing depth and ACR of the “FD Multi-SNP Mixture Kit”. (A) Coverage depth. (B) Distribution of loci across different ACR ranges.

The average ACR of the panel was 0.92 ± 0.05, which is excellent according to the current knowledge (Figure 1B). The forensic community defines a DNA mixture as a biological sample derived from two or more contributors (Bieber et al., 2016). Thus, a good intralocus balance is more favorable for mixture deconvolution (Pang et al., 2020). Due to the characteristics of multi-SNP markers, a cold case from several decades ago was resolved (Chen et al., 2024). The case involved a sample from a campstool, which was suspected to be the potential weapon in a homicide. DNA extraction and STR profiling of the campstool revealed that the DNA was a mixture. However, the results from CE-STR were fraught with uncertainty and controversy. In contrast, analysis using NGS with the “FD Multi-SNP Mixture Kit” indicated that the probability of the suspect’s DNA being present on the campstool was as high as 99.999159%. This result contradicted the suspect’s statement and was instrumental in solving the case, demonstrating the effectiveness of the FD Multi-SNP Mixture Kit for degraded and trace DNA mixtures in forensic investigations.

3.2 Marker detection performance of the “FD Multi-SNP mixture kit”

For approximately 30 years, STR markers with 3-4 bp repeats have been considered the conventional marker in forensic genetics. However, the stochastic effects of this method would likely to reduce the success rate of low-input DNA analysis (Nwawuba Stanley et al., 2020). To investigate whether multi-SNP markers are more suitable for mixture analysis, a back-to-back comparison was made between the standard CE-STR method and the NGS of “FD Multi-SNP Mixture Kit”. The result showed higher detection rates (93.76%–100%) with the NGS of “FD multi-SNP Mixture Kit”, while lower detection rates (33.93%–100%) were observed with the conventional approach (Figure 2A). The higher the resolution of detection, the more reliable the mixture deconvolution. In this study, the NGS method of “FD Multi-SNP Mixture Kit” accurately detected more alleles, suggesting that it may provide more opportunities for further interpretation of DNA mixture. In addition, the success rates of minor contributors were lower than those of major contributors, regardless of the typing method (Figure 2B, Additional file: Supplementary Table S3), further demonstrating that absolute amounts of DNA play a critical role in detection rates and sequencing performance of the “FD Multi-SNP Mixture Kit”.

Figure 2
www.frontiersin.org

Figure 2. Comparison of detection rates between the two methods for 37 mixed samples. (A) Comparison of the overall detection rates of the CE typing methods. (B) Typing accuracy of each contributor in the mixture samples (P1, P2, P3, P4, P5 and P6 represent six unrelated individuals).

3.3 Sensitivity performance of the “FD Multi-SNP mixture kit”

To further investigate the sensitivity of this panel, a series of diluted DNA samples were tested. The fluctuating trends were observed in the counts of loci detected from 5 ng, 2.5 ng, 1.25 ng and 0.625 ng of DNA input. Among these, approximately 95%–97% of the loci were correctly typed. The low standard deviations within the same DNA inputs showed comparable results in the triplicates. The number of loci detected decreased when the input amount was less than 0.625 ng. At 0.3125 ng DNA, more than 90% of the profiles could be typed correctly. However, many allelic dropouts began at 0.15625 ng DNA, and only 15% of the loci could be detected when the input was reduced to 0.00976563 ng. The abundance of loci in this panel provide more opportunities for DNA identification (Bruijns et al., 2018). For this panel, over 70 loci were detected with an input of 0.00976563 ng DNA (Figure 3A). It has been reported that 40 to 60 SNPs have a comparable power of 13–15 STRs (Butler et al., 2007), and the multi-SNP markers have more alleles di-allelic SNPs. Therefore, the remaining markers should have retained a comparable discrimination power of 20–23 STR loci.

Figure 3
www.frontiersin.org

Figure 3. Sensitivity study of the “FD Multi-SNP Mixture Kit”. (A) Counts of loci detected at different input amounts. (B) Sequencing depths and average ACR for different DNA input amounts.

In this study, sequencing depths of over 1,000× were measured when the input amount was above 0.3125 ng, and the ACR remained above 0.7 even when the DNA input was as low as 0.15625 ng. A dramatic decrease in depth and ACR occurred when the input was reduced to 0.0390625 ng. Under this circumstance, the average depth was no more than 300×, and the average ACR was only 0.5. To fully utilize the flow cells and ensure consistent sequencing depth for all samples, diluted libraries of variable concentration should be volumetrically pooled (Bruinsma et al., 2018). Based on the observed depth variations, the amount of 5ng to 0.625 ng of DNA could potentially be attributed to pipetting errors during library normalization. Compared to sequencing depth, the ACR values were less sensitive to the amount of DNA input. High quality ACR values (>0.7) could be maintained even when the input DNA was reduced to 0.15625 ng (Figure 3B). Balanced ACR is advantageous for distinguishing components from a DNA mixture (Pang et al., 2020), and the results demonstrated its potential applications in the detection and deconvolution of low-amount DNA.

3.4 Mixture detection performance of the “FD Multi-SNP mixture kit”

For single-source DNA samples, we observed a decrease in the number of loci detected as DNA input decreased (Figure 3A). The detection of loci is expected to correlate with the quantity of DNA input rather than the number of contributors (Figure 4A). The detection rate of loci was as high as 93.95% in a 10 ng DNA mixture, but decreased to 81.18% when the DNA input was reduced to 0.05 ng. This decrease aligns with the results of the sensitivity study, further emphasizing the direct correlation between the detection rate and the quantity of DNA input. These findings underscore that the number of contributors does not affect the detection rate, and that a higher DNA input can result in a larger number of detected loci.

Figure 4
www.frontiersin.org

Figure 4. Detection performance of the “FD multi-SNP Mixture Kit”. (A) Histogram of the detection rates in the mixtures with different numbers of contributors and different DNA inputs. The x-axis represents the number of individuals (left) and the total input amounts in the mixtures, while the y-axis shows the average detection rates of the samples. (B) Heatmap showing the detection rate for each of the DNA mixtures in the study (N = 75). (C) Detection rate of loci from contributors at different proportions and different input amounts. The x-axis represents the DNA amounts of the individuals, while the y-axis shows the average detection rates.

Higher input amounts have been shown to increase the detection rate, thereby providing more genetic markers for subsequent analysis. However, in mixture detection, the DNA amount of each contributor varies in proportion in each sample, leading to variations in typings even when the same DNA input is used. To minimize this difference, the contributor with the lowest DNA amount of was defined as the same person (individual B) in the unbalanced DNA mixtures. As expected, individual B consistently exhibited the lowest detection rate, except in Sample M37 (Figure 4B). Signals from low-contributing sources are prone to being masked by those from high contributors (van Oorschot et al., 2010). Similar to previous findings, more accurately typed loci were observed in the major contributors. Though researchers often describe detection rates in terms of component ratios, the significance of DNA input can be overlooked. Figure 4C demonstrates that while the ratio of components in mixtures remained constant, there was significant variation in detection rates. For example, individual B had a contribution proportion of 10% in 43 DNA mixtures, but the detection rate of the loci was inconsistent. These results suggest that, in addition to component proportions, another critical factor must influence contributor genotyping.

In the sensitivity analysis, the number of loci decreased with decreased DNA input (Figure 2B). The same trend in detection rate was observed in the mixture study. Within the range of 0.078 ng–0.156 ng DNA, approximately 99.68% of the contributor’s loci could be accurately typed, which was similar to the rate of 99.93% at 0.15625 ng DNA in the sensitivity study. For single source DNA samples, over 78.69% of the loci were correctly typed at an input of 0.00976563 ng. A comparable level (81.31%) was observed in the mixture study for contributors with DNA inputs between 0 and 0.001 ng (Figure 4C), indicating a similar detection limit between single- DNA and multi-sourced DNA samples. As detection success was closely linked to the DNA input of the minor contributor, we anticipated similar detection rates amongst the contributors in the balanced DNA mixtures. As expected, the contributors showed indistinguishable detection rates in the balanced 5- and 10-person DNA mixtures (Figure 4B).

3.5 Statistical parameters of the “FD multi-SNP mixture kit”

High-throughput SNP genotyping has been shown to be effective due to smaller amplicons, but the limited heterozygosity of SNP markers may make them less informative for certain purposes (Hares, 2015). Kidd et al. (Soundararajan et al., 2016) suggest that SNPs must contain multiple alleles to capture sufficient genetic diversity. This panel consists of 567 multi-SNPs, which is 30 times the size of the expanded Combined DNA Index System (CODIS) core loci, suggesting that the panel may have higher discrimination power. The 567 markers generated 1,855 unique allelotypes, and the corresponding allele frequencies are shown in Figure 1. Most of the loci had 2 to 4 unique allelotypes, with CHM020 being the most polymorphic locus with 16 unique genotypes, followed by CHM519 with 8 unique allelotypes. On the other hand, CHM462 and CHM085 showed poor polymorphism, with only one allelotype found among the 409 unrelated individuals. However, these loci accounted for only a small proportion of the panel (Additional file: Supplementary Table S1; Figure 5).

Figure 5
www.frontiersin.org

Figure 5. Allele frequencies of the 567 markers in the “FD multi-SNP Mixture Kit” within the Chinese Han population.

Other forensic statistical parameters, such as Hardy-Weinberg equilibrium (HWE), expected heterozygosity (He), observed heterozygosity (Ho), match probability (MP), discrimination power (DP), and power of exclusion (PE), were also determined. More than 98.06% (556/567) of the loci followed HWE, and 11 loci deviated from HWE after Bonferroni correction (p = 0.00008818). The average He and Ho were 0.55 ± 0.12 and 0.45 ± 0.12, respectively. The MP values ranged from 0.04 to 1, and DP values ranged from 0 to 0.96. The highest DP was found at CHM020 (0.96), and the highest PE was observed at CHM348 (1.00) (Additional file: Supplementary Table S2). Large panels are rich in genetic markers and are expected to provide more opportunities for genetic identification. However, large panels may also increase the likelihood of allelic dropouts (Shestak et al., 2021). In this study, the panel comprised 567 markers, with an average of 98.21% ± 1.44% of markers genotyped for each individual. Furthermore, 446 loci were stably detected in all 409 unrelated individuals, with combined power of exclusion (CPE), combined match probability (CMP), and total discrimination power (TDP) values of 1–0.00×10°, 0.00×10° and 1–0.00×10°, respectively, suggesting that this panel is powerful for the purpose of forensic identification.

4 Conclusion

In this study, our aim to develop novel genetic markers for the deconvolution of complex human DNA mixtures. We utilized the D value to evaluate the multi-SNP markers, which ultimately led to the creation of the “FD Multi-SNP Mixture Kit”. This kit features 567 multi-SNPs and is highly effective in DNA mixture detection. The high polymorphism and intra-locus balance across the loci may enhance the utility in mixture deconvolution. Our approach demonstrated exceptional sensitivity in analyzing both single-source and mixed DNA samples, underscoring its potential for identifying donors in low-yield samples.

Data availability statement

The datasets presented in this study can be found in NCBI Sequence Read Archive (SRA) repository, accession number PRJNA1163920 (https://www.ncbi.nlm.nih.gov/sra/PRJNA1163920).

Ethics statement

The studies involving humans were approved by Ethics Committee at the Academy of Forensic Science, Ministry of Justice, China. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

AC: Data curation, Formal Analysis, Investigation, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. LL: Data curation, Formal Analysis, Investigation, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. JZ: Investigation, Validation, Writing–review and editing. TL: Investigation, Validation, Writing–review and editing. CY: Data curation, Writing–review and editing. HP: Methodology, Project administration, Resources, Supervision, Writing–review and editing. CL: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Writing–review and editing. SZ: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by grants from the National Natural Science Fund of China (Nos 82072123 and 81930056).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1432378/full#supplementary-material

References

Aly, S. M., and Sabri, D. M. (2015). Next generation sequencing (NGS): a golden tool in forensic toolkit. Arch. Med. sadowej i Kryminol. 65, 260–271. doi:10.5114/amsik.2015.61029

PubMed Abstract | CrossRef Full Text | Google Scholar

Babol-Pokora, K., and Berent, J. (2008). SNP-minisequencing as an excellent tool for analysing degraded DNA recovered from archival tissues. Acta Biochim. Pol. 55, 815–819. doi:10.18388/abp.2008_3045

PubMed Abstract | CrossRef Full Text | Google Scholar

Bieber, F. R., Buckleton, J. S., Budowle, B., Butler, J. M., and Coble, M. D. (2016). Evaluation of forensic DNA mixture evidence: protocol for evaluation, interpretation, and statistical calculations using the combined probability of inclusion. BMC Genet. 17, 125. doi:10.1186/s12863-016-0429-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, H., Thompson, R., Murphy, G., Peters, D., La Rue, B., King, J., et al. (2017). Development and validation of a novel multiplexed DNA analysis system, InnoTyper® 21. Genetics 29, 80–99. doi:10.1016/j.fsigen.2017.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruijns, B., Tiggelaar, R., and Gardeniers, H. (2018). Massively parallel sequencing techniques for forensics: a review. Electrophoresis 39, 2642–2654. doi:10.1002/elps.201800082

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruinsma, S., Burgess, J., Schlingman, D., Czyz, A., Morrell, N., Ballenger, C., et al. (2018). Bead-linked transposomes enable a normalization-free workflow for NGS library preparation. BMC Genomics 19, 722. doi:10.1186/s12864-018-5096-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Butler, J. M. (2023). Recent advances in forensic biology and forensic DNA typing: INTERPOL review 2019-2022. Forensic Sci. Int. Synerg 6, 100311. doi:10.1016/j.fsisyn.2022.100311

PubMed Abstract | CrossRef Full Text | Google Scholar

Butler, J. M., Coble, M. D., and Vallone, P. M. (2007). STRs vs. SNPs: thoughts on the future of forensic DNA testing. Forensic Sci. Med. pathology 3, 200–205. doi:10.1007/s12024-007-0018-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, A., Xiong, L., Qu, Y., Xi, S., Tao, R., Li, C., et al. (2022). Opportunity of next-generation sequencing-based short tandem repeat system for tumor source identification. Front. Oncol. 12, 800028. doi:10.3389/fonc.2022.800028

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, A., Yang, Y., Yang, Q., Tao, R., Chen, C., Zhang, S., et al. (2021a). Analytical validation of an RI sample cartridge with the RapidHIT® ID system. Int. J. Leg. Med. 135, 1257–1265. doi:10.1007/s00414-021-02553-8

CrossRef Full Text | Google Scholar

Chen, A., Zhang, S., Li, J., Ji, C., Chen, J., and Li, C. (2020). Detecting genetic hypermutability of gastrointestinal tumor by using a forensic STR kit. Front. Med. 14, 101–111. doi:10.1007/s11684-019-0698-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, A., Zhang, S., Xiong, L., Xi, S., Tao, R., Chen, C., et al. (2021b). Investigation of an alternative marker for hypermutability evaluation in different tumors. Genes (Basel) 12, 197. doi:10.3390/genes12020197

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Chen, A., Tao, R., Zhu, R., Zhang, H., You, X., et al. (2024). Solution to a case involving the interpretation of trace degraded DNA mixtures. Int. J. Leg. Med. doi:10.1007/s00414-024-03302-3

CrossRef Full Text | Google Scholar

Chen, P., Deng, C., Li, Z., Pu, Y., Yang, J., Yu, Y., et al. (2019). A microhaplotypes panel for massively parallel sequencing analysis of DNA mixtures. Genetics 40, 140–149. doi:10.1016/j.fsigen.2019.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Crespillo, M., Barrio, P. A., Luque, J. A., Alves, C., Aler, M., Alessandrini, F., et al. (2014). GHEP-ISFG collaborative exercise on mixture profiles of autosomal STRs (GHEP-MIX01, GHEP-MIX02 and GHEP-MIX03): results and evaluation. Genetics 10, 64–72. doi:10.1016/j.fsigen.2014.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang., Z., Li, L., Zhou, J., You, A., Gao, L., Li, T., et al. (2021). Multiple nucleotide polymorphism DNA markers for the accurate evaluation of genetic variations. bioRxiv, 434561. doi:10.1101/2021.03.09.434561

CrossRef Full Text | Google Scholar

Freire-Aradas, A., Fondevila, M., Kriegel, A. K., Phillips, C., Gill, P., Prieto, L., et al. (2012). A new SNP assay for identification of highly degraded human DNA. Forensic Sci. Int. Genet. 6, 341–349. doi:10.1016/j.fsigen.2011.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, F., Yu, J., Zhang, L., and Li, J. (2017). Massively parallel sequencing of forensic STRs and SNPs using the Illumina® ForenSeq™ DNA signature Prep kit on the MiSeq FGx™ forensic genomics system. Forensic Sci. Int. Genet. 31, 135–148. doi:10.1016/j.fsigen.2017.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hares, D. R. (2015). Selection and implementation of expanded CODIS core loci in the United States. Genetics 17, 33–34. doi:10.1016/j.fsigen.2015.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Isaacson, J., Schwoebel, E., Shcherbina, A., Ricke, D., Harper, J., Petrovick, M., et al. (2015). Robust detection of individual forensic profiles in DNA mixtures. Forensic Sci. Int. Genet. 14, 31–37. doi:10.1016/j.fsigen.2014.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Jiang, Y., Wang, S., Zhang, Q., and Ding, X. (2019). Optimal sequencing depth design for whole genome re-sequencing in pigs. BMC Bioinforma. 20, 556. doi:10.1186/s12859-019-3164-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaur, N., Fonneløp, A. E., and Egeland, T. (2014). Regression models for DNA-mixtures. Forensic Sci. Int. Genet. 11, 105–110. doi:10.1016/j.fsigen.2014.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kidd, K. K., Evsanaa, B., Togtokh, A., Brissenden, J. E., Roscoe, J. M., Dogan, M., et al. (2022). North Asian population relationships in a global context. Sci. Rep. 12, 7214. doi:10.1038/s41598-022-10706-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kidd, K. K., and Speed, W. C. (2015). Criteria for selecting microhaplotypes: mixture detection and deconvolution. Investig. Genet. 6, 1. doi:10.1186/s13323-014-0018-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kidd, K. K., Speed, W. C., Pakstis, A. J., Podini, D. S., Lagacé, R., Chang, J., et al. (2017). Evaluating 130 microhaplotypes across a global set of 83 populations. Genetics 29, 29–37. doi:10.1016/j.fsigen.2017.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kureshi, A., Li, J., Wen, D., Sun, S., Yang, Z., and Zha, L. (2020). Construction and forensic application of 20 highly polymorphic microhaplotypes. R. Soc. Open Sci. 7, 191937. doi:10.1098/rsos.191937

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi:10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Nwawuba Stanley, U., Mohammed Khadija, A., Bukola, A. T., Omusi Precious, I., and Ayevbuomwan Davidson, E. (2020). Forensic DNA profiling: autosomal short tandem repeat as a prominent marker in crime investigation. Malays J. Med. Sci. 27, 22–35. doi:10.21315/mjms2020.27.4.3

CrossRef Full Text | Google Scholar

Pang, J. B., Rao, M., Chen, Q. F., Ji, A. Q., Zhang, C., Kang, K. L., et al. (2020). A 124-plex microhaplotype panel based on next-generation sequencing developed for forensic applications. Sci. Rep. 10, 1945. doi:10.1038/s41598-020-58980-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shestak, A. G., Bukaeva, A. A., Saber, S., and Zaklyazminskaya, E. V. (2021). Allelic dropout is a common phenomenon that reduces the diagnostic yield of PCR-based sequencing of targeted gene panels. Front. Genet. 12, 620337. doi:10.3389/fgene.2021.620337

PubMed Abstract | CrossRef Full Text | Google Scholar

Soundararajan, U., Yun, L., Shi, M., and Kidd, K. K. (2016). Minimal SNP overlap among multiple panels of ancestry informative markers argues for more international collaboration. Genetics 23, 25–32. doi:10.1016/j.fsigen.2016.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoler, N., and Nekrutenko, A. (2021). Sequencing error profiles of Illumina sequencing instruments. Nar. Genom Bioinform 3, lqab019. doi:10.1093/nargab/lqab019

PubMed Abstract | CrossRef Full Text | Google Scholar

Tao, R., Wang, S., Chen, A., Xia, R., Zhang, X., Yang, Q., et al. (2021). Parallel sequencing of 87 STR and 294 SNP markers using the prototype of the SifaMPS panel on the MiSeq FGx™ system. Forensic Sci. Int. Genet. 52, 102490. doi:10.1016/j.fsigen.2021.102490

PubMed Abstract | CrossRef Full Text | Google Scholar

Tao, R., Yang, Q., Xia, R., Zhang, X., Chen, A., Li, C., et al. (2022). A sequence-based 163plex microhaplotype assay for forensic DNA analysis. Front. Genet. 13, 988223. doi:10.3389/fgene.2022.988223

PubMed Abstract | CrossRef Full Text | Google Scholar

van Oorschot, R. A., Ballantyne, K. N., and Mitchell, R. J. (2010). Forensic trace DNA: a review. Investig. Genet. 1, 14. doi:10.1186/2041-2223-1-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Y., Zhu, Q., Wang, H., Cao, Y., Li, X., Zhang, X., et al. (2024). Pairwise kinship inference and pedigree reconstruction using 91 microhaplotypes. Forensic Sci. Int. Genet. 72, 103090. doi:10.1016/j.fsigen.2024.103090

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, L., Chu, X., Zheng, J., Xiao, C., Zhang, Z., Huang, G., et al. (2019). Targeted capture and sequencing of 1245 SNPs for forensic applications. Genetics 42, 227–234. doi:10.1016/j.fsigen.2019.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, R., Li, H., Li, R., Peng, D., Wang, N., Shen, X., et al. (2021). Identification and sequencing of 59 highly polymorphic microhaplotypes for analysis of DNA mixtures. Int. J. Leg. Med. 135, 1137–1149. doi:10.1007/s00414-020-02483-x

CrossRef Full Text | Google Scholar

Zheng, H., Tao, R., Zhang, J., Zhang, J., Wang, S., Yang, Z., et al. (2019). Development and validation of a novel SiFaSTR(TM) 23-plex system. Electrophoresis 40, 2644–2654. doi:10.1002/elps.201900045

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: genetic markers, multi-SNPs, individual identification, mixture detection, forensic

Citation: Chen A, Li L, Zhou J, Li T, Yuan C, Peng H, Li C and Zhang S (2024) Human complex mixture analysis by “FD Multi-SNP Mixture Kit”. Front. Genet. 15:1432378. doi: 10.3389/fgene.2024.1432378

Received: 14 May 2024; Accepted: 05 September 2024;
Published: 27 September 2024.

Edited by:

Hongyu Sun, Sun Yat-sen University, China

Reviewed by:

Guanglin He, Sichuan University, China
Bofeng Zhu, Southern Medical University, China

Copyright © 2024 Chen, Li, Zhou, Li, Yuan, Peng, Li and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hai Peng, penghai138@163.com; Chengtao Li, lichengtao@fudan.edu.cn; Suhua Zhang, zhang_suhua@fudan.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.