Evaluation of DNA Extraction Methods Developed for Forensic and Ancient DNA Applications Using Bone Samples of Different Age.

The efficient extraction of DNA from challenging samples, such as bones, is critical for the success of downstream genotyping analysis in molecular genetic disciplines. Even though the ancient DNA community has developed several protocols targeting small DNA fragments that are typically present in decomposed or old specimens, only recently forensic geneticists have started to adopt those protocols. Here, we compare an ancient DNA extraction protocol (Dabney) with a bone extraction method (Loreille) typically used in forensics. Real-time quantitative PCR and forensically representative typing methods including fragment size analysis and sequencing were used to assess protocol performance. We used four bone samples of different age in replicates to study the effects of both extraction methods. Our results confirm Loreille's overall increased gain of DNA when enough tissue is available and Dabney's improved efficiency for retrieving shorter DNA fragments that is beneficial when highly degraded DNA is present. The results suggest that the choice of extraction method needs to be based on available sample, degradation state, and targeted genotyping method. We modified the Dabney protocol by pooling parallel lysates prior to purification to study gain and performance in single tube typing assays and found that up to six parallel lysates lead to an almost linear gain of extracted DNA. These data are promising for further forensic investigations as the adapted Dabney protocol combines increased sensitivity for degraded DNA with necessary total DNA amount for forensic applications.

Performance of ancestry-informative SNP and microhaplotype markers.

The use of microhaplotypes (MHs) for ancestry inference has added to an increasing number of ancestry-informative markers (AIMs) for forensic application that includes autosomal single nucleotide polymorphisms (SNPs) and insertions/deletions (indels). This study compares bi-allelic and tri-allelic SNPs as well as MH markers for their ability to differentiate African, European, South Asian, East Asian, and American population groups from the 1000 Genomes Phase 3 database. A range of well-established metrics were applied to rank each marker according to the population differentiation potential they measured. These comprised: absolute allele frequency differences (δ); Rosenberg's informativeness for (ancestry) assignment (In); the fixation index (FST); and the effective number of alleles (Ae). A panel consisting of all three marker types resulted in the lowest mean divergence per population per individual (MDPI = 2.16%) when selected by In. However, when marker types were not mixed, MHs were the highest performing markers by most metrics (MDPI < 4%) for differentiation between the five continental populations.

SD quants-Sensitive detection tetraplex-system for nuclear and mitochondrial DNA quantification and degradation inference.

Measuring the quantity of DNA present in a forensic sample is relevant in a number of ways. First, it informs the analyst about the general DNA content to adjust the volume of DNA extract used for the genotyping assay to the optimal conditions (when possible). Second, quantification values can serve as plausibility checks for the performance of the DNA extraction method used as extraction positive and negative controls demand expected values. Third and relevant to highly compromised specimens, DNA quantification can inform about the degradation state of the DNA extracted from the unknown biological sample and aid the choice of downstream genotyping assays. While there are different, commercial products for the quantification of nuclear DNA available, commercial mitochondrial DNA (mtDNA) quantification systems are rare. Even more so, the simultaneous quantification of nuclear and mtDNA that is of relevance in highly degraded forensic specimens has rarely been described. We present here a novel real-time qPCR based tetraplex system termed SD quants that targets two different-sized mtDNA and a nuclear DNA region and includes an internal positive control to monitor potential inhibition. SD quants was compared to other existing quantification systems and subjected to analysis of severely degraded DNA present in ancient DNA and aged forensic specimens. This study complies with the MIQE (Bustin et al., 2009) guidelines (when applicable).

Resolving the matrilineal relationship of seven Late Bronze Age individuals from Stillfried, Austria.

In 1976 human remains of seven individuals were discovered in a storage pit located within the Late Bronze Age (9th century B.C.) settlement Stillfried an der March, Austria. In contrast to the common funeral rite of cremation typical for the Urnfield culture (1300-800 B.C.) the individuals' skeletal remains were found outstandingly preserved (Figure S1). As a result, the burial was subject to various investigations, including two conflicting genealogical pedigree reconstructions, one of which was favoured by later geological fingerprinting. We performed mitochondrial (mt)DNA testing in order to genetically characterize the remains and shed light into the matrilineal relationship of the seven individuals that were earlier anthropologically identified as three adults (two women and a man) and four subadults (one female and three males). MtDNA was analysed using Primer Extension Capture and Massively Parallel Sequencing. The results were by and large in conflict with both pedigree models but confirmed some of the details that were elaborated in previous studies. Whereas both pedigree models suggested that all children were related to one or both females, mtDNA analyses revealed that only one subadult male resulted in the same mitotype as one adult female. All other children yielded different mitotypes indicating that they were maternally unrelated to the two females and between each other.

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.

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.

Open source software EuroForMix can be used to analyse complex SNP mixtures.

A series of two- and three-person mixtures of varying dilutions were prepared and analysed with Life Technologies' HID-Ion AmpliSeq™ Identity Panel v2.2 using the Ion PGM™ massively parallel sequencing (MPS) system. From this panel we used 134 autosomal SNPs. Using the reference samples of three donors, we evaluated the strength of evidence with likelihood ratio (LR) calculations using the open-source quantitative EuroForMix program and compared the results with a previous study using a qualitative software (LRmix). SNP analysis is a special case of STRs, restricted to a maximum of two alleles per locus. We showed that simple two-person mixtures can be readily analysed with both LRmix and Euroformix, but the performance of three- or more person mixtures is generally inefficient with LRmix. Taking account of the "peak height" information, by substituting 'sequence read' coverage values from the MPS data for each SNP allele, greatly improves the discrimination between true and non-contributors. The higher the mixture proportion (Mx) of the person of interest is, the higher the LR. Simulation experiments (up to six contributors) showed that the strength of the evidence is dependent upon Mx, but relatively insensitive to the number of contributors. If a higher number of loci were multiplexed, the analysis of mixtures would be much improved, because the extra information would enable lower Mx values to be evaluated. In summary, incorporating the 'sequence read' (coverage) into the quantitative model shows a significant benefit over the qualitative approach. Calculations are quite fast (six seconds for three contributors).

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.

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.

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.

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.