Abstract

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Long-read whole genome sequencing is prohibitively expensive for mtDNA heteroplasmy detection and often does not recapitulate the full mtDNA length. Here we describe a method to target, multiplex and sequence full-length, native single-molecule the human mitochondrial genome utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks as barcodes with long-read sequencing, we implemented a protocol in an optimal setting for both high or low integrity genomic DNA to target the circular mitochondrial genome with extremely high coverage. Our analytical pipeline efficiently detects single nucleotide heteroplasmy, physically determines phase and can accurately disentangle complex deletion patterns. This workflow is a unique tool for studying mtDNA variation in health and disease, and will accelerate mitochondrial research.