Mitochondrial Genome

The mitochondrial genome is the complete set of genetic material found within mitochondria, the organelles responsible for cellular energy production. This compact, circular DNA molecule is inherited maternally in most animals and encodes 13 essential proteins for oxidative phosphorylation, along with tRNAs and rRNAs. This guide is intended for biologists, bioinformatics trainees, and clinicians who need a practical, source grounded understanding of mitochondrial genome structure, analysis workflows, and interpretation limits. For foundational context, see the NCBI Bookshelf entry on mitochondrial genetics NCBI Bookshelf.

At a Glance

Core Concepts

The mitochondrial genome sits in the matrix of each mitochondrion. Unlike nuclear DNA, it lacks introns, has a very high mutation rate (especially in the D loop), and uses a distinct genetic code (e.g., AUA codes for methionine instead of isoleucine in humans). Its replication and transcription are driven by nuclear encoded factors. Understanding these features is critical for assembly and annotation. The EMBL EBI training modules provide a comprehensive introduction to mitochondrial genomics EMBL EBI Training. Researchers studying plant mitochondrial genomes, such as the assembly of Ligusticum chuanxiong Assembly of the complete mitochondrial genome of Ligusticum chuanxiong and its evolutionary implications, often encounter larger circular molecules with frequent recombination and acquisition of foreign DNA. In contrast, animal mitochondrial genomes are more conserved in size and gene order.

Decision Points for Mitochondrial Genome Analysis

Before starting a mitochondrial genome project, consider these decision points to avoid wasted effort.

Sample type and storage. Buccal swabs stored on FTA cards can yield sufficient DNA for whole genome sequencing, but DNA quality declines over decades as noted by the authors of a study on DNA quantity and quality in buccal samples stored on FTA cards Nearly two decades on paper: DNA quantity and quality in buccal samples stored on FTA cards. If DNA is degraded, short read sequencing may suffice. For high quality DNA, long reads can resolve repetitive regions.

Sequencing depth. Because of high copy number, mitochondrial genomes can be enriched by PCR or captured. However, off target nuclear mitochondrial segments (NUMTs) can contaminate assemblies. A read depth of 50x after enrichment is typical. The Galaxy Training Network offers workflows that include coverage thresholds for reliable assembly Galaxy Training Network.

Assembly strategy. Use a reference guided approach for close relatives of well characterized species. For novel or highly divergent genomes, de novo assembly is required. Tools like NOVOPlasty (for circular genomes) or MitoFinder (for annotation) are common choices. The Bioconductor repository provides R packages such as MitochondrialGenetics for downstream analysis Bioconductor.

Validation. Always confirm circularization (the assembly should form a circle). Check for NUMT contamination by mapping reads to the assembled genome and looking for regions with abrupt coverage drop off.

Practical Workflow for Assembly and Annotation

The following steps assume you have raw sequencing reads (FASTQ) from whole genome sequencing or enriched libraries.

1. Quality control. Run FastQC or MultiQC to check base quality, adapter content, and duplication levels. Remove adapters and low quality bases with Trimmomatic or cutadapt. Raw reads can be retrieved from the NCBI Sequence Read Archive NCBI Sequence Read Archive.

2. Mitochondrial read extraction. Map all reads to a reference mitochondrial genome (e.g., human rCRS) using a sensitive mapper like BWA. Extract mapped reads. Alternatively, use baiting and iterative mapping with MITObim for circular genomes.

3. De novo assembly. Assemble extracted reads using SPAdes with the careful flag or use MitoZ. For long reads, Canu or Flye can produce a single contig. Check that the assembly forms a circle by inspecting the contig ends.

4. Circularization and correction. If the assembly is not circular, identify overlapping ends using BLAST against itself. Align to a related genome to confirm orientation. Polish the assembly with Pilon using the original reads.

5. Annotation. Use MITOS, MitoAnnotator, or GeSeq for gene annotation. Verify start codons and tRNAs. For plant genomes, use IPMGA. Compare annotation with closely related genomes. The assembly of Setaria viridis mitochondrial genome illustrates how annotation can reveal lineage specific gene losses and rearrangements Assembly and characterization of the complete mitochondrial genome of Setaria viridis (L.) Beauv..

6. Variant calling. To identify heteroplasmy or somatic mutations, call variants using GATK Mutect2 or LoFreq on the reassembled genome. Filter for strand bias and read position.

Quality Checks and Validation

After assembly and annotation, perform these checks.

Coverage uniformity. Create a coverage plot across the genome. Even coverage suggests no large NUMT contamination. A sharp drop may indicate a collapsed repeat.

Codon usage and genetic code. Ensure predicted protein sequences match the known mitochondrial genetic code. Compare translation with BLAST to related species.

Circular genome verification. Confirm that the assembly circularizes with no remaining linear ends. Use dot plots to look for repeat induced misassembly.

Functional completeness. Check that all core genes (13 in animals) are present and have proper reading frames. Missing genes may indicate assembly gaps or pseudogenization, common in some parasitic lineages as seen in the mitochondrial genomes of Klossia species Mitochondrial Genomes of Klossia Species: The First Complete Extrachromosomal Sequence Data From Parasites in the Monoxenous Apicomplexan Family Adeleidae.

Common Mistakes

Ignoring NUMTs. Nuclear pseudogenes of mitochondrial origin can be co assembled as part of the mitochondrial genome. This inflates genome size and introduces sequence errors. Filter by requiring consistent read depth and by mapping reads back to the full nuclear genome.

Assuming circularity without proof. Linear assemblies that are actually circular lead incorrect genome length and gene order. Always search for overlapping ends.

Using the wrong genetic code. Translating mitochondrial genes with the standard nuclear code produces nonsense proteins. Set translation table to Table 13 (for most animal mitochondria) or appropriate table for plants/fungi.

Overlooking heteroplasmy. Low frequency variants may represent true heteroplasmy, not sequencing errors. Validate with a second sequencing method or PCR.

Poorly annotated tRNAs. tRNA genes have cloverleaf structure and are often missed by automated annotation. Manually inspect with tRNAscan SE.

Limits and Interpretation

Mitochondrial genome analysis has several inherent limits.

Heteroplasmy detection power. Variants below 1% frequency may be difficult to detect with short read sequencing. For clinical applications, deep sequencing (1000x or more) is recommended. The clinical relevance of low level heteroplasmy remains an area of active research, as discussed in the context of mitochondrial ncRNAs and cardiovascular diseases Mitochondrial ncRNAs: From Pathological Regulation to Targeted Therapy in Cardiovascular Diseases.

Phylogenetic resolution. While mitochondrial genomes are widely used for phylogenetics, high mutation rates lead to homoplasy and saturation at deep evolutionary timescales. For recent divergences, they are excellent markers. For ancient splits, combine with nuclear markers.

Functional inference. A complete mitochondrial genome sequence does not reveal mitochondrial function. To assess oxidative phosphorylation capacity, complementary proteomics or respirometry is required. The case report on nevus lipomatosus with comprehensive genetic analysis illustrates that even a complete mitochondrial genome may not explain a phenotype without functional studies Generalized Nevus Lipomatosus Cutaneous Superficialis: A Case Report with Comprehensive Genetic Analysis.

NUMT interference in forensics. Mitochondrial DNA typing for identification can be confused by NUMTs. Laboratories should use primers that avoid NUMT amplification and sequence from multiple overlapping amplicons.

Frequently Asked Questions

1. Why is the mitochondrial genome used in forensics? Mitochondrial DNA is abundant (high copy number per cell) and persists in degraded samples such as hair shafts or old bones. Its hypervariable regions allow discrimination between individuals when nuclear DNA is insufficient.

2. Can mitochondrial DNA be inherited from the father? In rare cases, paternal leakage occurs in some species, including humans very rarely. However, the standard model is strict maternal inheritance in mammals. Specialized studies may need to consider this possibility.

3. How do I choose between reference guided and de novo assembly? If a high quality reference genome from a closely related species exists (e.g., within the same genus), reference guided assembly is faster and more accurate. For novel species or if you suspect large rearrangements, use de novo assembly.

4. What is a NUMT and how do I avoid it? NUMT stands for nuclear mitochondrial segment. It is a piece of mitochondrial DNA that has been inserted into the nuclear genome. To avoid including NUMTs in your assembly, cross map reads to both the nuclear and mitochondrial references. Use read depth analysis, NUMTs typically have lower coverage than the true mitochondrial genome.

References and Further Reading

• NCBI Bookshelf: Mitochondrial Genetics and Human Disease NCBI Bookshelf
• EMBL EBI Training: Introduction to Mitochondrial Genomics EMBL EBI Training
• Galaxy Training Network: Mitochondrial Genome Assembly Galaxy Training Network
• Bioconductor: MitochondrialGenetics Package Bioconductor
• NCBI Sequence Read Archive: Repository for Raw Sequencing Data NCBI Sequence Read Archive
• Assembly of the complete mitochondrial genome of Ligusticum chuanxiong and its evolutionary implications. BMC Genomics. PubMed
• Assembly and characterization of the complete mitochondrial genome of Setaria viridis (L.) Beauv. Front Plant Sci. PubMed
• Mitochondrial ncRNAs: From Pathological Regulation to Targeted Therapy in Cardiovascular Diseases. J Cardiovasc Transl Res. PubMed
• Nearly two decades on paper: DNA quantity and quality in buccal samples stored on FTA cards. J Forensic Sci. PubMed
• Generalized Nevus Lipomatosus Cutaneous Superficialis: A Case Report with Comprehensive Genetic Analysis. Appl Clin Genet. PubMed
• Mitochondrial Genomes of Klossia Species: The First Complete Extrachromosomal Sequence Data From Parasites in the Monoxenous Apicomplexan Family Adeleidae. J Eukaryot Microbiol. PubMed

Related Articles

• Protein Synthesis
• Incomplete Dominance Gene
• Cell Membrane Function Biology
• Dna Structure
• Protein Structure