Forensic mitochondrial coding region analysis for increased discrimination using pyrosequencing technology

Pyrosequencing: Current forensic methodology and future applications-a review

The recent development of small, single-amplicon-based benchtop systems for pyrosequencing has opened up a host of novel procedures for applications in forensic science. Pyrosequencing is a sequencing by synthesis technique, based on chemiluminescent inorganic pyrophosphate detection. This review explains the pyrosequencing workflow and illustrates the step-by-step chemistry, followed by a description of the assay design and factors to keep in mind for an exemplary assay. Existing and potential forensic applications are highlighted using this technology. Current applications include identifying species, identifying bodily fluids, and determining smoking status. We also review progress in potential applications for the future, including research on distinguishing monozygotic twins, detecting alcohol and drug abuse, and other phenotypic characteristics such as diet and body mass index. Overall, the versatility of the pyrosequencing technologies renders it a useful tool in forensic genomics.

Developmental validation of SpeID: A pyrosequencing-based assay for species identification

In crime scenes, biological exhibits are often human in origin, yet biological stains from other fauna may also be present at a crime scene, creating confusion during an investigation. Furthermore, identifying the source of a biological sample can be critical during an investigation. To identify the presence of biological material from non-human sources, it is common to use genetic markers within mitochondrial DNA such as cytochrome b, 16S rRNA, and 12S rRNA genes. This process usually requires DNA sequencing, a process that is neither quick nor easy. In general, a faster, more standardized method for species identification from tissue and body fluids is desirable.For this reason, we have developed a vertebrate specific real-time quantitation method that is followed by an automated pyrosequencing-based procedure that sequences a short fragment within the 12S rRNA gene. Using no more than 35 bases, the assay can distinguish between 32 different species commonly found in and around a household with a turnaround time of 6 h from extraction to sequencing. -Using this procedure, up to 48 samples can be run at a time without the need for expensive reagents or bioinformatic skills.

Forensic Analysis of Mitochondrial and Autosomal Markers Using Pyrosequencing®

Forensic casework analyses often face challenges, such as limited genetic material with or without fragmentation and damage. To compensate for low amounts and degradation, shorter amplicons are often applied in the analysis. Also, a change of markers might be necessary using mitochondrial instead of autosomal markers. In addition, forensic research often involves analysis of large number of samples for marker evaluation and population-database compilation. Therefore, a flexible, robust but also rapid method for the detection of variation is highly useful. Pyrosequencing(®) is a rapid, reliable, easy-to-use method for sequence analysis. The method is well suited for rapid forensic analysis of a few targets or analysis of a single target in many samples. It allows sequencing of very short amplicons, which facilitates analysis of degraded DNA. Here we present the use of Pyrosequencing, a robust method for sensitive forensic analysis of mitochondrial DNA, autosomal STRs, and Y-chromosome STRs and SNPs.

Helena, the hidden beauty: Resolving the most common West Eurasian mtDNA control region haplotype by massively parallel sequencing an Italian population sample

The analysis of mitochondrial (mt)DNA is a powerful tool in forensic genetics when nuclear markers fail to give results or maternal relatedness is investigated. The mtDNA control region (CR) contains highly condensed variation and is therefore routinely typed. Some samples exhibit an identical haplotype in this restricted range. Thus, they convey only weak evidence in forensic queries and limited phylogenetic information. However, a CR match does not imply that also the mtDNA coding regions are identical or samples belong to the same phylogenetic lineage. This is especially the case for the most frequent West Eurasian CR haplotype 263G 315.1C 16519C, which is observed in various clades within haplogroup H and occurs at a frequency of 3-4% in many European populations. In this study, we investigated the power of massively parallel complete mtGenome sequencing in 29 Italian samples displaying the most common West Eurasian CR haplotype - and found an unexpected high diversity. Twenty-eight different haplotypes falling into 19 described sub-clades of haplogroup H were revealed in the samples with identical CR sequences. This study demonstrates the benefit of complete mtGenome sequencing for forensic applications to enforce maximum discrimination, more comprehensive heteroplasmy detection, as well as highest phylogenetic resolution.

Development of a multiplex PCR system of 59 mitochondrial SNPs and genetic analysis in Chinese population

The analysis of SNPs located on the mitochondrial DNA can provide information on maternal genetics. In the present study, a set of 59 SNPs were detected simultaneously using three multiplex allele-specific PCR and subsequent CE. Allele-specific primers were designed with different sizes to allow for specifically amplified paired alleles in the same reaction. An allelic ladder based on reference alleles was also created to maintain high-quality analysis standard. Samples from 400 unrelated individuals (200 of Han population and 200 of Uyghur population, China) were successfully analyzed and assigned into 106 relevant haplotypes, resulting in a discrimination power of 98.5%. The haplotype diversity was 0.978 for Han and 0.972 for Uyghur, respectively. Pairwise comparison of haplotype frequency distributions showed significant difference across ethnicities. These results suggest that the 59-SNP PCR system is a reliable, rapid, and economical method for large-scale screening of mitochondrial DNA variation, adding a new aspect for forensic individual identification.

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.

Excap: maximization of haplotypic diversity of linked markers

Genetic markers, defined as variable regions of DNA, can be utilized for distinguishing individuals or populations. As long as markers are independent, it is easy to combine the information they provide. For nonrecombinant sequences like mtDNA, choosing the right set of markers for forensic applications can be difficult and requires careful consideration. In particular, one wants to maximize the utility of the markers. Until now, this has mainly been done by hand.

We propose an algorithm that finds the most informative subset of a set of markers. The algorithm uses a depth first search combined with a branch-and-bound approach. Since the worst case complexity is exponential, we also propose some data-reduction techniques and a heuristic.

We implemented the algorithm and applied it to two forensic caseworks using mitochondrial DNA, which resulted in marker sets with significantly improved haplotypic diversity compared to previous suggestions. Additionally, we evaluated the quality of the estimation with an artificial dataset of mtDNA. The heuristic is shown to provide extensive speedup at little cost in accuracy.

Second generation sequencing allows for mtDNA mixture deconvolution and high resolution detection of heteroplasmy

Aim

To use parallel array pyrosequencing to deconvolute mixtures of mitochondrial DNA (mtDNA) sequence and provide high resolution analysis of mtDNA heteroplasmy.

Methods

The hypervariable segment 1 (HV1) of the mtDNA control region was analyzed from 30 individuals using the 454 GS Junior instrument. Mock mixtures were used to evaluate the system’s ability to deconvolute mixtures and to reliably detect heteroplasmy, including heteroplasmic differences between 5 family members of the same maternal lineage. Amplicon sequencing was performed on polymerase chain reaction (PCR) products generated with primers that included multiplex identifiers (MID) and adaptors for pyrosequencing. Data analysis was performed using NextGENe^® software. The analysis of an autosomal short tandem repeat (STR) locus (D18S51) and a Y-STR locus (DYS389 I/II) was performed simultaneously with a portion of HV1 to illustrate that multiplexing can encompass different markers of forensic interest.

Results

Mixtures, including heteroplasmic variants, can be detected routinely down to a component ratio of 1:250 (20 minor variant copies with a coverage rate of 5000 sequences) and can be readily detected down to 1:1000 (0.1%) with expanded coverage. Amplicon sequences from D18S51, DYS389 I/II, and the second half of HV1 were successfully partitioned and analyzed.

Conclusions

The ability to routinely deconvolute mtDNA mixtures down to a level of 1:250 allows for high resolution analysis of mtDNA heteroplasmy, and for differentiation of individuals from the same maternal lineage. The pyrosequencing approach results in poor resolution of homopolymeric sequences, and PCR/sequencing artifacts require a filtering mechanism similar to that for STR stutter and spectral bleed through. In addition, chimeric sequences from jumping PCR must be addressed to make the method operational.

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.

Analysis of the putative remains of a European patron saint--St. Birgitta

Saint Birgitta (Saint Bridget of Sweden) lived between 1303 and 1373 and was designated one of Europe's six patron saints by the Pope in 1999. According to legend, the skulls of St. Birgitta and her daughter Katarina are maintained in a relic shrine in Vadstena abbey, mid Sweden. The origin of the two skulls was assessed first by analysis of mitochondrial DNA (mtDNA) to confirm a maternal relationship. The results of this analysis displayed several differences between the two individuals, thus supporting an interpretation of the two skulls not being individuals that are maternally related. Because the efficiency of PCR amplification and quantity of DNA suggested a different amount of degradation and possibly a very different age for each of the skulls, an orthogonal procedure, radiocarbon dating, was performed. The radiocarbon dating results suggest an age difference of at least 200 years and neither of the dating results coincides with the period St. Birgitta or her daughter Katarina lived. The relic, thought to originate from St. Birgitta, has an age corresponding to the 13^th century (1215–1270 cal AD, 2σ confidence), which is older than expected. Thus, the two different analyses are consistent in questioning the authenticity of either of the human skulls maintained in the Vadstena relic shrine being that of St. Birgitta. Of course there are limitations when interpreting the data of any ancient biological materials and these must be considered for a final decision on the authenticity of the remains.