Unearthing who and Y at Harewood Cemetery and inference of George Washington's Y-chromosomal haplotype.

An excavation conducted at Harewood Cemetery to identify the unmarked grave of Samuel Washington resulted in the discovery of burials presumably belonging to George Washington's paternal grandnephews and their mother, Lucy Payne. To confirm their identities this study examined Y-chromosomal, mitochondrial, and autosomal DNA from the burials and a living Washington descendant. The burial's Y-STR profile was compared to FamilyTreeDNA's database, which resulted in a one-step difference from the living descendant and an exact match to another Washington. A more complete Y-STR and Y-SNP profile from the descendant was inferred to be the Washington Y profile. Kinship comparisons performed in relation to the descendant, who is a 4th and 5th degree relative of the putative individuals, resulted in >37,000 overlapping autosomal SNPs and strong statistical support with likelihood ratios exceeding one billion. This study highlights the benefits of a multi-marker approach for kinship prediction and DNA-assisted identification of historical remains.

Complete Mitochondrial DNA Genome Variation in the Swedish Population.

The development of complete mitochondrial genome (mitogenome) reference data for inclusion in publicly available population databases is currently underway, and the generation of more high-quality mitogenomes will only enhance the statistical power of this forensically useful locus. To characterize mitogenome variation in Sweden, the mitochondrial DNA (mtDNA) reads from the SweGen whole genome sequencing (WGS) dataset were analyzed. To overcome the interference from low-frequency nuclear mtDNA segments (NUMTs), a 10% variant frequency threshold was applied for the analysis. In total, 934 forensic-quality mitogenome haplotypes were characterized. Almost 45% of the SweGen haplotypes belonged to haplogroup H. Nearly all mitogenome haplotypes (99.1%) were assigned to European haplogroups, which was expected based on previous mtDNA studies of the Swedish population. There were signature northern Swedish and Finnish haplogroups observed in the dataset (e.g., U5b1, W1a), consistent with the nuclear DNA analyses of the SweGen data. The complete mitogenome analysis resulted in high haplotype diversity (0.9996) with a random match probability of 0.15%. Overall, the SweGen mitogenomes provide a large mtDNA reference dataset for the Swedish population and also contribute to the effort to estimate global mitogenome haplotype frequencies.

Evaluating the Usefulness of Human DNA Quantification to Predict DNA Profiling Success of Historical Bone Samples.

This study assessed the usefulness of DNA quantification to predict the success of historical samples when analyzing SNPs, mtDNA, and STR targets. Thirty burials from six historical contexts were utilized, ranging in age from 80 to 800 years postmortem. Samples underwent library preparation and hybridization capture with two bait panels (FORCE and mitogenome), and STR typing (autosomal STR and Y-STR). All 30 samples generated small (~80 bp) autosomal DNA target qPCR results, despite mean mappable fragments ranging from 55-125 bp. The qPCR results were positively correlated with DNA profiling success. Samples with human DNA inputs as low as 100 pg resulted in ≥80% FORCE SNPs at 10X coverage. All 30 samples resulted in mitogenome coverage ≥100X despite low human DNA input (as low as 1 pg). With PowerPlex Fusion, ≥30 pg human DNA input resulted in >40% of auSTR loci. At least 59% of Y-STR loci were recovered with Y-target qPCR-based inputs of ≥24 pg. The results also indicate that human DNA quantity is a better predictor of success than the ratio of human to exogenous DNA. Accurate quantification with qPCR is feasible for historical bone samples, allowing for the screening of extracts to predict the success of DNA profiling.

The Value of Whole-Genome Sequencing for Mitochondrial DNA Population Studies: Strategies and Criteria for Extracting High-Quality Mitogenome Haplotypes.

Whole-genome sequencing (WGS) data present a readily available resource for mitochondrial genome (mitogenome) haplotypes that can be utilized for genetics research including population studies. However, the reconstruction of the mitogenome is complicated by nuclear mitochondrial DNA (mtDNA) segments (NUMTs) that co-align with the mtDNA sequences and mimic authentic heteroplasmy. Two minimum variant detection thresholds, 5% and 10%, were assessed for the ability to produce authentic mitogenome haplotypes from a previously generated WGS dataset. Variants associated with NUMTs were detected in the mtDNA alignments for 91 of 917 (~8%) Swedish samples when the 5% frequency threshold was applied. The 413 observed NUMT variants were predominantly detected in two regions (nps 12,612-13,105 and 16,390-16,527), which were consistent with previously documented NUMTs. The number of NUMT variants was reduced by ~97% (400) using a 10% frequency threshold. Furthermore, the 5% frequency data were inconsistent with a platinum-quality mitogenome dataset with respect to observed heteroplasmy. These analyses illustrate that a 10% variant detection threshold may be necessary to ensure the generation of reliable mitogenome haplotypes from WGS data resources.

Ancient DNA Methods Improve Forensic DNA Profiling of Korean War and World War II Unknowns.

The integration of massively parallel sequencing (MPS) technology into forensic casework has been of particular benefit to the identification of unknown military service members. However, highly degraded or chemically treated skeletal remains often fail to provide usable DNA profiles, even with sensitive mitochondrial (mt) DNA capture and MPS methods. In parallel, the ancient DNA field has developed workflows specifically for degraded DNA, resulting in the successful recovery of nuclear DNA and mtDNA from skeletal remains as well as sediment over 100,000 years old. In this study we use a set of disinterred skeletal remains from the Korean War and World War II to test if ancient DNA extraction and library preparation methods improve forensic DNA profiling. We identified an ancient DNA extraction protocol that resulted in the recovery of significantly more human mtDNA fragments than protocols previously used in casework. In addition, utilizing single-stranded rather than double-stranded library preparation resulted in increased attainment of reportable mtDNA profiles. This study emphasizes that the combination of ancient DNA extraction and library preparation methods evaluated here increases the success rate of DNA profiling, and likelihood of identifying historical remains.

Extended kinship analysis of historical remains using SNP capture.

DNA-assisted identification of historical remains requires the genetic analysis of highly degraded DNA, along with a comparison to DNA from known relatives. This can be achieved by targeting single nucleotide polymorphisms (SNPs) using a hybridization capture and next-generation sequencing approach suitable for degraded skeletal samples. In the present study, two SNP capture panels were designed to target ~ 25,000 (25 K) and ~ 95,000 (95 K) nuclear SNPs, respectively, to enable distant kinship estimation (up to 4th degree relatives). Low-coverage SNP data were successfully recovered from 14 skeletal elements 75 years postmortem using an Illumina MiSeq benchtop sequencer. All samples contained degraded DNA but were of varying quality with mean fragment lengths ranging from 32 bp to 170 bp across the 14 samples. SNP comparison with DNA from known family references was performed in the Parabon Fx Forensic Analysis Platform, which utilizes a likelihood approach for kinship prediction that was optimized for low-coverage sequencing data with cytosine deamination. The 25 K panel produced 15,000 SNPs on average, which allowed for accurate kinship prediction with strong statistical support in 16 of the 21 pairwise comparisons. The 95 K panel increased the average SNPs to 42,000 and resulted in an additional accurate kinship prediction with strong statistical support (17 of 21 pairwise comparisons). This study demonstrates that SNP capture combined with massively parallel sequencing on a benchtop platform can yield sufficient SNP recovery from compromised samples, enabling accurate, extended kinship predictions.

The FORCE Panel: An All-in-One SNP Marker Set for Confirming Investigative Genetic Genealogy Leads and for General Forensic Applications.

The FORensic Capture Enrichment (FORCE) panel is an all-in-one SNP panel for forensic applications. This panel of 5422 markers encompasses common, forensically relevant SNPs (identity, ancestry, phenotype, X- and Y-chromosomal SNPs), a novel set of 3931 autosomal SNPs for extended kinship analysis, and no clinically relevant/disease markers. The FORCE panel was developed as a custom hybridization capture assay utilizing ~20,000 baits to target the selected SNPs. Five non-probative, previously identified World War II (WWII) cases were used to assess the kinship panel. Each case included one bone sample and associated family reference DNA samples. Additionally, seven reference quality samples, two 200-year-old bone samples, and four control DNAs were processed for kit performance and concordance assessments. SNP recovery after capture resulted in a mean of ~99% SNPs exceeding 10X coverage for reference and control samples, and 44.4% SNPs for bone samples. The WWII case results showed that the FORCE panel could predict first to fifth degree relationships with strong statistical support (likelihood ratios over 10,000 and posterior probabilities over 99.99%). To conclude, SNPs will be important for further advances in forensic DNA analysis. The FORCE panel shows promising results and demonstrates the utility of a 5000 SNP panel for forensic applications.

Capture enrichment and massively parallel sequencing for human identification.

In the past decade, hybridization capture has gained attention within the forensic field for its possible use in human identification. One of the primary benefits to capture enrichment is its applicability to degraded DNA fragments that, due to their reduced size, are not amenable to traditional PCR enrichment techniques. Hybridization capture is typically introduced after genomic library preparation of extracted DNA templates for the subsequent enrichment of mitochondrial DNA or single nucleotide polymorphisms within the nuclear genome. The enriched molecules are then subjected to massively parallel sequencing (MPS) for sensitive and high-throughput DNA sequence generation. Bioinformatic analysis of capture product removes PCR duplicates that were introduced during the laboratory workflow in order to characterize the original DNA template molecules. In the case of aged and degraded skeletal remains, the fraction of endogenous human DNA may be very low; therefore low-coverage sequence analysis may be required. This review contains an overview of current capture methodologies and the primary literature on hybridization capture as evaluated for forensic applications.

Platinum-Quality Mitogenome Haplotypes from United States Populations.

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.

Pathogenic Variant Filtering for Mitochondrial Genome Haplotype Reporting.

Given the enhanced discriminatory power of the mitochondrial DNA (mtDNA) genome (mitogenome) over the commonly sequenced control region (CR) portion, the scientific merit of mitogenome sequencing is generally accepted. However, many laboratories remain beholden to CR sequencing due to privacy policies and legal requirements restricting the use of disease information or coding region (codR) information. In this report, we present an approach to obviate the reporting of sensitive codR data in forensic haplotypes. We consulted the MitoMap database to identify 92 mtDNA codR variants with confirmed pathogenicity. We determined the frequencies of these pathogenic variants in literature-quality and forensic-quality databases to be very low, at 1.2% and 0.36%, respectively. The observed effect of pathogenic variant filtering on random match statistics in 2488 forensic-quality mitogenome haplotypes from four populations was nil. We propose that pathogenic variant filtering should be incorporated into variant calling algorithms for mitogenome haplotype reporting to maximize the discriminatory power of the locus while minimizing the reveal of sensitive genetic information.

Next generation sequencing of STR artifacts produced from historical bone samples.

STR artifacts are commonly observed in electrophoretic data and can complicate interpretation of the profiles produced. Even when a consensus approach is applied, reproducible artifacts have the potential to convolute the analysis. DNA obtained from historical bone samples is often heavily degraded and damaged, requiring the use of more sensitive procedures to increase allele recovery. Additionally, skeletal remains exposed to environmental conditions may be afflicted with microbial DNA contamination that cross-reacts with the primers during short tandem repeat (STR) multiplex amplification. STR artifacts manifested as a result of these circumstances can be sourced and characterized using new sequencing technologies to potentially ease the analysis burden. For this study, PCR product from 17 low-quality bone samples exhibiting reproducible autosomal and Y-chromosomal STR (Y-STR) artifacts in capillary electrophoresis (CE) data were sequenced with next generation sequencing (NGS). Sequenced reads were bioinformatically sorted using STRait Razor to determine the authenticity of alleles and confirm the profile generated by CE. Sequence data from the PCR products and a subset of the associated extracts were further analyzed with Kaiju to classify the microbial species present and identify potential sources of artifact peaks. A suspected Y-STR artifact was similar in sequence to Pseudomonas sp. BAY1663, a species ubiquitously found in soil. Regions of homology were observed between the Pseudomonas genome and the presumed primer binding locations for Y-STRs included in the AmpFlSTR Y-Filer STR kit. Characterization of such supposed artifact peaks may aid in interpretation of CE data and ultimately lead to increased confidence in the reported results.

A Forensic Genomics Approach for the Identification of Sister Marija Crucifiksa Kozulic.

Sister Marija Krucifiksa Kozulic (1852-1922) was a Croatian nun who is in consideration for beatification by the Vatican, which is facilitated by the identification of her 20th-century remains. Sister Marija was buried in a tomb in Rijeka, Croatia, along with other nuns including her biological sister, Tereza Kozulic (1861-1933). When the remains were exhumed in 2011, they were found in a deteriorated state and commingled with several other sets of remains. Thus, mitochondrial genome sequencing of the long bones was performed to sort the remains by mitochondrial haplotype. Two similar but unique haplotypes belonging to haplogroup H1bu were identified, and samples from these bones were subjected to autosomal short tandem repeat (STR) and single nucleotide polymorphism (SNP) sequencing. Although only partial profiles were obtained, the data were sufficient for kinship analysis with the profile of a paternal niece of Sister Marija (Fides Kozulic). The data indicate that it is 574,195-fold more likely that the two sets of skeletal remains represent 2nd-degree relatives of Fides than sisters who are unrelated to Fides. Although it is impossible to discern which set of remains belongs to Marija and which belongs to Tereza, forensic genomics methods have enabled identification of the sisters.

Mitochondrial DNA haplogrouping to assist with the identification of unknown service members from the World War II Battle of Tarawa.

The World War II Battle of Tarawa, 1943, was a devastating conflict that resulted in losses of more than 1100 American and 4690 Japanese troops. The United States government aims to identify and repatriate the remains of all missing American service members through the Defense Prisoner of War/Missing in Action (POW/MIA) Accounting Agency (DPAA) and its partners such as the Armed Forces Medical Examiner System's Armed Forces DNA Identification Laboratory (AFMES-AFDIL). Remains associated with the Battle of Tarawa have been recovered from field excavations conducted by History Flight, a DPAA strategic partner, as well as from the National Memorial Cemetery of the Pacific (NMCP) in Hawaii where unknowns have been disinterred for identification. DNA testing at the AFMES-AFDIL has produced mitochondrial DNA (mtDNA) sequences from 1027 case samples to date. Haplogroup assignments indicate that more than one third (36.2 %) of field-collected samples are likely of Asian maternal ancestry. Therefore the field collections from the Tarawa battlefield comprise the remains of American service members but also those of foreign nationals from Asia. The mtDNA of the NMCP unknowns is similar in ancestry proportion to the family reference sample distribution. The DPAA uses the ancestry information gleaned from mtDNA sequence data in conjunction with anthropological evidence to make foreign national determinations. In this way, mtDNA haplogrouping is used to sort the commingled and fragmentary remains recovered from Tarawa between Americans and foreign nationals, which are then repatriated to their country of origin.

Impact of the sequencing method on the detection and interpretation of mitochondrial DNA length heteroplasmy.

Advancements in sequencing technologies allow for rapid and efficient analysis of mitochondrial DNA (mtDNA) in forensic laboratories, which is particularly beneficial for specimens with limited nuclear DNA. Next generation sequencing (NGS) offers higher throughput and sensitivity over traditional Sanger-type sequencing (STS) as well as the ability to quantitatively analyze the data. Changes in sample preparation, sequencing method and analysis required for NGS may alter the mtDNA haplotypes compared to previously generated STS data. Thus, the present study aimed to characterize the impact of different sequencing workflows on the detection and interpretation of length heteroplasmy (LHP), a particularly complicated aspect of mtDNA analysis. Whole mtDNA genome (mitogenome) data were generated for 16 high-quality samples using well-established Illumina and Ion methods, and the NGS data were compared to previously-generated STS mtDNA control region data. Although the mitogenome haplotypes were concordant with the exception of length and low-level variants (<30 % variant frequency), LHP in the hypervariable segment (HVS) polycytosine regions (C-tracts) differed across sequencing methods. Consistent with previous studies, LHP in HVS1 was observed in samples with nine or more consecutive cytosines (Cs) and eight Cs in the HVS2 region in the STS data. The Illumina data produced a similar pattern of LHP as the STS data, whereas the Ion data were noticeably different. More complex LHP (i.e. more length molecules) was observed in the Ion data, as length variation occurred in multiple homopolymer stretches within the targeted HVS regions. Further, the STS dominant or major molecule (MM) differed from the Ion MM in 11 (37 %) of the 30 regions evaluated and six instances (20 %) in Illumina data. This is of particular interest, as the MM is used by many forensic laboratories to report the HVS C-tract in the mtDNA haplotype. In general, the STS MMs were longer than the Illumina MMs, while the Ion MMs were the shortest. The higher rate of homopolymer indels in Ion data likely contributed to these differences. Supplemental analysis with alternative approaches demonstrated that the LHP pattern may also be altered by the bioinformatic tool and workflow used for data interpretation. The broader application of NGS in forensic laboratories will undoubtedly result in the use of varying sample preparation and sequencing methods. Based on these findings, minor LHP differences are expected across sequencing workflows, and it will be important that C-tract indels continue to be ignored for forensic queries and comparisons.

Resolving mitochondrial haplogroups B2 and B4 with next-generation mitogenome sequencing to distinguish Native American from Asian haplotypes.

Mitochondrial haplogroup information can be useful in forensic contexts that rely primarily on mitochondrial DNA (mtDNA) testing, which often involve limited or degraded DNA. Due to the phylogeographic patterning of mtDNA in human populations, mitochondrial haplogroups are indicative of maternal ancestry (as mtDNA is a maternally inherited marker). In certain circumstances, maternal ancestry inferred from mitochondrial haplogrouping could be beneficial to forensic investigations. For example, ancestry information could assist in the identification of unknown service members from past conflicts, such as the World War II Battle of Tarawa involving American and Japanese forces. In this context, it could be useful to distinguish Native American mtDNA from Asian mtDNA to bolster the anthropological and circumstantial evidence leading to an identification or foreign national determination. Although most of the founding Native American haplogroups contain diagnostic variants in the mitochondrial control region (CR), haplogroup B2 does not, and this makes it more difficult to distinguish B2 from the parental B4 and closely related B4b haplogroups found in Asia. In this paper, the amount of mtDNA information required to distinguish Native American haplotypes from Asian haplotypes within haplogroup B was examined. Fifty-six samples belonging to subtypes of B2 and B4 were sequenced for the entire mitogenome. Haplogroups were estimated from three ranges of mitochondrial DNA (HV1 and 2, CR, and full mitogenome). Half of the samples could not be precisely haplogrouped without full mitogenome data, although enough variants were often provided to make an accurate B2 versus B4 distinction. Native American B2 haplotypes were distinguishable using CR data alone in 82% of samples, though the remaining samples required full mitogenome data for haplogroup B2 designation. The use of full mitogenome data consistently enables accurate haplogroup determination, and opens the possibility for gaining information on maternal ancestry.

Mitochondrial DNA control region variation in Lebanon, Jordan, and Bahrain.

This study investigated the mitochondrial DNA (mtDNA) control region variation in Middle Eastern populations (610 individuals from Lebanon, Jordan and the Kingdom of Bahrain) for which population data are scarce. FST comparison among populations revealed that there are significant differences in mtDNA distributions between Bahrain and the two other populations, while Lebanon and Jordan showed no significant differences. This was also reflected by the distribution of the observed lineages that differed prominently between Bahrain and the other two investigated populations. Jordan and Lebanon fit the hitherto known genetic results of the Levant population. Data are available via EMPOP (https://empop.online) and GenBank.

Bioinformatic removal of NUMT-associated variants in mitotiling next-generation sequencing data from whole blood samples.

Nuclear mitochondrial DNA segments (NUMTs) have arisen because of the transposition of segments of the mitochondrial DNA genome (mitogenome) into the nuclear genome. When using a "mitotiling" strategy, NUMTs may be more readily amplified when targeting the entire mitogenome compared to the control region, as hundreds of primers are required for complete sequencing coverage. In samples with a high percentage of nuclear DNA copies per cell, such as whole blood, NUMT coenrichment may be exacerbated. The present study examined bioinformatic approaches for removing NUMTs and NUMT-associated variants (NAVs) from next-generation sequence data generated using two mitotiling kits (Precision ID and QIAseq). Across 16 samples with low mtDNA copy number, NUMT coenrichment produced 890 NAVs with >5% variant frequency. The use of the consensus sequence to eliminate NUMT reads proved to be effective for QIAseq data, and resulted in >85% NAV removal in Precision ID data. This method was bolstered by NAV filtering in Precision ID analysis. Alternative high stringency mapping to the revised Cambridge Reference Sequence (rCRS) and the human genome reference GRCh38 for the QIAseq data caused a reduction in mitogenome coverage without complete NUMT removal. These bioinformatic solutions facilitate mitotiling sequence data analysis for low-level variant detection.

Repair of DNA damage caused by cytosine deamination in mitochondrial DNA of forensic case samples.

DNA sequence damage from cytosine deamination is well documented in degraded samples, such as those from ancient and forensic contexts. This study examined the effect of a DNA repair treatment on mitochondrial DNA (mtDNA) from aged and degraded skeletal samples. DNA extracts from 21 non-probative, degraded skeletal samples (aged 50-70 years) were utilized for the analysis. A portion of each sample extract was subjected to DNA repair using a commercial repair kit, the New England BioLabs' NEBNext FFPE DNA Repair Kit (Ipswich, MA). MtDNA was enriched using PCR and targeted capture in a side-by-side experiment of untreated and repaired DNA. Sequencing was performed using both traditional (Sanger-type; STS) and next-generation sequencing (NGS) methods Although cytosine deamination was evident in the mtDNA sequence data, the observed level of damaged bases varied by sequencing method as well as by enrichment type. The STS PCR amplicon data did not show evidence of cytosine deamination that could be distinguished from background signal in either the untreated or repaired sample set. However, the same PCR amplicons showed 850 C → T/G → A substitutions consistent with cytosine deamination with variant frequencies (VFs) of up to 25% when sequenced using NGS methods The occurrence of base misincorporation due to cytosine deamination was reduced by 98% (to 10) in the NGS amplicon data after repair. The NGS capture data indicated low levels (1-2%) of cytosine deamination in mtDNA fragments that was effectively mitigated by DNA repair. The observed difference in the level of cytosine deamination between the PCR and capture enrichment methods can be attributed to the greater propensity for stochastic effects from the PCR enrichment technique employed (e.g., low template input, increased PCR cycles). Altogether these results indicate that DNA repair may be required when sequencing PCR-amplified DNA from degraded forensic case samples with NGS methods.

Developmental validation of a Nextera XT mitogenome Illumina MiSeq sequencing method for high-quality samples.

Generating mitochondrial genome (mitogenome) data from reference samples in a rapid and efficient manner is critical to harnessing the greater power of discrimination of the entire mitochondrial DNA (mtDNA) marker. The method of long-range target enrichment, Nextera XT library preparation, and Illumina sequencing on the MiSeq is a well-established technique for generating mitogenome data from high-quality samples. To this end, a validation was conducted for this mitogenome method processing up to 24 samples simultaneously along with analysis in the CLC Genomics Workbench and utilizing the AQME (AFDIL-QIAGEN mtDNA Expert) tool to generate forensic profiles. This validation followed the Federal Bureau of Investigation's Quality Assurance Standards (QAS) for forensic DNA testing laboratories and the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines. The evaluation of control DNA, non-probative samples, blank controls, mixtures, and nonhuman samples demonstrated the validity of this method. Specifically, the sensitivity was established at ≥25 pg of nuclear DNA input for accurate mitogenome profile generation. Unreproducible low-level variants were observed in samples with low amplicon yields. Further, variant quality was shown to be a useful metric for identifying sequencing error and crosstalk. Success of this method was demonstrated with a variety of reference sample substrates and extract types. These studies further demonstrate the advantages of using NGS techniques by highlighting the quantitative nature of heteroplasmy detection. The results presented herein from more than 175 samples processed in ten sequencing runs, show this mitogenome sequencing method and analysis strategy to be valid for the generation of reference data.

Performance evaluation of a mitogenome capture and Illumina sequencing protocol using non-probative, case-type skeletal samples: Implications for the use of a positive control in a next-generation sequencing procedure.

Next-generation ancient DNA technologies have the potential to assist in the analysis of degraded DNA extracted from forensic specimens. Mitochondrial genome (mitogenome) sequencing, specifically, may be of benefit to samples that fail to yield forensically relevant genetic information using conventional PCR-based techniques. This report summarizes the Armed Forces Medical Examiner System's Armed Forces DNA Identification Laboratory's (AFMES-AFDIL) performance evaluation of a Next-Generation Sequencing protocol for degraded and chemically treated past accounting samples. The procedure involves hybridization capture for targeted enrichment of mitochondrial DNA, massively parallel sequencing using Illumina chemistry, and an automated bioinformatic pipeline for forensic mtDNA profile generation. A total of 22 non-probative samples and associated controls were processed in the present study, spanning a range of DNA quantity and quality. Data were generated from over 100 DNA libraries by ten DNA analysts over the course of five months. The results show that the mitogenome sequencing procedure is reliable and robust, sensitive to low template (one ng control DNA) as well as degraded DNA, and specific to the analysis of the human mitogenome. Haplotypes were overall concordant between NGS replicates and with previously generated Sanger control region data. Due to the inherent risk for contamination when working with low-template, degraded DNA, a contamination assessment was performed. The consumables were shown to be void of human DNA contaminants and suitable for forensic use. Reagent blanks and negative controls were analyzed to determine the background signal of the procedure. This background signal was then used to set analytical and reporting thresholds, which were designated at 4.0X (limit of detection) and 10.0X (limit of quantiation) average coverage across the mitogenome, respectively. Nearly all human samples exceeded the reporting threshold, although coverage was reduced in chemically treated samples resulting in a ∼58% passing rate for these poor-quality samples. A concordance assessment demonstrated the reliability of the NGS data when compared to known Sanger profiles. One case sample was shown to be mixed with a co-processed sample and two reagent blanks indicated the presence of DNA above the analytical threshold. This contamination was attributed to sequencing crosstalk from simultaneously sequenced high-quality samples to include the positive control. Overall this study demonstrated that hybridization capture and Illumina sequencing provide a viable method for mitogenome sequencing of degraded and chemically treated skeletal DNA samples, yet may require alternative measures of quality control.