Factors affecting the detection and quantification of mitochondrial point heteroplasmy using Sanger sequencing and SNaPshot minisequencing

Estimating Copy-Number Proportions: The Comeback of Sanger Sequencing

Determination of the relative copy numbers of mixed molecular species in nucleic acid samples is often the objective of biological experiments, including Single-Nucleotide Polymorphism (SNP), indel and gene copy-number characterization, and quantification of CRISPR-Cas9 base editing, cytosine methylation, and RNA editing. Standard dye-terminator chromatograms are a widely accessible, cost-effective information source from which copy-number proportions can be inferred. However, the rate of incorporation of dye terminators is dependent on the dye type, the adjacent sequence string, and the secondary structure of the sequenced strand. These variable rates complicate inferences and have driven scientists to resort to complex and costly quantification methods. Because these complex methods introduce their own biases, researchers are rethinking whether rectifying distortions in sequencing trace files and using direct sequencing for quantification will enable comparable accurate assessment. Indeed, recent developments in software tools (e.g., TIDE, ICE, EditR, BEEP and BEAT) indicate that quantification based on direct Sanger sequencing is gaining in scientific acceptance. This commentary reviews the common obstacles in quantification and the latest insights and developments relevant to estimating copy-number proportions based on direct Sanger sequencing, concluding that bidirectional sequencing and sophisticated base calling are the keys to identifying and avoiding sequence distortions.

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two. Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

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.

Get it off, but keep it: Efficient cleaning of hair shafts with parallel DNA extraction of the surface stain

The analysis of hair samples is a common task in forensic investigations. Material transferred to the surface of a hair during a crime challenges the analysis as it has to be removed efficiently. However, the removal of the stain can also lead to a loss of information on stain contributors. DNA analysis of the stain itself might thus be helpful for the forensic investigation. The aim of this study was the examination of different methods to remove common biological surface stains completely from human hair shafts without hampering the parallel DNA extraction of the cleaned hair shaft and the isolated surface stain (blood, saliva, vaginal secretion, semen, and skin flocks). Four different methods of cleaning (water, lysis buffer, swabbing, NaClO) were compared to their cleaning efficiency as well as their success of mtDNA analysis of three hair donors and the original five stains on the hair. In order to test the suitability of this procedure for future analysis methods, a selection of samples were also sequenced with MPS. Additionally, nuclear DNA analysis of the stain DNA was performed using a screening STR assay to test the potential success for detection of a STR profile. The most efficient removal of the stain was achieved using NaClO, however compromising further analysis of the stain DNA. The best results for cleaning and parallel stain analysis were obtained using a swab moistened with 0.5 % SDS for surface cleaning. Especially water failed to remove stains efficiently, leading to a high amount of mixed mtDNA in the DNA extracts. MPS showed an increased sensitivity for detection of minute mixtures.

Base specific variation rates at mtDNA positions 16093 and 16183 in human hairs

Small variations between haplotypes detected in different tissues from the same individual have been previously described. These differences complicate the interpretation of mtDNA results in real forensic casework. mtDNA haplotypes detected in hair strands collected at the crime scene have to be frequently compared with haplotypes of reference samples (buccal swabs) from victims or suspects. Nucleotide position 16093 is a well-known hot spot where differences can accumulate between different tissues of the same individual. Intra individual variation was also detected at positions 16182 and 16183 in haplotypes showing an uninterrupted HV1 poly-C stretch (with 16189C). In order to better characterize the type of variation in these positions between buccal cells and hair strands from the same individual, we have performed Sanger sequencing in 25-28 hair strands (411 in total) from 15 individuals showing either an uninterrupted HV1 polyC-stretch (16189C) or 16093C/Y in their buccal cells. The results have been evaluated by also taking into account our previous results published in [19]. We have found that no variation among hair strands was detected in individuals showing T16093 in buccal cells, while variation in hair strands (T16093, 16093C and 16093Y) were detected in individuals showing 16093C or 16093Y in buccal cells. Regarding nucleotide positions 16182 and 16183 in combination with an uninterrupted polyC-stretch, no variation was detected in hairs from individuals showing A16182 16183C in their buccal cells. In contrast, individuals A16182 A16183 showed hair strands with A16182 16183 M and A16182 16183C. And finally, individuals with 16182C 16183C showed some variation in a small amount of their hair strands (some hairs with 16182 M 16183C). These results can be relevant for forensic practitioners when comparing reference samples with hair strands, which is the type of sample most tested by using mtDNA analysis in forensic casework.

Mitochondrial and nuclear disease panel (Mito-aND-Panel): Combined sequencing of mitochondrial and nuclear DNA by a cost-effective and sensitive NGS-based method

Background

The diagnosis of mitochondrial disorders is challenging because of the clinical variability and genetic heterogeneity of these conditions. Next‐Generation Sequencing (NGS) technology offers a robust high‐throughput platform for nuclear and mitochondrial DNA (mtDNA) analyses.

Method

We developed a custom Agilent SureSelect Mitochondrial and Nuclear Disease Panel (Mito‐aND‐Panel) capture kit that allows parallel enrichment for subsequent NGS‐based sequence analysis of nuclear mitochondrial disease‐related genes and the complete mtDNA genome. Sequencing of enriched mtDNA simultaneously with nuclear genes was compared with the separated sequencing of the mitochondrial genome and whole exome sequencing (WES).

Results

The Mito‐aND‐Panel permits accurate detection of low‐level mtDNA heteroplasmy due to a very high sequencing depth compared to standard diagnostic procedures using Sanger sequencing/SNaPshot and WES which is crucial to identify maternally inherited mitochondrial disorders.

Conclusion

We established a NGS‐based method with combined sequencing of the complete mtDNA and nuclear genes which enables a more sensitive heteroplasmy detection of mtDNA mutations compared to traditional methods. Because the method promotes the analysis of mtDNA variants in large cohorts, it is cost‐effective and simple to setup, we anticipate this is a highly relevant method for sequence‐based genetic diagnosis in clinical diagnostic applications.

EditR: A Method to Quantify Base Editing from Sanger Sequencing

CRISPR-Cas9-Cytidine deaminase fusion enzymes—termed “base editors”—allow targeted editing of genomic deoxycytidine to deoxythymidine (C:G→T:A) without the need for double-stranded break induction. Base editors represent a paradigm shift in gene editing technology due to their unprecedented efficiency to mediate targeted, single-base conversion. However, current analysis of base editing outcomes rely on methods that are either imprecise or expensive and time-consuming. To overcome these limitations, we developed a simple, cost-effective, and accurate program to measure base editing efficiency from fluorescence-based Sanger sequencing, termed “EditR.” We provide EditR as a free online tool or downloadable desktop application requiring a single Sanger sequencing file and guide RNA sequence. EditR is more accurate than enzymatic assays, and provides added insight to the position, type, and efficiency of base editing. Furthermore, EditR is likely amenable to quantify base editing from the recently developed adenosine deaminase base editors that act on either DNA (adenosine deaminase base editors [ABEs]) or RNA (REPAIRs) (catalyzes A:T→G:C). Collectively, we demonstrate that EditR is a robust, inexpensive tool that will facilitate the broad application of base editing technology, thereby fostering further innovation in this burgeoning field.

Recent Advances in Detecting Mitochondrial DNA Heteroplasmic Variations

The co-existence of wild-type and mutated mitochondrial DNA (mtDNA) molecules termed heteroplasmy becomes a research hot point of mitochondria. In this review, we listed several methods of mtDNA heteroplasmy research, including the enrichment of mtDNA and the way of calling heteroplasmic variations. At the present, while calling the novel ultra-low level heteroplasmy, high-throughput sequencing method is dominant while the detection limit of recorded mutations is accurate to 0.01% using the other quantitative approaches. In the future, the studies of mtDNA heteroplasmy may pay more attention to the single-cell level and focus on the linkage of mutations.

Machine Learned Replacement of N-Labels for Basecalled Sequences in DNA Barcoding

This study presents a machine learning method that increases the number of identified bases in Sanger Sequencing. The system post-processes a KB basecalled chromatogram. It selects a recoverable subset of N-labels in the KB-called chromatogram to replace with basecalls (A,C,G,T). An N-label correction is defined given an additional read of the same sequence, and a human finished sequence. Corrections are added to the dataset when an alignment determines the additional read and human agree on the identity of the N-label. KB must also rate the replacement with quality value of in the additional read. Corrections are only available during system training. Developing the system, nearly 850,000 N-labels are obtained from Barcode of Life Datasystems, the premier database of genetic markers called DNA Barcodes. Increasing the number of correct bases improves reference sequence reliability, increases sequence identification accuracy, and assures analysis correctness. Keeping with barcoding standards, our system maintains an error rate of percent. Our system only applies corrections when it estimates low rate of error. Tested on this data, our automation selects and recovers: 79 percent of N-labels from COI (animal barcode); 80 percent from matK and rbcL (plant barcodes); and 58 percent from non-protein-coding sequences (across eukaryotes).

Hairy matters: MtDNA quantity and sequence variation along and among human head hairs

Hairs from the same donor have been found to differ in mtDNA sequence within and among themselves and from other tissues, which impacts interpretation of results obtained in a forensic setting. However, little is known on the magnitude of this phenomenon and published data on systematic studies are scarce. We addressed this issue by generating mtDNA control region (CR) profiles of >450 hair fragments from 21 donors by Sanger-type sequencing (STS). To mirror forensic scenarios, we compared hair haplotypes from the same donors to each other, to the corresponding buccal swab reference haplotypes and analyzed several fragments of individual hairs. We also investigated the effects of hair color, donor sex and age, mtDNA haplogroup and chemical treatment on mtDNA quantity, amplification success and variation. We observed a wide range of individual CR sequence variation. The reference haplotype was the only or most common (≥75%) hair haplotype for most donors. However, in two individuals, the reference haplotype was only found in about a third of the investigated hairs, mainly due to differences at highly variable positions. Similarly, most hairs revealed the reference haplotype along their entire length, however, about a fifth of the hairs contained up to 71% of segments with deviant haplotypes, independent of the longitudinal position. Variation affected numerous positions, typically restricted to the individual hair and in most cases heteroplasmic, but also fixed (i.e. homoplasmic) substitutions were observed. While existing forensic mtDNA interpretation guidelines were found still sufficient for all comparisons to reference haplotypes, some comparisons between hairs from the same donor could yield false exclusions when those guidelines are strictly followed. This study pinpoints the special care required when interpreting mtDNA results from hair in forensic casework.

Evidence for frequent and tissue-specific sequence heteroplasmy in human mitochondrial DNA

Mitochondrial point heteroplasmy is a common event observed not only in patients with mitochondrial diseases but also in healthy individuals. We here report a comprehensive investigation of heteroplasmy occurrence in human including the whole mitochondrial control region from nine different tissue types of 100 individuals. Sanger sequencing was used as a standard method and results were supported by cloning, minisequencing, and massively parallel sequencing. Only 12% of all individuals showed no heteroplasmy, whereas 88% showed at least one heteroplasmic position within the investigated tissues. In 66% of individuals up to 8 positions were affected. The highest relative number of heteroplasmies was detected in muscle and liver (79%, 69%), followed by brain, hair, and heart (36.7%-30.2%). Lower percentages were observed in bone, blood, lung, and buccal cells (19.8%-16.2%). Accumulation of position-specific heteroplasmies was found in muscle (positions 64, 72, 73, 189, and 408), liver (position 72) and brain (partial deletion at position 71). Deeper analysis of these specific positions in muscle revealed a non-random appearance and position-specific dependency on age. MtDNA heteroplasmy frequency and its potential functional importance have been underestimated in the past and its occurrence is ubiquitous and dependent at least on age, tissue, and position-specific mutation rates.

Quantitative assessment of heteroplasmy of mitochondrial genome: perspectives in diagnostics and methodological pitfalls

The role of alterations of mitochondrial DNA (mtDNA) in the development of human pathologies is not understood well. Most of mitochondrial mutations are characterized by the phenomenon of heteroplasmy which is defined as the presence of a mixture of more than one type of an organellar genome within a cell or tissue. The level of heteroplasmy varies in wide range, and the expression of disease is dependent on the percent of alleles bearing mutations, thus allowing consumption that an upper threshold level may exist beyond which the mitochondrial function collapses. Recent findings have demonstrated that some mtDNA heteroplasmic mutations are associated with widely spread chronic diseases, including atherosclerosis and cancer. Actually, each etiological mtDNA mutation has its own heteroplasmy threshold that needs to be measured. Therefore, quantitative evaluation of a mutant allele of mitochondrial genome is an obvious methodological challenge, since it may be a keystone for diagnostics of individual genetic predisposition to the disease. This review provides a comprehensive comparison of methods applicable to the measurement of heteroplasmy level of mitochondrial mutations associated with the development of pathology, in particular, in atherosclerosis and its clinical manifestations.

Analysis of the whole mitochondrial genome: translation of the Ion Torrent Personal Genome Machine system to the diagnostic bench?

Next-generation sequencing (NGS), an innovative sequencing technology that enables the successful analysis of numerous gene sequences in a massive parallel sequencing approach, has revolutionized the field of molecular biology. Although NGS was introduced in a rather recent past, the technology has already demonstrated its potential and effectiveness in many research projects, and is now on the verge of being introduced into the diagnostic setting of routine laboratories to delineate the molecular basis of genetic disease in undiagnosed patient samples. We tested a benchtop device on retrospective genomic DNA (gDNA) samples of controls and patients with a clinical suspicion of a mitochondrial DNA disorder. This Ion Torrent Personal Genome Machine platform is a high-throughput sequencer with a fast turnaround time and reasonable running costs. We challenged the chemistry and technology with the analysis and processing of a mutational spectrum composed of samples with single-nucleotide substitutions, indels (insertions and deletions) and large single or multiple deletions, occasionally in heteroplasmy. The output data were compared with previously obtained conventional dideoxy sequencing results and the mitochondrial revised Cambridge Reference Sequence (rCRS). We were able to identify the majority of all nucleotide alterations, but three false-negative results were also encountered in the data set. At the same time, the poor performance of the PGM instrument in regions associated with homopolymeric stretches generated many false-positive miscalls demanding additional manual curation of the data.

Clinicopathological relevance of BRAF mutations in human cancer

BRAF represents one of the most frequently mutated protein kinase genes in human tumours. The mutation is commonly tested in pathology practice. BRAF mutation is seen in melanoma, papillary thyroid carcinoma (including papillary thyroid carcinoma arising from ovarian teratoma), ovarian serous tumours, colorectal carcinoma, gliomas, hepatobiliary carcinomas and hairy cell leukaemia. In these cancers, various genetic aberrations of the BRAF proto-oncogene, such as different point mutations and chromosomal rearrangements, have been reported. The most common mutation, BRAF V600E, can be detected by DNA sequencing and immunohistochemistry on formalin fixed, paraffin embedded tumour tissue. Detection of BRAF V600E mutation has the potential for clinical use as a diagnostic and prognostic marker. In addition, a great deal of research effort has been spent in strategies inhibiting its activity. Indeed, recent clinical trials involving BRAF selective inhibitors exhibited promising response rates in metastatic melanoma patients. Clinical trials are underway for other cancers. However, cutaneous side effects of treatment have been reported and therapeutic response to cancer is short-lived due to the emergence of several resistance mechanisms. In this review, we give an update on the clinical pathological relevance of BRAF mutation in cancer. It is hoped that the review will enhance the direction of future research and assist in more effective use of the knowledge of BRAF mutation in clinical practice.

Mass spectrometric base composition profiling: Implications for forensic mtDNA databasing

In forensic genetics mitochondrial DNA (mtDNA) is usually analyzed by direct Sanger-type sequencing (STS). This method is known to be laborious and sometimes prone to human error. Alternative methods have been proposed that lead to faster results. Among these are methods that involve mass-spectrometry resulting in base composition profiles that are, by definition, less informative than the full nucleotide sequence. Here, we applied a highly automated electrospray ionization mass spectrometry (ESI-MS) system (PLEX-ID) to an mtDNA population study to compare its performance with respect to throughput and concordance to STS. We found that the loss of information power was relatively low compared to the gain in speed and analytical standardization. The detection of point and length heteroplasmy turned out to be roughly comparable between the technologies with some individual differences related to the processes. We confirm that ESI-MS provides a valuable platform for analyzing mtDNA variation that can also be applied in the forensic context.

Diagnosis of mitochondrial disorders applying massive pyrosequencing

Mitochondrial disorders are a frequent cause of neurological disability affecting children and adults. Traditionally, molecular diagnosis of mitochondrial diseases was mostly accomplished by the use of Sanger sequencing and PCR-RFLP. However, there are particular drawbacks associated with the use of these methods. Recent multidisciplinary advances have led to new sequencing methods that may overcome these limitations. Our goal was to explore the use of a next generation sequencing platform in the molecular diagnosis of mitochondrial diseases reporting our findings in adult patients that present with a clinical-pathological diagnosis of a mitochondrial encephalomyopathy. Complete genomic sequences of mitochondrial DNA were obtained by 454 massive pyrosequencing from blood samples. The analysis of these sequences allowed us to identify two diagnostic pathogenic mutations and 74 homoplasmic polymorphisms, useful for obtaining high-resolution mitochondrial haplogroups. In summary, molecular diagnosis of mitochondrial disorders could be efficiently done from readily accessible samples, such as blood, with the use of a new sequencing platform.