Evaluating 130 microhaplotypes across a global set of 83 populations

Assessing the performance of multi-InDel panels for human identification among admixed Brazilians

Insertion/deletion polymorphisms, or InDels, are widely present in the human genome. They have been considered as potential markers for forensic analysis because they can be genotyped using the CE platform and compatible typing techniques used in forensic laboratories. Additionally, InDels have lower mutation rates and often short amplicon sizes, making them ideal for detecting degraded samples. However, most InDels are bi-allelic; therefore, their discrimination power is relatively low. A new set of genetic marker called multi-InDels was reported to improve InDel's informativeness. Multi-InDel markers are generally designated as microhaplotypes encompassing two or more InDels within a short distance, usually less than 200 bp. In this study, we evaluated the applicability of three previously proposed panels of multi-InDel markers, designed for Asian populations, for human identification in Brazil. We assessed all the multi-InDel markers using high-coverage whole-genome sequencing data from a census-based cohort of 1171 Brazilians residing in São Paulo, the largest Brazilian capital. The results showed that most markers are informative for Brazilian individuals since they present more than three frequent haplotypes with different sizes. However, most markers are prone to amplification/sequencing errors due to repetitive or low-complexity regions. Among the tested panels, the one from Huang et al. (2014) is the most promising for forensic use in Brazil, with a combined match probability and cumulative power of exclusion of 4.92 ×10-14 and 0.9991, respectively. Nevertheless, these values are low compared to the ones obtained with CODIS STRs (short tandem repeats) and larger SNP (single nucleotide polymorphisms) panels. Therefore, new attempts to scan the human genome for highly polymorphic multi-InDel markers are still necessary to obtain a suitable panel of multi-InDels for worldwide populations.

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.

Human complex mixture analysis by "FD Multi-SNP Mixture Kit"

Introduction

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

Methods

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

Results

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

Conclusion

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

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.

Performance comparison of a previously validated microhaplotype panel and a forensic STR panel for DNA mixture analysis

Short Tandem Repeats (STRs) are the most widespread markers in forensic genetics. However, STR stutter peaks can mask alleles from a minor contributor when analysing mixtures, hindering the interpretation of complex profiles. In this study we compared the performance of a previously described panel of microhaplotypes (MHs), an alternative type of forensic marker, against a standard STR kit. The parameters evaluated included: capability of determining the minimum number of contributors in the mixture; percentages of allele drop-outs and drop-ins; retrieval of alleles belonging to the minor contributor, and estimation of likelihood ratio (LR) values. In addition, the capacity of EuroForMix software to estimate each donor's percentage of contribution was tested, as well as the impact on results when using manually, or automatically prepared libraries. The MH panel showed better performance than STRs for the detection of 2-contributor mixtures, but the lower degree of polymorphism per MH marker hindered the task of deconvolution with multiple contributors. MHs presented higher drop-in rates and lower drop-out rates, a higher capability to recover the minor contributor's alleles and provided higher LR values than STRs, likely due to the much higher number of loci combined in the panel. Estimations of contributor ratios using EuroForMix showed promising results and marginal differences were found in these values between manually and automatically prepared libraries. Overall, results showed that the mixture detection performance of the MH panel was better or equal to the standard forensic autosomal STR panel, indicating microhaplotypes are informative markers for this purpose.

Genetic diversity of North African populations in the 17q21 genomic region

The demographic history of human populations in North Africa has been characterized by complex migration processes that have determined the current genetic structure of these populations. We examined the autosomal markers of eight sampled populations in northern Africa (Tunisia and Libya) to explore their genetic structure and to place them in a global context. We genotyped a set of 30 autosomal single-nucleotide polymorphisms (SNPs) extending 9.5 Mb and encompassing the 17q21 inversion region. Our data include 403 individuals from Tunisia and Libya. To put our populations in the global context, we analyzed our data in comparison with other populations, including those of the 1000 Genomes Project. To evaluate the data, we conducted genetic diversity, principal component, STRUCTURE, and haplotype analyses. The analysis of genetic composition revealed the genetic heterogeneity of North African populations. The principal component and STRUCTURE analyses converged and revealed the intermediate position of North Africans between Europeans and Asians. Haplotypic analysis demonstrated that the normal (H1) and inverted (H2) polymorphisms in the chromosome 17q21 region occur in North Africa at frequencies similar to those found in European and Southwest Asian populations. The results highlight the complex demographic history of North Africa, reflecting the influence of genetic flow from Europe and the Near East that dates to the prehistoric period. These gene flows added to demographic factors (inbreeding, endogamy), natural factors (topography, Sahara), and cultural factors that play a role in the emergence of the diverse and heterogeneous genetic structures of North African populations. This study contributes to a better understanding of the complex structure of North African populations.

MHBase: A comprehensive database of short microhaplotypes for advancing forensic genetic analysis

Microhaplotypes (MHs) were first recommended by Prof. Kidd for use in forensics because they can improve human identification, kinship analysis, mixture deconvolution, and ancestry prediction. Since their introduction, extensive research has demonstrated the advantages of MHs in forensic applications and provided useful data for different populations. Currently, two databases, ALFRED (ALlele FREquency Database) and MicroHapDB (MicroHaplotype DataBase), house the published MH information and population data. We previously constructed a single nucleotide polymorphism SNP-SNP MH database (D-SNPsDB) of MHs within 50 bp on the whole human genome for 26 populations integrating basic data such as physical genome positions, mapping of variant identifiers (rsIDs), allele frequencies, and basic variant information. Building upon the previous research, we further selected MHs containing at least two variants (SNPs and/or insertions/deletions [InDels]) within a short DNA fragment (≤ 50 bp) in 26 populations based on the 1000 Genomes Project dataset (Phase 3) to construct a more comprehensive database. Subsequently, we established a user-friendly website that allows users to search the MH database (MHBase) based on their research objectives and study population to find suitable loci and provides other functions such as querying reported loci, performing online calculations using the PHASE software, and calculating ancestral-related parameters. The loci in the database are classified as SNP-based MHs, which include only SNPs, and InDel-including MHs, which contain at least one InDel. Here, we provide a detailed overview of the MHBase and an analysis of shared loci at the global and continental levels, ancestral markers, the genetic distance within loci, and mapping with the genome annotation file. The website is an accessible and useful tool for researchers engaged in marker discovery, population studies, assay development, and panel design.

Identification of the body fluid donor in mixtures through target mRNA cSNP sequencing

In forensic practice, mixture stains containing various body fluids are common, presenting challenges for interpretation, particularly in multi-contributor mixtures. Traditional STR profiles face difficulties in such scenarios. Over recent years, RNA has emerged as a promising biomarker for body fluid identification, and mRNA polymorphism has shown excellent performance in identifying body fluid donors in previous studies. In this study, a massively parallel sequencing assay was developed, encompassing 202 coding region SNPs (cSNPs) from 45 body fluid/tissue-specific genes to identify both body fluid/tissue origin and the respective donors, including blood, saliva, semen, vaginal secretion, menstrual blood, and skin. The specificity was evaluated by examining the single-source body fluids/tissue and revealed that the same body fluid exhibited similar expression profiles and the tissue origin could be identified. For laboratory-generated mixtures containing 2-6 different components and mock case mixtures, the donor of each component could be successfully identified, except for the skin donor. The discriminatory power for all body fluids ranged from 0.997176329 (menstrual blood) to 0.99999999827 (blood). The concordance of DNA typing and mRNA typing for the cSNPs in this system was also validated. This cSNP typing system exhibits excellent performance in mixture deconvolution.

Investigating the effectiveness of forensic genetics and population genetic diversity using a multi-InDel system in Chinese Hezhou and Southern Shaanxi Han populations

INTRODUCTION

Multiple insertion-deletion (multi-InDel) has greater potential in forensic genetics than InDel, and its efficacy in kinship testing, individual identification, DNA mixture detection and ancestry inference remains to be explored.

METHODS

Consequently, we designed an efficient and robust system consisting of 41 multi-InDels to evaluate its efficacy in forensic applications in Chinese Hezhou Han (HZH) and Southern Shaanxi Han (SNH) populations and explore the genetic relationships between the SNH, HZH, and 26 reference populations.

RESULTS AND CONCLUSION

The obtained results showed that 38 out of the 41 multi-InDels had fairly high genetic variations. The the cumulative probability of discrimination and exclusion values of the multi-InDels (except MI38) in HZH and SNH populations both exceeded 1-e-25 and 1-e-6, correspondingly. The genetic compositions of HZH and SNH individuals were similar to that of East Asians and the Naive Bayes model could well distinguish East Asians, Africans and Americans. These results indicated that the multi-InDel systerm can serve as an effective tool to provide important evidence for the development of multi-InDels in forensic practice and better analyse the genetic background of the Han Chinese populations.

Biogeographical Ancestry Analyses Using the ForenSeqTM DNA Signature Prep Kit and Multiple Prediction Tools

The inference of biogeographical ancestry (BGA) can assist in police investigations of serious crime cases and help to identify missing people and victims of mass disasters. In this study, we evaluated the typing performance of 56 ancestry-informative SNPs in 177 samples using the ForenSeq™ DNA Signature Prep Kit on the MiSeq FGx system. Furthermore, we compared the prediction accuracy of the tools Universal Analysis Software v1.2 (UAS), the FROG-kb, and GenoGeographer when inferring the ancestry of 503 Europeans, 22 non-Europeans, and 5 individuals with co-ancestry. The kit was highly sensitive with complete aiSNP profiles in samples with as low as 250pg input DNA. However, in line with others, we observed low read depth and occasional drop-out in some SNPs. Therefore, we suggest not using less than the recommended 1ng of input DNA. FROG-kb and GenoGeographer accurately predicted both Europeans (99.6% and 91.8% correct, respectively) and non-Europeans (95.4% and 90.9% correct, respectively). The UAS was highly accurate when predicting Europeans (96.0% correct) but performed poorer when predicting non-Europeans (40.9% correct). None of the tools were able to correctly predict individuals with co-ancestry. Our study demonstrates that the use of multiple prediction tools will increase the prediction accuracy of BGA inference in forensic casework.

Developmental validation of a custom-designed Multi-InDel panel: A five-dye multiplex amplification system for challenging DNA samples

The continual investigation of novel genetic markers has yielded promising solutions for addressing the challenges encountered in forensic DNA analysis. In this study, we have introduced a custom-designed panel capable of simultaneously amplifying 41 novel Multi-insertion/deletion (Multi-InDel) markers and an amelogenin locus using the capillary electrophoresis platform. Through a developmental validation study conducted in accordance with guidelines recommended by the Scientific Working Group on DNA Analysis Methods, we demonstrated that the new Multi-InDel system exhibited the sensitivity to produce reliable genotyping profiles with as little as 62.5 pg of template DNA. Accurate and complete genotyping profiles could be obtained even in the presence of specific concentrations of PCR inhibitors. Furthermore, the maximum amplicon size for this system was limited to under 220 bp in the genotyping profile, resulting in its superior efficiency compared to commercially available short tandem repeat kits for both naturally and artificially degraded samples. In the context of mixed DNA analysis, the Multi-InDel system was proved informative in the identification of two-person DNA mixture, even when the template DNA of the minor contributor was as low as 50 pg. In conclusion, a series of performance evaluation studies have provided compelling evidence that the new Multi-InDel system holds promise as a valuable tool for forensic DNA 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.

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.

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.

Genetic investigation of 124 SNPs in a Myanmar population using the Precision ID Identity Panel and the Illumina MiSeq

Single nucleotide polymorphisms (SNPs) have become popular in forensic genetics as an alternative to short tandem repeats (STRs). The Precision ID Identity Panel (Thermo Fisher Scientific), consisting of 90 autosomal SNPs and 34 Y-chromosomal SNPs, enabled human identification studies on global populations through next-generation sequencing (NGS). However, most previous studies on the panel have used the Ion Torrent platform, and there are few reports on the Southeast Asian population. Here, a total of 96 unrelated males from Myanmar (Yangon) were analyzed with the Precision ID Identity Panel on a MiSeq (Illumina) using an in-house TruSeq compatible universal adapter and a custom variant caller, Visual SNP. The sequencing performance evaluated by locus balance and heterozygote balance was comparable to that of the Ion Torrent platform. For 90 autosomal SNPs, the combined match probability (CMP) was 6.994 × 10^-34, lower than that of 22 PowerPlex Fusion autosomal STRs (3.130 × 10^-26). For 34 Y-SNPs, 14 Y-haplogroups (mostly O2 and O1b) were observed. We found 51 cryptic variations (42 haplotypes) around target SNPs, of which haplotypes corresponding to 33 autosomal SNPs decreased CMP. Interpopulation analysis revealed that the Myanmar population is genetically closer to the East and Southeast Asian populations. In conclusion, the Precision ID Identity Panel can be successfully analyzed on the Illumina MiSeq and provides high discrimination power for human identification in the Myanmar population. This study broadened the accessibility of the NGS-based SNP panel by expanding the available NGS platforms and adopting a robust NGS data analysis tool.

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 DNA Phenotyping of appearance, ancestry and age

Forensic DNA Phenotyping (FDP) comprises the prediction of a person's externally visible characteristics regarding appearance, biogeographic ancestry and age from DNA of crime scene samples, to provide investigative leads to help find unknown perpetrators that cannot be identified with forensic STR-profiling. In recent years, FDP has advanced considerably in all of its three components, which we summarize in this review article. Appearance prediction from DNA has broadened beyond eye, hair and skin color to additionally comprise other traits such as eyebrow color, freckles, hair structure, hair loss in men, and tall stature. Biogeographic ancestry inference from DNA has progressed from continental ancestry to sub-continental ancestry detection and the resolving of co-ancestry patterns in genetically admixed individuals. Age estimation from DNA has widened beyond blood to more somatic tissues such as saliva and bones as well as new markers and tools for semen. Technological progress has allowed forensically suitable DNA technology with largely increased multiplex capacity for the simultaneous analysis of hundreds of DNA predictors with targeted massively parallel sequencing (MPS). Forensically validated MPS-based FDP tools for predicting from crime scene DNA i) several appearance traits, ii) multi-regional ancestry, iii) several appearance traits together with multi-regional ancestry, and iv) age from different tissue types, are already available. Despite recent advances that will likely increase the impact of FDP in criminal casework in the near future, moving reliable appearance, ancestry and age prediction from crime scene DNA to the level of detail and accuracy police investigators may desire, requires further intensified scientific research together with technical developments and forensic validations as well as the necessary funding.

Systematically exploring the performance of a self-developed Multi-InDel system in forensic identification, ancestry inference and genetic structure analysis of Chinese Manchu and Mongolian groups

The insertion/deletion (InDel) polymorphism has promising applications in forensic DNA analysis. However, the insufficient forensic efficiencies of the present InDel-based systems restrict their applications in parentage testing, due to the lower genetic polymorphism of the biallelic InDel locus and the limited number of InDel loci in a multiplex amplification system. Here, we introduced an in-house developed system which contained 41 polymorphic Multi-InDel markers (equivalent to 82 InDels in total), to serve as an efficient and reliable tool for different forensic applications in the Manchu and Mongolian groups. We demonstrated that the new system exhibited potential efficiencies for personal identification, parentage testing, two-person DNA mixture interpretation and ancestry inference of intercontinental populations. Meanwhile, we explored the genetic backgrounds of the Manchu and Mongolian groups by conducting a series of population genetic analyses. We showed that the Manchu and Mongolian groups shared closer genetic relationships with East Asian populations, especially Han Chinese populations in northern China. Moreover, more similar genetic compositions were detected between the Manchu group and the northern Han populations in this study, suggesting that the Manchu group had higher genetic affinities with northern Han populations than the Mongolian group. Overall. this study provided the necessary evidence that these Multi-InDel genetic markers could play an important role in forensic applications.

Development and evaluations of the ancestry informative markers of the VISAGE Enhanced Tool for Appearance and Ancestry

The VISAGE Enhanced Tool for Appearance and Ancestry (ET) has been designed to combine markers for the prediction of bio-geographical ancestry plus a range of externally visible characteristics into a single massively parallel sequencing (MPS) assay. We describe the development of the ancestry panel markers used in ET, and the enhanced analyses they provide compared to previous MPS-based forensic ancestry assays. As well as established autosomal single nucleotide polymorphisms (SNPs) that differentiate sub-Saharan African, European, East Asian, South Asian, Native American, and Oceanian populations, ET includes autosomal SNPs able to efficiently differentiate populations from Middle East regions. The ability of the ET autosomal ancestry SNPs to distinguish Middle East populations from other continentally defined population groups is such that characteristic patterns for this region can be discerned in genetic cluster analysis using STRUCTURE. Joint cluster membership estimates showing individual co-ancestry that signals North African or East African origins were detected, or cluster patterns were seen that indicate origins from central and Eastern regions of the Middle East. In addition to an augmented panel of autosomal SNPs, ET includes panels of 85 Y-SNPs, 16 X-SNPs and 21 autosomal Microhaplotypes. The Y- and X-SNPs provide a distinct method for obtaining extra detail about co-ancestry patterns identified in males with admixed backgrounds. This study used the 1000 Genomes admixed African and admixed American sample sets to fully explore these enhancements to the analysis of individual co-ancestry. Samples from urban and rural Brazil with contrasting distributions of African, European, and Native American co-ancestry were also studied to gauge the efficiency of combining Y- and X-SNP data for this purpose. The small panel of Microhaplotypes incorporated in ET were selected because they showed the highest levels of haplotype diversity amongst the seven population groups we sought to differentiate. Microhaplotype data was not formally combined with single-site SNP genotypes to analyse ancestry. However, the haplotype sequence reads obtained with ET from these loci creates an effective system for de-convoluting two-contributor mixed DNA. We made simple mixture experiments to demonstrate that when the contributors have different ancestries and the mixture ratios are imbalanced (i.e., not 1:1 mixtures) the ET Microhaplotype panel is an informative system to infer ancestry when this differs between the contributors.

Comparative evaluation of the MAPlex, Precision ID Ancestry Panel, and VISAGE Basic Tool for biogeographical ancestry inference

Biogeographical ancestry (BGA) inference from ancestry-informative markers (AIMs) has strong potential to support forensic investigations. Over the past two decades, several forensic panels composed of AIMs have been developed to predict ancestry at a continental scale. These panels typically comprise fewer than 200 AIMs and have been designed and tested with a limited set of populations. How well these panels recover patterns of genetic diversity relative to larger sets of markers, and how accurately they infer ancestry of individuals and populations not included in their design remains poorly understood. The lack of comparative studies addressing these aspects makes the selection of appropriate panels for forensic laboratories difficult. In this study, the model-based genetic clustering tool STRUCTURE was used to compare three popular forensic BGA panels: MAPlex, Precision ID Ancestry Panel (PIDAP), and VISAGE Basic Tool (VISAGE BT) relative to a genome-wide reference set of 10k SNPs. The genotypes for all these markers were obtained for a comprehensive set of 3957 individuals from 228 worldwide human populations. Our results indicate that at the broad continental scale (K=6) typically examined in forensic studies, all forensic panels produced similar genetic structure patterns compared to the reference set (G'≈90%) and had high classification performance across all regions (average AUC-PR > 97%). However, at K= 7 and K= 8, the forensic panels displayed some region-specific clustering deviations from the reference set, particularly in Europe and the region of East and South-East Asia, which may be attributed to differences in the design of the respective panels. Overall, the panel with the most consistent performance in all regions was VISAGE BT with an average weighted AUC̅_W score of 96.26% across the three scales of geographical resolution investigated.

Evaluation of a SNP-STR haplotype panel for forensic genotype imputation

Short tandem repeat polymorphism (STR)-based individual identification is a popular and reliable method in many forensic applications. However, STRs still frequently fail to find any matched records. In such cases, if known STRs could provide more information, it would be very helpful to solve specific problems. Genotype imputation has long been used in the study of single nucleotide polymorphisms (SNPs) and has recently been introduced into forensic fields. The idea is that, through a reference haplotype panel containing SNPs and STRs, we can obtain unknown genetic information through genotype imputation based on known STR or SNP genotypes. Several recent studies have already demonstrated this exciting idea, and a 1000 Genomes SNP-STR haplotype panel has also been released. To further study the performance of genotype imputation in forensic fields, we collected STR, microhaplotype (MH) and SNP array genotypes from Chinese Han population individuals and then performed genotype imputation analysis based on the released reference panel. As a result, the average locus imputation accuracy was ∼83 % (or ∼70 %) when SNPs in the SNP array (or MH SNPs) were imputed from STRs, and was ∼30 % when highly polymorphic markers (STRs and MHs) were imputed from each other. When STRs were imputed from SNP array, the average locus imputation accuracy increased to ∼48 %. After analyzing the match scores between real STRs and the STRs imputed from SNPs, ∼80 % of studied STR records can be connected to corresponding SNP records, which may help for individual identification. Our results indicate that genotype imputation has great potential for forensic applications.

A sequence-based 163plex microhaplotype assay for forensic DNA analysis

Novel genetic marker microhaplotype has led to an upsurge in forensic genetic research. This study established a 163 microhaplotype (MH) multiplex assay based on next-generation sequencing (NGS) and evaluated the assay's performance and applicability. Our results showed that the 163 MH assay was accurate, repeatable and reliable, and could distinguish between African, European-American, Southern Asia and Eastern Asia populations. Among the 163 MH makers, 48 MHs with Ae > 3.0 in China Eastern Han were selected and confirmed to be highly polymorphic, with a combined power of discrimination of 1-8.26 × 10-44 and the combined power of exclusion in duos and trios of 1-1.26 × 10-8 and 1-8.27 × 10-16, respectively. Moreover, the mixture study demonstrated the realizability of the MHs in deconvoluting mixtures with different proportions of two to five-person. In conclusion, our findings support the use of this MH assay for ancestry inference, human identification, paternity testing and mixture deconvolution in forensic research.

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.

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.

State of the Art for Microhaplotypes

In recent years, the number of publications on microhaplotypes has averaged more than a dozen papers annually. Many have contributed to a significant increase in the number of highly polymorphic microhaplotype loci. This increase allows microhaplotypes to be very informative in four main areas of forensic uses of DNA: individualization, ancestry inference, kinship analysis, and mixture deconvolution. The random match Probability (RMP) can be as small as 10−100 for a large panel of microhaplotypes. It is possible to measure the heterozygosity of an MH as the effective number of alleles (Ae). Ae > 7.5 exists for African populations and >4.5 exists for Native American populations for a smaller panel of two dozen selected microhaplotypes. Using STRUCTURE, at least 10 different ancestral clusters can be defined by microhaplotypes. The Ae for a locus is also identical to the Paternity Index (PI), the measure of how informative a locus will be in parentage testing. High Ae loci can also be useful in missing persons cases. Finally, high Ae microhaplotypes allow the near certainty of seeing multiple additional alleles in a mixture of two or more individuals in a DNA sample. In summary, a panel of higher Ae microhaplotypes can outperform the standard CODIS markers.

Set of 15 SNP-SNP Markers for Detection of Unbalanced Degraded DNA Mixtures and Noninvasive Prenatal Paternity Testing

Unbalanced and degraded mixtures (UDM) are very common in forensic DNA analysis. For example, DNA signals from criminal suspects are masked by a large amount of DNA from victims, or cell-free fetal DNA (cffDNA) in maternal plasma is masked by a high background of maternal DNA. Currently, detecting minor DNA in these mixtures is complex and challenging. We developed a new set of SNP-SNP microhaplotypes with short amplicons, and we successfully genotyped them using the new method of amplification-refractory mutation system PCR (ARMS-PCR) combined with SNaPshot technology based on a capillary electrophoresis (CE) platform. This panel reflects a high polymorphism in the Southwest Chinese Han population and thus has excellent potential for mixture studies. We evaluated the feasibility of this panel for UDM detection and noninvasive prenatal paternity testing (NIPPT). Fifteen SNP-SNPs detected minor DNA of homemade DNA mixtures, with a sensitivity of 0.025-0.05 ng and a specificity of 1:1,000. In addition, the panel successfully genotyped degraded DNA from single and mixed samples. Finally, 15 SNP-SNPs were applied to 26 trios. All samples displayed positive results with at least one marker to detect cffDNA. Besides, all fetal alleles in maternal plasma were confirmed by genotyping fetal genomic DNA from amniocentesis and paternal genomic DNA from peripheral blood. The results indicated that the SNP-SNP strategy based on the CE platform was useful for UDM detection and NIPPT.

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.

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

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

Enhanced mixture interpretation with macrohaplotypes based on long-read DNA sequencing

Deconvoluting mixture samples is one of the most challenging problems confronting DNA forensic laboratories. Efforts have been made to provide solutions regarding mixture interpretation. The probabilistic interpretation of Short Tandem Repeat (STR) profiles has increased the number of complex mixtures that can be analyzed. A portion of complex mixture profiles, particularly for mixtures with a high number of contributors, are still being deemed uninterpretable. Novel forensic markers, such as Single Nucleotide Variants (SNV) and microhaplotypes, also have been proposed to allow for better mixture interpretation. However, these markers have both a lower discrimination power compared with STRs and are not compatible with CODIS or other national DNA databanks worldwide. The short-read sequencing (SRS) technologies can facilitate mixture interpretation by identifying intra-allelic variations within STRs. Unfortunately, the short size of the amplicons containing STR markers and sequence reads limit the alleles that can be attained per STR. The latest long-read sequencing (LRS) technologies can overcome this limitation in some samples in which larger DNA fragments (including both STRs and SNVs) with definitive phasing are available. Based on the LRS technologies, this study developed a novel CODIS compatible forensic marker, called a macrohaplotype, which combines a CODIS STR and flanking variants to offer extremely high number of haplotypes and hence very high discrimination power per marker. The macrohaplotype will substantially improve mixture interpretation capabilities. Based on publicly accessible data, a panel of 20 macrohaplotypes with sizes of ~ 8 k bp and the maximum high discrimination powers were designed. The statistical evaluation demonstrates that these macrohaplotypes substantially outperform CODIS STRs for mixture interpretation, particularly for mixtures with a high number of contributors, as well as other forensic applications. Based on these results, efforts should be undertaken to build a complete workflow, both wet-lab and bioinformatics, to precisely call the variants and generate the macrohaplotypes based on the LRS technologies.

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.

Autosomal STR and SNP characterization of populations from the Northeastern Peruvian Andes with the ForenSeq™ DNA Signature Prep Kit

Autosomal DNA data from Peru for human identity testing purposes are scarce in the scientific literature, which hinders obtaining an appropriate portrait of the genetic variation of the resident populations. In this study we genetically characterize five populations from the Northeastern Peruvian Andes (Chachapoyas, Awajún, Wampís, Huancas and Cajamarca). Autosomal short tandem repeat (aSTR) and identity informative single nucleotide polymorphism (iiSNP) data from a total of 233 unrelated individuals are provided, and forensic genetic parameters are calculated for each population and for the combined set Northeastern Peruvian Andes. After correction for multiple testing in the whole dataset of the Northeastern Peruvian Andes, the only departure from Hardy-Weinberg equilibrium was observed in locus rs2111980. Twenty one out of 27 aSTR loci exhibited an increased number of alleles due to sequence variation in the repeat motif and flanking regions. For iiSNPs 33% of the loci displayed flanking region variation. The combined random match probability (RMP), assuming independence of all loci (aSTRs and iiSNPs), in the Chachapoyas, the population with the largest samples size (N = 172), was 8.14 × 10^-62 for length-based data while for sequence-based was 4.15 × 10^-67. In the merged dataset (Northeastern Peruvian Andes; N = 233), the combined RMP when including all markers were 2.96 × 10^-61 (length-based) and 3.21 × 10^-66 (sequence-based). These new data help to fill up some of the gaps in the genetic canvas of South America and provide essential length- and sequence-based background information for other forensic genetic studies in Peru.

Considering the flanking region variants of nonbinary SNP and phenotype-informative SNP to constitute 30 microhaplotype loci for increasing the discriminative ability of forensic applications

The flanking region variants of nonbinary SNPs and phenotype-informative SNPs (piSNPs) have been observed, which may greatly improve the discriminative ability after constituting microhaplotype. In this study, 30 microhaplotype loci based on the nonbinary SNPs and piSNPs (shown to be related to phenotypes such as hair and eye color) were selected. Genotyping were conducted on 100 unrelated northern Han Chinese, and the 26 populations from the 1000 Genome Project were also included for comparison of populations differentiation. The simulated study was conducted for evaluating the efficiency of kinship testing. These 30 microhaplotype loci we selected had good polymorphism, with a mean effective number of alleles (Ae) of 3.46. The average Ae increase was 1.27 compared with the target SNPs. The populations from the five regions worldwide could also be distinguished using these loci. The results of kinship testing showed that these microhaplotype loci had the similar ability as 15 STR loci of AmpFlSTR^R Identifiler^R PCR Amplification Kit to identify the biological parent and a stronger ability to exclude the nonbiological parents. So, these 30 microhaplotype loci may be multifunctional for forensic application, including the ability of personal identification and kinship testing equivalent to 15 STR loci, and the power of ancestry inference for distinguishing the main intercontinental population. Moreover, our selected phenotypic microhaplotype loci may theoretically have phenotype prediction capabilities. But the phenotype prediction efficiency of these phenotypic microhaplotype loci may be worse than that of piSNPs and the detailed prediction accuracy of different populations needs to be further studied.

Evaluation of microhaplotypes in forensic kinship analysis from a Swedish population perspective

The development of massively parallel sequencing (MPS) technology has enabled the discovery of several new types of forensic markers where microhaplotypes are one of these promising novel genetic markers. Microhaplotypes are, commonly, less than 300 nucleotides in length and consist of two or more closely linked single-nucleotide polymorphisms (SNPs). In this study, we have examined a custom-made QIAseq Microhaplotype panel (Qiagen), including 45 different microhaplotype loci. DNA libraries were prepared according to the GeneRead DNAseq Targeted Panels V2 library preparation workflow (Qiagen) and sequenced on a MiSeq FGx instrument (Verogen). We evaluated the performance of the panel based on 75 samples of Swedish origin and haplotype frequencies were established. We performed sensitivity studies and could detect haplotypes at input amounts down to 0.8 ng. We also studied mixture samples with two contributors for which haplotypes, for the minor contributor, were detectable down to the level of 1:100. Furthermore, we executed kinship simulations to evaluate the usefulness of this panel in kinship analysis. The results showed that both paternity and full sibling cases can clearly be solved. When simulating a half sibling versus unrelated case scenario, there were, however, some overlap of the likelihood ratio distributions potentially resulting in inconclusiveness. To conclude, the results of this initial study are promising for further implementation of this microhaplotype assay into the forensic field, although we noticed some primer design issues that could be optimized, which possibly would increase the power of the assay.

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.