Nanopore sequencing: An enrichment-free alternative to mitochondrial DNA sequencing

Estimating bloodstain age in the short term based on DNA fragment length using nanopore sequencer

We used a nanopore sequencer to quantify DNA fragments > 10,000 bp in size and then evaluated their relationship with short-term bloodstain age. Moreover, DNA degradation was investigated after bloodstains were wetted once with water. Bloodstain samples on cotton gauze were stored at room temperature and low humidity for up to 6 months. Bloodstains stored for 1 day were wetted with nuclease-free water, allowed to dry, and stored at room temperature and low humidity for up to 1 week. The proportion of fragments > 20,000 bp in dry bloodstains tended to decrease over time, particularly for fragments > 50,000 bp in size. This trend was modeled using a power approximation curve, with the highest R2 value (0.6475) noted for fragments > 50,000 bp in size; lower values were recorded for shorter fragments. The proportion of longer fragments was significantly reduced in bloodstains that were dried after being wetted once, and there was significant difference in fragments > 50,000 bp between dry conditions and once-wetted. This result suggests that even temporary exposure to water causes significant DNA fragmentation, but not extensive degradation. Thus, bloodstains that appear fresh but have a low proportion of long DNA fragments may have been wetted previously. Our results indicate that evaluating the proportion of long DNA fragments yields information on both bloodstain age and the environment in which they were stored.

Genetic characteristics of complete mtDNA genome sequence of Indonesian local rabbit (Oryctolagus cuniculus)

BACKGROUND

Indonesian local rabbit (Oryctolagus cuniculus) is a local breed in Indonesia. We reveal the mitochondrial genome sequence of the Indonesian local Rabbit for the first time. A better understanding of the mechanisms underlying these beneficial aspects of local breeds over imported ones requires detailed genetic investigations, of which mtDNA genome sequencing is of particular importance. Such an investigation will solve the major issues of misidentification with Javanese hares (Lepus nigricollis) and maternal lineage. In addition, this information will guide better statistics on the Indonesian local rabbit breed population and strategies for its conservation and breeding plans. This study aimed to identify and explore the characteristics of the mtDNA genomes of Indonesian local rabbits.

RESULT

This study observed that the length of the mtDNA genome is 17,469 bp, consisting of two ribosomal RNA (12S rRNA, 16S rRNA), 22 transfer RNA genes (trnR, trnG, trnK, trnD, trnS, trnY, trnC, trnN, trnA, trnW, trnM, trnQ, trnl, trnL, trnV, trnF, trnP, trnT, trnE, trnL, trnS, trnH), 13 protein-coding genes (PCGs) (ND4l, ND3, COX3, ATP6, ATP8, COX2, COX1, ND2, ND1, CYTB, ND6, ND5, ND4), a replication origin, and a noncoding control region (D-loop).

CONCLUSIONS

mtDNA genome of Indonesian local rabbit was the longest and had the most extended D-loop sequence among the other references of Oryctolagus cuniculus. Other specific differences were also found in the percentage of nucleotides and variation in most of the PCGs when they were aligned with Oryctolagus cuniculus references from GenBank. Indonesian local Rabbits strongly suspected brought from Europe during the colonial period in Indonesia.

Nanopore long-read next-generation sequencing for detection of mitochondrial DNA large-scale deletions

Primary mitochondrial diseases are progressive genetic disorders affecting multiple organs and characterized by mitochondrial dysfunction. These disorders can be caused by mutations in nuclear genes coding proteins with mitochondrial localization or by genetic defects in the mitochondrial genome (mtDNA). The latter include point pathogenic variants and large-scale deletions/rearrangements. MtDNA molecules with the wild type or a variant sequence can exist together in a single cell, a condition known as mtDNA heteroplasmy. MtDNA single point mutations are typically detected by means of Next-Generation Sequencing (NGS) based on short reads which, however, are limited for the identification of structural mtDNA alterations. Recently, new NGS technologies based on long reads have been released, allowing to obtain sequences of several kilobases in length; this approach is suitable for detection of structural alterations affecting the mitochondrial genome. In the present work we illustrate the optimization of two sequencing protocols based on long-read Oxford Nanopore Technology to detect mtDNA structural alterations. This approach presents strong advantages in the analysis of mtDNA compared to both short-read NGS and traditional techniques, potentially becoming the method of choice for genetic studies on mtDNA.

Enhancing Molecular Testing for Effective Delivery of Actionable Gene Diagnostics

There is a deep need to navigate within our genomic data to find, understand and pave the way for disease-specific treatments, as the clinical diagnostic journey provides only limited guidance. The human genome is enclosed in every nucleated cell, and yet at the single-cell resolution many unanswered questions remain, as most of the sequencing techniques use a bulk approach. Therefore, heterogeneity, mosaicism and many complex structural variants remain partially uncovered. As a conceptual approach, nanopore-based sequencing holds the promise of being a single-molecule-based, long-read and high-resolution technique, with the ability of uncovering the nucleic acid sequence and methylation almost in real time. A key limiting factor of current clinical genetics is the deciphering of key disease-causing genomic sequences. As the technological revolution is expanding regarding genetic data, the interpretation of genotype-phenotype correlations should be made with fine caution, as more and more evidence points toward the presence of more than one pathogenic variant acting together as a result of intergenic interplay in the background of a certain phenotype observed in a patient. This is in conjunction with the observation that many inheritable disorders manifest in a phenotypic spectrum, even in an intra-familial way. In the present review, we summarized the relevant data on nanopore sequencing regarding clinical genomics as well as highlighted the importance and content of pre-test and post-test genetic counselling, yielding a complex approach to phenotype-driven molecular diagnosis. This should significantly lower the time-to-right diagnosis as well lower the time required to complete a currently incomplete genotype-phenotype axis, which will boost the chance of establishing a new actionable diagnosis followed by therapeutical approach.

Next generation biosecurity: Towards genome based identification to prevent spread of agronomic pests and pathogens using nanopore sequencing

The unintentional movement of agronomic pests and pathogens is steadily increasing due to the intensification of global trade. Being able to identify accurately and rapidly early stages of an invasion is critical for developing successful eradication or management strategies. For most invasive organisms, molecular diagnostics is today the method of choice for species identification. However, the currently implemented tools are often developed for certain taxa and need to be adapted for new species, making them ill-suited to cope with the current constant increase in new invasive species. To alleviate this impediment, we developed a fast and accurate sequencing tool allowing to modularly obtain genetic information at different taxonomical levels. Using whole genome amplification (WGA) followed by Oxford nanopore MinION sequencing, our workflow does not require any a priori knowledge on the investigated species and its classification. While mainly focusing on harmful plant pathogenic insects, we also demonstrate the suitability of our workflow for the molecular identification of bacteria (Erwinia amylovora and Escherichia coli), fungi (Cladosporium herbarum, Colletotrichum salicis, Neofabraea alba) and nematodes (Globodera rostochiensis). On average, the pairwise identity between the generated consensus sequences and best GenBank BLAST matches was 99.6 ± 0.6%. Additionally, assessing the generated insect genomic dataset, the potential power of the workflow to detect pesticide resistance genes, as well as arthropod-infecting viruses and endosymbiotic bacteria is demonstrated.

Benchmarking Low-Frequency Variant Calling With Long-Read Data on Mitochondrial DNA

Background: Sequencing quality has improved over the last decade for long-reads, allowing for more accurate detection of somatic low-frequency variants. In this study, we used mixtures of mitochondrial samples with different haplogroups (i.e., a specific set of mitochondrial variants) to investigate the applicability of nanopore sequencing for low-frequency single nucleotide variant detection.

Methods: We investigated the impact of base-calling, alignment/mapping, quality control steps, and variant calling by comparing the results to a previously derived short-read gold standard generated on the Illumina NextSeq. For nanopore sequencing, six mixtures of four different haplotypes were prepared, allowing us to reliably check for expected variants at the predefined 5%, 2%, and 1% mixture levels. We used two different versions of Guppy for base-calling, two aligners (i.e., Minimap2 and Ngmlr), and three variant callers (i.e., Mutserve2, Freebayes, and Nanopanel2) to compare low-frequency variants. We used F₁ score measurements to assess the performance of variant calling.

Results: We observed a mean read length of 11 kb and a mean overall read quality of 15. Ngmlr showed not only higher F₁ scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F₁ score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F₁ score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F₁ scores (5% level: F₁ score >0.99, 2% level: F₁ score >0.54, and 1% level: F₁ score >0.70) across all callers and mixture levels.

Conclusion: We here present the benchmarking for low-frequency variant calling with nanopore sequencing by identifying current limitations.

Translocation Behaviors of Synthetic Polyelectrolytes through Alpha-Hemolysin (α-HL) and Mycobacterium smegmatis Porin A (MspA) Nanopores

DNAs have been used as probes for nanopore sensing of noncharged biomacromolecules due to its negative phosphate backbone. Inspired by this, we explored the potential of diblock synthetic polyelectrolytes as more flexible and inexpensive nanopore sensing probes by investigating translocation behaviors of PEO-b-PSS and PEO-b-PVBTMA through commonly used alpha-hemolysin (α-HL) and Mycobacterium smegmatis porin A (MspA) nanopores. Translocation recordings in different configurations of pore orientation and testing voltage indicated efficient PEO-b-PSS translocations through α-HL and PEO-b-PVBTMA translocations through MspA. This work provides insight into synthetic polyelectrolyte-based probes to expand probe selection and flexibility for nanopore sensing.

Estimating individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq

We tried to estimate individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq. The BAM files produced by MiSeq were viewed using Integrative Genomics Viewer (IGV) to verify mixed bases. By sorting the reads according to base type for each mixed base, partial haplotypes were determined. Then, the BAM files produced by MinKNOW were viewed using IGV. To determine haplotypes with IGV, only mixed bases determined by MiSeq were used as target bases. By sorting the reads according to base type for each target base, each contributor's haplotype was estimated. In mixed samples from two contributors, even a haplotype with a minor contribution of 5% could be distinguished from the haplotype of the major contributor. In mixed samples of three contributors (mixture ratios of 1:1:1 and 4:2:1), each haplotype could also be distinguished. Sequences of C-stretches were determined very inaccurately in the MinION analysis. Although the analysis method was simple, each haplotype was correctly detected in all mixed samples with two or three contributors in various mixture ratios by combining MinION and MiSeq. This should be useful for identifying contributors to mixed samples.

Accurate profiling of forensic autosomal STRs using the Oxford Nanopore Technologies MinION device

The high variability characteristic of short tandem repeat (STR) markers is harnessed for human identification in forensic genetic analyses. Despite the power and reliability of current typing techniques, sequence-level information both within and around STRs are masked in the length-based profiles generated. Forensic STR typing using next generation sequencing (NGS) has therefore gained attention as an alternative to traditional capillary electrophoresis (CE) approaches. In this proof-of-principle study, we evaluate the forensic applicability of the newest and smallest NGS platform available - the Oxford Nanopore Technologies (ONT) MinION device. Although nanopore sequencing on the handheld MinION offers numerous advantages, including low startup cost and on-site sample processing, the relatively high error rate and lack of forensic-specific analysis software has prevented accurate profiling across STR panels in previous studies. Here we present STRspy, a streamlined method capable of producing length- and sequence-based STR allele designations from noisy, error-prone third generation sequencing reads. To assess the capabilities of STRspy, seven reference samples (female: n = 2; male: n = 5) were amplified at 15 and 30 PCR cycles using the Promega PowerSeq 46GY System and sequenced on the ONT MinION device in triplicate. Basecalled reads were then processed with STRspy using a custom database containing alleles reported in the STRSeq BioProject NIST 1036 dataset. Resultant STR allele designations and flanking region single nucleotide polymorphism (SNP) calls were compared to the manufacturer-validated genotypes for each sample. STRspy generated robust and reliable genotypes across all autosomal STR loci amplified with 30 PCR cycles, achieving 100% concordance based on both length and sequence. Furthermore, we were able to identify flanking region SNPs in the 15-cycle dataset with > 90% accuracy. These results demonstrate that when analyzed with STRspy ONT reads can reveal additional variation in and around STR loci depending on read coverage. As the first and only third generation sequencing platform-specific method to successfully profile the entire panel of autosomal STRs amplified by a commercially available multiplex, STRspy significantly increases the feasibility of nanopore sequencing in forensic applications.

Nanopore Single-Molecule Sequencing for Mitochondrial DNA Methylation Analysis: Investigating Parkin-Associated Parkinsonism as a Proof of Concept

Objective: To establish a workflow for mitochondrial DNA (mtDNA) CpG methylation using Nanopore whole-genome sequencing and perform first pilot experiments on affected Parkin biallelic mutation carriers (Parkin-PD) and healthy controls. Background: Mitochondria, including mtDNA, are established key players in Parkinson's disease (PD) pathogenesis. Mutations in Parkin, essential for degradation of damaged mitochondria, cause early-onset PD. However, mtDNA methylation and its implication in PD is understudied. Herein, we establish a workflow using Nanopore sequencing to directly detect mtDNA CpG methylation and compare mtDNA methylation between Parkin-related PD and healthy individuals. Methods: To obtain mtDNA, whole-genome Nanopore sequencing was performed on blood-derived from five Parkin-PD and three control subjects. In addition, induced pluripotent stem cell (iPSC)-derived midbrain neurons from four of these patients with PD and the three control subjects were investigated. The workflow was validated, using methylated and unmethylated 897 bp synthetic DNA samples at different dilution ratios (0, 50, 100% methylation) and mtDNA without methylation. MtDNA CpG methylation frequency (MF) was detected using Nanopolish and Megalodon. Results: Across all blood-derived samples, we obtained a mean coverage of 250.3X (SD ± 80.5X) and across all neuron-derived samples 830X (SD ± 465X) of the mitochondrial genome. We detected overall low-level CpG methylation from the blood-derived DNA (mean MF ± SD = 0.029 ± 0.041) and neuron-derived DNA (mean MF ± SD = 0.019 ± 0.035). Validation of the workflow, using synthetic DNA samples showed that highly methylated DNA molecules were prone to lower Guppy Phred quality scores and thereby more likely to fail Guppy base-calling. CpG methylation in blood- and neuron-derived DNA was significantly lower in Parkin-PD compared to controls (Mann-Whitney U-test p < 0.05). Conclusion: Nanopore sequencing is a useful method to investigate mtDNA methylation architecture, including Guppy-failed reads is of importance when investigating highly methylated sites. We present a mtDNA methylation workflow and suggest methylation variability across different tissues and between Parkin-PD patients and controls as an initial model to investigate.

NanoCaller for accurate detection of SNPs and indels in difficult-to-map regions from long-read sequencing by haplotype-aware deep neural networks

Long-read sequencing enables variant detection in genomic regions that are considered difficult-to-map by short-read sequencing. To fully exploit the benefits of longer reads, here we present a deep learning method NanoCaller, which detects SNPs using long-range haplotype information, then phases long reads with called SNPs and calls indels with local realignment. Evaluation on 8 human genomes demonstrates that NanoCaller generally achieves better performance than competing approaches. We experimentally validate 41 novel variants in a widely used benchmarking genome, which could not be reliably detected previously. In summary, NanoCaller facilitates the discovery of novel variants in complex genomic regions from long-read sequencing.

DNA-nanopore technology: a human perspective

The purpose of this article is to give a brief overview of the current state of nanopore sequencing in relation to forensic science with a brief outline of where it stands in relation to current methods, its potential uses in forensic science and factors which may influence acceptance of this technology by forensic practitioners, the judiciary and law enforcement. Perhaps most importantly consideration is also given to concerns which may influence the acceptance of the technology by the general public.

The Isolation and Deep Sequencing of Mitochondrial DNA

In recent years, next-generation sequencing (NGS) has become a powerful tool for studying both inherited and somatic heteroplasmic mitochondrial DNA (mtDNA) variation. NGS has proved particularly powerful when combined with single-cell isolation techniques, allowing the investigation of low-level heteroplasmic variants both between cells and within tissues. Nevertheless, there remain significant challenges, especially around the selective enrichment of mtDNA from total cellular DNA and the avoidance of nuclear pseudogenes. This chapter summarizes the techniques needed to enrich, amplify, sequence, and analyse mtDNA using NGS .

Applying genomic and transcriptomic advances to mitochondrial medicine

Next-generation sequencing (NGS) has increased our understanding of the molecular basis of many primary mitochondrial diseases (PMDs). Despite this progress, many patients with suspected PMD remain without a genetic diagnosis, which restricts their access to in-depth genetic counselling, reproductive options and clinical trials, in addition to hampering efforts to understand the underlying disease mechanisms. Although they represent a considerable improvement over their predecessors, current methods for sequencing the mitochondrial and nuclear genomes have important limitations, and molecular diagnostic techniques are often manual and time consuming. However, recent advances in genomics and transcriptomics offer realistic solutions to these challenges. In this Review, we discuss the current genetic testing approach for PMDs and the opportunities that exist for increased use of whole-genome NGS of nuclear and mitochondrial DNA (mtDNA) in the clinical environment. We consider the possible role for long-read approaches in sequencing of mtDNA and in the identification of novel nuclear genomic causes of PMDs. We examine the expanding applications of RNA sequencing, including the detection of cryptic variants that affect splicing and gene expression and the interpretation of rare and novel mitochondrial transfer RNA variants.

The 287,403 bp Mitochondrial Genome of Ectomycorrhizal Fungus Tuber calosporum Reveals Intron Expansion, tRNA Loss, and Gene Rearrangement

In the present study, the mitogenome of Tuber calosporum was assembled and analyzed. The mitogenome of T. calosporum comprises 15 conserved protein-coding genes, two rRNA genes, and 14 tRNAs, with a total size of 287,403 bp. Fifty-eight introns with 170 intronic open reading frames were detected in the T. calosporum mitogenome. The intronic region occupied 69.41% of the T. calosporum mitogenome, which contributed to the T. calosporum mitogenome significantly expand relative to most fungal species. Comparative mitogenomic analysis revealed large-scale gene rearrangements occurred in the mitogenome of T. calosporum, involving gene relocations and position exchanges. The mitogenome of T. calosporum was found to have lost several tRNA genes encoding for cysteine, aspartate, histidine, etc. In addition, a pair of fragments with a total length of 32.91 kb in both the nuclear and mitochondrial genomes of T. calosporum was detected, indicating possible gene transfer events. A total of 12.83% intragenomic duplications were detected in the T. calosporum mitogenome. Phylogenetic analysis based on mitochondrial gene datasets obtained well-supported tree topologies, indicating that mitochondrial genes could be reliable molecular markers for phylogenetic analyses of Ascomycota. This study served as the first report on mitogenome in the family Tuberaceae, thereby laying the groundwork for our understanding of the evolution, phylogeny, and population genetics of these important ectomycorrhizal fungi.

Nanopore Sequencing of the Fungal Intergenic Spacer Sequence as a Potential Rapid Diagnostic Assay

Fungal infections are being caused by a broadening spectrum of fungi, yet in many cases, identification to the species level is required for proper antifungal selection. We investigated the fungal intergenic spacer (IGS) sequence in combination with nanopore sequencing for fungal identification. We sequenced isolates from two Cryptococcus species complexes, C. gattii and C. neoformans, which are the main pathogenic members of this genus, using the Oxford Nanopore Technologies MinION device and Sanger sequencing. There is enough variation within the two complexes to argue for further resolution into separate species, which we wanted to see if nanopore sequencing could detect. Using the R9.4.1 flow cell, IGS sequence identities averaged 99.57% compared to Sanger sequences of the same region. When the newer R10.3 flow cell was used, accuracy increased to 99.83% identity compared to the same Sanger sequences. Nanopore sequencing errors were predominantly in regions of homopolymers, with G homopolymers displaying the largest number of errors and C homopolymers displaying the least. Phylogenetic analysis of the nanopore- and Sanger-derived sequences resulted in indistinguishable trees. Comparison of average percent identities between the C. gattii and C. neoformans species complexes resulted in only a 74 to 77% identity between the two complexes. Sequencing using the nanopore platform could be completed in less than an hour, and samples could be multiplexed in groups as large as 24 sequences in a single run. These results suggest that sequencing the IGS region using nanopore sequencing could be a potential new molecular diagnostic strategy.

A new method for long-read sequencing of animal mitochondrial genomes: application to the identification of equine mitochondrial DNA variants

BACKGROUND

Mitochondrial DNA is remarkably polymorphic. This is why animal geneticists survey mitochondrial genomes variations for fundamental and applied purposes. We present here an approach to sequence whole mitochondrial genomes using nanopore long-read sequencing. Our method relies on the selective elimination of nuclear DNA using an exonuclease treatment and on the amplification of circular mitochondrial DNA using a multiple displacement amplification step.

RESULTS

We optimized each preparative step to obtain a 100 million-fold enrichment of horse mitochondrial DNA relative to nuclear DNA. We sequenced these amplified mitochondrial DNA using nanopore sequencing technology and obtained mitochondrial DNA reads that represented up to half of the sequencing output. The sequence reads were 2.3 kb of mean length and provided an even coverage of the mitochondrial genome. Long-reads spanning half or more of the whole mtDNA provided a coverage that varied between 118X and 488X. We evaluated SNPs identified using these long-reads by Sanger sequencing as ground truth and found a precision of 100.0%; a recall of 93.1% and a F1-score of 0.964 using the Twilight horse mtDNA reference. The choice of the mtDNA reference impacted variant calling efficiency with F1-scores varying between 0.947 and 0.964.

CONCLUSIONS

Our method to amplify mtDNA and to sequence it using the nanopore technology is usable for mitochondrial DNA variant analysis. With minor modifications, this approach could easily be applied to other large circular DNA molecules.

Approaches to Whole Mitochondrial Genome Sequencing on the Oxford Nanopore MinION

Traditional approaches for interrogating the mitochondrial genome often involve laborious extraction and enrichment protocols followed by Sanger sequencing. Although preparation techniques are still demanding, the advent of next-generation or massively parallel sequencing has made it possible to routinely obtain nucleotide-level data with relative ease. These short-read sequencing platforms offer deep coverage with unparalleled read accuracy in high-complexity genomic regions but encounter numerous difficulties in the low-complexity homopolymeric sequences characteristic of the mitochondrial genome. The inability to discern identical units within monomeric repeats and resolve copy-number variations for heteroplasmy detection results in suboptimal genome assemblies that ultimately complicate downstream data analysis and interpretation of biological significance. Oxford Nanopore Technologies offers the ability to generate long-read sequencing data on a pocket-sized device known as the MinION. Nanopore-based sequencing is scalable, portable, and theoretically capable of sequencing the entire mitochondrial genome in a single contig. Furthermore, the recent development of a nanopore protein with dual reader heads allows for clear identification of nucleotides within homopolymeric stretches, significantly increasing resolution throughout these regions. The unrestricted read lengths, superior homopolymeric resolution, and affordability of the MinION device make it an attractive alternative to the labor-intensive, time-consuming, and costly mainstay deep-sequencing platforms. This article describes three approaches to extract, prepare, and sequence mitochondrial DNA on the Oxford Nanopore MinION device. Two of the workflows include enrichment of mitochondrial DNA prior to sequencing, whereas the other relies on direct sequencing of native genomic DNA to allow for simultaneous assessment of the nuclear and mitochondrial genomes. © 2019 by John Wiley & Sons, Inc. Basic Protocol: Enrichment-free mitochondrial DNA sequencing Alternate Protocol 1: Mitochondrial DNA sequencing following enrichment with polymerase chain reaction (PCR) Alternate Protocol 2: Mitochondrial DNA sequencing following enrichment with PCR-free hybridization capture Support Protocol 1: DNA quantification and quality assessment using the Agilent 4200 TapeStation System Support Protocol 2: AMPure XP bead clean-up Support Protocol 3: Suggested data analysis pipeline.

Ion torrent high throughput mitochondrial genome sequencing (HTMGS)

The implementation and popularity of next generation sequencing (NGS) has led to the development of various rapid whole mitochondrial genome sequencing techniques. We summarise an efficient and cost-effective NGS approach for mitochondrial genomic DNA in humans using the Ion Torrent platform, and further discuss our bioinformatics pipeline for streamlined variant calling. Ion 316 chips were utilised with the Ion Torrent semi-conductor platform Personal Genome Machine (PGM) to perform tandem sequencing of mitochondrial genomes from the core pedigree (n = 315) of the Norfolk Island Health Study. Key improvements from commercial methods focus on the initial PCR step, which currently requires extensive optimisation to ensure the accurate and reproducible elongation of each section of the complete mitochondrial genome. Dual-platform barcodes were incorporated into our protocol thereby extending its potential application onto Illumina-based systems. Our bioinformatics pipeline consists of a modified version of GATK best practices tailored for mitochondrial genomic data. When compared with current commercial methods, our method, termed high throughput mitochondrial genome sequencing (HTMGS), allows high multiplexing of samples and the use of alternate library preparation reagents at a lower cost per sample (~1.7 times) when compared to current commercial methodologies. Our HTMGS methodology also provides robust mitochondrial sequencing data (>450X average coverage) that can be applied and modified to suit various study designs. On average, we were able to identify ~30 variants per sample with 572 variants observed across 315 samples. We have developed a high throughput sequencing and analysis method targeting complete mitochondrial genomes; with the potential to be platform agnostic with analysis options that adhere to current best practices.