Sequencing-induced artefacts in NGS STR data.

Significant progress has been made in recent years in the development of techniques for Next Generation Sequencing (NGS), or Massively Parallel Sequencing (MPS), of forensically relevant short tandem repeat (STR) loci. However, as these technologies are investigated and adopted by forensic laboratories, new challenges unfold that require further scrutiny. In the analysis of DNA profiles generated using the MiSeq FGx sequencing system, we have observed noise sequences with relatively high readcounts that are challenging to distinguish from genuine alleles. These high read count noise sequences appear as allele sequences with one or a few substituted bases compared to a known allele sequence within the profile. An examination of ForenSeq DNA Signature Prep Kit STR noise sequences revealed that the substituted base of a parent allele can align to the same position on the sequence across noise sequences. This suggests that these substitution events occur at specific positions within the amplicon, resulting in multiple noise reads with substitutions at the same position. Mapping of the noise events onto the original raw read positions revealed a high number of events, or "noise spikes", occurring at specific positions within a given sequencing run. These noise spikes affected reads across the entire run, agnostic of locus or sample, while the position, occurrence, and amplitude of the spikes differed across runs. The majority of noise sequences with high read counts in a DNA profile were generated from base changes at these spike positions, and could be classified as "noise spike artefacts". In this paper we present evidence of the noise spike artefacts and their genesis during the sequencing process in the sequencing-by-synthesis (SBS) cycles, as well as the methods developed to detect them. The information and methods will assist laboratories with detecting noise spikes in MiSeq FGx sequencing runs, differentiating authentic allele sequences from noise spike artefacts, and developing protocols for analyst review and handling of MiSeq FGx data.

MPSproto: An extension of EuroForMix to evaluate MPS-STR mixtures.

We have developed MPSproto as an extension of EuroForMix to improve handling of stutter artefacts and other typing errors that commonly occur in MPS-STR data. MPSproto implements two models for read depth: gamma and negative binomial. It differs from EuroForMix in that calibration is required before mixtures are interpreted. In this study a mixture dataset (2-4 persons) was revisited, where EuroForMix interpretation of MPS-STR mixtures using the LUS+ format was first described; the performance of this model was compared to the MPSproto models. Results indicated that, overall, the MPSproto models performed better than the conventional EuroForMix model, and the gamma model implemented in MPSproto performed best. Differences were highlighted and further investigated to establish causality. Goodness of fit tests showed that the MPSproto models were adequate for the sequence reads when a low analytical threshold was applied.

Modeling allelic analyte signals for aSTRs in NGS DNA profiles.

We describe an adaption of Bright et al.'s work modeling peak height variability in CE-DNA profiles to the modeling of allelic aSTR (autosomal short tandem repeats) read counts from NGS-DNA profiles, specifically for profiles generated from the ForenSeq™ DNA Signature Prep Kit, DNA Primer Mix B. Bright et al.'s model consists of three key components within the estimation of total allelic product-template, locus-specific amplification efficiencies, and degradation. In this work, we investigated the two mass parameters-template and locus-specific amplification efficiencies-and used MLE (maximum likelihood estimation) and MCMC (Markov chain Monte Carlo) methods to obtain point estimates to calculate the total allelic product. The expected read counts for alleles were then calculated after proportioning some of the expected stutter product from the total allelic product. Due to preferential amplicon selection introduced by the sample purification beads, degradation is difficult to model from the aSTR outputs alone. Improved modeling of the locus-specific amplification efficiencies may mask the effects of degradation. Whilst this model could be improved by introducing locus specific variances in addition to locus specific priors, our results demonstrate the suitability of adapting Bright et al.'s allele peak height model for NGS-DNA profiles. This model could be incorporated into continuous probabilistic interpretation approaches for mixed DNA profiles.

Variability and additivity of read counts for aSTRs in NGS DNA profiles.

There has been an increase in the number of laboratories and researchers adopting new sequencing technologies, known as next-generation sequencing (NGS). An understanding of the behaviour of NGS DNA profiles is needed to enable for the development of probabilistic genotyping methods for the interpretation of such profiles. In this work, we investigate NGS analyte signal variation, specifically heterozygous balance and stutter variability from profiles generated using the ForenSeq™ DNA Signature Prep Kit, DNA Primer Mix B. We also investigate additivity of analyte signals in NGS profiles for overlapping allelic and stutter signals originating from the same or different contributors. We describe models that can be used to inform a continuous method for the interpretation of DNA profiling data.

An examination of STR nomenclatures, filters and models for MPS mixture interpretation.

The increased interest in the use of Massively Parallel Sequencing (MPS) technologies to type traditional autosomal STR markers raises multiple questions regarding interpretation of the results via probabilistic genotyping. To begin to address some of those questions, we examined the effects of using differing degrees of sequence information, pre-filtering, and data modeling to interpret complex MPS-STR mixtures in a probabilistic genotyping software. Sixty ForenSeq typing results for mixtures of from two to four contributors were: 1) represented using three separate formats that captured different degrees of sequence information, and 2) were analyzed using three different filtering approaches prior to probabilistic interpretation. All mixtures for the different format and filtering variants were subsequently interpreted with respect to ten reference profiles, using both qualitative (LRmix) and quantitative (EuroForMix) models to calculate the likelihood ratio (LR). The LR results indicated moderate information gain when the STR nomenclature was based upon the longest uninterrupted stretch (LUS) compared with conventional capillary electrophoresis repeat units (RU), whereas additional gains were very small when the complete sequence information was utilised. Use of a static analytical threshold for data pre-filtering improved LRs compared to a dynamic (percentage-based) threshold, as the static threshold prevented excessive filtering of alleles originating from minor contributors. For interpretations performed using a quantitative model, a small improvement in performance was observed if a stutter model was employed instead of using stutter thresholds to pre-filter the data, whereas - as expected - performance worsened considerably under the qualitative model when stutter was not pre-filtered. Given the empirical and theoretical findings in this study we discuss the value of utilizing sequence-level information and potential paths forward to increase information gain using MPS systems.

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.

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.

A closer look at Verogen's Forenseq™ DNA Signature Prep kit autosomal and Y-STR data for streamlined analysis of routine reference samples.

Massively parallel sequencing (MPS) provides forensic DNA laboratories an option to overcome the limitations associated with CE and current STR assays. Verogen's MPS ForenseqTM DNA Signature kit concomitantly amplifies 27 autosomal, 7 X-, and 24 Y-STRs. In addition, 94 identity, 56 ancestry, and 22 phenotypic-informative SNPs are included for a total of over 200 markers in one multiplex. An internal validation of this platform was conducted using reference samples to investigate whether the Forenseq™ DNA Signature Prep kit, specifically primer panel B, has the capability to provide consistent and accurate typing/sequencing data. The data presented in this report is limited to that corresponding to autosomal and Y-STRs. Results suggest that the system can consistently generate accurate genotyping data when up to 40 high-quality, high-quantity (i.e. 1 ng) single source samples are pooled into a sequencing reaction. The results generated were used to determine appropriate analysis parameters and thresholds for streamlined data interpretation of reference samples.

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.

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.

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.

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.

Short tandem repeat typing on the 454 platform: strategies and considerations for targeted sequencing of common forensic markers.

To investigate the feasibility of next generation sequencing technology (NGS) for the multiplex detection and sequence production of short tandem repeats (STRs) from degraded and low DNA quantity samples, standard polymerase chain reaction amplification methods were used to enrich for commonly employed STR markers. Samples were amplified with two multiplexing strategies: a multiplex containing thirteen miniSTR markers and a series of multiplexes containing four miniSTR markers each. Each sample multiplex was barcoded with a sample-specific multiplex identifier for subsequent parallel tagged sequencing on the GS Junior System (454 Life Sciences, a Roche company, Branford, CT). Sequencing results from over fifty DNA extracts representing both pristine samples and low-quality evidentiary specimens reflected known genotypes and were consistent across multiple extracts and/or amplifications of the same sample. Furthermore, the NGS data revealed sequence information not available with standard capillary electrophoresis-based detection alone. For the population samples tested, a total of 152 single nucleotide polymorphisms or insertions/deletions were identified in over 935 recovered alleles, averaging one polymorphism for every six recovered alleles. For three of the loci, the sequence information doubled the number of alleles detected via traditional STR typing by fragment analysis. In addition, twenty-eight of these variants were only seen once within our dataset, highlighting the potential for discrimination among individuals. These additional data are likely to be particularly valuable in missing persons and disaster victim identification cases for which only partial profiles may be recovered and/or only distant kin are available as references. And, considering the opportunity to target only small amplicons with NGS, this type of STR typing will allow for greater information recovery from challenging casework samples. While our results highlight the potential of new technologies for recovering discriminatory genetic information from evidentiary specimens, our data also reveal the complexities of NGS-based STR typing, both in terms of the laboratory assays themselves as well as the downstream data processing and analysis.

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.

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.

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.

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.

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.