A high-throughput Sanger strategy for human mitochondrial genome sequencing.

BACKGROUND: A population reference database of complete human mitochondrial genome (mtGenome) sequences is needed to enable the use of mitochondrial DNA (mtDNA) coding region data in forensic casework applications. However, the development of entire mtGenome haplotypes to forensic data quality standards is difficult and laborious. A Sanger-based amplification and sequencing strategy that is designed for automated processing, yet routinely produces high quality sequences, is needed to facilitate high-volume production of these mtGenome data sets. RESULTS: We developed a robust 8-amplicon Sanger sequencing strategy that regularly produces complete, forensic-quality mtGenome haplotypes in the first pass of data generation. The protocol works equally well on samples representing diverse mtDNA haplogroups and DNA input quantities ranging from 50 pg to 1 ng, and can be applied to specimens of varying DNA quality. The complete workflow was specifically designed for implementation on robotic instrumentation, which increases throughput and reduces both the opportunities for error inherent to manual processing and the cost of generating full mtGenome sequences. CONCLUSIONS: The described strategy will assist efforts to generate complete mtGenome haplotypes which meet the highest data quality expectations for forensic genetic and other applications. Additionally, high-quality data produced using this protocol can be used to assess mtDNA data developed using newer technologies and chemistries. Further, the amplification strategy can be used to enrich for mtDNA as a first step in sample preparation for targeted next-generation sequencing.

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory.

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

Sequence-based autosomal STR characterization in four US populations using PowerSeq™ Auto/Y system.

The forensic science community is poised to utilize modern advances in massively parallel sequencing (MPS) technologies to better characterize biological samples with higher resolution. A critical component towards the advancement of forensic DNA analysis with these technologies is a comprehensive understanding of the diversity and population distribution of sequence-based short tandem repeat (STR) alleles. Here we analyzed 786 samples of individuals from different population groups, including four of the mostly commonly encountered in forensic casework in the USA. DNA samples were amplified with the PowerSeq™ Auto/Y System Prototype Kit (Promega Corp.), and sequencing was performed on an Illumina® MiSeq instrument. Sequence data were analyzed using a bioinformatics processing tool, Altius. For additional data analysis and profile comparison, capillary electrophoresis (CE) size-based STR genotypes were generated for a subset of individuals, and where possible, also with a second commercially available MPS STR assay. Autosomal STR loci were analyzed and frequencies were calculated based on sequence composition. Also, population genetics studies were performed, with Hardy-Weinberg equilibrium, polymorphic information content (PIC), and observed and expected heterozygosity all assessed. Overall, sequence-based allelic variants of the repeat region were observed in 20 out of 22 different STR loci commonly used in forensic DNA genotyping, with the highest number of sequence variation observed at locus D12S391. The highest increase in allelic diversity and in PIC through sequence-based genotyping was observed at loci D3S1358 and D8S1179. Such detailed sequence analysis, as the one performed in the present study, is important to help understand the diversity of sequence-based STR alleles across different populations and to demonstrate how such allelic variation can improve statistics used for forensic casework.

Use of the LUS in sequence allele designations to facilitate probabilistic genotyping of NGS-based STR typing results.

Some of the expected advantages of next generation sequencing (NGS) for short tandem repeat (STR) typing include enhanced mixture detection and genotype resolution via sequence variation among non-homologous alleles of the same length. However, at the same time that NGS methods for forensic DNA typing have advanced in recent years, many caseworking laboratories have implemented or are transitioning to probabilistic genotyping to assist the interpretation of complex autosomal STR typing results. Current probabilistic software programs are designed for length-based data, and were not intended to accommodate sequence strings as the product input. Yet to leverage the benefits of NGS for enhanced genotyping and mixture deconvolution, the sequence variation among same-length products must be utilized in some form. Here, we propose use of the longest uninterrupted stretch (LUS) in allele designations as a simple method to represent sequence variation within the STR repeat regions and facilitate - in the nearterm - probabilistic interpretation of NGS-based typing results. An examination of published population data indicated that a reference LUS region is straightforward to define for most autosomal STR loci, and that using repeat unit plus LUS length as the allele designator can represent greater than 80% of the alleles detected by sequencing. A proof of concept study performed using a freely available probabilistic software demonstrated that the LUS length can be used in allele designations when a program does not require alleles to be integers, and that utilizing sequence information improves interpretation of both single-source and mixed contributor STR typing results as compared to using repeat unit information alone. The LUS concept for allele designation maintains the repeat-based allele nomenclature that will permit backward compatibility to extant STR databases, and the LUS lengths themselves will be concordant regardless of the NGS assay or analysis tools employed. Further, these biologically based, easy-to-derive designations uphold clear relationships between parent alleles and their stutter products, enabling analysis in fully continuous probabilistic programs that model stutter while avoiding the algorithmic complexities that come with string based searches. Though using repeat unit plus LUS length as the allele designator does not capture variation that occurs outside of the core repeat regions, this straightforward approach would permit the large majority of known STR sequence variation to be used for mixture deconvolution and, in turn, result in more informative mixture statistics in the near term. Ultimately, the method could bridge the gap from current length-based probabilistic systems to facilitate broader adoption of NGS by forensic DNA testing laboratories.

DNA identification of "Earthquake McGoon" 50 years postmortem.

This report describes the genetic identification of James "Earthquake McGoon" McGovern, a WWII fighter ace who perished in Laos while providing supplies to French troops during the French Indochina war. Because reference samples were unavailable for all of the potential casualties, testing of the entire mitochondrial genome, autosomal STRs and Y-chromosomal STRs was performed to increase the genetic information available for analysis. Kinship analyses performed on the evidentiary data and numerous indirect family references for McGovern excluded other possible casualties and definitively established McGovern's identity. This particular case demonstrates the practical utility of novel research technologies and aggressive genetic typing protocols in the identification of aged, degraded remains.

Development of a multiplex single base extension assay for mitochondrial DNA haplogroup typing.

AIM: To provide a screening tool to reduce time and sample consumption when attempting mitochondrial DNA (mtDNA) haplogroup typing. METHODS: A single base primer extension assay was developed to enable typing, in a single reaction, of twelve mtDNA haplogroup specific polymorphisms. For validation purposes a total of 147 samples were tested including 73 samples successfully haplogroup typed using mtDNA control region (CR) sequence data, 20 samples inconclusively haplogroup typed by CR sequence data, 21 samples previously haplogroup typed using RFLP analysis, and 31 samples of known ancestral origin without previous haplogroup typing. Additionally, two highly degraded human bones embalmed and buried in the early 1950s were analyzed using the SNP multiplex. RESULTS: When the SNP multiplex was used to type the 96 previously CR sequenced specimens, an increase in haplogroup or macrohaplogroup assignment relative to conventional CR sequence analysis was observed. The single base extension assay was also successfully used to assign a haplogroup to decades-old, embalmed skeletal remains dating to World War II. CONCLUSION: The SNP multiplex was successfully used to obtain haplogroup status of highly degraded human bones, and demonstrated the ability to eliminate possible contributors. The SNP multiplex provides a low-cost, high throughput method for typing of mtDNA haplogroups A, B, C, D, E, F, G, H, L1/L2, L3, M, and N that could be useful for screening purposes for human identification efforts and anthropological studies.

Performance and concordance of the ForenSeq™ system for autosomal and Y chromosome short tandem repeat sequencing of reference-type specimens.

Though the utility of next-generation sequencing (NGS) technologies for forensic short tandem repeat (STR) typing has been evident for several years, commercially available assays and software solutions developed specifically to meet forensic needs have only recently become available. One of these, the ForenSeq™ DNA Signature Prep Kit (Illumina, Inc.) sequences 27 autosomal STR (aSTR) and 24 Y chromosome STR (Y-STR) loci (concurrent with additional nuclear markers) per multiplexed sample, with automated secondary and tertiary analyses of the data accomplished via the associated ForenSeq™ Universal Analysis Software (UAS). In this study we investigated the performance of the ForenSeq system for aSTR and Y-STR typing by examination of 151 sample libraries developed from high quality DNAs amplified at the target 1ng template. Utilizing PCR Primer Mix B, greater than 99.5% of aSTR loci and 97.0% of Y-STR loci were recovered when 42 or fewer sample libraries were pooled for sequencing. A direct comparison of UAS developed fragment length results to capillary electrophoresis (CE) based data identified only two allele call discrepancies when no UAS quality flag was triggered. Review of the ForenSeq data indicated that most samples with total sequence read counts exceeding 40,000 could be interpreted to develop nearly complete aSTR genotypes or Y-STR haplotypes. However, markers D22S1045 and DYS392 produced poor or inconsistent results even when sample read counts were greater than 85,000. Excluding these two loci, analyst-interpreted aSTR and Y-STR ForenSeq profiles were 99.96% and 100% concordant, respectively, with CE data. In addition to demonstrating concordance on par with other CE kit to kit comparisons, the results from this study will assist laboratories seeking to develop workflows for high volume processing and analysis of aSTRs and Y-STRs from reference-type specimens using the ForenSeq system.

Internal validation of STRmix™ for the interpretation of single source and mixed DNA profiles.

The interpretation of DNA evidence can entail analysis of challenging STR typing results. Genotypes inferred from low quality or quantity specimens, or mixed DNA samples originating from multiple contributors, can result in weak or inconclusive match probabilities when a binary interpretation method and necessary thresholds (such as a stochastic threshold) are employed. Probabilistic genotyping approaches, such as fully continuous methods that incorporate empirically determined biological parameter models, enable usage of more of the profile information and reduce subjectivity in interpretation. As a result, software-based probabilistic analyses tend to produce more consistent and more informative results regarding potential contributors to DNA evidence. Studies to assess and internally validate the probabilistic genotyping software STRmix™ for casework usage at the Federal Bureau of Investigation Laboratory were conducted using lab-specific parameters and more than 300 single-source and mixed contributor profiles. Simulated forensic specimens, including constructed mixtures that included DNA from two to five donors across a broad range of template amounts and contributor proportions, were used to examine the sensitivity and specificity of the system via more than 60,000 tests comparing hundreds of known contributors and non-contributors to the specimens. Conditioned analyses, concurrent interpretation of amplification replicates, and application of an incorrect contributor number were also performed to further investigate software performance and probe the limitations of the system. In addition, the results from manual and probabilistic interpretation of both prepared and evidentiary mixtures were compared. The findings support that STRmix™ is sufficiently robust for implementation in forensic laboratories, offering numerous advantages over historical methods of DNA profile analysis and greater statistical power for the estimation of evidentiary weight, and can be used reliably in human identification testing. With few exceptions, likelihood ratio results reflected intuitively correct estimates of the weight of the genotype possibilities and known contributor genotypes. This comprehensive evaluation provides a model in accordance with SWGDAM recommendations for internal validation of a probabilistic genotyping system for DNA evidence interpretation.

Genetic variation for 15 autosomal STR loci (PowerPlex 16) in a population sample from northern Greece.

Fifteen autosomal STR loci were typed in a population sample of 318 unrelated individuals from northern Greece. Allele frequencies were determined, and Hardy-Weinberg equilibrium tested, for the following loci: D5S818, D7S820, TH01, TPOX, vWA, CSF1PO, D16S539, D13S317, D8S1179, D3S1358, FGA, PENTA D, PENTA E, D21S11, D18S51.

Mitochondrial DNA heteroplasmy in the emerging field of massively parallel sequencing.

Long an important and useful tool in forensic genetic investigations, mitochondrial DNA (mtDNA) typing continues to mature. Research in the last few years has demonstrated both that data from the entire molecule will have practical benefits in forensic DNA casework, and that massively parallel sequencing (MPS) methods will make full mitochondrial genome (mtGenome) sequencing of forensic specimens feasible and cost-effective. A spate of recent studies has employed these new technologies to assess intraindividual mtDNA variation. However, in several instances, contamination and other sources of mixed mtDNA data have been erroneously identified as heteroplasmy. Well vetted mtGenome datasets based on both Sanger and MPS sequences have found authentic point heteroplasmy in approximately 25% of individuals when minor component detection thresholds are in the range of 10-20%, along with positional distribution patterns in the coding region that differ from patterns of point heteroplasmy in the well-studied control region. A few recent studies that examined very low-level heteroplasmy are concordant with these observations when the data are examined at a common level of resolution. In this review we provide an overview of considerations related to the use of MPS technologies to detect mtDNA heteroplasmy. In addition, we examine published reports on point heteroplasmy to characterize features of the data that will assist in the evaluation of future mtGenome data developed by any typing method.

Full mtGenome reference data: development and characterization of 588 forensic-quality haplotypes representing three U.S. populations.

Though investigations into the use of massively parallel sequencing technologies for the generation of complete mitochondrial genome (mtGenome) profiles from difficult forensic specimens are well underway in multiple laboratories, the high quality population reference data necessary to support full mtGenome typing in the forensic context are lacking. To address this deficiency, we have developed 588 complete mtGenome haplotypes, spanning three U.S. population groups (African American, Caucasian and Hispanic) from anonymized, randomly-sampled specimens. Data production utilized an 8-amplicon, 135 sequencing reaction Sanger-based protocol, performed in semi-automated fashion on robotic instrumentation. Data review followed an intensive multi-step strategy that included a minimum of three independent reviews of the raw data at two laboratories; repeat screenings of all insertions, deletions, heteroplasmies, transversions and any additional private mutations; and a check for phylogenetic feasibility. For all three populations, nearly complete resolution of the haplotypes was achieved with full mtGenome sequences: 90.3-98.8% of haplotypes were unique per population, an improvement of 7.7-29.2% over control region sequencing alone, and zero haplotypes overlapped between populations. Inferred maternal biogeographic ancestry frequencies for each population and heteroplasmy rates in the control region were generally consistent with published datasets. In the coding region, nearly 90% of individuals exhibited length heteroplasmy in the 12418-12425 adenine homopolymer; and despite a relatively high rate of point heteroplasmy (23.8% of individuals across the entire molecule), coding region point heteroplasmies shared by more than one individual were notably absent, and transversion-type heteroplasmies were extremely rare. The ratio of nonsynonymous to synonymous changes among point heteroplasmies in the protein-coding genes (1:1.3) and average pathogenicity scores in comparison to data reported for complete substitutions in previous studies seem to provide some additional support for the role of purifying selection in the evolution of the human mtGenome. Overall, these thoroughly vetted full mtGenome population reference data can serve as a standard against which the quality and features of future mtGenome datasets (especially those developed via massively parallel sequencing) may be evaluated, and will provide a solid foundation for the generation of complete mtGenome haplotype frequency estimates for forensic applications.

Toward increased utility of mtDNA in forensic identifications.

The utility of mtDNA in forensic identifications is limited by its low power of discrimination and the absence of high quality mtDNA databases. Single nucleotide polymorphisms (SNPs) in the control region outside of hypervariable regions I and II (HVI/HVII), and in the coding region of the mtDNA genome, can provide additional discrimination in mtDNA testing. We have identified particularly useful SNP sites via high throughput sequencing of the entire mtDNA genome. We report here two cases in which an 11-plex SNP assay (panel "A") targeting the most common HVI/HVII type successfully resolved two cases in which identifications could not be made on the basis of HVI/HVII sequencing. Additionally, we established a database of 286 samples for SNP panel "A" generated with robotic protocols. We have addressed the need for high quality mtDNA control region (CR) databases by developing robotic protocols for lab processing, and a carefully devised electronic data review process. A large-scale databasing effort targeting several populations underrepresented in current mtDNA databases is underway.

Development of forensic-quality full mtGenome haplotypes: success rates with low template specimens.

Forensic mitochondrial DNA (mtDNA) testing requires appropriate, high quality reference population data for estimating the rarity of questioned haplotypes and, in turn, the strength of the mtDNA evidence. Available reference databases (SWGDAM, EMPOP) currently include information from the mtDNA control region; however, novel methods that quickly and easily recover mtDNA coding region data are becoming increasingly available. Though these assays promise to both facilitate the acquisition of mitochondrial genome (mtGenome) data and maximize the general utility of mtDNA testing in forensics, the appropriate reference data and database tools required for their routine application in forensic casework are lacking. To address this deficiency, we have undertaken an effort to: (1) increase the large-scale availability of high-quality entire mtGenome reference population data, and (2) improve the information technology infrastructure required to access/search mtGenome data and employ them in forensic casework. Here, we describe the application of a data generation and analysis workflow to the development of more than 400 complete, forensic-quality mtGenomes from low DNA quantity blood serum specimens as part of a U.S. National Institute of Justice funded reference population databasing initiative. We discuss the minor modifications made to a published mtGenome Sanger sequencing protocol to maintain a high rate of throughput while minimizing manual reprocessing with these low template samples. The successful use of this semi-automated strategy on forensic-like samples provides practical insight into the feasibility of producing complete mtGenome data in a routine casework environment, and demonstrates that large (>2kb) mtDNA fragments can regularly be recovered from high quality but very low DNA quantity specimens. Further, the detailed empirical data we provide on the amplification success rates across a range of DNA input quantities will be useful moving forward as PCR-based strategies for mtDNA enrichment are considered for targeted next-generation sequencing workflows.

Characterization of a modified amplification approach for improved STR recovery from severely degraded skeletal elements.

Degraded skeletal remains generally contain limited quantities of genetic material and thus DNA-based identification efforts often target the mitochondrial DNA (mtDNA) control region due to the relative abundance of intact mtDNA as compared to nuclear DNA. In many missing person cases, however, the discriminatory power of mtDNA is inadequate to permit identification when associated anthropological, odontological, or contextual evidence is also limited, and/or the event involves a large number of individuals. In situations such as these, more aggressive amplification protocols which can permit recovery of STR data are badly needed as they may represent the last hope for conclusive identification. We have previously demonstrated the potential of a modified Promega PowerPlex 16 amplification strategy for the recovery of autosomal STR data from severely degraded skeletal elements. Here, we further characterize the results obtained under these modified parameters on a variety of sample types including pristine control DNA and representative case work specimens. Not only is the amplification approach evaluated here sensitive to extremely low authentic DNA input quantities (6 pg), but when the method was applied to thirty-one challenging casework specimens, nine or more alleles were reproducibly recovered from 69% of the samples tested. Moreover, when we independently considered bone samples extracted with a protocol that includes complete demineralization of the bone matrix, the percentage of samples yielding nine or more reproducible alleles increased to 95% with the modified amplification parameters. Overall, direct comparisons between the modified amplification protocol and the standard amplification protocol demonstrated that allele recovery was significantly greater using the aggressive parameters, with only a minimal associated increase in artifactual data.

mtGenome reference population databases and the future of forensic mtDNA analysis.

Mitochondrial DNA (mtDNA) testing in the forensic context requires appropriate, high quality population databases for estimating the rarity of questioned haplotypes. Currently, however, available forensic mtDNA reference databases only include information from the mtDNA control region. While this information is obviously strengthening the foundation upon which current mtDNA identification efforts are based, these data do not adequately prepare the field for recent and rapid advancements in mtDNA typing technologies. Novel tools that quickly and easily permit access to mtDNA coding region data for increased discrimination are now available in the form of single nucleotide polymorphism assays, sequence specific oligonucleotide probes, mass spectrometry instrumentation and next generation sequencing technologies. However, the randomly sampled entire mtGenome reference population data required for statistical interpretation of coding region data are lacking. As a result, in the near future, it seems that routine use of mtDNA coding region data in forensic case work will depend more upon the availability of high-quality entire mtGenome population reference data than the ease with which coding region data can be generated from evidence specimens. Until mtGenome reference databases are available, the utility of novel mtDNA typing technologies and the benefits of recovering mtDNA coding region information from forensic specimens will be limited. Thus, future mtDNA databasing efforts are needed for the development of entire mtDNA genome reference population data suitable for forensic comparisons.

Titanic's unknown child: the critical role of the mitochondrial DNA coding region in a re-identification effort.

This report describes a re-examination of the remains of a young male child recovered in the Northwest Atlantic following the loss of the Royal Mail Ship Titanic in 1912 and buried as an unknown in Halifax, Nova Scotia shortly thereafter. Following exhumation of the grave in 2001, mitochondrial DNA (mtDNA) hypervariable region 1 sequencing and odontological examination of the extremely limited skeletal remains resulted in the identification of the child as Eino Viljami Panula, a 13-month-old Finnish boy. This paper details recent and more extensive mitochondrial genome analyses that indicate the remains are instead most likely those of an English child, Sidney Leslie Goodwin. The case demonstrates the benefit of targeted mtDNA coding region typing in difficult forensic cases, and highlights the need for entire mtDNA sequence databases appropriate for forensic use.

Mitochondrial control region sequences from an African American population sample.

Entire mitochondrial control region data were generated for 248 African American individuals, which had been previously typed for 15 autosomal STRs [J.M. Butler, R. Schoske, P.M. Vallone, J.W. Redman, M.C. Kline, Allele frequencies for 15 autosomal STR loci on U.S. Caucasian, African American, and Hispanic populations, J. Forensic Sci. 48 (2003) 908-911].

The use of mitochondrial DNA single nucleotide polymorphisms to assist in the resolution of three challenging forensic cases.

Mitochondrial DNA (mtDNA) single nucleotide polymorphisms (SNPs) in an 11-plex assay were typed in three missing person cases involving highly degraded human remains. Unlike the traditional forensic approach to analyzing mtDNA which focuses on sequencing portions of the noncoding Control Region, this assay targets discriminatory SNPs that reside principally in the coding region. In two of the cases, the SNP typing successfully excluded one of two reference families that could not be excluded on the basis of mtDNA hypervariable region sequencing alone, and resulted in the final resolution of both decades-old cases. In a third case, SNP typing confirmed the sorting and reassociation of multiple commingled skeletal elements. The application of a specific mtDNA SNP assay in these cases demonstrates its utility in distinguishing samples when the most common Caucasian hypervariable region type is encountered in forensic casework.

Mitochondrial control region sequences from a U.S. "Hispanic" population sample.

Entire mitochondrial control region data was generated for 128 "Hispanics" from the United States. These samples have been previously typed for 15 autosomal STRs [J.M. Butler, R. Schoske, P.M. Vallone, J.W. Redman, M.C. Kline, Allele frequencies for 15 autosomal STR loci on U.S. Caucasian, African American, and Hispanic populations, J. Forensic Sci. 48 (2003) 908-911]. High-throughput robotics, a redundant sequencing approach, and several quality control checks were implemented to generate a high-quality database. The data presented here will augment Hispanic reference data available for forensic mtDNA comparisons.