Concordance and reproducibility of a next generation mtGenome sequencing method for high-quality samples using the Illumina MiSeq

The effect of library preparation protocol on the efficiency of heteroplasmy detection in mitochondrial DNA using two massively parallel sequencing Illumina systems

Massively parallel sequencing (MPS) technology has become the gold standard in mitochondrial DNA research due to its high sensitivity in detecting mtDNA heteroplasmy, a prognostic marker in various medical applications. Various MPS technologies and platforms used for mtDNA analysis exist. Obtaining reliable and sensitive results requires deep and uniform coverage of the entire mtDNA sequence, which is heavily influenced by the choice of library preparation method and sequencing platform. Here, we present a comparison of the sequencing coverage and the ability to heteroplasmy detection using two library preparation protocols (Nextera XT DNA Library Preparation Kit and Nextera DNA Flex Library Preparation Kit) and two different (MiSeq FGx and ISeq 100) Illumina MPS platforms. Our study indicates that the Nextera DNA Flex Library protocol provides a more balanced coverage along the mitogenome and a reliable heteroplasmy detection with both MiSeq and iSeq Illumina MPS systems.

Capability of the iSeq 100 sequencing system from Illumina to detect low-level substitutions in the human mitochondrial genome

The significance of mtDNA heteroplasmy in forensic and medical genetics has increased recently because massively parallel sequencing (MPS) technologies enable more accurate and precise detection of minority nucleotide variants. Recent reports have shown that detection of low-level substitutions may depend on library preparation or sequencing protocol, and can vary for different MPS platforms. The MiSeq (Illumina) and Ion S5 (Thermo Fisher Scientific) are mainly used for heteroplasmy detection, but no data are available regarding the iSeq 100, an Illumina platform of the smallest throughput. Notably, unlike the other systems, the machine utilizes sequencing by synthesis one-channel chemistry to determine DNA sequences. Thus, it is important to verify the capability of the iSeq 100 system to determine mitochondrial haplotypes and detect heteroplasmic substitutions. In this study, previously determined entire mitochondrial genomes were sequenced with the iSeq 100 system. Each mitogenome was sequenced twice, giving approximately 2000x and 10,000x coverage. All homoplasmic mutations and minority variants above the 19 % level detected with the iSeq 100 system were also observed after dideoxy sequencing. Moreover, all heteroplasmic substitutions above the 2 % level were consistently detected with SBS one-channel chemistry. However, detection of low-level mtDNA variants may require additional, confirmatory experiments. In summary, the iSeq 100 system enables reproducible and accurate sequencing of human mitochondrial genomes. Detection of mtDNA minority variants depends on the laboratory protocol and sequencing platform used, but homoplasmic mutations and heteroplasmy above the 2 % level can be correctly detected with the iSeq 100 system.

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

Low-level point heteroplasmy detection in human mitogenomes amplified with different polymerases and sequenced on MiSeq FGx platform

Introduction

Massively parallel sequencing of mitogenomes usually requires prior amplification. The PCR step may influence the quality of the data obtained, especially when low-level heteroplasmy detection is applied.

Aim

The aim of this study was to compare the reliability of two different DNA polymerases in detecting homoplasmic and heteroplasmic substitutions in human mitogenomes.

Material and methods

Mitogenomes of five samples were amplified with Long PCR Enzyme Mix from Fermentas or TaKaRa LA Taq DNA Polymerase from TaKaRa. Then, NexteraTM XT DNA libraries were sequenced on MiSeq FGx platform (Illumina). mtDNA substitutions were called for alternative variants above the 1% level.

Results

All homoplasmic substitutions detected in amplicons generated with polymerases studied here and sequenced on MiSeq FGx system were consistently identified as homoplasmies with alternative sequencing methods. TaKaRa LA Taq DNA Polymerase was found to be less accurate in low-level heteroplasmy detection than Long PCR Enzyme Mix enzyme as more false negative and false positive results were observed for minority variants called above the 1% level. Nevertheless, both PCR systems studied can be successfully used to detect authentic mtDNA substitutions, for which minority variants exceed the 3.61% level assuming at least 10,000x coverage and sequencing Nextera XT DNA libraries on MiSeq FGx machine.

Conclusions

The accuracy and sensitivity of point heteroplasmy detection with the MiSeq FGx instrument varies on polymerase used for mtDNA amplification. Therefore, it is recommended to validate the laboratory protocols used for mtDNA substitution detection prior to their implementation for the forensic or medical genetics purposes.

Plasticity of Mature B Cells Between Follicular and Classic Hodgkin Lymphomas: A Series of 22 Cases Expanding the Spectrum of Transdifferentiation

Follicular lymphoma and classic Hodgkin lymphoma can be associated in composite and/or sequential lymphomas. Common IGH and BCL2 rearrangements have already been identified between both contingents of these entities, but mutation profiles have not yet been investigated. The main objective of this study was to analyze the transdifferentiation process that may occur between Hodgkin and follicular contingents in sequential and composite lymphomas to better characterize these entities. From 2004 to 2020, a retrospective multicentric study was performed, including 9 composite and 13 sequential lymphomas. Clinical data were retrospectively collected. Fluorescent in situ hybridization of BCL2 and BCL6 rearrangements, polymerase chain reaction of IGH and IGK rearrangements, next-generation sequencing of IGK rearrangement, and targeted next-generation sequencing (TNGS) on a panel of genes frequently mutated in lymphomas were performed on each contingent of composite and sequential lymphomas. For TNGS, each contingent was isolated by laser capture microdissection. Clinical presentation and evolution were more aggressive in sequential than composite lymphomas. By fluorescent in situ hybridization, common rearrangements of BCL6 and BCL2 were identified between both contingents. Similarly, a common clonal relationship was established by evaluating IGH and IGK rearrangement by polymerase chain reaction or next-generation sequencing. By TNGS, the same pathogenic variants were identified in both contingents in the following genes: CREBBP, KMT2D, BCL2, EP300, SF3B1, SOCS1, ARID1A, and BCOR. Specific pathogenic variants for each contingent were also identified: XPO1 for Hodgkin lymphoma contingent and FOXO1, TNFRSF14 for follicular lymphoma contingent. This study reinforces the hypothesis of a transdifferentiation process between Hodgkin and follicular contingent of sequential/composite lymphomas.

Human Mitochondrial Control Region and mtGenome: Design and Forensic Validation of NGS Multiplexes, Sequencing and Analytical Software

Forensic mitochondrial DNA (mtDNA) analysis conducted using next-generation sequencing (NGS), also known as massively parallel sequencing (MPS), as compared to Sanger-type sequencing brings modern advantages, such as deep coverage per base (herein referred to as read depth per base pair (bp)), simultaneous sequencing of multiple samples (libraries) and increased operational efficiencies. This report describes the design and developmental validation, according to forensic quality assurance standards, of end-to-end workflows for two multiplexes, comprised of ForenSeq mtDNA control region and mtDNA whole-genome kits the MiSeq FGx^TM instrument and ForenSeq universal analysis software (UAS) 2.0/2.1. Polymerase chain reaction (PCR) enrichment and a tiled amplicon approach target small, overlapping amplicons (60–150 bp and 60–209 bp for the control region and mtGenome, respectively). The system provides convenient access to data files that can be used outside of the UAS if desired. Studies assessed a range of environmental and situational variables, including but not limited to buccal samples, rootless hairs, dental and skeletal remains, concordance of control region typing between the two multiplexes and as compared to orthogonal data, assorted sensitivity studies, two-person DNA mixtures and PCR-based performance testing. Limitations of the system and implementation considerations are discussed. Data indicated that the two mtDNA multiplexes, MiSeq FGx and ForenSeq software, meet or exceed forensic DNA quality assurance (QA) guidelines with robust, reproducible performance on samples of various quantities and qualities.

Assessment of Illumina® Human mtDNA Genome assay: workflow evaluation with development of analysis and interpretation guidelines

Mitochondrial DNA (mtDNA) is a small but significant part of the human genome, whose applicability potential has gradually increased with the advent of massively parallel sequencing (MPS) technology. Knowledge of the particular workflow, equipment, and reagents used, along with extensive usage of negative controls to monitor all preparation steps constitute the prerequisites for confident reporting of results. In this study, we performed an assessment of Illumina® Human mtDNA Genome assay on MiSeq FGx™ instrument. Through analysis of several types of negative controls, as well as mtDNA positive controls, we established thresholds for data analysis and interpretation, consisting of several components: minimum read depth (220 reads), minimum quality score (41), percentage of minor allele sufficient for analysis (3.0%), percentage of minor allele sufficient for interpretation (6.0%), and percentage of major allele sufficient for homoplasmic variant call (97.0%). Based on these criteria, we defined internal guidelines for analysis and interpretation of mtDNA results obtained by MPS. Our study shows that the whole mtDNA assay on MiSeq FGx™ produces repeatable and reproducible results, independent of the analyst, which are also concordant with Sanger-type sequencing results for mtDNA control region, as well as with MPS results produced by NextSeq®. Overall, established thresholds and interpretation guidelines were successfully applied for the sequencing of complete mitochondrial genomes from high-quality samples. The underlying principles and proposed methodology on the definition of internal laboratory guidelines for analysis and interpretation of MPS results may be applicable to similar MPS workflows, e.g. targeting good-quality samples in forensic genetics and molecular diagnostics.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

Massive parallel sequencing in forensics: advantages, issues, technicalities, and prospects

In the last decade, next-generation sequencing (NGS) technology, alternatively massive parallel sequencing (MPS), was applied to all fields of biological research. Its introduction to the field of forensics was slower, mainly due to lack of accredited sequencers, kits, and relatively higher sequencing error rates as compared with standardized Sanger sequencing. Currently, a majority of the problematic issues have been solved, which is proven by the body of reports in the literature. Here, we discuss the utility of NGS sequencing in forensics, emphasizing the advantages, issues, the technical aspects of the experiments, commercial solutions, and the potentially interesting applications of MPS.

Characterization of background noise in MiSeq MPS data when sequencing human mitochondrial DNA from various sample sources and library preparation methods

Improved resolution of massively parallel sequencing (MPS) allows for the characterization of mitochondrial (mt) DNA heteroplasmy to levels previously unattainable with traditional sequencing approaches. An essential criterion for the reporting of heteroplasmy is the ability of the MPS method to distinguish minor sequence variants (MSVs) from system noise, or error. Therefore, an assessment of the background noise in the MPS method is desirable to identify the point at which reliable data can be reported. Substitution and sequence specific error (SSE) was evaluated for a variety of sample types and two library preparations. Substitution error rates ranged from 0.18 to 0.49 per 100 nucleotides with C positions generally having the highest rate of misincorporation. Comparison of error rates across sample types indicated a significant increase for samples with damaged DNA. The positions of error were varied across datasets (pairwise concordance 0-68%), but had greater consistency within the damaged samples (80-96%). The most commonly observed motif preceding error in forward reads was CCG, while GGT was most common in reverse reads, both consistent with previous findings. The findings illustrate that for datasets containing samples with damaged DNA, reporting thresholds for heteroplasmy may have to be modified and individual sites with error levels exceeding thresholds should be scrutinized. Collectively, the shifting error profiles observed across the various sample types and library preparation methods demonstrates the need for an assessment of error under these varying circumstances. Characterization of the applicable background noise will help to ensure that thresholds are reliably set for detection of true MSVs.

Validation of a 52-mtSNP minisequencing panel for haplogroup classification of forensic DNA samples

Mitochondrial DNA (mtDNA) is a useful tool in forensic investigation as it provides information about the matrilineal ancestry of individuals. In addition, mtDNA can be analyzed when the analysis of other nuclear markers is underperforming. Recently, we developed a minisequencing panel for the simultaneous analysis of 52 mtDNA SNPs to classify maternal lineages into the main haplogroups and their phylogeographic origin. In order to make this panel suitable for forensic genetics laboratories, a validation study has been performed in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines, including species specificity, reproducibility, sensitivity, and stability tests. The results demonstrate that the panel of 52 mtDNA SNPs is highly sensitive, since it enables to obtain complete genetic profiles of samples containing minimal amounts of DNA (1 pg). Furthermore, it provides sufficient genetic information to detect the matrilineal biogeographical origin of highly degraded samples, i.e., ancient dating skeletal remains, and samples with the presence of inhibitors, such as hematin and humic acid. In addition, this panel can detect mixtures in samples whose mtDNA haplogroups of contributors are different. Overall, the results of this study demonstrate the suitability of this minisequencing panel of 52 mtDNA SNPs to be used in forensic cases, with samples of low amount or degraded DNA.

Comparison of whole mitochondrial genome variants between hair shafts and reference samples using massively parallel sequencing

Hair shafts are one of the most common types of evidence at crime scenes, and mitochondrial DNA (mtDNA) has been analyzed as a valuable genetic marker for hair shafts in forensic casework. However, the mtDNA analysis strategy may vary according to the quantity and quality of DNA extracted from a forensic sample and the available massively parallel sequencing (MPS) platform in laboratories. Forensic practitioners often have to interpret mtDNA sequences exhibiting point heteroplasmy (PHP) that are analyzed using different analytical methods. In the present study, the whole mitochondrial genome (mtGenome) variants of hair shaft samples obtained from 20 donors, which were sampled in duplicate and stored at room temperature for > 1 year, were analyzed using the Precision ID mtDNA Whole Genome Panel and Ion S5 system. The whole mtGenome variants of 20 blood and 20 buccal swab samples (reference samples) from the hair shaft donors were analyzed using the Nextera XT DNA Library Prep Kit and MiSeq System. A total of 20 unique mtGenome haplotypes were observed, and 56 PHP variants were identified across the 4 sets of tissue. When the major nucleotide of PHP was considered, 16 of 20 haplotypes of the hair shaft samples matched those of the corresponding blood and buccal swab samples. In four donors, the major nucleotide of PHP was inverted at one nucleotide position between the hair shaft and reference samples. However, the data obtained on MPS, showing high PHP resolution, provided substantial information to avoid false exclusion when comparing two haplotypes containing PHP with inverted major nucleotides. In conclusion, the present study demonstrates the utility of MPS in forensic casework in the comparative analysis of mtGenome variants containing PHP.

Estimation of the number of contributors to mixed samples of DNA by mitochondrial DNA analyses using massively parallel sequencing

We evaluated whether the number of contributors to mixed DNA samples can be estimated by analyzing the D-loop of mitochondrial DNA using massively parallel sequencing. The A- (positions 16,209-16,400) and B- (positions 30-284) amplicons in hypervariable regions 1 and 2, respectively, were sequenced using MiSeq with 2 × 251 cycles. Sequence extraction and trimming were performed using CLC Genomics Workbench 11 and the number of observed haplotypes was counted for each amplicon type using Microsoft Excel. The haplotype ratios were calculated by dividing the number of counted reads of the corresponding haplotype by the total number of sequence reads. Haplotypes that were over the threshold (5%) were defined as positive haplotypes. The number of larger positive haplotypes in either of the two amplicon types was defined as the number of contributors. Samples were collected from seven individuals. Seventeen mixed samples were prepared by mixing DNA from two to five contributors at various ratios. The number of contributors was correctly estimated from almost all of the mixed samples containing equal amounts of DNA from two to five people. In mixed samples of two or three people, the minor components were detected down to a ratio of 20:1 or 8:2:1. However, heteroplasmy, base deletions, and sharing of the same haplotypes caused incorrect estimations of the number of contributors. Although this method still has room for improvement, it may be useful for estimating the number of contributors in a mixed sample, as it does not rely on forensic mathematics.

Identification of the perpetrator among identical twins using next-generation sequencing technology: A case report

BACKGROUND

Monozygotic (MZ) twins, considered genetically identical, cannot be distinguished using regular short tandem repeats (STR) typing, thus presenting a challenge for forensic geneticists. In paternity testing, single nucleotide polymorphisms (SNPs) in nuclear DNA can help distinguish MZ twins. However, the unique features of the mitochondrial genome, such as high copy number, small genome size, and high substitution rate, make it a promising source for applications in forensic science.

METHODS

Whole-genome sequencing (WGS) was performed on blood samples, and bioinformatic analysis was used to distinguish between MZ twins. Amplification refractory mutation system polymerase chain reaction (ARMS-PCR) was used to confirm the WGS results. This methodology was further applied to forensic samples from criminal cases. Amplicon sequencing was also performed to further exclude the innocent twin.

RESULTS

The monozygosity of the twins was confirmed using STR typing. Only one potential somatic mutation, m.6903 T > C (2.6%), in the mitochondrial genome of one of the twins was verified when the sequence depth was set to 2000-fold, while no other distinguishing locus in the nuclear genome was identified. By dividing the number of C-reads by total reads, WGS data confirmed the amount of the minor component C to be 2.6%, which was further confirmed by ARMS-PCR. In addition, the heterogeneous locus was used to identify samples obtained from four criminal cases for forensic testing. Two heterogeneous loci in the sperm DNA of the other twin were identified by amplicon sequencing, and the amount of minor component T in m.6935C > T and m.6938C > T was estimated to be 17.91% and 18.79%, respectively.

CONCLUSION

The biological samples taken from the MZ twins were distinguished using a combination of WGS, allele-specific PCR, and deep-amplicon sequencing. Compared with nuclear DNA, mitochondrial DNA exhibited a higher potential for distinguishing between the MZ twins. The distinguishing feature of the mitochondrial DNA was the heterogeneous SNPs that occurred in only one twin. One SNP was further verified in the samples from the criminal cases and helped identify the perpetrator in case 1 and case 2. Furthermore, two heterogeneous SNPs found by amplicon sequencing helped to exclude the innocent twin in all four cases. Our findings demonstrated that a combination of deep sequencing and molecular analysis can be an effective way to distinguish between identical twins and can be used to analyze samples from criminal cases.

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.

Mitochondrial DNA in human identification: a review

Mitochondrial DNA (mtDNA) presents several characteristics useful for forensic studies, especially related to the lack of recombination, to a high copy number, and to matrilineal inheritance. mtDNA typing based on sequences of the control region or full genomic sequences analysis is used to analyze a variety of forensic samples such as old bones, teeth and hair, as well as other biological samples where the DNA content is low. Evaluation and reporting of the results requires careful consideration of biological issues as well as other issues such as nomenclature and reference population databases. In this work we review mitochondrial DNA profiling methods used for human identification and present their use in the main cases of humanidentification focusing on the most relevant issues for forensics.

Recovery of mtDNA from unfired metallic ammunition components with an assessment of sequence profile quality and DNA damage through MPS analysis

Recovery of suitable amounts of quality DNA from copper and brass surfaces, like those encountered in ammunition, has been a challenge for the forensic community. The ability of copper ions to rapidly facilitate oxidative damage leading to fragmentation of DNA significantly reduces the pool of templates for PCR amplification. We compared two methods for recovering mitochondrial (mt) DNA from the surface of unfired copper projectiles, brass casings, and aluminum casings, and found that using a cotton swab moistened with 0.5M EDTA was the favored approach, especially when the metallic surface was etched. Degradation was significantly higher for DNA samples recovered from copper and brass surfaces, when compared to aluminum. Massively parallel sequencing (MPS) of the control region, using the PowerSeq™ CRM Nested System kit and the Illumina MiSeq instrument, produced full haplotypes for aluminum samples regardless of the method used to deposit or collect DNA, while less than 60% of the copper and brass samples produced partial or full profile information. Touch DNA collected from copper and brass samples produced higher rates of partial or full MPS profile information (∼88-96%), while collection with 0.5M EDTA produced better results than when collection was performed with water; average of ∼70% versus ∼47%. While MPS data was not impacted by noise in the sequencing process, a higher than expected rate of noise was observed, potentially due to an increase in low-level damage lesions. Noise patterns were strikingly different when compared to control data, suggesting that noisy sites may be predictable when testing samples with high levels of oxidative damage. Library preparation was a poor predictor of MPS data quality, as a large percentage of reads did not align with the reference genome. This may impact the number of samples that can be run when a deep-coverage MPS approach is being considered for analysis of mtDNA heteroplasmy. Overall, when applying an MPS approach to the analysis of mtDNA recovered from ammunition, results are expected from touch DNA, will be limited for copper and brass components when the DNA is exposed to an aqueous environment, and DNA degradation will be accelerated when DNA comes in contact with copper or brass surfaces. Practitioners should consider collecting DNA from metallic surfaces with 0.5M EDTA, as this will maximize yield and mitigate degradation. The results of this study directly impact MPS analysis of minor mtDNA sequence variants from metallic surfaces, and are particularly relevant to forensic investigations.

Development of a triplex mtDNA qPCR assay to assess quantification, degradation, inhibition, and amplification target copy numbers

A hybrid absolute/relative qPCR assay which provides information regarding the condition of mitochondrial DNA (mtDNA) in a DNA sample is described. MtDNA concentration (copy number/μL) is determined via absolute quantification using a standard curve of a synthetic duplex DNA previously described (Kavlick et al., 2011). The state of mtDNA degradation is determined via the relative quantification of a mtDNA target found within the 16 s rRNA gene which is 3× longer than that of the short target in the former duplex assay, using the delta, delta Ct (ΔΔCt) method. The presence or absence of PCR inhibitors in the sample is qualitatively determined using a custom internal positive control (IPC) system which targets a unique and non-naturally occurring duplex DNA sequence. This IPC effectively detected inhibition by humic acid, tannic acid, melanin, and EDTA. All three assay components utilize sensitive and specific hydrolysis probes. The utility of ΔΔCt method was demonstrated in a series of experiments involving laboratory-fragmented DNA. Also described is a method for estimating copy number of any mtDNA target longer than the two targets amplified. The described triplex assay works well for intact and for fragmented or degraded mtDNA and therefore may be useful in forensic and ancient DNA disciplines as well as in biomedical research or practice.

Recovery of whole mitochondrial genome from compromised samples via multiplex PCR and massively parallel sequencing

In forensic casework, compromised samples often possess limited or degraded nuclear DNA, rendering mitochondrial DNA a more feasible option for forensic DNA analyses. The emergence of massively parallel sequencing (MPS) has enabled the recovery of extensive sequence information from very low quantities of DNA. We have developed a multiplex PCR method that amplifies the complete mitochondrial genome in a range of forensically relevant samples including single cells, cremated remains, bone, maggot and hairs isolated from dust bunnies. Following library preparation, MPS yields complete or nearly complete mitochondrial genome coverage for all samples. To confirm concordance between sample types and between sequencing platforms, we compared sequencing results from hair and buccal swabs from two references. Low initial DNA input into the multiplex PCR allows for conservation of precious DNA while MPS maximizes recovery of genetic information.

Compromised samples are routine in forensic casework, and evidence DNA is often limited and/or degraded. Mitochondrial DNA is often a suitable option for forensic analysis as it is available in multiple copies per cell. We have established a method that amplifies the complete mitochondrial genome in a range of compromised samples and utilizes new high-throughput sequencing technologies to recover maximal genetic information from a small amount of sample. Combination of the aforementioned amplification method and high-throughput sequencing allows the analyst to conserve DNA in forensic cases where DNA is often limited.

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.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.

Evaluation of GeneMarker® HTS for improved alignment of mtDNA MPS data, haplotype determination, and heteroplasmy assessment

Existing software has not allowed for effective alignment of mitochondrial (mt) DNA sequence data generated using a massively parallel sequencing (MPS) approach, combined with the ability to perform a detailed assessment of the data. The regions of sequence that are typically difficult to align are homopolymeric stretches, isolated patterns of SNPs (single nucleotide polymorphisms), and INDELs (insertions/deletions). A custom software solution, GeneMarker^® HTS, was developed and evaluated to address these limitations, and to provide a user-friendly interface for forensic practitioners and others interested in mtDNA analysis of MPS data. GeneMarker^® HTS generates an exportable consensus mtDNA sequence that produces phylogenetically correct SNP and INDEL calls using a customizable motif-based alignment algorithm. Sequence data from 500 individuals, with various alignment asymmetries and levels of heteroplasmy, were used to assess the software. Accuracy in producing mtDNA haplotypes, the ability to correctly identify low-level heteroplasmic sequence variants, and the user-based features of the software were evaluated. Analyzed sequences yielded correct mtDNA haplotypes, and heteroplasmic variants were properly identified with minimal manual interpretation. The software offers numerous user-defined parameters for filtering the data that address the interests of researchers and practitioners, and provides multiple options for viewing and navigating through the data.

MPS analysis of the mtDNA hypervariable regions on the MiSeq with improved enrichment

The non-coding displacement (D) loop of the human mitochondrial (mt) genome contains two hypervariable regions known as HVR1 and HVR2 that are most often analyzed by forensic DNA laboratories. The massively parallel sequencing (MPS) protocol from Illumina (Human mtDNA D-Loop Hypervariable Region protocol) utilizes four sets of established PCR primer pairs for the initial amplification (enrichment) step that span the hypervariable regions. Transposase adapted (TA) sequences are attached to the 5'-end of each primer, allowing for effective library preparation prior to analysis on the MiSeq, and AmpliTaq Gold DNA polymerase is the enzyme recommended for amplification. The amplification conditions were modified by replacing AmpliTaq Gold with TaKaRa Ex Taq® HS, along with an enhanced PCR buffer system. The resulting method was compared to the recommended protocol and to a conventional non-MPS approach used in an operating forensic DNA laboratory. The modified amplification conditions gave equivalent or improved results, including when amplifying low amounts of DNA template from hair shafts which are a routine evidence type in forensic mtDNA cases. Amplification products were successfully sequenced using an MPS approach, addressing sensitivity of library preparation, evaluation of precision and accuracy through repeatability and reproducibility, and mixture studies. These findings provide forensic laboratories with a robust and improved enrichment method as they begin to implement the D-loop protocol from Illumina. Given that Ex Taq® HS is a proofreading enzyme, using this approach should allow for improved analysis of low-level mtDNA heteroplasmy.