State of the Art for Microhaplotypes

MhGeneS: An Analytical Pipeline to Allow for Robust Microhaplotype Genotyping

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

Screening a new set of microhaplotypes in exonic regions for sample identity testing and paternity testing during whole exome sequencing analysis

Whole exome sequencing (WES) is widely used in clinical diagnosis. Before obtaining an accurate diagnosis, it is essential to conduct sample identity testing and paternity testing on trio samples. Currently, there is a lack of optimal genetic markers for these purposes, with limited literature available in this area. Microhaplotypes (MHs) are promising genetic markers due to their high polymorphism, low mutation rate, short amplified fragments, absence of stutter and amplification bias. These characteristics make them suitable for sample tracking and paternity testing during WES analysis. In this study, we screened out a set of polymorphic MHs in exonic regions for the above purposes. The results showed that the power of discrimination (PD) and probability of exclusion (PE) of this set of markers ranged from 0.2682 to 0.8878 and 0.0178 to 0.4583, respectively. Both the cumulative power of discrimination (CPD) and cumulative probability of exclusion (CPE) exceeded 0.999999, indicating the great value of these markers in paternity testing and individual identification in the study population. However, these markers had the effective number of alleles (Ae) values ranging from 1.1784 to 3.8727 (average 2.1805) and informativeness (In) values ranging from 0.0151 to 0.2209 (average 0.0766), showing limited value in DNA mixture analysis and biogeographical ancestry inference.

A free lunch: microhaplotype discovery in an existing amplicon panel improves parentage assignment for the highly polymorphic Pacific oyster

Amplicon panels using genotyping by sequencing methods are now common, but have focused on characterizing SNP markers. We investigate how microhaplotype (MH) discovery within a recently developed Pacific oyster (Magallana gigas) amplicon panel could increase the statistical power for relationship assignment. Trios (offspring and two parents) from three populations in a newly established breeding program were genotyped on a 592 locus panel. After processing, 92% of retained amplicons contained polymorphic MH variants and 85% of monomorphic SNP markers contained MH variation. The increased allelic richness resulted in substantially improved power for relationship assignment with much lower estimated false positive rates. No substantive differences in assignment accuracy occurred between SNP and MH datasets, but using MHs increased the separation in log-likelihood values between true parents and highly related potential parents (aunts and uncles). A high number of Mendelian incompatibilities among trios were observed, likely due to null alleles. Further development of a MH panel, including removing loci with high rates of null alleles, would enable high-throughput genotyping by reducing panel size and therefore cost for Pacific oyster research and breeding programs.

Development and validation of a multiplex panel with 232 microhaplotypes and software for forensic kinship analysis

In this study, we developed and validated a novel microhaplotype (MH) panel, the FGID Microhaplotype Kit, which contains 232 loci and was specifically designed for forensic kinship analysis. The performance of the panel was evaluated through rigorous testing that included sensitivity, species specificity, inhibitor resistance, uniformity, stability, accuracy and mixture deconvolution. The results showed that the kit is capable of reliably detecting all loci with minimal DNA input. It showed high species specificity for 12 non-human DNA samples and resistance to common inhibitors. In addition, forensic statistical analysis revealed a combined discriminatory power (cDP) of 1-1.68e-223 and superior combined exclusion power for duo and trio cases compared to standard STR panels. The panel was also tested for kinship analyzes with simulated and real pedigree samples and showed significantly higher likelihood ratios (LR) for detecting relationships between parents and offspring, full siblings, half siblings and first cousins, especially for more distant kinship types where conventional STR panels have difficulties. Using the FGID kinship software with the MH panel significantly improved the accuracy of kinship analysis, allowing even closely related individuals to be effectively discriminated while reducing the number of false negatives. In addition, principal component analysis (PCA) showed that the panel can distinguish the major world populations and East Asian subpopulations. Taken together, these results suggest that the FGID Microhaplotype Kit and associated software provide an efficient and accurate solution for forensic kinship analysis that offers better discriminatory power and reliability than traditional STR-based methods.

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

BACKGROUND

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

RESULTS AND DISCUSSION

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

CONCLUSIONS

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

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

Using simulated microhaplotype genotyping data to evaluate the value of machine learning algorithms for inferring DNA mixture contributor numbers

Inferring the number of contributors (NoC) is a crucial step in interpreting DNA mixtures, as it directly affects the accuracy of the likelihood ratio calculation and the assessment of evidence strength. However, obtaining the correct NoC in complex DNA mixtures remains challenging due to the high degree of allele sharing and dropout. This study aimed to analyze the impact of allele sharing and dropout on NoC inference in complex DNA mixtures when using microhaplotypes (MH). The effectiveness and value of highly polymorphic MH for NoC inference in complex DNA mixtures were evaluated through comparing the performance of three NoC inference methods, including maximum allele count (MAC) method, maximum likelihood estimation (MLE) method, and random forest classification (RFC) algorithm. In this study, we selected the top 100 most polymorphic MH from the Southern Han Chinese (CHS) population, and simulated over 40 million complex DNA mixture profiles with the NoC ranging from 2 to 8. These profiles involve unrelated individuals (RM type) and related pairs of individuals, including parent-offspring pairs (PO type), full-sibling pairs (FS type), and second-degree kinship pairs (SE type). Our results indicated that how the number of detected alleles in DNA mixture profiles varied with the markers' polymorphism, kinship's involvement, NoC, and dropout settings. Across different types of DNA mixtures, the MAC and MLE methods performed best in the RM type, followed by SE, FS, and PO types, while RFC models showed the best performance in the PO type, followed by RM, SE, and FS types. The recall of all three methods for NoC inference were decreased as the NoC and dropout levels increased. Furthermore, the MLE method performed better at low NoC, whereas RFC models excelled at high NoC and/or high dropout levels, regardless of the availability of a priori information about related pairs of individuals in DNA mixtures. However, the RFC models which considered the aforementioned priori information and were trained specifically on each type of DNA mixture profiles, outperformed RFC_ALL model that did not consider such information. Finally, we provided recommendations for model building when applying machine learning algorithms to NoC inference.

A proof-of-principle study: The potential application of MiniHap biomarkers in ancestry inference based on the QNome nanopore sequencing

Haplotyped SNPs convey forensic-related information, and microhaplotypes (MHs), as the most representative of this kind of marker, have proved the potential value for human forensics. In recent years, nanopore sequencing technology has developed rapidly, with its outstanding ability to sequence long continuous DNA fragments and obtain phase information, making the detection of longer haplotype marker possible. In this proof-of-principle study, we proposed a new type of forensic marker, MiniHap, based on five or more SNPs within a molecular distance less than 800 bp, and investigated the haplotype data of 56 selected MiniHaps in five Chinese populations using the QNome nanopore sequencing. The sequencing performance, allele (haplotype) frequencies, forensic parameters, effective number of alleles (Ae), and informativeness (In) were subsequently calculated. In addition, we performed principal component analysis (PCA), phylogenetic tree, and structure analysis to investigate the population genetic relationships and ancestry components among the five investigated populations and 26 worldwide populations. MiniHap-04 exhibited remarkable forensic efficacy, with 148 haplotypes reported and the Ae was 66.9268. In addition, the power of discrimination (PD) was 0.9934, the probability of exclusion (PE) was 0.9898, and the In value was 0.7893. Of the 56 loci, 85.71% had PD values above 0.85, 66.07% had PE values above 0.54, 67.86% had Ae values over 7.0%, and 55.36% were with In values above 0.2 across all samples, indicating that most of the MiniHaps are suitable for individual identification, paternity testing, mixture deconvolution, and ancestry inference. Moreover, the results of PCA, phylogenetic tree and structure analysis demonstrated that this MiniHap panel had the competency in continental population ancestry inference, but the differentiation within intracontinental/linguistically restricted subpopulations was not ideal. Such findings suggested that the QNome device for MiniHap detection was feasible and this novel marker has the potential in ancestry inference. Yet, the establishment of a more comprehensive database with sufficient reference population data remains necessary to screen more suitable MiniHaps.

Adapting an established Ampliseq microhaplotype panel to nanopore sequencing through direct PCR

We have adapted an established Ampliseq microhaplotype panel for nanopore sequencing with the Oxford Nanopore Technologies (ONT) system, as a cost-effective and highly scalable solution for forensic genetics applications. For this purpose, we designed a protocol combining direct PCR amplification from unextracted DNA with ONT library construction and sequencing using the MinION device and workflow. The analysis of reference samples at input amounts of 5-10 ng of DNA demonstrates stable coverage patterns, allele balance, and strand bias, reaching profile completeness and concordance rates of ∼95%. Similar levels were achieved when using direct-PCR from blood, buccal and semen swabs. Dilution series results indicate sensitivity is maintained down to 250 pg of input DNA, and informative profiles are produced down to 62.5 pg. Finally, we demonstrated the forensic utility of the nanopore workflow by analyzing two third degree pedigrees that showed low likelihood ratio values after the analysis of an extended panel of 38 STRs, achieving likelihood ratios 2-3 orders of magnitude higher when testing with the MinION-based haplotype data.

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.

Massively parallel sequencing of 74 microhaplotypes and forensic characteristics in three Chinese Sino-Tibetan populations

Microhaplotype (MH), as an emerging type of forensic genetic marker in recent years, has the potential to support multiple forensic applications, especially for mixture deconvolution and biogeographic ancestry inference. Herein, we investigated the genotype data of 74 MHs included in a novel MH panel, the Ion AmpliSeq MH-74 Plex Microhaplotype Research Panel, in three Chinese Sino-Tibetan populations (Han, Tibetan, and Yi) using the Ion Torrent semiconductor sequencing. The sequencing performance, allele frequencies, effective number of alleles (Ae), informativeness (In), and forensic parameters were subsequently estimated and calculated. In addition, principal component analysis (PCA) and structure analysis were performed to explore the population relationships among the three populations and the ancestry component distribution. Overall, this novel MH panel is robust and reliable, and has an excellent sequencing performance. The Ae values ranged from 1.0126 to 7.0855 across all samples, and 75.68 % of MHs had Ae values >2.0000. Allele frequencies at some loci varied considerably among the three studied populations, and the mean In value was 0.0195. Moreover, the genetic affinity between Tibetans and Yis was closer than that between Tibetans and Hans. The aforementioned results suggest that the Ion AmpliSeq MH-74 Plex Microhaplotype Research Panel is highly polymorphic in three investigated populations and could be used as an effective tool for human forensics. Although these 74 MHs have demonstrated the competency in continental population stratification, a higher resolution for distinguishing intracontinental subpopulations and a more comprehensive database with sufficient reference population data still remain to be accomplished.

Haplotype based testing for a better understanding of the selective architecture

BACKGROUND

The identification of genomic regions affected by selection is one of the most important goals in population genetics. If temporal data are available, allele frequency changes at SNP positions are often used for this purpose. Here we provide a new testing approach that uses haplotype frequencies instead of allele frequencies.

RESULTS

Using simulated data, we show that compared to SNP based test, our approach has higher power, especially when the number of candidate haplotypes is small or moderate. To improve power when the number of haplotypes is large, we investigate methods to combine them with a moderate number of haplotype subsets. Haplotype frequencies can often be recovered with less noise than SNP frequencies, especially under pool sequencing, giving our test an additional advantage. Furthermore, spurious outlier SNPs may lead to false positives, a problem usually not encountered when working with haplotypes. Post hoc tests for the number of selected haplotypes and for differences between their selection coefficients are also provided for a better understanding of the underlying selection dynamics. An application on a real data set further illustrates the performance benefits.

CONCLUSIONS

Due to less multiple testing correction and noise reduction, haplotype based testing is able to outperform SNP based tests in terms of power in most scenarios.

From chip to SNP: Rapid development and evaluation of a targeted capture genotyping-by-sequencing approach to support research and management of a plaguing rodent

The management of invasive species has been greatly enhanced by population genetic analyses of multilocus single-nucleotide polymorphism (SNP) datasets that provide critical information regarding pest population structure, invasion pathways, and reproductive biology. For many applications there is a need for protocols that offer rapid, robust and efficient genotyping on the order of hundreds to thousands of SNPs, that can be tailored to specific study populations and that are scalable for long-term monitoring schemes. Despite its status as a model laboratory species, there are few existing resources for studying wild populations of house mice (Mus musculus spp.) that strike this balance between data density and laboratory efficiency. Here we evaluate the utility of a custom targeted capture genotyping-by-sequencing approach to support research on plaguing house mouse populations in Australia. This approach utilizes 3,651 hybridization capture probes targeting genome-wide SNPs identified from a sample of mice collected in grain-producing regions of southeastern Australia genotyped using a commercially available microarray platform. To assess performance of the custom panel, we genotyped wild caught mice (N = 320) from two adjoining farms and demonstrate the ability to correctly assign individuals to source populations with high confidence (mean >95%), as well as robust kinship inference within sites. We discuss these results in the context of proposed applications for future genetic monitoring of house mice in Australia.

Exploring the forensic effectiveness and population genetic differentiation by self-constructed 41 multi-InDel panel in Yunnan Zhuang group

Yunnan is one of the main residences of the Zhuang group which is one of the 55 ethnic minorities in China. At present, there are relatively few researches on population genetics and forensic science of the Yunnan Zhuang group. Therefore, this study used a self-constructed panel containing 41 multi-InDel markers to analyze the genetic polymorphisms of 173 individuals from Yunnan Zhuang group. The results indicated that these 41 multi-InDels in Yunnan Zhuang group were highly polymorphic markers expect for three markers. The cumulative match probability and combined exclusion probability values of the 40 multi-InDels (MI38 marker was excluded) were 8.0671E-26 and 0.9999995959, respectively. In addition, population genetic analyses were performed on genotyping data of 41 multi-InDel markers among the Yunnan Zhuang and 26 reference populations, revealing that the Yunnan Zhuang group was genetically close to the five populations in East Asia. According to the STRUCTURE analysis, the Yunnan Zhuang group presented similar ancestral compositions to the five populations from East Asia, and when the K value was three, the five intercontinental populations showed their different genetic structures. In conclusion, the 41 multi-InDel markers could be used as an effective tool for individual identification and paternity testing of the Zhuang group in Yunnan province, as well as for their ancestry information inference studies.

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.

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.