Analysis of mixtures using next generation sequencing of mitochondrial DNA hypervariable regions

Aim To apply massively parallel and clonal sequencing (next generation sequencing or NGS) to the analysis of forensic mixed samples. Methods A duplex polymerase chain reaction (PCR) assay targeting the mitochondrial DNA (mtDNA) hypervariable regions I/II (HVI/HVII) was developed for NGS analysis on the Roche 454 GS Junior instrument. Eight sets of multiplex identifier-tagged 454 fusion primers were used in a combinatorial approach for amplification and deep sequencing of up to 64 samples in parallel. Results This assay was shown to be highly sensitive for sequencing limited DNA amounts ( ~ 100 mtDNA copies) and analyzing contrived and biological mixtures with low level variants ( ~ 1%) as well as “complex” mixtures (≥3 contributors). PCR artifact “hybrid” sequences generated by jumping PCR or template switching were observed at a low level (<2%) in the analysis of mixed samples but could be eliminated by reducing the PCR cycle number. Conclusion This study demonstrates the power of NGS technologies targeting the mtDNA HVI/HVII regions for analysis of challenging forensic samples, such as mixtures and specimens with limited DNA

DNA sequencing by MALDI-TOF MS using alkali cleavage of RNA/DNA chimeras

Approaches developed for sequencing DNA with detection by mass spectrometry use strategies that deviate from the Sanger-type methods. Procedures demonstrated so far used the sequence specificity of RNA endonucleases, as unfortunately equivalent enzymes for DNA do not exist and therefore require transcription of DNA into RNA prior to fragmentation

The Frequency of Heteroplasmy in the HVII Region of mtDNA Differs across Tissue Types and Increases with Age

An immobilized sequence-specific oligonucleotide (SSO) probe system consisting of 16 SSO probes that detect sequence polymorphisms within five regions of the mtDNA control region was used to investigate the frequency of heteroplasmy in human mtDNA. Five regions of hypervariable region II (HVII) of the control region were studied in blood-, muscle-, heart-, and brain-tissue samples collected from 43 individuals during autopsy. An initial search for heteroplasmy was conducted by use of the SSO probe system. Samples in which multiple probe signals were detected within a region were sequenced for the HVII region, to verify the typing-strip results. The frequency of heteroplasmy was 5 of 43 individuals, or 11.6%. The frequency of heteroplasmy differed across tissue types, being higher in muscle tissue. The difference in the frequency of heteroplasmy across different age groups was statistically significant, which suggests that heteroplasmy increases with age. As a test for contamination and to confirm heteroplasmy, the samples were sequenced for the HVI region and were typed by use of a panel of five polymorphic nuclear markers. Portions of the tissues that appeared to be heteroplasmic were extracted at least one additional time; all gave identical results. The results from these tests indicate that the multiple sequences present in individual samples result from heteroplasmy and not from contamination

Correction: Shelly Y. Shih; et al.; Applications of Probe Capture Enrichment Next Generation Sequencing for Whole Mitochondrial Genome and 426 Nuclear SNPs for Forensically Challenging Samples. Genes 2018, 9, 49

The authors wish to make the following change to their paper [1][...

Applications of Probe Capture Enrichment Next Generation Sequencing for Whole Mitochondrial Genome and 426 Nuclear SNPs for Forensically Challenging Samples

The application of next generation sequencing (NGS) for the analysis of mitochondrial (mt) DNA, short tandem repeats (STRs), and single nucleotide polymorphism (SNPs) has demonstrated great promise for challenging forensic specimens, such as degraded, limited, and mixed samples. Target enrichment using probe capture rather than PCR amplification offers advantages for analysis of degraded DNA since two intact PCR primer sites in the template DNA molecule are not required. Furthermore, NGS software programs can help remove PCR duplicates to determine initial template copy numbers of a shotgun library. Moreover, the same shotgun library prepared from a limited DNA source can be enriched for mtDNA as well as nuclear markers by hybrid capture with the relevant probe panels. Here, we demonstrate the use of this strategy in the analysis of limited and mock degraded samples using our custom probe capture panels for massively parallel sequencing of the whole mtgenome and 426 SNP markers. We also applied the mtgenome capture panel in a mixed sample and analyzed using both phylogenetic and variant frequency based bioinformatics tools to resolve the minor and major contributors. Finally, the results obtained on individual telogen hairs demonstrate the potential of probe capture NGS analysis for both mtDNA and nuclear SNPs for challenging forensic specimens

A twin study of mitochondrial DNA polymorphisms shows that heteroplasmy at multiple sites is associated with mtDNA variant 16093 but not with zygosity.

The mitochondrial theory of ageing proposes that damage to mitochondria and diminished mitochondrial DNA (mtDNA) repair are major contributors to cellular dysfunction and age-related diseases. We investigate the prevalence of heteroplasmy in the mtDNA control region in buccal swab and blood derived samples for 178 women from the TwinsUK cohort (41 DZ pair 39 MZ pairs, 18 singletons, mean age 57.5 range 28-82) and its relationship to age, BMI and fasting insulin and glucose serum levels. The overall estimated prevalence of heteroplasmy for both tissues in the control region measured for 37 sites was 17%. The prevalence of heteroplasmy was higher among the older half of the study subjects than in the younger half (23% vs 10% p<0.03), primarily reflecting the increase in the prevalence of a heteroplasmic dinucleotide CA repeat in variable region II (VRII) with age. The VRII 523-524 heteroplasmic site (heteroplasmic in 25 subjects) was also associated with a decrease in BMI. In addition, concordance rates for common heteroplasmy were observed to be near complete for both dizygotic (DZ = 94%) and monozygotic twin pairs (MZ = 100%), consistent with previous reports that suggest variation in heteroplasmy rates between generations are determined by bottlenecks in maternal transmission of mitochondria. Differences in the prevalence of heteroplasmy were observed overall between samples derived from buccal swabs (19%) and blood (15%, p<0.04). These were particularly marked at position 16093 of hypervariable region I (HVI, 7% vs 0%, respectively, p<4×10(-11)). The presence of the C allele at position 16093 in blood was associated with the presence of heteroplasmy in buccal swabs at this position (p = 3.5×10(-14)) and also at VRII (p = 2×10(-4)) suggesting a possible predisposing role for this site in the accumulation of heteroplasmy. Our data indicate that BMI is potentially associated with control region heteroplasmy

Spectra of the 7-mer region of the reverse GTP trace of HVI of sample C004 (a) and C018 (b) Cambridge Reference Sequence

<p><b>Copyright information:</b></p><p>Taken from "DNA sequencing by MALDI-TOF MS using alkali cleavage of RNA/DNA chimeras"</p><p></p><p>Nucleic Acids Research 2007;35(8):e62-e62.</p><p>Published online 10 Apr 2007</p><p>PMCID:PMC1885642.</p><p>© 2007 The Author(s)</p> A fragment at 2191.4 Da which corresponds to TACTATG disappeared. These spectra were recorded on a Bruker Ultraflex II in reflectron mode which achieves isotopic resolution in this mass range

Spectra of the 9-mer region of the reverse CTP trace of HVI of sample C004 (a) and C018 (b) Cambridge Reference Sequence

<p><b>Copyright information:</b></p><p>Taken from "DNA sequencing by MALDI-TOF MS using alkali cleavage of RNA/DNA chimeras"</p><p></p><p>Nucleic Acids Research 2007;35(8):e62-e62.</p><p>Published online 10 Apr 2007</p><p>PMCID:PMC1885642.</p><p>© 2007 The Author(s)</p> The polymorphism 16311 of sample C004 results in a shift of the peak at 2799.8 Da (TTTATGTAC) to 2815.8 Da (TTTATGTGC) due to an A to G base change in the fragment. These spectra were recorded on a Bruker Autoflex in linear mode