Mayra Eduardoff

Inter-laboratory evaluation of SNP-based forensic identification by massively parallel sequencing using the Ion PGM™.

Next generation sequencing (NGS) offers the opportunity to analyse forensic DNA samples and obtain massively parallel coverage of targeted short sequences with the variants they carry. We evaluated the levels of sequence coverage, genotyping precision, sensitivity and mixed DNA patterns of a prototype version of the first commercial forensic NGS kit: the HID-Ion AmpliSeq™ Identity Panel with 169-markers designed for the Ion PGM™ system. Evaluations were made between three laboratories following closely matched Ion PGM™ protocols and a simple validation framework of shared DNA controls. The sequence coverage obtained was extensive for the bulk of SNPs targeted by the HID-Ion AmpliSeq™ Identity Panel. Sensitivity studies showed 90-95% of SNP genotypes could be obtained from 25 to 100pg of input DNA. Genotyping concordance tests included Coriell cell-line control DNA analyses checked against whole-genome sequencing data from 1000 Genomes and Complete Genomics, indicating a very high concordance rate of 99.8%. Discordant genotypes detected in rs1979255, rs1004357, rs938283, rs2032597 and rs2399332 indicate these loci should be excluded from the panel. Therefore, the HID-Ion AmpliSeq™ Identity Panel and Ion PGM™ system provide a sensitive and accurate forensic SNP genotyping assay. However, low-level DNA produced much more varied sequence coverage and in forensic use the Ion PGM™ system will require careful calibration of the total samples loaded per chip to preserve the genotyping reliability seen in routine forensic DNA. Furthermore, assessments of mixed DNA indicate the users control of sequence analysis parameter settings is necessary to ensure mixtures are detected robustly. Given the sensitivity of Ion PGM™, this aspect of forensic genotyping requires further optimisation before massively parallel sequencing is applied to routine casework.

Building a forensic ancestry panel from the ground up: The EUROFORGEN Global AIM-SNP set

Emerging next-generation sequencing technologies will enable DNA analyses to add pigmentation predictive and ancestry informative (AIM) SNPs to the range of markers detectable from a single PCR test. This prompted us to re-appraise current forensic and genomics AIM-SNPs and from the best sets, to identify the most divergent markers for a five population group differentiation of Africans, Europeans, East Asians, Native Americans and Oceanians by using our own online genome variation browsers. We prioritized careful balancing of population differentiation across the five group comparisons in order to minimize bias when estimating co-ancestry proportions in individuals with admixed ancestries. The differentiation of European from Middle East or South Asian ancestries was not chosen as a characteristic in order to concentrate on introducing Oceanian differentiation for the first time in a forensic AIM set. We describe a complete set of 128 AIM-SNPs that have near identical population-specific divergence across five continentally defined population groups. The full set can be systematically reduced in size, while preserving the most informative markers and the balance of population-specific divergence in at least four groups. We describe subsets of 88, 55, 28, 20 and 12 AIMs, enabling both new and existing SNP genotyping technologies to exploit the best markers identified for forensic ancestry analysis.

Massively parallel sequencing of complete mitochondrial genomes from hair shaft samples

Though shed hairs are one of the most commonly encountered evidence types, they are among the most limited in terms of DNA quantity and quality. As a result, DNA testing has historically focused on the recovery of just about 600 base pairs of the mitochondrial DNA control region. Here, we describe our success in recovering complete mitochondrial genome (mtGenome) data (∼16,569bp) from single shed hairs. By employing massively parallel sequencing (MPS), we demonstrate that particular hair samples yield DNA sufficient in quantity and quality to produce 2-3kb mtGenome amplicons and that entire mtGenome data can be recovered from hair extracts even without PCR enrichment. Most importantly, we describe a small amplicon multiplex assay comprised of sixty-two primer sets that can be routinely applied to the compromised hair samples typically encountered in forensic casework. In all samples tested here, the MPS data recovered using any one of the three methods were consistent with the control Sanger sequence data developed from high quality known specimens. Given the recently demonstrated value of complete mtGenome data in terms of discrimination power among randomly sampled individuals, the possibility of recovering mtGenome data from the most compromised and limited evidentiary material is likely to vastly increase the utility of mtDNA testing for hair evidence.

Evaluation of DNA Variants Associated with Androgenetic Alopecia and Their Potential to Predict Male Pattern Baldness

Androgenetic alopecia, known in men as male pattern baldness (MPB), is a very conspicuous condition that is particularly frequent among European men and thus contributes markedly to variation in physical appearance traits amongst Europeans. Recent studies have revealed multiple genes and polymorphisms to be associated with susceptibility to MPB. In this study, 50 candidate SNPs for androgenetic alopecia were analyzed in order to verify their potential to predict MPB. Significant associations were confirmed for 29 SNPs from chromosomes X, 1, 5, 7, 18 and 20. A simple 5-SNP prediction model and an extended 20-SNP model were developed based on a discovery panel of 305 males from various European populations fitting one of two distinct phenotype categories. The first category consisted of men below 50 years of age with significant baldness and the second; men aged 50 years or older lacking baldness. The simple model comprised the five best predictors: rs5919324 near AR, rs1998076 in the 20p11 region, rs929626 in EBF1, rs12565727 in TARDBP and rs756853 in HDAC9. The extended prediction model added 15 SNPs from five genomic regions that improved overall prevalence-adjusted predictive accuracy measured by area under the receiver characteristic operating curve (AUC). Both models were evaluated for predictive accuracy using a test set of 300 males reflecting the general European population. Applying a 65% probability threshold, high prediction sensitivity of 87.1% but low specificity of 42.4% was obtained in men aged <50 years. In men aged ≥50, prediction sensitivity was slightly lower at 67.7% while specificity reached 90%. Overall, the AUC=0.761 calculated for men at or above 50 years of age indicates these SNPs offer considerable potential for the application of genetic tests to predict MPB patterns, adding a highly informative predictive system to the emerging field of forensic analysis of externally visible characteristics.

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

The EUROFORGEN Global ancestry-informative SNP (AIM-SNPs) panel is a forensic multiplex of 128 markers designed to differentiate an individuals ancestry from amongst the five continental population groups of Africa, Europe, East Asia, Native America, and Oceania. A custom multiplex of AmpliSeq™ PCR primers was designed for the Global AIM-SNPs to perform massively parallel sequencing using the Ion PGM™ system. This study assessed individual SNP genotyping precision using the Ion PGM™, the forensic sensitivity of the multiplex using dilution series, degraded DNA plus simple mixtures, and the ancestry differentiation power of the final panel design, which required substitution of three original ancestry-informative SNPs with alternatives. Fourteen populations that had not been previously analyzed were genotyped using the custom multiplex and these studies allowed assessment of genotyping performance by comparison of data across five laboratories. Results indicate a low level of genotyping error can still occur from sequence misalignment caused by homopolymeric tracts close to the target SNP, despite careful scrutiny of candidate SNPs at the design stage. Such sequence misalignment required the exclusion of component SNP rs2080161 from the Global AIM-SNPs panel. However, the overall genotyping precision and sensitivity of this custom multiplex indicates the Ion PGM™ assay for the Global AIM-SNPs is highly suitable for forensic ancestry analysis with massively parallel sequencing.

Evaluation of the predictive capacity of DNA variants associated with straight hair in Europeans.

DNA-based prediction of hair morphology, defined as straight, curly or wavy hair, could contribute to an improved description of an unknown offender and allow more accurate forensic reconstructions of physical appearance in the field of forensic DNA phenotyping. Differences in scalp hair morphology are significant at the worldwide scale and within Europe. The only genome-wide association study made to date revealed the Trichohyalin gene (TCHH) to be significantly associated with hair morphology in Europeans and reported weaker associations for WNT10A and FRAS1 genes. We conducted a study that centered on six SNPs located in these three genes with a sample of 528 individuals from Poland. The predictive capacity of the candidate DNA variants was evaluated using logistic regression; classification and regression trees; and neural networks, by applying a 10-fold cross validation procedure. Additionally, an independent test set of 142 males from six European populations was used to verify performance of the developed prediction models. Our study confirmed association of rs11803731 (TCHH), rs7349332 (WNT10A) and rs1268789 (FRAS1) SNPs with hair morphology. The combined genotype risk score for straight hair had an odds ratio of 2.7 and these predictors explained ∼ 8.2% of the total variance. The selected three SNPs were found to predict straight hair with a high sensitivity but low specificity when a 10-fold cross validation procedure was applied and the best results were obtained using the neural networks approach (AUC=0.688, sensitivity=91.2%, specificity=23.0%). Application of the neural networks model with 65% probability threshold on an additional test set gave high sensitivity (81.4%) and improved specificity (50.0%) with a total of 78.7% correct calls, but a high non-classification rate (66.9%). The combined TTGGGG SNP genotype for rs11803731, rs7349332, rs1268789 (European frequency=4.5%) of all six straight hair-associated alleles was identified as the best predictor, giving >80% probability of straight hair. Finally, association testing of 44 SNPs previously identified to be associated with male pattern baldness revealed a suggestive association with hair morphology for rs4679955 on 3q25.1. The study results reported provide the starting point for the development of a predictive test for hair morphology in Europeans. More studies are now needed to discover additional determinants of hair morphology to improve the predictive accuracy of this trait in forensic analysis.

Optimized mtDNA Control Region Primer Extension Capture Analysis for Forensically Relevant Samples and Highly Compromised mtDNA of Different Age and Origin

The analysis of mitochondrial DNA (mtDNA) has proven useful in forensic genetics and ancient DNA (aDNA) studies, where specimens are often highly compromised and DNA quality and quantity are low. In forensic genetics, the mtDNA control region (CR) is commonly sequenced using established Sanger-type Sequencing (STS) protocols involving fragment sizes down to approximately 150 base pairs (bp). Recent developments include Massively Parallel Sequencing (MPS) of (multiplex) PCR-generated libraries using the same amplicon sizes. Molecular genetic studies on archaeological remains that harbor more degraded aDNA have pioneered alternative approaches to target mtDNA, such as capture hybridization and primer extension capture (PEC) methods followed by MPS. These assays target smaller mtDNA fragment sizes (down to 50 bp or less), and have proven to be substantially more successful in obtaining useful mtDNA sequences from these samples compared to electrophoretic methods. Here, we present the modification and optimization of a PEC method, earlier developed for sequencing the Neanderthal mitochondrial genome, with forensic applications in mind. Our approach was designed for a more sensitive enrichment of the mtDNA CR in a single tube assay and short laboratory turnaround times, thus complying with forensic practices. We characterized the method using sheared, high quantity mtDNA (six samples), and tested challenging forensic samples (n = 2) as well as compromised solid tissue samples (n = 15) up to 8 kyrs of age. The PEC MPS method produced reliable and plausible mtDNA haplotypes that were useful in the forensic context. It yielded plausible data in samples that did not provide results with STS and other MPS techniques. We addressed the issue of contamination by including four generations of negative controls, and discuss the results in the forensic context. We finally offer perspectives for future research to enable the validation and accreditation of the PEC MPS method for final implementation in forensic genetic laboratories.

Evaluation of the precision ID whole MtDNA genome panel for forensic analyses

Mitochondrial DNA (mtDNA) amplification and Massively Parallel Sequencing (MPS) using an early access version of the Precision ID Whole MtDNA Genome Panel (Thermo Fisher Scientific) and the Ion Personal Genome Machine (PGM) were evaluated using 15 forensically relevant samples. Samples were selected to represent typical forensic specimens for mtDNA analysis including hairs, hair shafts, swabs and ancient solid tissue samples (bones and teeth) that were stored in the freezer for up to several years after having been typed with conventional Sanger-type Sequencing and Capillary Electrophoresis. The MPS haplotypes confirmed the earlier results in all samples and provided additional sequence information that improved discrimination power and haplogroup estimation. The results raised the appetite for further experiments to validate and apply the new technology in forensic practice.

Mass spectrometric base composition profiling: Implications for forensic mtDNA databasing

In forensic genetics mitochondrial DNA (mtDNA) is usually analyzed by direct Sanger-type sequencing (STS). This method is known to be laborious and sometimes prone to human error. Alternative methods have been proposed that lead to faster results. Among these are methods that involve mass-spectrometry resulting in base composition profiles that are, by definition, less informative than the full nucleotide sequence. Here, we applied a highly automated electrospray ionization mass spectrometry (ESI-MS) system (PLEX-ID) to an mtDNA population study to compare its performance with respect to throughput and concordance to STS. We found that the loss of information power was relatively low compared to the gain in speed and analytical standardization. The detection of point and length heteroplasmy turned out to be roughly comparable between the technologies with some individual differences related to the processes. We confirm that ESI-MS provides a valuable platform for analyzing mtDNA variation that can also be applied in the forensic context.