Solution to a case involving the interpretation of trace degraded DNA mixtures.

DNA mixture analysis poses a significant challenge in forensic genetics, particularly when dealing with degraded and trace amount DNA samples. Multi-SNPs (MNPs) are genetic markers similar to microhaplotypes but with smaller molecular sizes (< 75 bp), making them theoretically more suitable for analyzing degraded and trace amount samples. In this case report, we investigated a cold case involving a campstool stored for over a decade, aiming to detect and locate the suspect's DNA. We employed both conventional capillary electrophoresis-based short tandem repeat (CE-STR) analysis and next-generation sequencing-based multi-SNP (NGS-MNP) analysis. The typing results and deconvolution of the mixed CE-STR profiles were inconclusive regarding the presence of the suspect's DNA in the mixed samples. However, through NGS-MNP analysis and presence probability calculations, we determined that the suspect's DNA was present in the samples from Sect. 4-1 with a probability of 1-8.41 × 10- 6 (99.999159%). This evidence contradicted the suspect's statement and aided in resolving the case. Our findings demonstrate the significant potential of MNP analysis for examining degraded and trace amount DNA mixtures in forensic investigations.

A Public Mid-Density Genotyping Platform for Hexaploid Sweetpotato (Ipomoea batatas [L.] Lam).

Small public breeding programs focused on specialty crops have many barriers to adopting technology, particularly creating and using genetic marker panels for genomic-based decisions in selection. Here, we report the creation of a DArTag panel of 3120 loci distributed across the sweetpotato (Ipomoea batatas [L.] Lam) genome for molecular-marker-assisted breeding and genomic prediction. The creation of this marker panel has the potential to bring cost-effective and rapid genotyping capabilities to sweetpotato breeding programs worldwide. The open access provided by this platform will allow the genetic datasets generated on the marker panel to be compared and joined across projects, institutions, and countries. This genotyping resource has the power to make routine genotyping a reality for any breeder of sweetpotato.

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.

Unique molecular identifier-based amplicon sequencing of microhaplotypes for background noise mitigation.

Microhaplotypes (MHs), comprising two or more single-nucleotide polymorphisms in a short fragment, are promising forensic markers owing to their remarkable polymorphic nature. Several studies have demonstrated the utility of MHs through massively parallel sequencing (MPS). Nevertheless, the background noise level associated with MHs in MPS, which imposes a practical detection limit for the system, remains uninvestigated. Currently, unique molecular identifier (UMI) systems are known to effectively mitigate background noise by tracking original DNA molecules and facilitating PCR and MPS error corrections. Hence, this study aimed to design a UMI-based amplicon sequencing system, designated MH-UMIseq, which can amplify 46 MHs simultaneously and generate MPS libraries in four steps: barcoding PCR, nuclease reaction, boosting PCR, and indexing PCR. The performance of the MH-UMIseq system was evaluated using the Illumina NextSeq 550 and MiniSeq systems with 31 sets for 5 ng, 1 ng, and 200 pg of input DNA. The fgbio toolkit was used in conjunction with STRait Razor 3.0 and Visual Microhap to analyze the UMI data on MHs. The corresponding average not suppressed noise proportion of MH-UMIseq were 0.1 %, 0.3 %, and 0.7 % for 5 ng, 1 ng, and 200 pg of DNA, respectively, which substantially suppressed the background noise for more than 1 ng of DNA. Interestingly, the proportion of not suppressed noise in MH-UMIseq notably decreased as the amount of input DNA increased. The number of UMI families was proportional to the copy number of the template DNA and closely correlated with the system resolution. Therefore, the resolution of MH-UMIseq system is expected to be higher than that of conventional MPS for the deconvolution of mixtures containing more than 1 ng of DNA.

Establishment and Application of a 42-plex Microhaplotype Assay in Forensic Genetics.

OBJECTIVES: To establish and forensically verify a 42 microhaplotypes (mircohaps, MHs) multiplex assay system based on next-generation sequencing (NGS), and to explore the application value of this system in the practice of forensic genetics. METHODS: A total of 42 highly polymorphic MHs were selected from previous studies, and sequenced by the MiSeq FGxTM platform to verify the repeata-bility, sensitivity, specificity, stability, and mixture analysis ability of the detection system. Through population genetic investigation of 102 unrelated Chinese Han individuals in Liyang City, Jiangsu Province, China, the application value of this system in forensic genetics was evaluated. RESULTS: The sequencing repeatability of the 42-plex MHs assay was 100% and the sensitivity was as low as 0.062 5 ng. The system had the ability to withstand the interference of indigo (≤2 500 ng/µL), humic acid (≤9 ng/µL), hemoglobin(≤20 µmol), and urea (≤200 ng/µL) and to detect mixtures of 2 people (1∶19), 3 people (1∶1∶9) and 4 people (1∶1∶1∶9). Based on 102 individual data, the combined power of discrimination and the combined power of exclusion were 1-3.45×10-30 and 1-3.77×10-11, respectively, and the average effect value of alleles was 2.899. CONCLUSIONS: The 42-plex MHs assay was successfully established in this study and this system has high repeatability and sensitivity, good anti-jamming ability and mixture analysis ability. The 42 MHs are highly polymorphism and have good application value in individual identification and paternity testing.

A preliminary study on identification of the blood donor in a body fluid mixture using a novel compound genetic marker blood-specific methylation-microhaplotype.

Blood-containing mixtures are frequently encountered at crime scenes involving violence and murder. However, the presence of blood, and the association of blood with a specific donor within these mixtures present significant challenges in forensic analysis. In light of these challenges, this study sought to address these issues by leveraging blood-specific methylation sites and closely linked microhaplotype sites, proposing a novel composite genetic marker known as "blood-specific methylation-microhaplotype". This marker was designed to the detection of blood and the determination of blood donor within blood-containing mixtures. According to the selection criteria mentioned in the Materials and Methods section, we selected 10 blood-specific methylation-microhaplotype loci for inclusion in this study. Among these loci, eight exhibited blood-specific hypomethylation, while the remaining two displayed blood-specific hypermethylation. Based on data obtained from 124 individual samples in our study, the combined discrimination power (CPD) of these 10 successfully sequenced loci was 0.999999298. The sample allele methylation rate (Ram) was obtained from massive parallel sequencing (MPS), which was defined as the proportion of methylated reads to the total clustered reads that were genotyped to a specific allele. To develop an allele type classification model capable of identifying the presence of blood and the blood donor, we used the Random Forest algorithm. This model was trained and evaluated using the Ram distribution of individual samples and the Ram distribution of simulated shared alleles. Subsequently, we applied the developed allele type classification model to predict alleles within actual mixtures, trying to exclude non-blood-specific alleles, ultimately allowing us to identify the presence of blood and the blood donor in the blood-containing mixtures. Our findings demonstrate that these blood-specific methylation-microhaplotype loci have the capability to not only detect the presence of blood but also accurately associate blood with the true donor in blood-containing mixtures with the mixing ratios of 1:29, 1:19, 1:9, 1:4, 1:2, 2:1, 7:1, 8:1, 31:1 and 36:1 (blood:non-blood) by DNA mixture interpretation methods. In addition, the presence of blood and the true blood donor could be identified in a mixture containing four body fluids (blood:vaginal fluid:semen:saliva = 1:1:1:1). It is important to note that while these loci exhibit great potential, the impact of allele dropouts and alleles misidentification must be considered when interpreting the results. This is a preliminary study utilising blood-specific methylation-microhaplotype as a complementary tool to other well-established genetic markers (STR, SNP, microhaplotype, etc.) for the analysis in blood-containing mixtures.

Performance of a 74-Microhaplotype Assay in Kinship Analyses.

Microhaplotypes (MHs) consisting of multiple SNPs and indels on short stretches of DNA are new and interesting loci for forensic genetic investigations. In this study, we analysed 74 previously defined MHs in two of the populations that our laboratory provides with forensic genetic services, Danes and Greenlanders. In addition to the 229 SNPs that originally made up the 74 MHs, 66 SNPs and 3 indels were identified in the two populations, and 45 of these variants were included in new definitions of the MHs, whereas 24 SNPs were considered rare and of little value for case work. The average effective number of alleles (Ae) was 3.2, 3.0, and 2.6 in Danes, West Greenlanders, and East Greenlanders, respectively. High levels of linkage disequilibrium were observed in East Greenlanders, which reflects the characteristics of this population that has a small size, and signs of admixture and substructure. Pairwise kinship simulations of full siblings, half-siblings, first cousins, and unrelated individuals were performed using allele frequencies from MHs, STRs and SNPs from Danish and Greenlandic populations. The MH panel outperformed the currently used STR and SNP marker sets and was able to differentiate siblings from unrelated individuals with a 0% false positive rate and a 1.1% false negative rate using an LR threshold of 10,000 in the Danish population. However, the panel was not able to differentiate half-siblings or first cousins from unrelated individuals. The results generated in this study will be used to implement MHs as investigative markers for relationship testing in our laboratory.

Improving DNA mixtures analysis using compound markers composed of InDels and SNPs screened from the whole genome with next-generation sequencing.

Next-generation sequencing (NGS) allows for better identification of insertion and deletion polymorphisms (InDels) and their combination with adjacent single nucleotide polymorphisms (SNPs) to form compound markers. These markers can improve the polymorphism of microhaplotypes (MHs) within the same length range, and thus, boost the efficiency of DNA mixture analysis. In this study, we screened InDels and SNPs across the whole genome and selected highly polymorphic markers composed of InDels and/or SNPs within 300 bp. Further, we successfully developed and evaluated an NGS-based panel comprising 55 loci, of which 24 were composed of both SNPs and InDels. Analysis of 124 unrelated Southern Han Chinese revealed an average effective number of alleles (Ae ) of 7.52 for this panel. The cumulative power of discrimination and cumulative probability of exclusion values of the 55 loci were 1-2.37 × 10-73 and 1-1.19 × 10-28 , respectively. Additionally, this panel exhibited high allele detection rates of over 97% in each of the 21 artificial mixtures involving from two to six contributors at different mixing ratios. We used EuroForMix to calculate the likelihood ratio (LR) and evaluate the evidence strength provided by this panel, and it could assess evidence strength with LR, distinguishing real and noncontributors. In conclusion, our panel holds great potential for detecting and analyzing DNA mixtures in forensic applications, with the capability to enhance routine mixture analysis.

Integrative lncRNA, circRNA, and mRNA analysis reveals expression profiles of six forensic body fluids/tissue.

RNAs have attracted much attention in forensic body fluid/tissue identification (BFID) due to their tissue-specific expression characteristics. Among RNAs, long RNAs (e.g., mRNA) have a higher probability of containing more polymorphic sites that can be used to assign the specific donor of the body fluid/tissue. However, few studies have characterized their overall profiles in forensic science. In this study, we sequenced the transcriptomes of 30 samples from venous blood, menstrual blood, semen, saliva, vaginal secretion, and skin tissue, obtaining a comprehensive picture of mRNA, lncRNA, and circRNA profiles. A total of 90,305 mRNAs, 102,906 lncRNAs (including 19,549 novel lncRNAs), and 40,204 circRNAs were detected. RNA type distribution, length distribution, and expression distribution were presented according to their annotation and expression level, and many novel body fluid/tissue-specific RNA markers were identified. Furthermore, the cognate relations among the three RNAs were analyzed according to gene annotations. Finally, SNPs and InDels from RNA transcripts were genotyped, and 21,611 multi-SNP and 4,471 multi-InDel transcriptomic microhaplotypes (tMHs) were identified. These results provide a comprehensive understanding of transcriptome profiles, which could provide new avenues for tracing the origin of the body fluid/tissue and identifying an individual.

Estimating microhaplotype allele frequencies from low-coverage or pooled sequencing data.

BACKGROUND: Microhaplotypes have the potential to be more cost-effective than SNPs for applications that require genetic panels of highly variable loci. However, development of microhaplotype panels is hindered by a lack of methods for estimating microhaplotype allele frequency from low-coverage whole genome sequencing or pooled sequencing (pool-seq) data. RESULTS: We developed new methods for estimating microhaplotype allele frequency from low-coverage whole genome sequence and pool-seq data. We validated these methods using datasets from three non-model organisms. These methods allowed estimation of allele frequency and expected heterozygosity at depths routinely achieved from pooled sequencing. CONCLUSIONS: These new methods will allow microhaplotype panels to be designed using low-coverage WGS and pool-seq data to discover and evaluate candidate loci. The python script implementing the two methods and documentation are available at https://www.github.com/delomast/mhFromLowDepSeq .

Construction and evaluation of a novel set of 90 microhaplotypes for forensic applications using NGS technology.

Microhaplotypes (MHs), small sets of linked single nucleotide polymorphisms (SNPs), are becoming a valuable tool for paternity testing, personal identification and other different forensic purposes due to their advantages of both short tandem repeats (STRs) and SNPs. However, only a small part of MHs with small segments have been developed and reported so far. And the current population genetic data of MHs are still insufficient. MHs with small segments possess unique advantages in mixture deconvolution, degradation material identification, noninvasive prenatal paternity testing and even medical tumor diagnostic applications. In the present study, a set of 90 autosomal MHs whose PCR amplicon lengths are from 90-150 bp, of which 58 MHs are less than or equal to 100 bp are selected, and assembled into an amplification multiplex system optimized for Ion S5™ System for forensic application. Genetic diversity study of 90 MHs in the populations from different intercontinental regions shows that the polymorphism information content (PIC) values of 83 MHs are greater than 0.4 in populations from East Asia (EAS), and the average PIC value of 90 MHs is greater than 0.5. A total of EAS populations shows the highest cumulative match probability (CMP) and cumulative probability of exclusion (CPE) values in five intercontinental populations. The CMP and CPE values of 90 MHs in EAS are 1.1688 × 10-54 and 0.999999999998954. The informativeness for assignment (In) values of the 90 MHs are calculated based on data from five intercontinental populations, and the In values of 20 MHs have greater than 0.1, indicating that the 20 MHs are high effectiveness in distinguishing different intercontinental populations, which can be used as candidate ancestry informative markers. Further, we have studied the polymorphisms of the 90 MHs based on 224 unrelated individuals of Henan Han population, China, and obtained the frequency data of the 90 MHs. In the Henan Han population, the effective number of alleles (Ae) of the 90 MHs ranges from 1.7649 (MH45) to 3.9792 (MH50), and the Ae values of 10 MHs reach to 3.0; the Ae values of 80 MHs are greater than 2, and the average Ae value for these MHs is 2.422. The average expected heterozygosity, observed heterozygosity, PIC, matching probability, discrimination power and probability of exclusion values of 90 MHs in the Henan Han population are 0.5788, 0.5851, 0.5039, 0.2608, 0.7392 and 0.2806, respectively. The CMP value of 90 MHs in the study population is less than 10-54, and their CPE value reaches 0.999999999999999923. Moreover, the results of the depth of coverage, allele coverage ratio and noise level indicate that the 90 MH amplification system has well sequencing performance, and the sequencing results are reliable. The results indicate the 90 MHs show higher polymorphisms in the study population. The present panel can be well used in paternity testing and individual identification in the study population and even the populations from EAS.

Forensic genetic analysis of single-nucleotide polymorphisms and microhaplotypes in Koreans through next-generation sequencing using precision ID identity panel.

BACKGROUND: Forensic DNA analysis has seen remarkable advancements with the advent of Next Generation Sequencing (NGS). In particular, NGS analysis of single nucleotide polymorphisms (SNPs) offers significant advantages in the analysis of challenging samples compared to conventional STR analysis. OBJECTIVE: This study aimed to investigate the SNPs of the Precision ID Identity Panel, a commercially available NGS panel for personal identification, by generating genetic profiles of 298 Koreans and comparing them with other global populations. METHODS: A total of 124 SNPs, including 90 autosomal and 34 Y-SNPs, were analyzed using the Precision ID Identity Panel, and forensic parameters, microhaplotypes, and population differences were investigated. RESULTS: The NGS data were successfully obtained from 298 Koreans. The analysis of forensic parameters exhibited a low combined match probability of 1.532 × 10- 34, which is comparable to that obtained from commonly used STR analysis. Additionally, the microhaplotype analysis revealed that the use of 16 microhaplotypes provided higher discriminatory power compared to single target SNPs. Furthermore, the adoption of microhaplotype data resulted in an increase of over 20% in expected heterozygosity at five loci. Inter-population analysis showed a close genetic relationship between Koreans and individuals from China and Myanmar in East and Southeast Asia, which are geographically adjacent to Korea. CONCLUSIONS: The results of this study show that the Precision ID Identity panel can be a useful alternative where traditional STR typing is not feasible. Also, the data from our study will be useful as a reference for Koreans in forensic investigations and the prosecution of criminal justice.

Identification Strategy of Biological Half Sibling Relationship.

OBJECTIVES: To compare the application value of the likelihood ratio (LR) method and identity by state (IBS) method in the identification involving half sibling relationships, and to provide a reference for the setting of relevant standards for identification of half sibling relationship. METHODS: (1) Based on the same genetic marker combinations, the reliability of computer simulation method was verified by comparing the distributions of cumulated identity by state score (CIBS) and combined full sibling index in actual cases with the distributions in simulated cases. (2) In different numbers of three genetic marker combinations, the simulation of full sibling, half sibling and unrelated individual pairs, each 1 million pairs, was obtained; the CIBS, as well as the corresponding types of cumulative LR parameters, were calculated. (3) The application value of LR method was compared with that of IBS method, by comparing the best system efficiency provided by LR method and IBS method when genetic markers in different amounts and of different types and accuracy were applied to distinguish the above three relational individual pairs. (4) According to the existing simulation data, the minimum number of genetic markers required to distinguish half siblings from the other two relationships using different types of genetic markers was estimated by curve fitting. RESULTS: (1) After the rank sum test, under the premise that the real relationship and the genetic marker combination tested were the same, there was no significant difference between the simulation method and the results obtained in the actual case. (2) In most cases, under the same conditions, the system effectiveness obtained by LR method was greater than that by IBS method. (3) According to the existing data, the number of genetic markers required for full-half siblings and half sibling identification could be obtained by curve fitting when the system effectiveness reached 0.95 or 0.99. CONCLUSIONS: When distinguishing half sibling from full sibling pairs or unrelated pairs, it is recommended to give preference to the LR method, and estimate the required number of markers according to the identification types and the population data, to ensure the identification effect.

Application of Microhaplotypes in Sibling Kinship Testing.

OBJECTIVES: To investigate the efficacy of different numbers of microhaplotype (MH) loci and the introduction of different reference samples on the identification of full sibling, half sibling and differentiation between full sibling and half sibling kinships, and to explore the effect of changing mutation rate on sibling testing. METHODS: First, a family map involving three generations was established, and four full sibling identification models, five half sibling identification models and five models distinguishing full and half siblings were constructed for different reference samples introduced. Based on the results of the previous study, two sets of nonbinary SNP-MH containing 34 and 54 loci were selected. Based on the above MH loci, 100 000 pairs of full sibling vs. unrelated individuals, 100 000 pairs of half sibling vs. unrelated individuals and 100 000 pairs of full sibling vs. half sibling were simulated based on the corresponding sibling kinship testing models, and the efficacy of each sibling kinship testing model was analyzed by the likelihood ratio algorithm under different thresholds. The mutant rate of 54 MH loci was changed to analyze the effect of mutation rate on sibling identification. RESULTS: In the same relationship testing model, the systematic efficacy of sibling testing was positively correlated with the number of MH loci detected. With the same number of MH loci, the efficacy of full sibling testing was better than that of uncle or grandfather when the reference sample introduced was a full sibling of A, but there was no significant difference in the identification efficacy of the four reference samples introduced for full sibling and half sibling differentiation testing. In addition, the mutation rate had a slight effect on the efficacy of sibling kinship testing. CONCLUSIONS: Increasing the number of MH loci and introducing reference samples of known relatives can increase the efficacy of full sibling testing, half sibling testing, and differentiation between full and half sibling kinships. The level of mutation rate in sibling testing by likelihood ratio method has a slight but insignificant effect on the efficacy.

Toward absolute abundance for conservation applications: Estimating the number of contributors via microhaplotype genotyping of mixed-DNA samples.

Molecular methods including metabarcoding and quantitative polymerase chain reaction have shown promise for estimating species abundance by quantifying the concentration of genetic material in field samples. However, the relationship between specimen abundance and detectable concentrations of genetic material is often variable in practice. DNA mixture analysis represents an alternative approach to quantify specimen abundance based on the presence of unique alleles in a sample. The DNA mixture approach provides novel opportunities to inform ecology and conservation by estimating the absolute abundance of target taxa through molecular methods; yet, the challenges associated with genotyping many highly variable markers in mixed-DNA samples have prevented its widespread use. To advance molecular approaches for abundance estimation, we explored the utility of microhaplotypes for DNA mixture analysis by applying a 125-marker panel to 1179 Chinook salmon (Oncorhynchus tshawytscha) smolts from the Sacramento-San Joaquin Delta, California, USA. We assessed the accuracy of DNA mixture analysis through a combination of mock mixtures containing DNA from up to 20 smolts and a trophic ecological application enumerating smolts in predator diets. Mock DNA mixtures of up to 10 smolts could reliably be resolved using microhaplotypes, and increasing the panel size would likely facilitate the identification of more individuals. However, while analysis of predator gastrointestinal tract contents indicated DNA mixture analysis could discern the presence of multiple prey items, poor and variable DNA quality prevented accurate genotyping and abundance estimation. Our results indicate that DNA mixture analysis can perform well with high-quality DNA, but methodological improvements in genotyping degraded DNA are necessary before this approach can be used on marginal-quality samples.

US Population Data for 94 Identity-Informative SNP Loci.

The US National Institute of Standards and Technology (NIST) analyzed a set of 1036 samples representing four major US population groups (African American, Asian American, Caucasian, and Hispanic) with 94 single nucleotide polymorphisms (SNPs) used for individual identification (iiSNPs). The compact size of iiSNP amplicons compared to short tandem repeat (STR) markers increases the likelihood of successful amplification with degraded DNA samples. Allele frequencies and relevant forensic statistics were calculated for each population group as well as the aggregate population sample. Examination of sequence data in the regions flanking the targeted SNPs identified additional variants, which can be combined with the target SNPs to form microhaplotypes (multiple phased SNPs within a short-read sequence). Comparison of iiSNP performance with and without flanking SNP variation identified four amplicons containing microhaplotypes with observed heterozygosity increases of greater than 15% over the targeted SNP alone. For this set of 1036 samples, comparison of average match probabilities from iiSNPs with the 20 CODIS core STR markers yielded an estimate of 1.7 × 10-38 for iiSNPs (assuming independence between all 94 SNPs), which was four orders of magnitude lower (more discriminating) than STRs where internal sequence variation was considered, and 10 orders of magnitude lower than STRs using established capillary electrophoresis length-based genotypes.

Evaluation of microhaplotype panels for complex kinship analysis using massively parallel sequencing.

In recent years, microhaplotypes (MHs) have become a research hotspot within the field of forensic genetics. Traditional MHs contain only SNPs that are closely linked within short fragments. Herein, we broaden the concept of general MHs to include short InDels. Complex kinship identification plays an important role in disaster victim identification and criminal investigations. For distant relatives (e.g., 3rd-degree), many genetic markers are required to enhance power of kinship testing. We performed genome-wide screening for new MH markers composed of two or more variants (InDel or SNP) within 220 bp based on the Chinese Southern Han from the 1000 Genomes Project. An NGS-based 67plex MH panel (Panel B) was successfully developed, and 124 unrelated individual samples were sequenced to obtain population genetic data, including alleles and allele frequencies. Of the 67 genetic markers, 65 MHs were, as far as we know, newly discovered, and 32 MHs had effective number of allele (Ae) values greater than 5.0. The average Ae and heterozygosity of the panel were 5.34 and 0.7352, respectively. Next, 53 MHs from a previous study were collected as Panel A (average Ae of 7.43), and Panel C with 87 MHs (average Ae of 7.02) was formed by combining Panels A and B. We investigated the utility of these three panels in kinship analysis (parent-child, full siblings, 2nd-degree, 3rd-degree, 4th-degree, and 5th-degree relatives), with Panel C exhibiting better performance than the two other panels. Panel C was able to separate parent-child, full-sibling, and 2nd-degree relative duos from unrelated controls in real pedigree data, with a small false testing level (FTL) of 0.11% in simulated 2nd-degree duos. For more distant relationships, the FTL was much higher: 8.99% for 3rd-degree, 35.46% for 4th-degree, and 61.55% for 5th-degree. When a carefully chosen extra relative was known, this may enhance the testing power for distant kinship analysis. Two twins from the Q family (2-5 and 2-7) and W family (3-18 and 3-19) shared the same genotypes in all tested MHs, which led to the incorrect conclusion that an uncle-nephew duo was classified as a parent-child duo. In addition, Panel C showed great capacity for excluding close relatives (2nd-degree and 3rd-degree relatives) during paternity tests. Among 18,246 real and 10,000 simulated unrelated pairs, none were misinterpreted as a relative within 2nd-degree at a log10(LR) cutoff of 4. The panels presented herein could provide supplementary power for the analysis of complex kinship.

Evaluation of large-scale highly polymorphic microhaplotypes in complex DNA mixtures analysis using RMNE method.

DNA mixture interpretation is one of the most challenging problems in forensics. Complex DNA mixtures are more difficult to analyze when there are more than two contributors or related contributors. Microhaplotypes (MHs) are polymorphic genetic markers recently discovered and employed in DNA mixture analysis. However, the evidentiary interpretation of the MH genotyping data needs more debate. The Random Man Not Excluded (RMNE) method analyzes DNA mixtures without using allelic peak height data or the number of contributors (NoC) assumptions. This study aimed to assess how well RMNE interpreted mixed MH genotyping data. We classified the MH loci from the 1000 Genomes Project database into groups based on their Ae values. Then we performed simulations of DNA mixtures with 2-10 unrelated contributors and DNA mixtures with a pair of sibling contributors. For each simulated DNA mixture, incorrectly included ratios were estimated for three types of non-contributors: random men, parents of contributors, and siblings of contributors. Meanwhile, RMNE probability was calculated for contributors and three types of non-contributors, allowing loci mismatch. The results showed that the MH number, the MH Ae values, and the NoC affected the RMNE probability of the mixture and the incorrectly included ratio of non-contributors. When there were more MHs, MHs with higher Ae values, and a mixture with less NoC, the RMNE probability, and the incorrectly included ratio decreased. The existence of kinship in mixtures complicated the mixture interpretation. Contributors' relatives as non-contributors and related contributors in the mixture increased the demands on the genetic markers to identify the contributors correctly. When 500 highly polymorphic MHs with Ae values higher than 5 were used, the four individual types could be distinguished according to the RMNE probabilities. This study reveals the promising potential of MH as a genetic marker for mixed DNA interpretation and the broadening of RMNE as a parameter indicating the relationship of a specific individual with a DNA mixture in the DNA database search.

An MPS-Based 50plex Microhaplotype Assay for Forensic DNA Analysis.

Microhaplotypes (MHs) are widely accepted as powerful markers in forensic studies. They have the advantage of both short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs), with no stutter and amplification bias, short fragments and amplicons, low mutation and recombination rates, and high polymorphisms. In this study, we constructed a panel of 50 MHs that are distributed on 21 chromosomes and analyzed them using the Multiseq multiple polymerase chain reaction (multi-PCR) targeted capture sequencing protocol based on the massively parallel sequencing (MPS) platform. The sizes of markers and amplicons ranged between 11-81 bp and 123-198 bp, respectively. The sensitivity was 0.25 ng, and the calling results were consistent with Sanger sequencing and the Integrative Genomics Viewer (IGV). It showed measurable polymorphism among sequenced 137 Southwest Chinese Han individuals. No significant deviations in the Hardy-Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were found at all MHs after Bonferroni correction. Furthermore, the specificity was 1:40 for simulated two-person mixtures, and the detection rates of highly degraded single samples and mixtures were 100% and 93-100%, respectively. Moreover, animal DNA testing was incomplete and low depth. Overall, our MPS-based 50-plex MH panel is a powerful forensic tool that provides a strong supplement and enhancement for some existing panels.

Forensic identity SNPs: Characterisation of flanking region variation using massively parallel sequencing.

Single nucleotide polymorphisms (SNPs) can be analysed for identity or kinship applications in forensic genetics to either provide an adjunct to traditional STR typing or as a stand-alone approach. The advent of massively parallel sequencing technology (MPS) has provided a useful opportunity to more easily deploy SNP typing in a forensic context, given the ability to simultaneously amplify a large number of markers. Furthermore, MPS also provides valuable sequence data for the targeted regions, which enables the detection of any additional variation seen in the flanking regions of amplicons. In this study we genotyped 977 samples across five UK-relevant population groups (White British, East Asian, South Asian, North-East African and West African) for 94 identity-informative SNP markers using the ForenSeq DNA Signature Prep Kit. Examination of flanking region variation allowed for the identification of 158 additional alleles across all populations studied. Here we present allele frequencies for all 94 identity-informative SNPs, both including and excluding the flanking region sequence of these markers. We also present information on the configuration of these SNPs in the ForenSeq DNA Signature Prep Kit, including performance metrics for the markers and investigation of bioinformatic and chemistry-based discordances. Overall, the inclusion of flanking region variation in the analysing workflow for these markers reduced the average combined match probability 2175 times across all populations, with a maximum reduction of 675,000-fold in the West African population. The gain due to flanking region-based discrimination increased the heterozygosity of some loci above that of some of the least useful forensic STR loci; thus demonstrating the benefit of enhanced analysis of currently targeted SNP markers for forensic applications.