Evaluation of a microhaplotypes panel for forensic genetics using massive parallel sequencing technology

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.

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.

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 biogeographical ancestry inference: recent insights and current trends

BACKGROUND

As a powerful complement to the paradigmatic DNA profiling strategy, biogeographical ancestry inference (BGAI) plays a significant part in human forensic investigation especially when a database hit or eyewitness testimony are not available. It indicates one's biogeographical profile based on known population-specific genetic variations, and thus is crucial for guiding authority investigations to find unknown individuals. Forensic biogeographical ancestry testing exploits much of the recent advances in the understanding of human genomic variation and improving of molecular biology.

OBJECTIVE

In this review, recent development of prospective ancestry informative markers (AIMs) and the statistical approaches of inferring biogeographic ancestry from AIMs are elucidated and discussed.

METHODS

We highlight the research progress of three potential AIMs (i.e., single nucleotide polymorphisms, microhaplotypes, and Y or mtDNA uniparental markers) and discuss the prospects and challenges of two methods that are commonly used in BGAI.

CONCLUSION

While BGAI for forensic purposes has been thriving in recent years, important challenges, such as ethics and responsibilities, data completeness, and ununified standards for evaluation, remain for the use of biogeographical ancestry information in human forensic investigations. To address these issues and fully realize the value of BGAI in forensic investigation, efforts should be made not only by labs/institutions around the world independently, but also by inter-lab/institution collaborations.

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 (A_e) values greater than 5.0. The average A_e 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 A_e of 7.43), and Panel C with 87 MHs (average A_e 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 log₁₀(LR) cutoff of 4. The panels presented herein could provide supplementary power for the analysis of complex kinship.

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.

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.

Application of a custom haplotype caller to analyze sequence-based data of 56 microhaplotypes

Microhaplotypes (microhaps) are recently introduced markers that aim to complement the limitations of conventional forensic markers such as short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs). With the potential of microhaps in forensics becoming clearer through massively parallel sequencing (MPS), MPS-based studies on microhaps are being actively reported. However, simpler workflow schemes for the generation and analysis of MPS data are still required to facilitate the practical application of MPS in forensics. In this study, we developed an in-house MPS panel that simultaneously amplifies 56 microhaps and a custom haplotype caller, Visual Microhap. The developed tool works on a web browser and provides four analysis options to extract SNP-based haplotypes from sequence-based data obtained by STRait Razor 3.0. To demonstrate the utility of the MPS panel and data analysis workflow scheme, we also analyzed 56 microhaps of 286 samples from four populations (African-American, Caucasian, Hispanic, and Korean). The average effective number of alleles (A_e) for the four groups was 3.45, ranging from 1.74 to 6.98. Forensic statistical parameters showed that this microhap panel is more powerful than conventional autosomal STRs for human identification. Meanwhile, the 56-plex panel mostly comprised microhaps with high Ae; however, the four populations were grossly distinguishable from each other by cluster analysis. Consequently, the developed in-house MPS panel for 56 microhaps and the adopted workflow using open-source tools can increase the utility of microhap MPS in forensic research and practice.

A multipurpose panel of microhaplotypes for use with STR markers in casework

A small panel of highly informative loci that can be genotyped on the same equipment as the standard CODIS short tandem repeat (STR) markers has strong potential for application in forensic casework. Single nucleotide polymorphisms (SNPs) can be typed by a couple of methods on capillary electrophoresis (CE) machines and on sequencers, but the amount of information relative to the laboratory effort has hindered use of SNPs in actual casework. Insertion-deletion markers (InDels) suffer from similar problems. Microhaplotypes (MHs) are much more informative per locus but have similar technical difficulties unless they are typed by massively parallel sequencing (MPS). As forensic labs are acquiring sequencing machines, MHs become more likely to be used in casework, especially if multiplexed with STRs. Here we present the details of a multipurpose panel of 24 MHs with the highest effective number of alleles (Ae) from previous work. An augmented STR panel of 24 loci (20 CODIS markers plus four commonly typed STRs) is also considered. The Ae and ancestry informativeness (In) distributions of these two datasets are compared. The MH panel is shown to have better individualization and population distinction than the augmented CODIS STRs. We note that the 24 MHs should be better for mixture analyses than the STRs. Finally, we suggest that a commercial kit including both the standard CODIS markers and this set of 24 MH would greatly improve the discrimination power over that of current commercial assays.

Evaluation of the MHSeqTyper47 kit for forensically challenging DNA samples

Microhaplotypes have been highly regarded for forensic mixture DNA deconvolution because they do not experience interference from stutters in the same way as short tandem repeat markers, and they tend to be more polymorphic than single nucleotide polymorphism markers. However, forensic microhaplotype kits have not been reported. The MHSeqTyper47 kit genotypes 47 microhaplotype loci. In this study, MiSeq FGx sequencing metrics for MHSeqTyper47 were presented, and the genotyping accuracy of this kit was examined. The sensitivity of MHSeqTyper47 reached 62.5 pg, and full genotyping results were obtained from degraded DNA samples with degradation indexes ≤ 3.00. Full genotypes were obtained in the presence of 100 ng/μL tannin, 50 μM heme, 25 ng/μL humic acid, and 1.25 μg/μL indigo dye. In DNA mixture studies, a minimum of 31 loci of the minor contributor were correctly genotyped at 1:99 or 99:1 mixing ratios, with the cumulative random matching probability of these loci reaching 4.54 × 10^-25. Mixing ratios could be reliably predicted from two-donor DNA mixtures based on the loci with four called alleles. Taken together, these data showed that the MHSeqTyper47 kit was effective for forensically challenging DNA analysis.

Forensic nanopore sequencing of microhaplotype markers using QitanTech's QNome

In recent years, extraordinary progress has been made in genome sequencing technologies, which has led to a decrease in cost and an increase in the diversity of sequenced genomes. Nanopore sequencing is one of the latest genome sequencing technologies. It aims to sequence longer contiguous pieces of DNA, which are essential for resolving structurally complex regions, and provides a new approach for forensic genetics to detect longer markers in real time. To date, multiple studies have been conducted to sequence forensic markers using MinION from Oxford Nanopore Technologies (ONT), and the results indicate that nanopore sequencing holds promise for forensic applications. Qitan Technology (QitanTech) recently launched its first commercial nanopore genome sequencer, QNome. It could achieve a read length of more than 150 kbp, and could generate approximately 500 Mb of data in 8 h. In this pilot study, we explored and validated this alternative nanopore sequencing device for microhaplotype (MH) profiling using a custom set of 15 MH loci. Seventy single-contributor samples were divided into 7 batches, each of which included 10 samples and control DNA 9947A and was sequenced by QNome. MH genotypes generated from QNome were compared to those from Ion Torrent sequencing (Ion S5XL system) to evaluate the accuracy and stability. Twelve samples randomly selected from the last three batches and Control DNA 9947A were also subjected to ONT MinION sequencing (with R9.4 flow cell) for parallel comparison. Based on MHtyper, a bioinformatics workflow developed for automated MH designation, all MH loci can be genotyped and reliably phased using the QNome data, with an overall accuracy of 99.83% (4 errors among 2310 genotypes). Three occurred near or in the region of homopolymer sequences, and one existed within 50 bp of the start of the sequencing reaction. In the last 15 samples (12 individual samples and 3 replicates of control DNA 9947A), two SNPs located at 4-mer homopolymers failed to obtain reliable genotypes on the MinION data. This study shows the potential of state-of-the-art nanopore sequencing methods to analyze forensic MH markers. Given the rapid pace of change, sporadic and nonrepetitive errors presented in this study are expected to be resolved by further developments of nanopore technologies and analysis tools.

Research progress on application of microhaplotype in forensic genetics

As a novel genetic marker, microhaplotype can be applied in the field of forensic genetics. Microhaplotype has the advantages of high polymorphism, low mutation rate, no stutter products and short amplification fragments. Microhaplotype can effectively detect mixture, and quantitatively analyze the contributors of mixture. DNA with severe fragmentation can be successfully genotyped by microhaplotype. It can be used as ancestry informative marker to effectively divide the global continental population according to genetic structure. Microhaplotype system can provide more information than traditional short tandem repeat and help to identify complex relationships. It can provide new ideas for tumor source identification, cell line identification and prenatal paternity testing. Here we review the applications of microhaplotype, intending to provide references for forensic practice.

The population genetics characteristics of a 90 locus panel of microhaplotypes

Single-nucleotide polymorphisms (SNPs) and small genomic regions with multiple SNPs (microhaplotypes, MHs) are rapidly emerging as novel forensic investigative tools to assist in individual identification, kinship analyses, ancestry inference, and deconvolution of DNA mixtures. Here, we analyzed information for 90 microhaplotype loci in 4009 individuals from 79 world populations in 6 major biogeographic regions. The study included multiplex microhaplotype sequencing (mMHseq) data analyzed for 524 individuals from 16 populations and genotype data for 3485 individuals from 63 populations curated from public repositories. Analyses of the 79 populations revealed excellent characteristics for this 90-plex MH panel for various forensic applications achieving an overall average effective number of allele values (A_e) of 4.55 (range 1.04–19.27) for individualization and mixture deconvolution. Population-specific random match probabilities ranged from a low of 10^–115 to a maximum of 10^–66. Mean informativeness (I_n) for ancestry inference was 0.355 (range 0.117–0.883). 65 novel SNPs were detected in 39 of the MHs using mMHseq. Of the 3018 different microhaplotype alleles identified, 1337 occurred at frequencies > 5% in at least one of the populations studied. The 90-plex MH panel enables effective differentiation of population groupings for major biogeographic regions as well as delineation of distinct subgroupings within regions. Open-source, web-based software is available to support validation of this technology for forensic case work analysis and to tailor MH analysis for specific geographical regions.

A highly polymorphic panel of 40-plex microhaplotypes for the Chinese Han population and its application in estimating the number of contributors in DNA mixtures

Microhaplotypes (MHs) have great potential in multiple forensic applications and have proven to be promising markers in complex DNA mixture analysis. In this study, we developed a multiplex panel of 40 highly polymorphic MHs for the Chinese Han population, evaluated its forensic values, and explored its application in predicting the number of contributors (NOCs) in DNA mixtures. The panel consisted of 20 newly proposed loci and 20 previously reported loci with lengths spanning less than 120 bp. The average effective number of alleles (A_e) was 3.77, and the cumulative matching probability (CMP) and the cumulative power of exclusion (CPE) reached 1.2E-37 and 1-2.1E-12, respectively, in the Chinese Han population from the 1000 Genomes Project. Further validation on 150 Chinese Han individuals showed that A_e ranged from 2.62 to 4.41 with a mean value of 3.61, and CMP and CPE were 3.61E-36 and 1-1.84E-12, respectively, indicating that this panel was informative for personal identification and paternity testing in the studied population. To estimate NOC in DNA mixtures, we developed a machine learning model based on this panel. As a result, the accuracies in artificial DNA mixtures reached 95.24% for 2- to 4-person mixtures and 83.33% for 2- to 6-person mixtures. Furthermore, the NOC estimation on simulated profiles with allele dropout showed that this panel was still robust under slight dropout. In conclusion, this panel has value for forensic identification and NOC estimation of DNA mixtures.

Pairwise kinship testing with microhaplotypes: Can advancements be made in kinship inference with these markers?

Kinship testing based on genetic relatedness is one of the major tasks in forensic genetics. Although short tandem repeats (STRs) are the "gold standard" biomarkers for relationship testing, microhaplotypes (MHs) have also emerged as viable options for kinship elucidation. In this work, the kinship testing efficiency of 54 highly polymorphic MHs was studied in two extended families consisting of parent-offspring, full siblings, grandparent-grandchildren, uncle/aunt-nephew/nieces, and first cousins. In addition, ten-thousand pairs of different degrees of relationships were simulated using various datasets including 54 MHs, 27 STRs plus 94 single nucleotide polymorphisms (SNPs) that were included in the ForenSeq DNA Signature Prep Kit (ForenSeq), 54 MHs plus loci in ForenSeq, and different subsets of 417-published MHs. The panels' system effectiveness in the kinship analysis were accessed by likelihood ratio distributions. The results showed that 54 MHs could be used in first-degree relationship testing with high reliability. The effectiveness of 54 MHs was slightly lower than ForenSeq but only by a narrow margin. Both 54 MHs and ForenSeq were not sufficient for distant relationship testing, and approximately 200 microhaplotypes with an average expected heterozygosity (He) = 0.79 were enough to determine second-degree relationships, but a panel of 417 MHs with an average He = 0.72 was not sufficient to first cousins testing according to the simulation analysis. In conclusion, 54 MHs could be used to serve as supplement markers for kinship testing; and well-established STR markers plus well-performing microhaplotype markers may become collective tools in forensic applications, though an enlarged pool of forensic markers is needed for distant relationship testing.

Genetic relationships of Southwest Asian and Mediterranean populations

The Southwest Asian, circum-Mediterranean, and Southern European populations (collectively, SWAMSE) together with Northern European populations form one of five "continental" groups of global populations in many analyses of population relationships. This region is of great anthropologic and forensic interest but relationships of large numbers of populations within the region have not been able to be cleanly resolved with autosomal genetic markers. To examine the genetic boundaries to the SWAMSE region and whether internal structure can be detected we have assembled data for a total of 151 separate autosomal genetic markers on populations in this region and other parts of the world for a global set of 95 populations. The markers include 83 ancestry informative SNPs as singletons and 68 microhaplotype loci defined by 204 SNPs. The 151 loci are ancestry informative on a global scale, identifying at least five biogeographic clusters. One of those clusters is a clear grouping of 37 populations containing the SWAMSE plus northern European populations to the exclusion of populations in South Central Asia and populations from farther East. A refined analysis of the 37 populations shows the northern European populations clustering separately from the SWAMSE populations. Within Southwest Asia the Samaritans and Shabaks are distinct outliers. The Yemenite Jews, Saudi, Kuwaiti, Palestinian Arabs, and Southern Tunisians cluster together loosely while the remaining populations from Northern Iraq, Mediterranean Europe, the Caucasus region, and Iran cluster in a more complex graded fashion. The majority of the SWAMSE populations from the mainland of Southwest Asia form a cluster with little internal structure reflecting a very complex history of endogamy and migrations. The set of 151 DNA polymorphisms not only distinguishes major geographical regions globally but can distinguish ancestry to a small degree within geographical regions such as SWAMSE. We discuss forensic characteristics of the polymorphisms and also identify those that rank highest by Rosenberg's I_n measure for the SWAMSE region populations and for the global set of populations analyzed. DATA AVAILABILITY: Genotypes on all 151 markers for all 3790 individuals typed in the Kidd Lab on the 72 Kidd lab populations have been deposited in the Zenodo archive and can be freely accessed at https://doi.org/10.5281/zenodo.4658892. Some of the data has been made public previously as supplemental files appended to publications. Data for the additional individuals included in the analyses was taken from already public datasets as indicated in the text.

Identification and sequencing of 59 highly polymorphic microhaplotypes for analysis of DNA mixtures

Mixture detection remains one of the major challenges within a forensic science context. In recent years, microhaplotypes were proposed to have great potential in mixture detection, although many of them are not as polymorphic as widely used short tandem repeat (STR) markers. In this study, 59 new highly polymorphic microhaplotypes were identified and sequenced with the NextSeq 500 Sequencer. Based on the whole 1000 Genomes Project dataset, the average effective number of alleles (Ae) of the 59 microhaplotypes was 5.44, and the Ae values of 36 of these microhaplotypes were > 5.00. Their genetic variations in 187 Han Chinese individuals were evaluated. The average allele coverage ratio (ACR) of heterozygotes across all loci was 0.96 ± 0.05. The number of observed alleles varied from 4 to 23, with an average of 8.8 alleles per microhaplotype locus. The average observed heterozygosity (Ho) of 59 loci was 0.77 ± 0.05, and the Ho values of 15 of these loci were > 0.80. All loci showed high polymorphisms with a discrimination power (DP) ranging from 0.80 to 0.97, and the average DP was 0.92 ± 0.03. The analysis of simulated mixtures demonstrated that the microhaplotypes reported here were highly polymorphic and performed well in forensic DNA mixture analysis. This study not only demonstrated the applicability of microhaplotypes in mixture analysis but also provided new choices for highly polymorphic microhaplotypes because after adding the markers identified here, the number of microhaplotypes with Ae values of > 4.00 will increase from ~ 50 to ~ 110.

Evaluation of three microhaplotypes in individual identification and ancestry inference

Microhaplotype as an emerging genetic marker has attracted more attention in forensic field. The purpose of this study was to evaluate the potential of microhaplotypes in individual identification and ancestry inference in Chinese Hainan Li and 26 1000 G populations. Three microhaplotypes were genotyped from 100 Li individuals using Agena MassARRAY. Moreover, 2504 individuals from 26 populations (1000 Genomes Project database) were enrolled. The genotypes frequencies of microhaplotypes in each population were calculated by the Plink software. We used Structure, Arlequin, and MEGA6 software to analyze the genetic structure, differentiation and genetic background difference, respectively. The forensic parameters of these microhaplotypes were calculated using Modified Powerstats software. The distribution of genotypes frequencies of three microhaplotypes elaborated the high diversities among the Li and 26 1000 G populations. Li population had a close genetic relationship with EAS populations using structure analysis. No differentiation was observed between Li and CHS population by Fst analysis. The NJ tree showed that the genetic background of Li and CHS is most similar. The average heterozygosity (HE), probability of match (PM), power of discrimination (PD), probability of exclusion (PE) and polymorphism information content (PIC) values for the three microhaplotypes in 27 populations were 0.535, 0.497, 0.465, 0.325, and 0.481, respectively. In conclusion, our results revealed three microhaplotypes as individual identification and ancestry inference genetic markers among Li population and 26 1000 G populations. Future studies are needed to confirm our results with larger samples and select much higher forensic efficacy microhaplotypes.

Multi-Indel: A Microhaplotype Marker Can Be Typed Using Capillary Electrophoresis Platforms

Since the concept of microhaplotypes was proposed by Kidd in 2013, various microhaplotype markers have been investigated for various forensic purposes, such as individual identification, deconvolution of DNA mixtures, or forensic ancestry inference. In our opinion, various compound markers are also regarded as generalized microhaplotypes, encompassing two or more variants in a short segment of DNA (e.g., 200 bp). That is, a set of variants (referred to herein as multi-variants) within a certain length includes single nucleotide polymorphisms (SNP), insertion/deletion polymorphisms (Indels), or short tandem repeat polymorphisms (STRs). At present, multi-variant is mainly aimed at multi-SNPs. However, the haplotype genotyping of multi-variants relies on single-strand analysis, mainly using massively parallel sequencing (MPS). Here, we describe a method based on a capillary electrophoresis (CE) platform that can directly obtain haplotypes of individuals. Several microhaplotypes consisting of three or more Indels with different insertion or deletion lengths in the range of less than 200 bp were screened out, each of which had at least three haplotypes. As a result, the haplotype of an individual was reflected by the length of its polymorphism. Finally, we established a multiplex amplification system containing 18 multi-Indel markers that could identify haplotypes on each chromosome of an individual. The combined power of discrimination (CPD) and the cumulative probability of exclusion (CPE) were 0.999999999997234 and 0.9984, respectively.

Broadening the Applicability of a Custom Multi-Platform Panel of Microhaplotypes: Bio-Geographical Ancestry Inference and Expanded Reference Data

The development of microhaplotype (MH) panels for massively parallel sequencing (MPS) platforms is gaining increasing relevance for forensic analysis. Here, we expand the applicability of a 102 autosomal and 11 X-chromosome panel of MHs, previously validated with both MiSeq and Ion S5 MPS platforms and designed for identification purposes. We have broadened reference population data for identification purposes, including data from 240 HGDP-CEPH individuals of native populations from North Africa, the Middle East, Oceania and America. Using the enhanced population data, the panel was evaluated as a marker set for bio-geographical ancestry (BGA) inference, providing a clear differentiation of the five main continental groups of Africa, Europe, East Asia, Native America, and Oceania. An informative degree of differentiation was also achieved for the population variation encompassing North Africa, Middle East, Europe, South Asia, and East Asia. In addition, we explored the potential for individual BGA inference from simple mixed DNA, by simulation of mixed profiles followed by deconvolution of mixture components.

Population genetic data of 74 microhaplotypes in four major U.S. population groups

Microhaplotypes (microhaps or MHs) are novel forensically relevant genetic markers that demand large and appropriate allele frequency datasets for their implementation in casework. In this study we report on the allele frequency data of 74 microhap loci (230 SNPs) included in a newly developed 74-plex assay. The panel was tested on the Ion S5 system on a total of 347 samples from four main U.S. population groups of African, European, East Asian and Southwest Hispanic descent. Overall, frequencies of individual alleles at each locus varied considerably among the different population groups. An increase in the average value of gene diversity was also observed as the number of SNPs per locus increased. Most microhap markers showed no significant deviation from Hardy-Weinberg ratios within any of the individual population samples displaying an average power of discrimination between 0.74 and 0.81 and an average probability of exclusion between 0.32 and 0.39. Moreover, the four population groups had no clear genetic affinities with the exception of U.S. European and U.S. Southwest Hispanic populations, which showed the lowest F_ST value. STRUCTURE and principal component analyses (PCA) analysis resulted in effective clustering of the four populations with the U.S. European and Southwest Hispanic showing some overlap. These results support the potential use of this sequence-based 74plex-microhaplotype assay for ancestry inference in addition to previously reported human identification and mixture deconvolution capabilities.

A Highly Polymorphic Panel Consisting of Microhaplotypes and Compound Markers with the NGS and Its Forensic Efficiency Evaluations in Chinese Two Groups

Novel genetic markers like microhaplotypes and compound markers show promising potential in forensic research. Based on previously reported single nucleotide polymorphism (SNP) and insertion/deletion (InDel) polymorphism loci, 29 genetic markers including 22 microhaplotypes and seven compound markers were identified. Genetic distributions of the 29 loci in five continental populations, Kazak and Mongolian groups in China were investigated. We found that the expected heterozygosity values of these 29 loci were >0.4 in these populations, indicating these loci were relatively high polymorphisms. Population genetic analyses of five continental populations showed that five loci displayed relatively high genetic variations among these continental populations and could be useful markers for ancestry analysis. In summary, the 29 loci displayed relatively high genetic diversities in continental populations and Chinese two groups and could be informative loci for forensic research.

Validation of novel forensic DNA markers using multiplex microhaplotype sequencing

Microhaplotypes (MH) are comprised of multiple single nucleotide polymorphisms (SNPs) that are located within 300 bases of genomic sequence. Improved tools are needed to facilitate broader application of microhaplotypes in a diverse range of populations and forensic settings. We designed an assay for multiplex sequencing of 90 microhaplotypes (mMHseq) that include 46 MH loci with high Effective Number of Alleles (A_e) from previous studies [1], and 44 high A_e MH loci containing between four to fourteen SNPs that were identified from the 1000 Genomes (1KG) Project. The unique design of mMHseq integrates a novel method for multiplex amplification from small DNA amounts, and multiplex sequencing of 48 samples in a single MiSeq run to detect all relevant MH variation. Assay performance was evaluated in a cohort of 156 individuals from seven different world populations from Africa, Asia, and Europe. Three of those populations from East Africa (Chagga, Sandawe, and Zaramo) and one from Eastern Europe (Adygei) had sufficient individuals sequenced by the assay to be included in statistical analyses with the 26 1KG populations. For those 30 populations the mean global average A_e was 5.08 (range: 2.7-11.54) and mean informativeness for biogeographic variation (I_n) was 0.30 (range: 0.08-0.70). Eighty-five novel SNPs were detected in 58 of the 90 microhaplotypes. Open-source, web-based software was developed to visualize haplotype phase data for each microhaplotype and individual. Our approach for multiplex microhaplotype sequencing can be customized and expanded as novel loci are being discovered.

DIP-microhaplotypes: new markers for detection of unbalanced DNA mixtures

The identification of a suspect in a degraded and unbalanced DNA mixture has been a challenge for the standard short tandem repeat polymorphisms (STR) typing. Several methods have been introduced to solve this problem, such as DIP-STR, DIP-SNP, and SNP-STR markers. In this study, we proposed DIP-microhaplotype (deletion/insertion linked a chain of SNPs) as a kind of new genetic marker to type the unbalanced and degraded DNA mixture. We established the detection method with ten DIP-microhaplotype markers including 26 SNPs using allele-specific multiplex PCR followed by SNaPshot assay. This novel compound marker allows us to detect the minor DNA with a sensitivity of 1:100 to 1:1000 in a DNA mixture of any gender. Most of the DIP-microhaplotype markers had a relatively high probability of informative alleles with an average informative value (I value) of 0.308. In all, we proposed DIP-microhaplotype as a novel type of DNA marker for the detection of minor contributor from unbalanced DNA mixtures. Due to their inherent shorter length, higher polymorphism, and sensitivity, DIP-microhaplotypes are promising markers for the examination of the degraded and unbalanced mixtures in forensic stains or clinical chimeras.

Construction and forensic application of 20 highly polymorphic microhaplotypes

Microhaplotype markers have become an important research focus in forensic genetics. However, many reported microhaplotype markers have limited polymorphisms. In this study, we developed a set of highly polymorphic microhaplotype markers based on tri-allelic single-nucleotide polymorphisms. Eleven newly discovered microhaplotypes along with nine previously identified in our laboratory were studied. The microhaplotype genotypes of unrelated individuals and familial samples were generated on the MiSeq PE300 platform. These 20 loci have an average greater than 3.5 effective number of alleles. Over the whole set, the cumulative power of discrimination was 1-3.3 × 10^-18, the cumulative power of exclusion was 1-1.928 × 10^-7 and the theoretical probability of detecting a mixture was 1-1.427 × 10^-6. Differentiation comparisons of 26 populations from the 1000 Genomes Project distinguished among East Asian, South Asian, African and European populations. Overall, these markers enrich the current microhaplotype marker databases and can be applied for individual identification, paternity testing and biogeographic ancestry distinction.

A 124-plex Microhaplotype Panel Based on Next-generation Sequencing Developed for Forensic Applications

Microhaplotypes are an emerging type of forensic genetic marker that are expected to support multiple forensic applications. Here, we developed a 124-plex panel for microhaplotype genotyping based on next-generation sequencing (NGS). The panel yielded intralocus and interlocus balanced sequencing data with a high percentage of effective reads. A full genotype was determined with as little as 0.1 ng of input DNA. Parallel mixture experiments and in-depth comparative analyses were performed with capillary-electrophoresis-based short tandem repeat (STR) and NGS-based microhaplotype genotyping, and demonstrated that microhaplotypes are far superior to STRs for mixture deconvolution. DNA from Han Chinese individuals (n = 256) was sequenced with the 124-plex panel. In total, 514 alleles were observed, and the forensic genetic parameters were calculated. A comparison of the forensic parameters for the 20 microhaplotypes with the top Ae values in the 124-plex panel and 20 commonly used forensic STRs showed that these microhaplotypes were as effective as STRs in identifying individuals. A linkage disequilibrium analysis showed that 106 of the 124 microhaplotypes were independently hereditary, and the combined match probability for these 106 microhaplotypes was 5.23 × 10-66. We conclude that this 124-plex microhaplotype panel is a powerful tool for forensic applications.

A microhap panel for kinship analysis through massively parallel sequencing technology

It is widely recognized that microhaps are powerful markers for different forensic purposes, mainly due to their advantages of both short tandem repeats and single nucleotide polymorphisms, including multiple alleles, low mutation rate, and absence of stutter peaks. In the present study, a panel of 60 microhap loci was developed and utilized in forensic kinship analysis as a preliminary study. Genotyping of microhap was performed by massively parallel sequencing and haplotypes were directly achieved from sequence reads of 73 samples from Chinese Han population. We observed that 49 out of 60 loci have effective number of alleles greater than 3.0 and 10 out of 60 have values above 4.0, with an average value of 3.5598. The heterozygosity values were in a range from 0.5840 to 0.8546 with an average of 0.7268 and the cumulative power of exclusion value of the 60 loci is equal to 1-4.78 × 10^-18 . Moreover, we demonstrated the applicability of this method by different relationship inference problems, including identification of single parent-offspring, full-sibling, and second-degree relative. The results indicated that the assembled microhap panel provided more power for relationship inference, than commonly used short tandem repeats or single nucleotide polymorphism system.

Building a custom large-scale panel of novel microhaplotypes for forensic identification using MiSeq and Ion S5 massively parallel sequencing systems

A large number of new microhaplotype loci were identified in the human genome by applying a directed search with selection criteria emphasizing short haplotype length (<120 nucleotides) and maximum levels of polymorphism in the composite SNPs. From these searches, 107 autosomal microhaplotypes and 11 X chromosome microhaplotypes were selected, with well-spaced autosomal positions to ensure their independence in relationship tests. The 118 microhaplotypes were assembled into a single multiplex assay for the analysis of forensic DNA with massively parallel sequencing (MPS). A single AmpliSeq-adapted primer set was made for Illumina MiSeq and Thermo Fisher Ion S5 MPS platforms and the performance of the assay was comprehensively evaluated in both systems. Five microhaplotypes showed critical sequencing failures in both MPS platforms and were removed, while a further 13 required manual checks and the application of sequence quality thresholds in one or both systems to ensure the successful analysis of low-level DNA in these loci. The targeting of short microhaplotype spans during marker selection, with an average length of 51 nucleotides in the 118 loci, led to a high level of sensitivity for the panel when sequencing the very degraded DNA typically encountered in forensic casework and the identification of missing persons.

Human Leukocyte Antigen alleles as an aid to STR in complex forensic DNA samples

Human biological samples with multiple contributors remain one of the most challenging aspects of DNA typing within a forensic science context. With the increasing sensitivity of commercially available kits allowing detection of low template DNA, complex mixtures are now a standard component of forensic DNA evidence. Over the years, various methods and techniques have been developed to try to resolve the issue of mixed profiles. However, forensic DNA analysis has relied on the same markers to generate DNA profiles for the past 30 years causing considerable challenges in the deconvolution of complex mixed samples. The future of resolving complicated DNA mixtures may rely on utilising markers that have been previously applied to gene typing of non-forensic relevance. With Massively Parallel Sequencing (MPS), techniques becoming more popular and accessible even epigenetic markers have become a source of interest for forensic scientists. The aim of this review is to consider the potential of alleles from the Human Leukocyte Antigen (HLA) complex as effective forensic markers. While Massively Parallel Sequencing of HLA is routinely used in clinical laboratories in fields such as transplantation, pharmacology or population studies, there have not been any studies testing its suitability for forensic casework samples.