Comparison of bone demineralisation procedures for DNA recovery from burned remains

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

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

A comparison of five DNA extraction methods from degraded human skeletal remains

Extracting DNA from degraded human remains poses a challenge for any forensic genetics laboratory, as it requires efficient high-throughput methods. While little research has compared different techniques, silica in suspension has been identified in the literature as the best method for recovering small fragments, which are often present in these types of samples. In this study, we tested five DNA extraction protocols on 25 different degraded skeletal remains. Including the humerus, ulna, tibia, femur, and petrous bone. The five protocols were organic extraction by phenol/chloroform/isoamyl alcohol, silica in suspension, High Pure Nucleic Acid Large Volume silica columns (Roche), InnoXtract™ Bone (InnoGenomics), and PrepFiler™ BTA with AutoMate™ Express robot (ThermoFisher). We analysed five DNA quantification parameters (small human target quantity, large human target quantity, human male target quantity, degradation index, and internal PCR control threshold), and five DNA profile parameters (number of alleles with peak height higher than analytic and stochastic threshold, average relative fluorescence units (RFU), heterozygous balance, and number of reportable loci) were analysed. Our results suggest that organic extraction by phenol/chloroform/isoamyl alcohol was the best performing method in terms of both quantification and DNA profile results. However, Roche silica columns were found to be the most efficient method.

Development and validation of a SYBR green-based mitochondrial DNA quantification method by following the MIQE and other guidelines

In forensic mitochondrial DNA (mtDNA) analysis, quantitative PCR (qPCR) is usually performed to obtain high-quality sequence data for subsequent Sanger or massively parallel sequencing. Unlike methods for nuclear DNA quantification using qPCR, a calibrator is necessary to obtain mtDNA concentrations (i.e., copies/µL). Herein, we developed and validated a mtDNA quantification method based on a SYBR Green assay by following MIQE [Bustin et al., Clin. Chem. 55 (2009) 611-22] and other guidelines. Primers were designed to amplify nucleotide positions 16,190-16,420 in hypervariable region 1 for qPCR using PowerUp SYBR Green and QuantStudio 5. The optimized conditions were 0.3 µM each primer and an annealing temperature of 60 °C under a 2-step cycling protocol. K562 DNA at 100 pg/µL was converted into a mtDNA concentration of 16,400 copies/µL using linearized plasmid DNA. This mtDNA calibrator was obtained by cloning the synthesized DNA fragments of mtDNA (positions 16,140-16,470) containing a 100-bp inversion. The linear dynamic range of the K562 standard curve was 10,000-0.1 pg/µL (r² ≥ 0.999). The accuracy was examined using NIST SRM 2372a, and its components A, B, and C were quantified with differences of -29.4%, -35.0%, and -22.0%, respectively, against the mtDNA concentrations calculated from published NIST data. We also examined the specificity of the primers, stability of the reaction mix, precision, tolerance against PCR inhibitors, and cross-reactivity against DNA from various animal taxa. Our newly developed mtDNA quantification method is expected to be useful for forensic mtDNA analysis.