Inter-laboratory evaluation of the EUROFORGEN Global ancestry-informative SNP panel by massively parallel sequencing using the Ion PGM™

Phenotype predictions of two-person mixture using single cell analysis

Over a decade after the publication of the first forensic DNA phenotyping (FDP) studies, DNA-based appearance predictions are now becoming a reality in routine crime scene investigations. The significant number of publications dedicated to the subject of FDP clearly demonstrates a sustained interest and a strong need for further method development. However, the implementation of FDP in routine work still encounters obstacles, and one of these challenges is making phenotype predictions from DNA mixtures. In this study, we examined single-cell sequencing as a potential tool to enable reliable phenotyping of contributors within mixtures. Two mock mixtures, each containing two contributors with similar and different physical appearances, were analyzed using two different workflows. In the first workflow, the mixtures were sequenced using the Ion AmpliSeq™ PhenoTrivium Panel, which includes 41 HIrisPlex-S (HPS) markers. Subsequently, the genotypes were analyzed using the HPS Deconvolution Tool to predict the phenotypes of both contributors. The second workflow involved the introduction of single-cell separation and collection using the DEPArray™ PLUS System. Two different PhenoTrivium amplification protocols were tested, and the phenotype predictions from single cells were compared with the results obtained using the HPS Tool. Our results suggest that the approach presented here allows for the obtainment of nearly complete HIrisPlex-S profiles with accurate genotypes and reliable phenotype predictions from single cells. This method proves successful in deconvoluting mixtures submitted to forensic DNA phenotyping.

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.

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.

A custom hybridisation enrichment forensic intelligence panel to infer biogeographic ancestry, hair and eye colour, and Y chromosome lineage

Massively parallel sequencing can provide genetic data for hundreds to thousands of loci in a single assay for various types of forensic testing. However, available commercial kits require an initial PCR amplification of short-to-medium sized targets which limits their application for highly degraded DNA. Development and optimisation of large PCR multiplexes also prevents creation of custom panels that target different suites of markers for identity, biogeographic ancestry, phenotype, and lineage markers (Y-chromosome and mtDNA). Hybridisation enrichment, an alternative approach for target enrichment prior to sequencing, uses biotinylated probes to bind to target DNA and has proven successful on degraded and ancient DNA. We developed a customisable hybridisation capture method, that uses individually mixed baits to allow tailored and targeted enrichment to specific forensic questions of interest. To allow collection of forensic intelligence data, we assembled and tested a custom panel of hybridisation baits to infer biogeographic ancestry, hair and eye colour, and paternal lineage (and sex) on modern male and female samples with a range of self-declared ancestries and hair/eye colour combinations. The panel correctly estimated biogeographic ancestry in 9/12 samples (75%) but detected European admixture in three individuals from regions with admixed demographic history. Hair and eye colour were predicted correctly in 83% and 92% of samples respectively, where intermediate eye colour and blond hair were problematic to predict. Analysis of Y-chromosome SNPs correctly assigned sex and paternal haplogroups, the latter complementing and supporting biogeographic ancestry predictions. Overall, we demonstrate the utility of this hybridisation enrichment approach to forensic intelligence testing using a combined suite of biogeographic ancestry, phenotype, and Y-chromosome SNPs for comprehensive biological profiling.

Systematic selections and forensic application evaluations of 111 individual identification SNPs in the Chinese Inner Mongolia Manchu group

Single nucleotide polymorphism (SNP) possesses a promising application in forensic individual identification due to its wide distribution in the human genome and the ability to carry out the genotyping of degraded biological samples by designing short amplicons. Some commonly used individual identification SNPs are less polymorphic in East Asian populations. In order to improve the individual identification efficiencies in East Asian populations, SNP genetic markers with relatively higher polymorphisms were selected from the 1,000 Genome Project phase III database in East Asian populations. A total of 111 individual identification SNPs (II-SNPs) with the observed heterozygosity values greater than 0.4 were screened in East Asian populations, and then, the forensic efficiencies of these selected SNPs were also evaluated in Chinese Inner Mongolia Manchu group. The observed heterozygosity and power of discrimination values at 111 II-SNPs in the Inner Mongolia Manchu group ranged from 0.4011 to 0.7005, and 0.5620 to 0.8025, respectively, and the average value of polymorphism information content was greater than 0.3978. The cumulative match probability and combined probability of exclusion values at II-SNPs were 7.447E^-51 and 1-4.17E^-12 in the Inner Mongolia Manchu group, respectively. The accumulative efficiency results indicated that the set of II-SNPs could be used as a potential tool for forensic individual identification and parentage testing in the Manchu group. The sequencing depths ranged from 781× to 12374×. And the mean allele count ratio and noise level were 0.8672 and 0.0041, respectively. The sequencing results indicated that the SNP genetic marker detection based on the massively parallel sequencing technology for SNP genetic markers had high sequencing performance and could meet the sequencing requirements of II-SNPs in the studied group.

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.

Multivariate statistical approach and machine learning for the evaluation of biogeographical ancestry inference in the forensic field

The biogeographical ancestry (BGA) of a trace or a person/skeleton refers to the component of ethnicity, constituted of biological and cultural elements, that is biologically determined. Nowadays, many individuals are interested in exploring their genealogy, and the capability to distinguish biogeographic information about population groups and subgroups via DNA analysis plays an essential role in several fields such as in forensics. In fact, for investigative and intelligence purposes, it is beneficial to inference the biogeographical origins of perpetrators of crimes or victims of unsolved cold cases when no reference profile from perpetrators or database hits for comparative purposes are available. Current approaches for biogeographical ancestry estimation using SNPs data are usually based on PCA and Structure software. The present study provides an alternative method that involves multivariate data analysis and machine learning strategies to evaluate BGA discriminating power of unknown samples using different commercial panels. Starting from 1000 Genomes project, Simons Genome Diversity Project and Human Genome Diversity Project datasets involving African, American, Asian, European and Oceania individuals, and moving towards further and more geographically restricted populations, powerful multivariate techniques such as Partial Least Squares-Discriminant Analysis (PLS-DA) and machine learning techniques such as XGBoost were employed, and their discriminating power was compared. PLS-DA method provided more robust classifications than XGBoost method, showing that the adopted approach might be an interesting tool for forensic experts to infer BGA information from the DNA profile of unknown individuals, but also highlighting that the commercial forensic panels could be inadequate to discriminate populations at intra-continental level.

Inference of recent admixture using genotype data

The inference of biogeographic ancestry (BGA) has become a focus of forensic genetics. Misinference of BGA can have profound unwanted consequences for investigations and society. We show that recent admixture can lead to misclassification and erroneous inference of ancestry proportions, using state of the art analysis tools with (i) simulations, (ii) 1000 genomes project data, and (iii) two individuals analyzed using the ForenSeq DNA Signature Prep Kit. Subsequently, we extend existing tools for estimation of individual ancestry (IA) by allowing for different IA in both parents, leading to estimates of parental individual ancestry (PIA), and a statistical test for recent admixture. Estimation of PIA outperforms IA in most scenarios of recent admixture. Furthermore, additional information about parental ancestry can be acquired with PIA that may guide casework.

Pushing the Boundaries: Forensic DNA Phenotyping Challenged by Single-Cell Sequencing

Single-cell sequencing is a fast developing and very promising field; however, it is not commonly used in forensics. The main motivation behind introducing this technology into forensics is to improve mixture deconvolution, especially when a trace consists of the same cell type. Successful studies demonstrate the ability to analyze a mixture by separating single cells and obtaining CE-based STR profiles. This indicates a potential use of the method in other forensic investigations, like forensic DNA phenotyping, in which using mixed traces is not fully recommended. For this study, we collected single-source autopsy blood from which the white cells were first stained and later separated with the DEPArray™ N×T System. Groups of 20, 10, and 5 cells, as well as 20 single cells, were collected and submitted for DNA extraction. Libraries were prepared using the Ion AmpliSeq™ PhenoTrivium Panel, which includes both phenotype (HIrisPlex-S: eye, hair, and skin color) and ancestry-associated SNP-markers. Prior to sequencing, half of the single-cell-based libraries were additionally amplified and purified in order to improve the library concentrations. Ancestry and phenotype analysis resulted in nearly full consensus profiles resulting in correct predictions not only for the cells groups but also for the ten re-amplified single-cell libraries. Our results suggest that sequencing of single cells can be a promising tool used to deconvolute mixed traces submitted for forensic DNA phenotyping.

Evaluation of a custom GeneRead™ massively parallel sequencing assay with 210 ancestry informative SNPs using the Ion S5™ and MiSeq platforms

A custom GeneRead DNAseq SNP panel with 210 markers was evaluated using the Ion S5 and MiSeq sequencing platforms. Sensitivity, PCR cycle number, and the use of half volume of reagents for target enrichment and library preparation were tested. Furthermore, genotype concordance between results obtained with the different sequencing platforms and with known profiles generated using other sequencing assays was analysed. The GeneRead DNASeq SNP assay gave reproducible results with an input of 200 pg DNA on both platforms. A total of 204 loci were successfully sequenced. Three loci failed completely in the PCR amplification, and three additional loci displayed frequent locus drop-outs due to low read depth or high heterozygote imbalance. Overall, the read depth across the loci was more well-balanced with the MiSeq, while the heterozygote balance was less variable with the Ion S5. Noise levels were low on both platforms (median< 0.2 %). Two simple criteria for genotyping were applied: A minimum threshold of 45 reads and an acceptable heterozygote balance range of 0.3-3.0. Complete concordance between platforms was observed except for three genotypes in one of the poorly performing loci, rs1470637. This locus had relatively low read depths on both platforms, skewed heterozygote balance, and frequent locus drop-outs. There was also full genotype concordance between the results from the GeneRead assay and known profiles generated with the QIAseq and Ion AmpliSeq assays. The few discordant results were either due to locus drop-outs in the poorly performing loci or allele drop-outs in the QIAseq assay. Profiles with a minimum of 179 SNPs were obtained from four challenging case work samples (blood swabs, bone, or blood from a corpse). Overall, the GeneRead DNASeq assay showed considerable potential and could provide a reliable method for SNP genotyping in cases involving identification of individuals, prediction of phenotypic traits, and ancestry inference.

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

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

Ancestry and phenotype predictions from touch DNA using massively parallel sequencing

Direct PCR can be used to successfully generate full STR profiles from DNA present on the surface of objects. STR profiles are only of use in cases where a potential donor profile is available for comparison, and DNA is of sufficient DNA quality and quantity to generate a reliable profile. Often, no donor information is available and only trace DNA is present on items. As a result, alternative techniques are required to generate genetic data that can provide investigative leads. Massively parallel sequencing (MPS) offers the ability to detect trace levels of DNA and improve DNA analysis success from touched items. Here, we present the first application of direct PCR coupled with MPS to generate forensic intelligence SNP data from latent DNA. The panels assessed are (1) the HIrisplex System that targets 24 SNPs to simultaneously predict hair and eye, and (2) the Precision ID Ancestry Panel that targets 165 autosomal SNPs indicative of biogeographic ancestry. For each panel, we analysed 60 touched samples across five individuals and four substrates (glass slide, fuse, zip-lock bag and wire) using Ion AmpliSeq Library Preparation Kit on the automated Ion Chef System and Ion Torrent PGM. We examine the SNP recovery, concordance with reference samples and the genotype reproducibility from different substrates and donors. The results demonstrate the application of this approach for obtaining informative genetic from trace amounts of DNA.

Forensic evaluation of the Asia Pacific ancestry-informative MAPlex assay

DNA intelligence, and particularly the inference of biogeographical ancestry (BGA) is increasing in interest, and relevance within the forensic genetics community. The majority of current MPS-based forensic ancestry-informative assays focus on the differentiation of major global populations. The recently published MAPlex (Multiplex for the Asia Pacific) panel contains 144 SNPs and 20 microhaplotypes and aims to improve the differentiation of populations in the Asia Pacific region. This study reports the first forensic evaluation of the MAPlex panel using AmpliSeq technology and Ion S5 sequencing. This study reports on the overall performance of MAPlex including the assay's sequence coverage distribution and stability, baseline noise and description of problematic SNPs. Dilution series, artificially degraded and mixed DNA samples were also analysed to evaluate the sensitivity of the panel with challenging or compromised forensic samples. As the first panel to combine biallelic SNPs, multiple-allele SNPs and microhaplotypes, the MAPlex assay demonstrated an enhanced capacity for mixture detection, not easily performed with common binary SNPs. This performance evaluation indicates that MAPlex is a robust, stable and highly sensitive assay that is applicable to forensic casework for the prediction of BGA.

Massive parallel sequencing in forensics: advantages, issues, technicalities, and prospects

In the last decade, next-generation sequencing (NGS) technology, alternatively massive parallel sequencing (MPS), was applied to all fields of biological research. Its introduction to the field of forensics was slower, mainly due to lack of accredited sequencers, kits, and relatively higher sequencing error rates as compared with standardized Sanger sequencing. Currently, a majority of the problematic issues have been solved, which is proven by the body of reports in the literature. Here, we discuss the utility of NGS sequencing in forensics, emphasizing the advantages, issues, the technical aspects of the experiments, commercial solutions, and the potentially interesting applications of MPS.

Evaluation of the VISAGE Basic Tool for Appearance and Ancestry Prediction Using PowerSeq Chemistry on the MiSeq FGx System

The study of DNA to predict externally visible characteristics (EVCs) and the biogeographical ancestry (BGA) from unknown samples is gaining relevance in forensic genetics. Technical developments in Massively Parallel Sequencing (MPS) enable the simultaneous analysis of hundreds of DNA markers, which improves successful Forensic DNA Phenotyping (FDP). The EU-funded VISAGE (VISible Attributes through GEnomics) Consortium has developed various targeted MPS-based lab tools to apply FDP in routine forensic analyses. Here, we present an evaluation of the VISAGE Basic tool for appearance and ancestry prediction based on PowerSeq chemistry (Promega) on a MiSeq FGx System (Illumina). The panel consists of 153 single nucleotide polymorphisms (SNPs) that provide information about EVCs (41 SNPs for eye, hair and skin color from HIrisPlex-S) and continental BGA (115 SNPs; three overlap with the EVCs SNP set). The assay was evaluated for sensitivity, repeatability and genotyping concordance, as well as its performance with casework-type samples. This targeted MPS assay provided complete genotypes at all 153 SNPs down to 125 pg of input DNA and 99.67% correct genotypes at 50 pg. It was robust in terms of repeatability and concordance and provided useful results with casework-type samples. The results suggest that this MPS assay is a useful tool for basic appearance and ancestry prediction in forensic genetics for users interested in applying PowerSeq chemistry and MiSeq for this purpose.

The Use of Forensic DNA Phenotyping in Predicting Appearance and Biogeographic Ancestry

BACKGROUND

Persons whose identifying DNA profile (STR profile) is not yet known to the ingvestigating authorities cannot be identified by standard forensic DNA analysis (STR profiling) as it is now practiced. In view of the current public debate, particularly in Germany, on the legalization of so-called forensic DNA phenotyping, we present its scientific basis, societal aspects, and forensic applications and describe the analytic techniques that are now available.

METHODS

This review is based on pertinent publications that were retrieved by a selective search in PubMed and in public media, and on the authors' own research.

RESULTS

Forensically validated DNA test systems are available for the categorization of eye, hair, and skin color and the inference of continental biogeographic ancestry. As for statistical measures of test accuracy, the AUC (area under the curve) values lie in the range 0.74-0.99 for eye color, 0.64-0.94 for hair color, and 0.72-0.99 for skin color, depending on the predictive model and color category used.The corre- sponding positive predictive values (PPV) are lower. Empirical social-scientific research on forensic DNA phenotyping has shown that preserving privacy and protecting against discrimination are major ethical and regulatory considerations.

CONCLUSION

All three methods of forensic DNA phenotyping-the predition of exter- nally visible characteristics, biogeographic ancestry, and the estimation of age from crime scene DNA-require a proper regulatory framework and should be used in conjunction with each other. Before forensic DNA phenotyping can be implemented in forensic practice, steps must be taken to minimize the risks of violation of privacy scrimination and to ensure that these methods are used transpar- ently and proportionately.

Differentiation of Hispanic biogeographic ancestry with 80 ancestry informative markers

Ancestry informative single nucleotide polymorphisms (SNPs) can identify biogeographic ancestry (BGA); however, population substructure and relatively recent admixture can make differentiation difficult in heterogeneous Hispanic populations. Utilizing unrelated individuals from the Genomic Origins and Admixture in Latinos dataset (GOAL, n = 160), we designed an 80 SNP panel (Setser80) that accurately depicts BGA through STRUCTURE and PCA. We compared our Setser80 to the Seldin and Kidd panels via resampling simulations, which models data based on allele frequencies. We incorporated Admixed American 1000 Genomes populations (1000 G, n = 347), into a combined populations dataset to determine robustness. Using multinomial logistic regression (MLR), we compared the 3 panels on the combined dataset and found overall MLR classification accuracies: 93.2% Setser80, 87.9% Seldin panel, 71.4% Kidd panel. Naïve Bayesian classification had similar results on the combined dataset: 91.5% Setser80, 84.7% Seldin panel, 71.1% Kidd panel. Although Peru and Mexico were absent from panel design, we achieved high classification accuracy on the combined populations for Peru (MLR = 100%, naïve Bayes = 98%), and Mexico (MLR = 90%, naïve Bayes = 83.4%) as evidence of the portability of the Setser80. Our results indicate the Setser80 SNP panel can reliably classify BGA for individuals of presumed Hispanic origin.

Characterization of sequence variation at 30 autosomal STRs in Chinese Han and Tibetan populations

Massively parallel sequencing (MPS) technologies have the ability to reveal sequence variations within STR alleles as well as their nominal allele lengths, which have traditionally been detected by CE instruments. Recently, Thermo Fisher Scientific has updated the MPS-STR panel, named the Precision ID GlobalFiler next-generation sequencing (NGS) STR Panel version 2, with primers redesigned to add two pentanucleotide tandem repeat loci and profile interpretation supported by the Converge software. Using the Ion Chef System, the Ion S5XL System, and the Converge software, genetic variations were characterized within STR repeat and flanking regions of 30 autosomal STR markers in 115 unrelated individuals from two Chinese population groups (58 Tibetans and 57 Hans). Nineteen STRs demonstrated a relative increase in diversity with the variant sequence alleles compared with those of traditional nominal length alleles. In total, 390 alleles were identified by their sequences compared with 258 alleles that were identified by length. Of these 92 sequence variants found within the STR repeat regions, 40 variants were located in STR flanking regions. Additionally, the agreement of the results with CE data was evaluated, as was the ability of this new MPS panel to analyze case-type (11 samples) and artificially degraded samples (seven samples in triplicate). The results generated from this study illustrate that extensive sequence variation exists in commonly used STR markers in the selected population samples and indicate that this NGS STR panel has the potential to be used as an effective tool for human forensics.

The MASTiFF panel-a versatile multiple-allele SNP test for forensics

Forensic identification tests often need recourse to markers that can successfully type highly degraded DNA, and binary single nucleotide polymorphisms (SNPs) have become the variants of choice for such analyses because of their short amplified fragment lengths. The two main drawbacks of SNPs are their reduced power of discrimination per marker compared with mainstream forensic STRs and an inability to robustly detect mixed DNA-particularly using capillary electrophoresis genotyping systems such as SNaPshot™, where the dye signals are much more imbalanced than those of STR profiles. This study compiled a compact set of multiple-allele SNPs consisting of loci that had three or four nucleotide variants at the same site in order to address the lack of mixture detection capability with binary SNP tests, as well as improving levels of polymorphism per SNP by transitioning to a maximum of six or ten genotypes per locus. We report the development and optimisation of a SNaPshot-based forensic test comprising 27 tri-allelic and 2 tetra-allelic SNPs, which we named MASTiFF: a multiple-allele SNP test for forensics. Assessments of the MASTiFF panel's levels of discrimination power in the five main population groups indicate random match probabilities ranging from 10^-15 down to 10^-20-improving the levels possible from an equivalent number of binary SNPs. The SNaPshot test was able to detect simple mixtures successfully with more than two alleles observed in 30% of SNPs. From allele frequency data, it is estimated that more than two alleles will be present in at least one MASTiFF SNP in 99.8% of two-person mixtures, making this panel an ideal supplementary test when SNPs are chosen for the analysis of degraded forensic DNA.

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

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

Inter-laboratory study on standardized MPS libraries: evaluation of performance, concordance, and sensitivity using mixtures and degraded DNA

We present results from an inter-laboratory massively parallel sequencing (MPS) study in the framework of the SeqForSTRs project to evaluate forensically relevant parameters, such as performance, concordance, and sensitivity, using a standardized sequencing library including reference material, mixtures, and ancient DNA samples. The standardized library was prepared using the ForenSeq DNA Signature Prep Kit (primer mix A). The library was shared between eight European laboratories located in Austria, France, Germany, The Netherlands, and Sweden to perform MPS on their particular MiSeq FGx sequencers. Despite variation in performance between sequencing runs, all laboratories obtained quality metrics that fell within the manufacturer’s recommended ranges. Furthermore, differences in locus coverage did not inevitably adversely affect heterozygous balance. Inter-laboratory concordance showed 100% concordant genotypes for the included autosomal and Y-STRs, and still, X-STR concordance exceeded 83%. The exclusive reasons for X-STR discordances were drop-outs at DXS10103. Sensitivity experiments demonstrated that correct allele calling varied between sequencing instruments in particular for lower DNA amounts (≤ 125 pg). The analysis of compromised DNA samples showed the drop-out of one sample (FA10013B01A) while for the remaining three degraded DNA samples MPS was able to successfully type ≥ 87% of all aSTRs, ≥ 78% of all Y-STRs, ≥ 68% of all X-STRs, and ≥ 92% of all iSNPs demonstrating that MPS is a promising tool for human identity testing, which in return, has to undergo rigorous in-house validation before it can be implemented into forensic routine casework.

A single nucleotide polymorphism panel for individual identification and ancestry assignment in Caucasians and four East and Southeast Asian populations using a machine learning classifier

Single nucleotide polymorphism (SNP) profiling is an effective means of individual identification and ancestry inferences in forensic genetics. This study established a SNP panel for the simultaneous individual identification and ancestry assignment of Caucasian and four East and Southeast Asian populations. We analyzed 220 SNPs (125 autosomal, 17 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs) of the DNA samples from 563 unrelated individuals of five populations (89 Caucasian, 234 Taiwanese Han, 90 Filipino, 79 Indonesian and 71 Vietnamese) and 18 degraded DNA samples. Informativeness for assignment (In) was used to select ancestry informative SNPs (AISNPs). A machine learning classifier, support vector machine (SVM), was used for ancestry assignment. Of the 220 SNPs, 62 were individual identification SNPs (IISNPs) (51 autosomal and 11 X-chromosomal SNPs) and 191 were AISNPs (100 autosomal, 13 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs). The 51 autosomal IISNPs offered cumulative random match probabilities (cRMPs) ranging from 1.56 × 10^-21 to 3.16 × 10^-22 among these five populations. Using AISNPs with the SVM, the overall accuracy rate of ancestry inference achieved in the testing dataset between Caucasian, Taiwanese Han, and Filipino populations was 88.9%, whereas it was 70.0% between Caucasians and each of the four East and Southeast Asian populations. For the 18 degraded DNA samples with incomplete profiling, the accuracy rate of ancestry assignment was 94.4%. We have developed a 220-SNP panel for simultaneous individual identification and ethnic origin differentiation between Caucasian and the four East and Southeast Asian populations. This SNP panel may assist with DNA analysis of forensic casework.

Age Estimation with DNA: From Forensic DNA Fingerprinting to Forensic (Epi)Genomics: A Mini-Review

Forensic genetics developed from protein-based techniques a quarter of a century ago and became famous as “DNA fingerprinting,” this being based on restriction fragment length polymorphisms (RFLPs) of high-molecular-weight DNA. The amplification of much smaller short tandem repeat (STR) sequences using the polymerase chain reaction soon replaced RFLP analysis and advanced to become the gold standard in genetic identification. Meanwhile, STR multiplexes have been developed and made commercially available which simultaneously amplify up to 30 STR loci from as little as 15 cells or fewer. The enormous information content that comes with the large variety of observed STR genotypes allows for genetic individualisation (with the exception of identical twins). Carefully selected core STR loci form the basis of intelligence-led DNA databases that provide investigative leads by linking unsolved crime scenes and criminals through their matched STR profiles. Nevertheless, the success of modern DNA fingerprinting depends on the availability of reference material from suspects. In order to provide new investigative leads in cases where such reference samples are absent, forensic scientists started to explore the prediction of phenotypic traits from the DNA of the evidentiary sample. This paradigm change now uses DNA and epigenetic markers to forecast characteristics that are useful to triage further investigative work. So far, the best investigated externally visible characteristics are eye, hair and skin colour, as well as geographic ancestry and age. Information on the chronological age of a stain donor (or any sample donor) is elemental for forensic investigations in a number of aspects and has, therefore, been explored by researchers in some detail. Among different methodological approaches tested to date, the methylation-sensitive analysis of carefully selected DNA markers (CpG sites) has brought the most promising results by providing prediction accuracies of ±3–4 years, which can be comparable to, or even surpass those from, eyewitness reports. This mini-review puts recent developments in age estimation via (epi)genetic methods in the context of the requirements and goals of forensic genetics and highlights paths to follow in the future of forensic genomics.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

A GHEP-ISFG collaborative study on the genetic variation of 38 autosomal indels for human identification in different continental populations

A collaborative effort was carried out by the Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) to promote knowledge exchange between associate laboratories interested in the implementation of indel-based methodologies and build allele frequency databases of 38 indels for forensic applications. These databases include populations from different countries that are relevant for identification and kinship investigations undertaken by the participating laboratories. Before compiling population data, participants were asked to type the 38 indels in blind samples from annual GHEP-ISFG proficiency tests, using an amplification protocol previously described. Only laboratories that reported correct results contributed with population data to this study. A total of 5839 samples were genotyped from 45 different populations from Africa, America, East Asia, Europe and Middle East. Population differentiation analysis showed significant differences between most populations studied from Africa and America, as well as between two Asian populations from China and East Timor. Low F_ST values were detected among most European populations. Overall diversities and parameters of forensic efficiency were high in populations from all continents.

Population and performance analyses of four major populations with Illumina's FGx Forensic Genomics System

The MiSeq FGx Forensic Genomics System (Illumina) enables amplification and massively parallel sequencing of 59 STRs, 94 identity informative SNPs, 54 ancestry informative SNPs, and 24 phenotypic informative SNPs. Allele frequency and population statistics data were generated for the 172 SNP loci included in this panel on four major population groups (Chinese, African Americans, US Caucasians, and Southwest Hispanics). Single-locus and combined random match probability values were generated for the identity informative SNPs. The average combined STR and identity informative SNP random match probabilities (assuming independence) across all four populations were 1.75E-67 and 2.30E-71 with length-based and sequence-based STR alleles, respectively. Ancestry and phenotype predictions were obtained using the ForenSeq™ Universal Analysis System (UAS; Illumina) based on the ancestry informative and phenotype informative SNP profiles generated for each sample. Additionally, performance metrics, including profile completeness, read depth, relative locus performance, and allele coverage ratios, were evaluated and detailed for the 725 samples included in this study. While some genetic markers included in this panel performed notably better than others, performance across populations was generally consistent. The performance and population data included in this study support that accurate and reliable profiles were generated and provide valuable background information for laboratories considering internal validation studies and implementation.

STRait Razor v2s: Advancing sequence-based STR allele reporting and beyond to other marker systems

STRait Razor has provided the forensic community a free-to-use, open-source tool for short tandem repeat (STR) analysis of massively parallel sequencing (MPS) data. STRait Razor v2s (SRv2s) allows users to capture physically phased haplotypes within the full amplicon of both commercial (ForenSeq) and "early access" panels (PowerSeq, Mixture ID). STRait Razor v2s may be run in batch mode to facilitate population-level analysis and is supported by all Unix distributions (including MAC OS). Data are reported in tables in string (haplotype), length-based (e.g., vWA allele 14), and International Society of Forensic Genetics (ISFG)-recommended (vWA [CE 14]-GRCh38-chr12:5983950-5984049 (TAGA)₁₀ (CAGA)₃ TAGA) formats. STRait Razor v2s currently contains a database of ∼2500 unique sequences. This database is used by SRv2s to match strings to the appropriate allele in ISFG-recommended format. In addition to STRs, SRv2s has configuration files necessary to capture and report haplotypes from all marker types included in these multiplexes (e.g., SNPs, InDels, and microhaplotypes). To facilitate mixture interpretation, data may be displayed from all markers in a format similar to that of electropherograms displayed by traditional forensic software. The download package for SRv2s may be found at https://www.unthsc.edu/graduate-school-of-biomedical-sciences/molecular-and-medical-genetics/laboratory-faculty-and-staff/strait-razor.

Massively parallel sequencing of 231 autosomal SNPs with a custom panel: a SNP typing assay developed for human identification with Ion Torrent PGM

The custom-designed single nucleotide polymorphism (SNP) panel amplified 231 autosomal SNPs in one PCR reaction and subsequently sequenced with massively parallel sequencing (MPS) technology and Ion Torrent personal genome machine (PGM). SNPs were chosen from SNPforID, IISNP, HapMap, dbSNP, and related published literatures. Full concordance was obtained between available MPS calling and Sanger sequencing with 9947A and 9948 controls. Ten SNPs (rs4606077, rs334355, rs430046, rs2920816, rs4530059, rs1478829, rs1498553, rs7141285, rs12714757 and rs2189011) with low coverage or heterozygote imbalance should be optimized or excluded from the panel. Sequence data had sufficiently high coverage and gave reliable SNP calling for the remaining 221 loci with the custom MPS–SNP panel. A default DNA input amount of 10 ng per reaction was recommended by Ampliseq technology but sensitivity testing revealed positive results from as little as 1 ng input DNA. Mixture testing with this panel is possible through analysis of the F_MAR (frequency of major allele reads) values at most loci with enough high coverage depth and low level of sequencing noise. These results indicate the potential advantage of the custom MPS–SNP assays and Ion Torrent PGM platform for forensic study.