Assessment of impact of DNA extraction methods on analysis of human remain samples on massively parallel sequencing success

Assessment of DNA quality for whole genome library preparation

In recent years, more sophisticated DNA technologies for genotyping have enabled considerable progress in various fields such as clinical genetics, archaeogenetics and forensic genetics. DNA samples previously rejected as too challenging to analyze due to low amounts of degraded DNA can now provide useful information. To increase the chances of success with the new methodologies, it is crucial to know the fragment size of the template DNA molecules, and whether the DNA in a sample is mostly single or double stranded. With this knowledge, an appropriate library preparation method can be chosen, and the DNA shearing parameters of the protocol can be adjusted to the DNA fragment size in the sample. In this study, we first developed and evaluated a user-friendly fluorometry-based protocol for estimation of DNA strandedness. We also evaluated different capillary electrophoresis methods for estimation of DNA fragmentation levels. Next, we applied the developed methodologies to a broad variety of DNA samples processed with different DNA extraction protocols. Our findings show that both the applied DNA extraction method and the sample type affect the DNA strandedness and fragmentation. The established protocols and the gained knowledge will be applicable for future sequencing-based high-density SNP genotyping in various fields.

Developmental validation of the ForenSeq® Kintelligence kit, MiSeq FGx® sequencing system and ForenSeq Universal Analysis Software

Forensic Investigative Genetic Genealogy, a recent sub discipline of forensic genomics, leverages the high throughput and sensitivity of detection of next generation sequencing and established genetic and genealogical approaches to support the identification of human remains from missing persons investigations and investigative lead generation in violent crimes. To facilitate forensic DNA evidence analysis, the ForenSeq® Kintelligence multiplex, consisting of 10,230 SNPs, was developed. Design of the ForenSeq Kintelligence Kit, the MiSeq FGx® Sequencing System and the ForenSeq Universal Analysis Software is described. Developmental validation in accordance with SWGDAM guidelines and forensic quality assurance standards, using single source samples, is reported for the end-to-end workflow from library preparation to data interpretation. Performance metrics support the conclusion that more genetic information can be obtained from challenging samples compared to other commercially available forensic targeted DNA assays developed for capillary electrophoresis (CE) or other current next generation sequencing (NGS) kits due to the higher number of markers, the overall shorter amplicon sizes (97.8% <150 bp), and kit design. Data indicate that the multiplex is robust and fit for purpose for a wide range of quantity and quality samples. The ForenSeq Kintelligence Kit and the Universal Analysis Software allow transfer of the genetic component of forensic investigative genetic genealogy to the operational forensic laboratory.

High-resolution genotyping of 58 STRs in 635 Northern Han Chinese with MiSeq FGx ® Forensic Genomics System

Sequence polymorphisms were characterized at 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), and 24 Y chromosomal STRs (Y-STRs) in 635 Northern Han Chinese with the ForenSeq DNA Signature Prep Kit on the MiSeq FGx Forensic Genomics System. Since repeat region (RR) and flanking region (FR) variation can be detected by massively parallel sequencing (MPS), the increase in the number of unique alleles and the average of gene diversity was 78.18% and 3.51% between sequence and length, respectively. A total of 74 novel RR variants were identified at 33 STRs compared with STRSeq and previous studies, and 13 FR variants (rs1770275883, rs2053373277, rs2082557941, rs1925525766, rs1926380862, rs1569322793, rs2051848492, rs2051848696, rs2016239814, rs2053269960, rs2044518192, rs2044536444, and rs2089968964) were first submitted to dbSNP. Also, 99.94% of alleles were concordant between the ForenSeq DNA Signature Prep Kit and commercial CE kits. Discordance resulted from the low performance at D22S1045 and occasionally at DYS392, flanking region deletions at D7S820 and DXS10074, and the strict alignment algorithm at DXS7132. Null alleles at DYS505 and DYS448 and multialleles at DYS387S1a/b, DYS385a/b, DYS448, DYS505, DXS7132, and HPRTB were validated with other MPS and CE kits. Thus, a high-resolution sequence-based (SB) and length-based (LB) allele frequencies dataset from Northern Han Chinese has been established already. As expected, forensic parameters increased significantly on combined power of discrimination (PD) and combined power of exclusion (PE) at A-STRs, mildly on combined PD and combined mean exclusion chance (MEC) at X-STRs, and barely on discrimination capacity (DC) at Y-STRs. Additionally, MiSeq FGx quality metrics and MPS performance were evaluated in this study, which presented the high-quality of the dataset at 20 consecutive runs, such as ≥ 60% bases with a quality score of 20 or higher (%≥ Q20), > 60% of effective reads, > 2000 × of depth of coverage (DoC), ≥ 60% of allele coverage ratio (ACR) or heterozygote balance, ≥ 70% of inter-locus balance, and ≤ 0.4 of the absolute value of observed minus expected heterozygosity (|H_exp - H_obs|). In conclusion, MiSeq FGx can help us generate a high-resolution and high-quality dataset for human identification and population genetic studies.

Developmental validation of the ForenSeq MainstAY kit, MiSeq FGx sequencing system and ForenSeq Universal Analysis Software

For human identification purposes, forensic genetics has primarily relied upon a core set of autosomal (and to a lesser extent Y chromosome) short tandem repeat (STR) markers that are enriched by amplification using the polymerase chain reaction (PCR) that are subsequently separated and detected using capillary electrophoresis (CE). While STR typing conducted in this manner is well-developed and robust, advances in molecular biology that have occurred over the last 15 years, in particular massively parallel sequencing (MPS) [1-7], offer certain advantages as compared to CE-based typing. First and foremost is the high throughput capacity of MPS. Current bench top high throughput sequencers enable larger batteries of markers to be multiplexed and multiple samples to be sequenced simultaneously (e.g., millions to billions of nucleotides can be sequenced in one run). Second, compared to the length-based CE approach, sequencing STRs increases discrimination power, enhances sensitivity of detection, reduces noise due to instrumentation, and improves mixture interpretation [4,8-23]. Third, since detection of STRs is based on sequence and not fluorescence, amplicons can be designed that are shorter in length and of similar lengths among loci, where possible, which can improve amplification efficiency and analysis of degraded samples. Lastly, MPS offers a single format approach that can be applied to analysis of a wide variety of genetic markers of forensic interest (e.g., STRs, mitochondrial DNA, single nucleotide polymorphisms, insertion/deletions). These features make MPS a desirable technology for casework [14,15,24,25-48]. The developmental validation of the ForenSeq MainstAY library preparation kit with the MiSeq FGx Sequencing System and ForenSeq Universal Software is reported here to assist with validation of this MPS system for casework [49]. The results show that the system is sensitive, accurate and precise, specific, and performs well with mixtures and mock case-type samples.

Optimization of InnoXtract™ extraction and purification system for DNA extraction from skeletal samples

The InnoXtract™ extraction and purification system is a purification method designed for DNA extraction from low-template samples, specifically rootless hair shafts. Its ability to successfully capture highly fragmented DNA suggests its suitability for use with other challenging sample types, including skeletal remains. However, the lysis and digestion parameters required modifications to successfully optimize the method for this sample type. A two-part digestion was developed utilizing a homebrew digestion buffer (0.5 M EDTA, 0.05% Tween 20, and 100 mM NaCl) and a supplemental lysis with the Hair Digestion Buffer included in the InnoXtract™ kit. Additionally, the magnetic bead volume was modified to improve DNA recovery from these challenging samples. With the altered protocol, the quality and quantity of DNA recovered from InnoXtract™ extracts were comparable to another commercial skeletal extraction method (PrepFiler™ BTA). This modified extraction method successfully purified sufficient amounts of quality DNA from a variety of skeletal samples to produce complete STR profiles. Successful STR typing from surface decomposition, burned, cremated, buried, and embalmed remains indicates the potential of this new method for challenging human identification and missing-person cases.

High-throughput DNA sequencing of environmentally insulted latent fingerprints after visualization with nanoscale columnar-thin-film technique

The goals of this study were to (a) ascertain human identity capabilities of DNA obtained from latent fingerprints that have been first environmentally insulted and then developed by the deposition of a columnar thin film (CTF), and (b) to determine if the CTF process and material are detrimental to single nucleotide polymorphism (SNP) analysis. Fingerprints were deposited on five different types of substrates and aged for one day, 7 days or 30 days while being environmentally insulted under one of the four conditions: 16.6 °C and 60% relative humidity (RH) (Condition A), 24.5 °C and 60% RH (Condition B), 35 °C and 67% RH (Condition C) and a cold condition (Condition D). Then CTF technique was then on 59% of these fingerprints. DNA samples from 805 fingerprints were extracted, quantified, subjected to manual library preparation using the Precision ID Identity Panel, and underwent high-throughput sequencing. The Ion S5™ platform was employed to sequence 124 SNP amplicons. SNPs were successfully sequenced from 802/805 samples. Total read depth was consistent across environmental conditions, and majority of samples had 100% profile completeness and 100% concordance. Anecdotally, libraries that were amplified with a higher cycle number had more 'Major Allele Frequency' flags compared to samples amplified with 23 cycle numbers, possibly due to stochastic effects. Neither the substrates nor the CTF process and materials inhibit downstream DNA analysis. DNA of low quality and quantity from the chosen samples can be sequenced using the Precision ID Identity Panel on the Ion S5™ platform which performed well, however, a different approach may be needed if spurious alleles are suspected.

Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Mitochondrial DNA (mtDNA) analysis using Sanger sequencing has been a routine practice for the identification of human skeletal remains. However, this process is usually challenging since DNA from the remains is highly degraded and at low concentration. Recently, the advent and implementation of massively parallel sequencing (MPS) have been offered the ability to improve mtDNA sequence data for forensic analysis. To assess the utility of the Ion S5™ system - an MPS platform for mtDNA analysis in challenging samples, we sequenced the mitochondrial control region of 52 age-old skeletal remains. Using the Precision ID mtDNA Control Region Panel, 50 full and two partial control region haplotypes at relatively high mean coverage of 2494 × were achieved for variant calling. Further variant analysis at 10% threshold for point heteroplasmy showed high degradation degree in terms of DNA damage in our bone samples. A higher point heteroplasmy threshold of 20% was required to diminish most of background noise caused by the damage. The results from this study indicated the potential application of the Ion S5™ system in sequencing degraded samples in Vietnam and provided valuable data sources for forensic analyses in the future.

The forensic landscape and the population genetic analyses of Hainan Li based on massively parallel sequencing DNA profiling

Due to the formation of the Qiongzhou Strait by climate change and marine transition, Hainan island was isolated from the mainland southern China during the Last Glacial Maximum. Hainan island, located at the southernmost part of China and separated from the Leizhou Peninsula by the Qiongzhou Strait, laid on one of the modern human northward migration routes from Southeast Asia to East Asia. The Hlai language-speaking Li minority, the second largest population after Han Chinese in Hainan island, is the direct descendants of the initial migrants in Hainan island and has unique ethnic properties and derived characteristics; however, the forensic-associated studies on Hainan Li population are still insufficient. Hence, 136 Hainan Li individuals were genotyped in this study using the MPS-based ForenSeq™ DNA Signature Prep Kit (DNA Primer Set A, DPMA) to characterize the forensic genetic polymorphism landscape, and DNA profiles were obtained from 152 different molecular genetic markers (27 autosomal STRs, 24 Y-STRs, 7 X-STRs, and 94 iiSNPs). A total of 419 distinct length variants and 586 repeat sequence sub-variants, with 31 novel alleles (at 17 loci), were identified across the 58 STR loci from the DNA profiles of Hainan Li population. We evaluated the forensic characteristics and efficiencies of DPMA, demonstrating that the STRs and iiSNPs in DPMA were highly polymorphic in Hainan Li population and could be employed in forensic applications. In addition, we set up three datasets, which included the genetic data of (i) iiSNPs (27 populations, 2640 individuals), (ii) Y-STRs (42 populations, 8281 individuals), and (iii) Y haplogroups (123 populations, 4837 individuals) along with the population ancestries and language families, to perform population genetic analyses separately from different perspectives. In conclusion, the phylogenetic analyses indicated that Hainan Li, with a southern East Asia origin and Tai-Kadai language-speaking language, is an isolated population relatively. But the genetic pool of Hainan Li influenced by the limited gene flows from other Tai-Kadai populations and Hainan populations. Furthermore, the establishment of isolated population models will be beneficial to clarify the exquisite population structures and develop specific genetic markers for subpopulations in forensic genetic fields.

Large fragment demineralization: an alternative pretreatment for forensic DNA typing of bones

In forensic genetics, the analysis of postmortem bones is one of the most challenging due to the low quantity of degraded endogenous DNA. The most widely used approach for sample preparation, in those cases, is pulverizing the bone. However, processing pulverized bone is extremely delicate, requiring strict laboratory conditions and operating procedures. In fact, several recent publications have focused on non-powder approaches. The objectives of this study were, thus, to validate a non-powder protocol for DNA extraction from forensic bones and an alternative pretreatment, large fragment demineralization (LFD). Thirty human femurs and tibiae received by the Legal Medicine Institute of Brescia, Italy, were included in the study. Bone powder and one transversal section of the diaphysis were sampled from each bone. DNA extraction from the powder was carried out using PrepFiler BTA (BTA), while the transversal section was submitted to the alternative demineralizing pretreatment (LFD) followed by DNA extraction using the QIAamp DNA Investigator. DNA extracts were assessed for human DNA quantity and degradation by means of a validated in-house qPCR assay and amplified with commercial kits. Inhibition assessment was carried out through Quality Sensor analysis using 24plex QS Kit. The differences in quantity, quality of human DNA, and number of alleles detected between both methods were comparable and not statistically significant. We propose the use of the LFD protocol as a complementary approach capable of confirming the genotypes or detect alleles not observed using BTA, without the need for pulverization.

Massively parallel sequencing and STR analysis from partial bloody fingerprints enhanced with columnar thin films

Fingerprint enhancement often includes either physical or chemical approaches, such as fingerprint powder or cyanoacrylate fuming, to improve the quality of a fingerprint for visualization and analysis. However, these methods become more complex when fingerprints are partial bloody, and these procedures may interfere with downstream DNA analysis. Columnar thin film (CTF) deposition is a type of nanotechnology that utilizes an evaporant material to enhance a fingerprint under low-pressure conditions. Short tandem repeat (STR) analysis is the traditional method employed in crime laboratories. When DNA is of poor quality and quantity, like that often obtained from fingerprints, little to no genetic information may be obtained. Single nucleotide polymorphisms (SNPs) may be used to glean additional information when STR analysis fails. In this pilot study, 100 partial bloody fingerprints were collected from two donors and deposited on five different crime scene substrates, in which half were enhanced with CTFs and were graded for quality by an IAI-certified latent fingerprint examiner. CTF-developed fingerprints, on average, had higher grades compared to non-developed partial bloody fingerprints. STR analysis using Fusion 6C was performed to assess inhibition from the evaporant materials, in which no inhibition was observed. Sequencing of SNPs using the Precision ID Identity Panel was also employed, in which genetic information that could not be obtained from STRs was acquired with SNPs. Various sample types (i.e. pristine, low quality, and contaminated) utilized in this project demonstrated the acceptable performance of the Precision ID Identity Panel.

The challenge of predicting human pigmentation traits in degraded bone samples with the MPS-based HIrisPlex-S system

Identification of human remains is an important part of human DNA analysis studies. STR and mitochondrial DNA markers are well suited for the analysis of degraded biological samples including bone material. However, these DNA markers may be useless when reference material is not available. In these cases, predictive DNA analysis can support the process of human identification by providing investigative leads. Forensic DNA phenotyping has progressed significantly by offering new methods based on massively parallel sequencing technology, but the frequent degradation processes observed in skeletal remains can make analysis of such samples challenging. In this study, we demonstrate the usefulness of a recently established Ion AmpliSeq^TM HIrisPlex-S panel using Ion Torrent technology for analyzing bone samples that show different levels of DNA degradation. In total, 63 bone samples at post-mortem intervals up to almost 80 years were genotyped and eye, hair and skin colour predictions were performed using the HIrisPlex-S models. Following the recommended coverage thresholds, it was possible to establish full DNA profiles comprising of 41 DNA variants for 35 samples (55.6%). For 5 samples (7.9%) no DNA profiles were generated. The remaining 23 samples (36.5%) produced partial profiles and showed a clear underperformance of 3 HIrisPlex-S SNPs - rs1545397 (OCA2), rs1470608 (OCA2) and rs10756819 (BNC2), all used for skin colour prediction only. None of the 23 samples gave complete genotypes needed for skin colour prediction was obtained, and in 7 of them (25.9%) the 3 underperformed SNPs were the cause. At the same time, the prediction of eye and hair colour using complete IrisPlex and HIrisPlex profiles could be made for these 23 samples in 20 (87.0%) and 12 cases (52.2%), respectively. Complete HIrisPlex-S profiles were generated from as little as 49 pg of template DNA. Five samples for which the HIrisPlex-S analysis failed, consistently failed in standard STR analysis. Importantly, the 3 underperforming SNPs produced significantly lower number of reads in good quality samples. Nonetheless, the AUC loss resulting from missing data for these 3 SNPs is not considered large (≤0.004) and the prediction of pigmentation from partial profiles is also available in the current HPS tool. The study shows that DNA degradation and the resulting loss of data are the most serious challenge to DNA phenotyping of skeletal remains. Although the newly developed HIrisPlex-S panel has been successfully validated in the current research, primer redesign for the 3 underperforming SNPs in the MPS design should be considered in the future.

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.

Utility of the Ion S5™ and MiSeq FGx™ sequencing platforms to characterize challenging human remains

Often in missing persons' and mass disaster cases, the samples remaining for analysis are hard tissues such as bones, teeth, nails, and hair. These remains may have been exposed to harsh environmental conditions, which pose challenges for downstream genotyping. Short tandem repeat analysis (STR) via capillary electrophoresis (CE) is still the gold standard for DNA typing; however, a newer technology known as massively parallel sequencing (MPS) could improve upon our current techniques by typing different and more markers in a single analysis, and consequently improving the power of discrimination. In this study, bone and tooth samples exposed to a variety of DNA insults (cremation, embalming, decomposition, thermal degradation, and fire) were assessed and sequenced using the Precision ID chemistry and a custom AmpliSeq™ STR and iiSNP panel on the Ion S5™ System, and the ForenSeq DNA Signature Prep Kit on the MiSeq FGx™ system, as well as the GlobalFiler™ PCR Amplification Kit on the 3500™ Genetic Analyzer. The results demonstrated that using traditional CE-based genotyping performed as expected, producing a partial or full DNA profile for all samples, and that both sequencing chemistries and platforms were able to recover sufficient STR and SNP information from a majority of the same challenging samples. Run metrics including profile completeness and mean read depth produced good results with each system, considering the degree of damage of some samples. Most sample insults (except decomposed) produced similar numbers of alleles for both MPS systems. Comparable markers produced full concordance between the two platforms.