Mitochondrial Genome Assembly for Parasite Taxonomy: From Sequencing to Species Identification and Drug Discovery

Abstract

This article provides a comprehensive resource for researchers and scientists utilizing mitochondrial genome assembly in parasite taxonomy and drug discovery. It explores the foundational principles of parasite mitochondrial genomics, detailing unique structural features and their phylogenetic significance. The content covers advanced methodological workflows from sample collection to genome annotation, addresses key troubleshooting strategies for complex genomes, and establishes rigorous validation and comparative analysis frameworks. By integrating current case studies and technological advances, this guide serves to enhance accurate species identification and support the development of novel therapeutic targets for parasitic diseases.

The Foundation of Parasite Mitochondrial Genomics: Structure, Variation, and Phylogenetic Power

Mitochondrial genomes in parasites exhibit remarkable structural diversity that deviates significantly from the standard circular model observed in most metazoans. These organellar genomes have evolved into various forms, including linear monomers, concatemers, and fragmented minichromosomes, providing valuable insights into evolutionary biology and serving as crucial molecular markers for taxonomic classification. The phylum Apicomplexa, which contains medically important parasites such as Plasmodium (malaria), Babesia, and Theileria, demonstrates particularly fascinating variations in mitochondrial architecture [1] [2]. Unlike typical animal mitochondrial genomes that range from 15-20 kb and contain 37 genes, parasitic protists often have significantly reduced mitochondrial genomes, sometimes as small as 6 kb, encoding only a handful of proteins alongside fragmented ribosomal RNA genes [1] [3]. This structural diversity reflects the complex evolutionary pathways and adaptive strategies these parasites have undergone, making mitochondrial genomics an invaluable tool for understanding parasite biology, evolution, and taxonomy.

Structural Variations Across Parasite Lineages

Linear Mitochondrial Genomes in Apicomplexan Parasites

Babesia and Theileria species, which cause piroplasmosis in animals, possess monomeric linear mitochondrial genomes ranging from 6.6 kb to 11.1 kb [1] [2]. These linear molecules consistently encode three protein-coding genes (cox1, cox3, and cob) and multiple fragmented large subunit (LSU) ribosomal RNA genes [1]. A defining characteristic of these linear genomes is the presence of terminal inverted repeats (TIRs) at both ends, which play crucial roles in replication and stability [1] [3].

Table 1: Characteristics of Linear Mitochondrial Genomes in Selected Apicomplexan Parasites

Parasite	Genome Size (bp)	Structure	Protein Genes	Terminal Repeats	Special Features
Babesia microti	11,100	Linear monomer	cox1, cob, cox3	Dual flip-flop inversion system (IR-A, IR-B)	Four distinct genome structures via inversion [2]
Babesia rodhaini	6,900	Linear monomer	cox1, cob, cox3	Dual flip-flop inversion system (IR-A, IR-B)	Four distinct genome structures via inversion [2]
Theileria velifera	6,125	Linear monomer	cox1, cob, cox3	Terminal inverted repeats (TIRs)	5 LSU rRNA fragments [3]
Theileria equi	8,200	Linear monomer	cox1, cob, cox3	Unusually long TIRs	Largest and most divergent Theileria mt genome [1]

The evolutionary significance of these linear genomes becomes apparent when compared to closely related parasites. Plasmodium species, despite being phylogenetically close to Babesia and Theileria, possess 6-kb concatenated linear mitochondrial genomes with different gene arrangements and transcriptional directions [1]. This structural divergence suggests distinct evolutionary pathways in these parasite lineages. Furthermore, the archaeopiroplasmid lineage, which branched off earlier from Babesia/Theileria, reveals intermediate forms, such as the novel dual flip-flop inversion system in Babesia microti and B. rodhaini that generates four distinct genome structures through inversions between two pairs of unique inverted repeats (IR-A and IR-B) [2].

Circular and Other Mitochondrial Genome Forms

While linear mitochondrial genomes dominate in certain parasite groups, circular and other unconventional forms also exist. Trypanosoma brucei, a kinetoplastid parasite, possesses a complex mitochondrial genome known as kinetoplast DNA (kDNA) organized as a catenated network of thousands of mini- and maxicircles [4] [5]. Recent research has revealed an intriguing phenomenon in T. brucei mitochondria: the presence of circular mRNAs [4]. These covalently closed circular RNAs represent a distinct subpopulation of mitochondrial mRNAs with different tail characteristics and UTR lengths compared to their linear counterparts, potentially representing a novel regulatory mechanism or degradation pathway [4].

Table 2: Diversity of Mitochondrial Genome Structures Across Parasite Taxa

Parasite Group	Genome Structure	Size Range	Gene Content	Notable Features
Plasmodium spp.	Concatenated linear	~6 kb	3 PCGs, 27 rRNA fragments	Tandemly repeated [1]
Babesia/Theileria	Linear monomer	6.1-11.1 kb	3 PCGs, rRNA fragments	Terminal inverted repeats [1] [3]
Eimeria spp.	Concatemeric	~6.2 kb	3 PCGs, 20 rRNA fragments	Similar to Plasmodium [2]
Trypanosomatids	Networked circles	Variable	18 PCGs, 2 rRNAs	Catenated mini/maxicircles [4]
Sucking lice	Minichromosomes	Small segments	37 genes total	18 separate chromosomes [5]

The mitochondrial genome of Trypanosoma brucei illustrates another fascinating aspect of mitochondrial diversity. This parasite possesses a complex kinetoplast DNA (kDNA) consisting of a catenated network of thousands of mini- and maxicircles [4]. Unlike the highly reduced mitochondrial genomes of apicomplexan parasites, kDNA contains 20 genes: 2 rRNAs and 18 mRNAs that mostly code for proteins of mitochondrial electron transport chain complexes [4]. The discovery of circular mRNAs in T. brucei mitochondria adds another layer of complexity to our understanding of mitochondrial genome expression and regulation in parasites [4].

Applications in Parasite Taxonomy and Phylogenetics

Mitochondrial genomes have become indispensable tools for parasite classification and evolutionary studies due to their unique characteristics, including maternal inheritance, high mutation rates, and conserved gene content [6]. The application of mitochondrial DNA in taxonomic identification, particularly through DNA barcoding using the cytochrome c oxidase I (COI) gene, has revolutionized species identification and enabled rapid classification across diverse parasitic taxa [6].

The use of complete mitochondrial genomes significantly enhances phylogenetic resolution compared to single-gene approaches. For example, in haemosporidian parasites, the standard 480-bp cytb barcode fragment has limitations in resolving mixed infections and co-infections, which are common in wildlife [7]. The complete mitochondrial genome (~6 kb) provides substantially more informative sites, yielding well-supported phylogenies and enabling more accurate species delimitation [7]. This approach has been successfully applied to various parasite groups, including the identification and classification of novel Didymozoidae parasites infecting yellowfin tuna [8] and the resolution of evolutionary relationships among Theileria species [3].

Mitochondrial genomes also provide insights into adaptive evolution and functional constraints across parasite lineages. Comparative analyses of evolutionary rates in protein-coding genes have revealed patterns of selection pressure related to parasite life history strategies [9]. Additionally, the presence or absence of specific mitochondrial genes, such as those encoding canonical Complex I of the electron transport chain, provides valuable phylogenetic markers that trace major evolutionary transitions [9].

Experimental Protocols for Mitochondrial Genome Analysis

PacBio HiFi Protocol for Haemosporidian Mitochondrial Genomes

Principle: This protocol utilizes long-read PacBio HiFi sequencing to generate high-fidelity mitochondrial genome sequences from haemosporidian parasites, enabling accurate detection of mixed infections and co-infections [7].

Procedure:

• Primer Design: Design barcoded oligonucleotides targeting approximately 6 kb of the haemosporidian mitochondrial genome using modified AE170 (forward) and AE171 (reverse) primers [7].
• Library Preparation:
  • Amplify mitochondrial targets from different specimens using barcoded primers
  • Multiplex amplicons in an SMRTbell library preparation
  • Sequence on PacBio HiFi platform (10-25 kb read length, ~99.5% accuracy)
• Data Analysis:
  • Process reads using HmtG-PacBio Pipeline
  • Apply machine learning algorithm (modified variational autoencoders) and clustering methods to identify mitochondrial haplotypes/species
  • Validate haplotypes against custom database of known mitochondrial sequences

Key Considerations:

• Minimum recommended coverage: 30X per haplotype based on error rates (average 0.2% per read)
• Enables detection of lineages present in very low parasitaemia
• Effectively resolves mixed infections and co-infections that standard cytb barcoding cannot detect [7]

Mitochondrial Genome Assembly from Short-Read Data

Principle: This approach extracts mitochondrial genome sequences from existing short-read genome sequence datasets, significantly expanding mitochondrial genome coverage across diverse taxa [9].

Procedure:

• Data Mining:
  • Identify mitochondrial sequences in existing nuclear genome assemblies based on size and elevated coverage
  • Extract mitochondrial contigs using customized pipelines (e.g., https://github.com/JFWolters/IdentifyMitoContigs)
• De Novo Assembly:
  • Reassemble mitochondrial genomes from raw reads using specialized assemblers (plasmidSPAdes or NOVOPlasty)
  • Assess assembly quality based on contiguity, completeness, and gene content
• Annotation and Validation:
  • Annotate protein-coding genes, rRNAs, and other features using MITOS web server
  • Validate assembly through comparison with known mitochondrial references

Applications: This method has been successfully applied to greatly expand the number of available yeast mitochondrial genomes, facilitating comparative studies of genome evolution across the subphylum Saccharomycotina [9].

Figure 1: Workflow for Mitochondrial Genome Analysis in Parasite Taxonomy

Oxford Nanopore Technologies for Parasite Genomics

Principle: This approach utilizes Oxford Nanopore Technologies (ONT) MinION sequencing to generate reference-quality mitochondrial genomes from parasitic nematodes using only long-read data [10].

Procedure:

• DNA Extraction: Modified protocols allow whole genome sequencing from single parasitic nematodes
• Library Preparation: Prepare libraries using ONT kits without fragmentation
• Sequencing: Sequence on MinION flow cells
• Assembly and Polishing:
  • Perform de novo assembly using only MinION data
  • Optionally polish with Illumina data for improved gene-level accuracy

Performance: Assemblies generated using only MinION data show similar or superior contiguity, completeness, and gene content compared to references, with 88.9-97.6% of complete coding sequences identical to those predicted in assemblies polished with Illumina data [10].

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Key Research Reagents for Mitochondrial Genome Studies in Parasites

Reagent/Material	Function	Application Examples
AE170/AE171 Primers	Amplify ~6 kb mitochondrial genome	Haemosporidian parasite mitochondrial genome amplification [7]
SMRTbell Library Prep Kit	PacBio long-read library preparation	High-fidelity mitochondrial genome sequencing [7]
TIANGEN Marine Animal Tissue DNA Kit	DNA extraction from parasite tissues	Mitochondrial DNA isolation from fish parasites [8]
NEBNext Ultra DNA Library Prep Kit	Illumina short-read library preparation	Mitochondrial genome sequencing from various parasites [8]
MITOS Web Server	Mitochondrial genome annotation	Automated annotation of mitochondrial genes and features [8] [3]
Trimmomatic	Quality control of raw sequencing data	Preprocessing of mitochondrial genome sequencing reads [8]
SPAdes/NOVOPlasty	Genome assembly from sequencing reads	De novo mitochondrial genome assembly [9]
IDBA Software	Assembly of mitochondrial genomes	Construction of complete mitochondrial sequences [3]
PROTAC BRD9 Degrader-6	PROTAC BRD9 Degrader-6, MF:C47H56N8O6, MW:829.0 g/mol	Chemical Reagent
Asn-pro-val-pabc-mmae tfa	Asn-pro-val-pabc-mmae tfa, MF:C63H97F3N10O15, MW:1291.5 g/mol	Chemical Reagent

The study of mitochondrial genome diversity in parasites continues to evolve with technological advancements and expanding taxonomic sampling. Future research directions should emphasize multi-genome studies that integrate mitochondrial and nuclear genomic data to provide comprehensive views of species relationships and evolutionary patterns [6]. Such integrative approaches are particularly important for resolving complex taxonomic relationships and understanding the mechanisms of species divergence in parasites.

Methodological innovations in long-read sequencing technologies, such as PacBio HiFi and Oxford Nanopore, are revolutionizing our ability to resolve complex mitochondrial genome structures and detect mixed infections that were previously challenging to characterize [7] [10]. These technologies, combined with advanced computational methods including machine learning approaches for haplotype identification, will enable more accurate and comprehensive characterization of parasite mitochondrial diversity [7].

There remains a critical need to expand mitochondrial genome sequencing to understudied parasite taxa and ecosystems to fully capture the extent of mitochondrial diversity [6]. Current sampling is heavily biased toward medically and veterinarily important species, leaving significant gaps in our understanding of mitochondrial evolution across the full spectrum of parasitic diversity. Filling these gaps will provide crucial insights into the evolutionary origins of the remarkable structural diversity observed in parasite mitochondrial genomes.

In conclusion, the diversity of mitochondrial genome architectures in parasites—from circular molecules to linear monomers, concatemers, and fragmented minichromosomes—provides a rich source of information for taxonomic classification, phylogenetic reconstruction, and understanding evolutionary processes. The continuous development of specialized protocols and reagents for parasite mitochondrial genomics will ensure that this field remains at the forefront of parasitology research, with important applications in disease diagnosis, surveillance, and control.

Within the context of mitochondrial genome assembly for parasite taxonomy, understanding the core set of mitochondrial genes is fundamental. The mitochondrial genome, while variable in size and structure across eukaryotes, maintains a conserved core of genes essential for oxidative phosphorylation and protein translation [11]. This application note details the conserved protein-coding genes and ribosomal RNA (rRNA) fragments found in mitochondrial genomes, with a specific emphasis on features relevant to apicomplexan and nematode parasites. We provide a standardized protocol for the identification and annotation of these core genetic elements, which serve as critical markers for phylogenetic studies and molecular detection assays.

Core Mitochondrial Gene Content

The core mitochondrial gene content consists of a suite of protein-coding genes and rRNA components, though their organization and structure can vary significantly between taxonomic groups.

Conserved Protein-Coding Genes

The essential protein-coding genes (PCGs) in mitochondrial genomes are predominantly subunits of the oxidative phosphorylation (OXPHOS) system complexes [12]. Table 1 summarizes the conserved PCGs and their functions across major parasitic groups.

Table 1: Conserved Mitochondrial Protein-Coding Genes in Parasites

Gene	Protein Complex	Function	Presence in Apicomplexa	Presence in Nematodes
cox1	Cytochrome c oxidase (Complex IV)	Terminal electron acceptor, proton pumping	Yes [3]	Yes [13]
cox3	Cytochrome c oxidase (Complex IV)	Proton channel formation	Yes [3]	Yes [13]
cob (cytb)	Cytochrome bc1 complex (Complex III)	Electron transfer from ubiquinol to cytochrome c	Yes [3]	Yes [13]
nad1	NADH dehydrogenase (Complex I)	Electron entry point, NADH oxidation	No (typically absent)	Yes [13]
nad4	NADH dehydrogenase (Complex I)	Proton translocation	No (typically absent)	Yes [13]
nad5	NADH dehydrogenase (Complex I)	Proton translocation	No (typically absent)	Yes [13]
atp6	ATP synthase (Complex V)	F0 subunit, proton channel	No (typically absent)	Yes [13]
atp8	ATP synthase (Complex V)	F0 subunit, function not fully defined	No (typically absent)	No (often absent) [13]

In apicomplexan parasites, like Theileria velifera, the mitochondrial genome is notably minimal, typically encoding only three PCGs: cox1, cox3, and cob [3]. In contrast, nematode mitochondria contain a larger complement of genes, often including at least 12 PCGs such as cox1-3, cob, nad1-6 and nad4L, and atp6 [13]. The gene atp8 is frequently absent from nematode mitochondrial genomes [13].

Ribosomal RNA (rRNA) Gene Organization

Mitochondrial ribosomal RNAs are crucial for the translation of the aforementioned PCGs within the organelle. A distinctive feature in many protist parasites, including apicomplexans and some green algae, is the fragmentation of the rRNA genes.

• Fragmented rRNAs in Apicomplexa: The mitochondrial genomes of apicomplexan parasites do not contain full-length large subunit (LSU) and small subunit (SSU) rRNA genes. Instead, these are broken into multiple short, fragmented coding modules that are scattered across the genome and are often transcribed individually [3] [14]. For example, the mitochondrial genome of Theileria velifera contains five fragments of the LSU rRNA gene, designated as LSU1, LSU3, LSU4, LSU5, and LSU6 [3].
• Fragmented and Scrambled rRNAs in Green Algae: This phenomenon is also well-documented in chlorophycean green algae like Polytomella parva and Chlamydomonas reinhardtii, where the mitochondrial SSU and LSU rRNA genes are broken into multiple pieces (e.g., 4 SSU and 8 LSU fragments in P. parva) that are scrambled in order and located on both DNA strands [14].
• Standard rRNA Genes in Nematodes: In contrast, nematode mitochondrial genomes typically contain standard, full-length rRNA genes for the small subunit (SSU or rrnS) and large subunit (LSU or rrnL) [13].

Table 2: Mitochondrial rRNA Gene Structure Across Organisms

Organism Group	rRNA Structure	Typical Number of Fragments	Example
Apicomplexan Parasites	Fragmented LSU rRNA	5-8 fragments [3]	Theileria velifera
Chlorophycean Green Algae	Fragmented and scrambled SSU & LSU rRNAs	4 SSU, 8 LSU fragments [14]	Polytomella parva
Nematode Parasites	Conventional, full-length SSU & LSU rRNAs	2 (1 SSU, 1 LSU) [13]	Ascaridia galli
Animals (Bilaterian)	Conventional, full-length SSU & LSU rRNAs	2 (1 SSU, 1 LSU) [12]	Homo sapiens

Experimental Protocol: Identification and Annotation of Core Mitochondrial Genes

This protocol describes a standardized workflow for identifying and annotating the core protein-coding and rRNA genes in a newly sequenced mitochondrial genome from a parasitic organism, using a combination of reference-based and ab initio tools.

The following diagram illustrates the complete experimental and computational workflow for mitochondrial gene annotation.

Step-by-Step Procedure

Materials & Reagents:

• Hardware: High-performance computing server.
• Software: MITOS2 Web Server [15] or local installation; BLAST+ suite; tRNAscan-SE [15]; MFannot; Geneious or Apollo software for manual curation.
• Data: Assembled mitochondrial genome sequence in FASTA format. Custom database of known fragmented mitochondrial rRNAs (compiled from sources like [3] [14]).

Procedure:

• Approximate Gene Finding and Initial Annotation

  • Submit the assembled mitochondrial genome FASTA file to the MITOS2 web server (http://mitos2.bioinf.uni-leipzig.de/) for initial annotation [15].
  • Alternatively, use the MFannot (http://megasun.bch.umontreal.ca/cgi-bin/mfannot/mfannotInterface.pl) tool for a complementary annotation, particularly for non-metazoan parasites [16].
  • These tools use hidden Markov models (HMMs) and BLAST searches against curated reference sets to identify the approximate locations of PCGs and rRNAs [15].

• Precise Boundary Prediction for Protein-Coding Genes

  • The initial HMM/BLAST predictions often lack precise start and stop codons. MITOS2 employs an improved probabilistic method to determine exact gene boundaries.
  • The algorithm considers:
    • Distance (δ): The proximity of a candidate start/stop codon to the approximate position from the HMM.
    • Codon Probability (Ï•): The empirical probability of a codon being used as a start or stop codon, based on frequencies from reference annotations (e.g., RefSeq) for the specific genetic code.
    • Gene Length (λ): The empirical probability of the resulting gene length, based on reference data [15].
  • The final boundaries are chosen by maximizing the product of these three factors: argmax(δ × ϕ × λ) [15].

• Identification of Fragmented rRNA Genes

  • For apicomplexans and related parasites, fragmented rRNA genes may not be fully annotated by standard tools.
  • Perform a BLASTN search of the mitogenome against a custom database of known fragmented mitochondrial rRNA sequences from related organisms (e.g., using sequences from Theileria parva or Polytomella parva as queries) [3] [14].
  • Manually inspect and annotate regions with significant sequence similarity to these rRNA fragments. Confirm their identity by ensuring they correspond to specific domains of the conventional LSU or SSU rRNA secondary structure [14].

• tRNA Gene Annotation

  • Execute tRNAscan-SE on the mitogenome sequence. Use the appropriate search mode (e.g., "Mito" for metazoans, "General" for other groups) to identify tRNA genes and predict their secondary structures and anticodons [15] [13].
  • Note that some parasitic groups, like apicomplexans, may entirely lack mitochondrial-encoded tRNAs, relying instead on nuclear-encoded tRNAs imported into the mitochondrion [3].

• Manual Curation and Validation

  • Visualize the annotation results using software like Geneious or Apollo [17].
  • Manually check and correct annotations, particularly for genes with non-canonical start codons (e.g., TTT, GTG) which are common in nematodes and other invertebrates [15] [13] [11].
  • Validate PCG annotations by examining their translated amino acid sequences for integrity and the absence of internal stop codons (assuming no RNA editing).

The Scientist's Toolkit

Table 3: Essential Research Reagents and Tools for Mitochondrial Gene Analysis

Tool / Reagent	Type	Primary Function	Key Application Note
MITOS2	Web Server / Software	Automated annotation of metazoan mitogenomes	Critical for precise PCG boundary prediction using its probabilistic model; also annotates rRNAs and tRNAs [15].
MFannot	Web Server	Automated annotation of protist mitogenomes	A valuable alternative for non-metazoan parasites, leveraging a different reference set and algorithm [16].
tRNAscan-SE	Software	Detection of tRNA genes	Accurately identifies tRNA genes and their secondary structures, including those with atypical features common in nematodes [15] [13].
BLAST+ Suite	Software	Sequence similarity search	Essential for identifying fragmented rRNA genes using custom databases and for initial homology-based gene finding [3].
Custom rRNA Fragment DB	Database	Curated collection of rRNA sequences	A necessary in-house resource for correctly annotating fragmented rRNAs in apicomplexans and other protists [3] [14].
PREP-Mt	Web Server	Prediction of RNA editing sites	Predicts C-to-U RNA editing sites in plant mitochondrial PCGs; less applicable to parasites but useful for understanding the phenomenon [18].
N1-(1,1-Difluoroethyl)pseudouridine	N1-(1,1-Difluoroethyl)pseudouridine, MF:C11H14F2N2O6, MW:308.24 g/mol	Chemical Reagent	Bench Chemicals
1-(alpha-L-Threofuranosyl)thymine	1-(alpha-L-Threofuranosyl)thymine, MF:C9H12N2O5, MW:228.20 g/mol	Chemical Reagent	Bench Chemicals

Application Note: Structural Variation in Parasite Mitochondrial Genomes

This application note details the critical role of structural variation analysis in mitochondrial genomes for parasite taxonomic classification and biological discovery. Characterizing variations such as gene order changes, gene absences, and the presence of terminal inverted repeats (TIRs) provides powerful insights into evolutionary relationships, genomic adaptation to parasitism, and potential drug targets.

Structural Variations as Taxonomic and Phylogenetic Markers

Mitochondrial (mt) genes are increasingly used as aids in phylogenetic and epidemiologic analyses of parasites due to their general lack of sexual recombination and uniparental inheritance [11]. The gene order and content of mitochondrial genomes can provide a strong phylogenetic signal, especially in complex evolutionary groups affected by rapid radiation or hybridization [19]. For instance, a study on three closely related oak species demonstrated that mitochondrial protein-coding genes (PCGs) could robustly resolve their phylogenetic relationships, forming a distinct clade separate from other species [19]. In parasites, comparative analysis of the linear mitochondrial genome of Theileria velifera (6,125 bp), which contains Terminal Inverted Repeats (TIRs) at both ends, helped clarify its evolutionary position and close relationship to T. annulata and T. parva [20]. Similarly, the highly unusual organization of kinetoplastid mtDNA, comprising catenated maxicircle and minicircle DNAs, is a defining characteristic of this protozoan group [11].

Genomic Adaptations in Parasitic Lineages

Adoptation to a parasitic lifestyle can leave distinct marks on the mitochondrial genome, most notably through gene loss. A comparative genomic study of the free-living red alga Gracilariopsis andersonii and its parasite Gracilariophila oryzoides revealed that the ATP8 and SDHC genes, which encode essential proteins, had become pseudogenes in the parasite's mitochondrial genome [21]. This finding indicates that these genes are no longer critical in the parasite's mitochondria, a conclusion supported by the observation that a parasite from a different class of red algae, Plocamiocolax puvinata, has lost the atp8 gene entirely [21]. Furthermore, nonadaptive processes like genetic drift, influenced by transmission mode, can impact genome architecture. Research on microsporidia has shown that vertical transmission is associated with larger genomes and a higher proportion of transposable elements (TEs), suggesting that population bottlenecks reduce the effectiveness of natural selection in purging mildly deleterious TE insertions [22].

Table 1: Documented Structural Variations in Parasite Mitochondrial Genomes

Parasite/Group	Variation Type	Specific Genomic Change	Functional/Taxonomic Implication	Citation
Theileria velifera (Apicomplexa)	TIRs	Linear mitochondrial genome with TIRs at both ends	Genome structure characteristic of many Apicomplexan parasites; used in phylogenetic analysis	[20]
Red Algal Adelphoparasites (e.g., Gracilariophila oryzoides)	Missing Genes / Pseudogenization	ATP8 and SDHC genes rendered pseudogenes	Loss of gene function as adaptation to parasitic lifestyle	[21]
Kinetoplastid Protozoa (e.g., Trypanosomatids)	Gene Order & Organization	Catenated maxicircle and minicircle DNA (kDNA)	Defining genomic feature requiring extraordinary gene expression mechanisms	[11]
Microsporidia	Transposable Element (TE) Abundance	Positive correlation between TE amount and genome size	Larger, TE-rich genomes associated with vertical transmission and genetic drift	[22]

Protocols for Analyzing Mitochondrial Structural Variations

Protocol 1: Assembling and Annotating Parasite Mitochondrial Genomes

This protocol outlines the steps for generating a high-quality mitochondrial genome assembly, which is a prerequisite for all downstream structural variation analyses. The method leverages third-generation sequencing technologies to overcome challenges posed by repetitive sequences and complex structural variants [19].

Materials and Equipment

• High-molecular-weight (HMW) genomic DNA from parasite tissue or cells.
• TRIzol Reagent or modified CTAB method for DNA extraction.
• PacBio HiFi or Oxford Nanopore long-read sequencing platform.
• High-fidelity DNA polymerase (e.g., Takara LA Taq Polymerase).
• Computational resources and software: Geneious Pro, MITOS web server, BLAST, Flye assembler.

Step-by-Step Procedure

• DNA Extraction: Extract HMW genomic DNA from flash-frozen parasite tissue using a method optimized for the organism (e.g., modified CTAB for plants, phenol/chloroform for other eukaryotes) [21] [19]. Assess DNA integrity via pulsed-field gel electrophoresis (fragments >50 kb desired) and purity via spectrophotometry (A260/A280 ≈ 1.8).
• Library Preparation and Sequencing: Prepare a SMRTbell library (e.g., using PacBio SMRTbell prep kit 3.0) from 5-7 μg of sheared HMW DNA (15-20 kb fragments). Sequence the library on a PacBio HiFi platform to generate long, accurate reads [19] [23].
• Genome Assembly: Assemble the mitochondrial genome from the sequencing reads using a long-read assembler like Flye, which is optimized for resolving complex repetitive structures [19]. Evaluate the completeness and continuity of the assembly by examining the uniformity of read coverage across the entire genome.
• Genome Annotation: Annotate the assembled mt-genome using the MITOS web server and BLAST against databases of known mitochondrial gene sequences [20]. Identify protein-coding genes (PCGs), ribosomal RNA (rRNA) gene fragments, and tRNA genes (using tRNAscan-SE). Manually verify and correct annotations.

Protocol 2: Detecting Gene Order Rearrangements and Missing Genes

This protocol describes a comparative genomics approach to identify structural variants by aligning a newly assembled mitochondrial genome to a reference.

Materials and Equipment

• Assembled and annotated mitochondrial genome (from Protocol 1).
• Reference mitochondrial genome(s) from related species.
• Computational software: MEGA 11.0, Clustal W, genome visualization tools (e.g., Gene Graphics).

Step-by-Step Procedure

• Data Collection: Obtain mitochondrial genome sequences of closely related parasite species from public databases (e.g., GenBank). These will serve as references for comparison.
• Synteny and Collinearity Analysis: Use genome visualization and alignment software to perform a global comparison between the newly assembled genome and the reference genomes. Look for large-scale inversions, translocations, and other rearrangements by assessing changes in gene order and orientation [23].
• Identification of Missing Genes: Compare the full complement of annotated genes from the target genome against the reference. Identify genes that are completely absent or appear as pseudogenes (e.g., due to frameshifts or premature stop codons), as was done for the atp8 and SDHC genes in red algal parasites [21].
• Phylogenetic Inference: Use nucleotide sequences of conserved mitochondrial PCGs (e.g., cox1, cob) from multiple species. Align the sequences using Clustal W in MEGA 11.0 software. Construct a phylogenetic tree (e.g., using Maximum Likelihood with 1,000 bootstrap replicates) to assess evolutionary relationships and correlate them with the observed structural variations [20].

Protocol 3: Identifying and Analyzing Terminal Inverted Repeats (TIRs)

This protocol focuses on the computational identification and characterization of TIRs, which are associated with certain linear mitochondrial genomes and TIR transposons.

Materials and Equipment

• Assembled mitochondrial genome contig(s).
• Computational tools: BLAST, custom scripts or pipelines for TIR identification (e.g., as developed for plant TIR transposon analysis [24]).

Step-by-Step Procedure

• Contig End Analysis: For each linear mitochondrial contig, extract the sequences from the 5' and 3' termini (e.g., first and last 500-1000 bp, depending on expected TIR size).
• Sequence Alignment and Inversion: Perform a local pairwise alignment of the terminal sequences. Also, create the reverse complement of one terminal sequence and align it with the other. Significant identity in the direct comparison indicates the presence of TIRs, as seen in the Theileria velifera mitochondrial genome [20].
• Characterization: Document the length and percentage identity of the TIRs. Investigate whether the TIRs are associated with known transposable elements, such as those from the Mutator-like element (MULE) superfamily, which is prevalent in plants [24].

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Reagents and Tools for Mitochondrial Structural Variation Analysis

Reagent / Tool	Function / Application	Example Use Case
PacBio HiFi Sequencing	Generates long, high-fidelity reads for assembling complex genomic regions.	Resolving repetitive structures and obtaining complete mt-genome assemblies without gaps [19].
MITOS Web Server	Automated annotation of mitochondrial genomes.	Rapid identification and annotation of protein-coding, rRNA, and tRNA genes in a new assembly [20].
MEGA 11.0 Software	Integrated tool for sequence alignment, evolutionary genetics, and phylogenetic reconstruction.	Aligning mitochondrial PCGs and constructing phylogenetic trees to place parasites in an evolutionary context [20].
Clustal W Algorithm	Multiple sequence alignment of nucleotide or protein sequences.	Creating accurate alignments of conserved genes (e.g., cox1, cob) for comparative and phylogenetic analysis [20].
Custom TIR Identification Pipeline	Systematically identifies active autonomous TIR transposons and TIR sequences in genomes.	Characterizing the structure of linear mitochondrial genomes and assessing the activity of TIR transposons [24].
[Lys8] Vasopressin Desglycinamide	[Lys8] Vasopressin Desglycinamide, MF:C44H63N11O12S2, MW:1002.2 g/mol	Chemical Reagent
Galectin-3 antagonist 2	Galectin-3 Antagonist 2	Galectin-3 Antagonist 2 is a high-affinity inhibitor for research into fibrosis, cancer, and metabolic diseases. This product is For Research Use Only. Not for human consumption.

Workflow and Pathway Visualizations

Mitochondrial Genome Assembly and Analysis Workflow

The diagram below outlines the comprehensive workflow for analyzing structural variations in parasite mitochondrial genomes, from sample preparation to biological insight.

Structural Variations in Parasite Mitochondrial Genomes

This diagram categorizes the primary types of structural variations discussed in this note and their relevance to parasite research.

The Role of Mitochondrial Genomes in Resolving Taxonomic Complexities and Evolutionary Relationships

Mitochondrial genomes (mtDNA) have become an invaluable tool in species classification and evolutionary studies, offering a powerful means to resolve complex taxonomic relationships [6]. Their unique characteristics, including maternal inheritance, high copy number per cell, and relatively high mutation rates compared to nuclear DNA, make them particularly suitable for phylogenetic analysis and species identification [6] [25]. In parasite taxonomy research, where morphological distinctions are often subtle or variable, mtDNA provides a robust molecular framework for delineating species boundaries and understanding evolutionary histories [6] [26].

The field of mitochondrial genomics has evolved from early biochemical studies to its current role in biodiversity assessment, with mitochondrial DNA barcoding—particularly utilizing the cytochrome c oxidase I (COI) gene—revolutionizing species identification across diverse taxa [6]. This approach enables rapid and accurate classification, which is crucial for understanding parasite systematics, host-parasite co-evolution, and for informing drug development strategies targeting taxonomically-defined groups.

Key Features of Mitochondrial Genomes for Taxonomic Resolution

Biological Properties and Advantages

Table 1: Characteristics of Mitochondrial Genomes Supporting Taxonomic Applications

Feature	Description	Taxonomic Utility
Maternal Inheritance	Generally uniparental inheritance without recombination [6]	Simplifies tracing of evolutionary lineages
High Copy Number	Hundreds to thousands of copies per cell [27]	Enables analysis from minimal or degraded samples
High Mutation Rate	5-10 times higher than nuclear genome [25]	Provides resolution for recently diverged taxa
Compact Structure	Small size with minimal non-coding DNA [25]	Facilitates efficient sequencing and analysis
Conserved Gene Content	13 protein-coding, 22 tRNA, and 2 rRNA genes in animals [25]	Enables consistent comparative analyses across diverse taxa

Mitochondrial genomes exhibit several biological properties that make them exceptionally useful for taxonomic studies. The lack of recombination and predominantly uniparental inheritance simplifies evolutionary analysis by reducing complexity in tracing lineages through evolutionary time [6] [25]. The higher mutation rate of mtDNA compared to nuclear DNA allows for the accumulation of genetic differences between recently diverged species, making it possible to distinguish even closely-related taxa [25].

The conserved gene content across metazoans provides a consistent framework for comparison, while variable regions offer characters for distinguishing taxa [25]. For parasite taxonomy, these properties are particularly valuable when working with small specimens, archived materials, or environmental samples where DNA quantity or quality may be limiting.

Mitochondrial DNA Barcoding and Beyond

Mitochondrial DNA barcoding using the COI gene has emerged as a standardized approach for species identification [6]. This method leverages the fact that intraspecific variation in COI sequences is generally low compared to interspecific divergence, creating a "barcode gap" that enables species discrimination. However, comprehensive mitochondrial genome analysis provides additional resolution through:

• Complete gene set analysis: Utilizing all 13 protein-coding genes for robust phylogenetic reconstruction
• Structural features: Gene order and genome rearrangements as taxonomic characters [26]
• RNA genes: Sequence variation in tRNA and rRNA genes [26]
• Non-coding regions: Including control region and intergenic spacers [25]

For parasitic taxa, complete mitogenome analysis has proven particularly valuable in resolving complexes of cryptic species—morphologically similar but genetically distinct organisms that may differ in host specificity, pathogenicity, or drug susceptibility [26].

Data Analysis and Interpretation in Taxonomic Contexts

Analytical Frameworks and Quantitative Measures

Table 2: Key Metrics and Thresholds in Mitochondrial Taxonomy

Metric	Calculation/Description	Interpretation Guidelines
Genetic Distance	Proportion of nucleotide differences between sequences	Higher values indicate greater evolutionary divergence
Nonsynonymous/Synonymous Substitution Ratio (dN/dS)	Ratio of amino acid-changing to silent substitutions	dN/dS < 1 suggests purifying selection; dN/dS > 1 suggests positive selection [25]
Neutrality Index	Measures direction of selection on amino acid variation [25]	Values >1 indicate excess amino acid polymorphisms relative to neutral expectations
Heteroplasmy Level	Proportion of variant mtDNA molecules in an individual [27]	May complicate species delimitation if high
Haplogroup Diversity	Group of similar haplotypes sharing a common ancestor	Defines evolutionary lineages within and between species

The interpretation of mitochondrial sequence data for taxonomic purposes requires careful consideration of evolutionary forces shaping genetic variation. The neutrality assumption—that most mtDNA variation is selectively neutral—has been challenged by research demonstrating complex interactions of selective pressures [25]. Analyses such as the McDonald-Kreitman test can distinguish neutral evolution from selection by comparing ratios of synonymous and nonsynonymous substitutions within and between species [25].

For species delimitation, multiple analytical approaches should be employed, including:

• Distance-based methods: Using genetic distance thresholds specific to the taxonomic group
• Tree-based methods: Identifying monophyletic groups in phylogenetic trees
• Character-based methods: Identifying fixed diagnostic characters
• Population genetic methods: Assessing gene flow and reproductive isolation

Addressing Analytical Challenges

Several challenges complicate the interpretation of mitochondrial data for taxonomy:

• Nuclear mitochondrial DNA segments (NUMTs): Nuclear copies of mitochondrial sequences can be co-amplified, leading to erroneous interpretations [6]. Careful bioinformatic filtering and validation are essential.
• Incomplete lineage sorting: Retention of ancestral polymorphism can cause discordance between gene trees and species trees [6].
• Hybridization and introgression: Transfer of mtDNA between species can obscure phylogenetic relationships [6].
• Heteroplasmy: Presence of multiple mtDNA variants within individuals [27] [25].
• Variable mutation rates: Differences in evolutionary rates across lineages can affect divergence time estimates [25].

To address these challenges, integrative taxonomy combining mitochondrial data with nuclear markers, morphological characters, ecological data, and other lines of evidence provides the most robust framework for species delimitation [26].

Experimental Protocol: Mitochondrial Genome Assembly for Parasite Taxonomy

Sample Preparation and DNA Extraction

Materials Required:

• Tissue samples (preserved in ethanol, RNAlater, or frozen)
• DNA extraction kit (e.g., E.Z.N.A. Mollusc DNA Kit for delicate specimens) [26]
• Proteinase K and appropriate digestion buffers
• RNase A for RNA removal
• Quantification system (e.g., Qubit fluorometer)

Procedure:

• Sample Lysis: Digest 10-25 mg tissue in lysis buffer with Proteinase K at 56°C for 4-24 hours (depending on tissue type)
• DNA Extraction: Follow manufacturer's protocol for DNA purification
• RNA Removal: Treat with RNase A (10-20 μg/mL) at room temperature for 10 minutes
• DNA Quantification: Measure DNA concentration using fluorometric methods
• Quality Assessment: Evaluate DNA integrity by agarose gel electrophoresis or fragment analyzer

Technical Notes:

• For ancient or degraded samples, specialized ancient DNA protocols should be used
• Minimize handling to reduce contamination risk
• Extract negative controls alongside samples to monitor contamination

Mitochondrial Genome Sequencing

Materials Required:

• PCR reagents for initial amplification (if using PCR-based approach)
• Library preparation kit appropriate for sequencing platform
• Next-generation sequencing platform (Illumina, PacBio, or Oxford Nanopore)
• Sanger sequencing reagents for gap filling and validation

Procedure for Shotgun Sequencing Approach:

• Library Preparation: Fragment DNA to appropriate size (300-800 bp for Illumina) and prepare sequencing library following manufacturer's protocol
• Enrichment (Optional): Use mitochondrial enrichment strategies (e.g., hybridization capture, PCR amplification) for samples with low mitochondrial content
• Sequencing: Perform sequencing on appropriate platform to achieve sufficient coverage (≥50× for mitochondrial genome)
• Validation: Design primers for gap filling and problematic regions for Sanger sequencing [26]

Procedure for Long-Read Sequencing:

• High Molecular Weight DNA Extraction: Use gentle extraction methods to preserve long DNA fragments
• Library Preparation: Prepare library without fragmentation for PacBio or Nanopore sequencing
• Sequencing: Perform sequencing to achieve sufficient coverage for assembly

Bioinformatics Analysis Pipeline

Materials Required:

• High-performance computing resources
• Bioinformatics software (see Research Toolkit section)
• Reference mitochondrial genomes

Procedure:

• Quality Control
  • Assess raw read quality using FastQC [26]
  • Trim adapters and low-quality bases using Cutadapt [26]

• Genome Assembly

  • For short reads: Perform de novo assembly using NOVOPlasty or other mitogenome-specific assemblers [26]
  • For long reads: Assemble using Canu, Flye, or other long-read assemblers
  • For hybrid approaches: Combine short and long read data

• Assembly Validation

  • Check for circularization of mitochondrial genome
  • Assess coverage uniformity
  • Verify absence of NUMTs by checking for unusual features (e.g., stop codons, indels)

• Genome Annotation

  • Annotate using MITOS2 or similar automated annotation tools [26]
  • Manually curate gene boundaries
  • Identify tRNA genes using tRNAscan-SE [26]

• Phylogenetic Analysis

  • Extract and align protein-coding genes
  • Construct phylogenetic trees using maximum likelihood or Bayesian methods
  • Perform molecular dating analyses if appropriate

Table 3: Research Reagent Solutions for Mitochondrial Genome Analysis

Category	Specific Products/Tools	Application Notes
DNA Extraction	E.Z.N.A. Mollusc DNA Kit [26], DNeasy Blood & Tissue Kit	Reliable yields from various sample types including difficult tissues
Amplification	Long-range PCR kits (e.g., LA Taq), Whole genome amplification	For enriching mitochondrial DNA from limited samples
Sequencing	Illumina kits (NovaSeq, MiSeq), PacBio SMRTbell, Oxford Nanopore ligation kits	Selection depends on required read length, accuracy, and budget
Assembly Software	NOVOPlasty [26], MITOS2 [26], Geneious, CLC Genomics Workbench	Specialized mitogenome assemblers outperform general tools
Annotation Tools	MITOS2 [26], tRNAscan-SE [26], ARWEN, DOGMA	Automated annotation with manual curation essential
Analysis Platforms	PhyloSuite [26], MEGA, BEAST, R with phylogenetic packages	Streamlined analysis workflows improve reproducibility

Mitochondrial genome analysis has transformed approaches to resolving taxonomic complexities and understanding evolutionary relationships in parasitic organisms. The protocols and analytical frameworks outlined here provide a roadmap for researchers engaged in parasite taxonomy, with applications spanning basic systematics to drug discovery pipelines.

Future developments in the field will likely focus on single-cell mitochondrial genomics, enabling analysis of individual parasites without cultivation; environmental DNA barcoding for detecting parasitic organisms in complex samples; and integration with multi-omics approaches to connect taxonomic identity with functional capacity. As sequencing technologies continue to advance and analytical methods become more sophisticated, mitochondrial genomics will remain a cornerstone of parasitology research, providing critical insights into the diversity, evolution, and biological characteristics of economically and medically significant parasites.

For researchers in drug development, accurately defining taxonomic boundaries through mitochondrial genomics provides the essential foundation for understanding distribution patterns, host specificity, and evolutionary trajectories of target species—all critical considerations for designing effective control strategies.

Within the broader scope of a thesis on mitochondrial genome assembly for parasite taxonomy, this application note details the experimental and bioinformatic protocols for characterizing the mitochondrial genome of Theileria velifera.

Theileria velifera is a tick-borne apicomplexan parasite that infects bovines, leading to economic losses in the livestock industry [3]. Precise parasite identification is crucial for disease control, yet traditional methods based on morphology can be subjective and limited [3]. Mitochondrial (mt) genomes, with their higher evolutionary rate compared to nuclear DNA, provide a powerful molecular marker for phylogenetic studies and taxonomic resolution [3] [28]. This case study demonstrates how the complete mt genome of T. velifera was sequenced, assembled, and analyzed to elucidate its phylogenetic placement among apicomplexan parasites.

Results and Data Analysis

General Features of theT. veliferaMitochondrial Genome

The complete mitochondrial genome of T. velifera was sequenced, assembled, and deposited in GenBank under accession number ON684327 [3]. Key characteristics are summarized below.

Table 1: General Features of the Theileria velifera Mitochondrial Genome

Feature	Characteristic
Genome Structure	Linear monomer [3]
Total Length	6,125 bp [3]
Protein-Coding Genes (PCGs)	3 genes: cox1, cob (cyt b), and cox3 [3]
rRNA Genes	5 large subunit (LSU) rRNA gene fragments (LSU1, LSU3, LSU4, LSU5, LSU6) [3]
Transfer RNA (tRNA) Genes	None identified [3]
Terminal Structures	Terminal Inverted Repeats (TIRs) at both ends [3]

Table 2: Nucleotide Composition and Skewness of the T. velifera Mitochondrial Genome

Parameter	Value/Calculation
AT-skew	(A - T) / (A + T) [3]
GC-skew	(G - C) / (G + C) [3]
The specific values for T. velifera were not explicitly detailed in the search results.

Table 3: Start and Stop Codon Usage in Theileria and Babesia spp.

Codon Type	Prevalent Codons
Start Codons	ATN, GTN, TTN [3]
Stop Codons	TAA, TAG, TGA [3]

Phylogenetic Placement

Phylogenetic analysis was conducted to resolve the evolutionary relationships of T. velifera within the apicomplexan parasites.

• Gene Selection: The analysis used concatenated amino acid sequences of two protein-coding genes, cox1 and cob [3]. The cox3 gene was omitted due to its significant sequence variation among Theileria and Babesia species [3].
• Method: The maximum likelihood (ML) method was employed with 1,000 bootstrap replicates, based on JTT and Freqs models [3].
• Finding: The analysis confirmed that T. velifera is closely related to T. annulata, T. parva, T. taurotragi, and T. lestoquardi [3].

The following workflow diagram illustrates the complete process from sample to phylogenetic insight:

Experimental Protocols

Sample Collection, DNA Extraction, and Sequencing

This protocol outlines the steps for obtaining high-quality genomic DNA from T. velifera-infected host blood for mitochondrial genome sequencing [3].

I. Sample Collection

• Collect venous blood from bovine hosts using sterile EDTA anticoagulant tubes [3].
• Store samples at -20°C prior to DNA extraction [3].
• Screen samples for T. velifera infection using established methods (e.g., amplification and sequencing of the 18S rRNA gene) to confirm positive status [3].

II. Genomic DNA Extraction

• Use a commercial genomic DNA extraction kit (e.g., TIANamp Genomic DNA Kit).
• Process 200 µL of EDTA-anticoagulated blood according to the manufacturer's instructions [3].
• Store eluted DNA at -20°C.

III. Library Preparation and Sequencing

• Prepare an Illumina paired-end library from the extracted DNA.
• Sequence the library on an Illumina Novoseq 6000 platform to generate high-throughput reads [3].

Genome Assembly and Annotation Protocol

This protocol describes the bioinformatic workflow for assembling the mitochondrial genome from raw sequencing reads and identifying its features [3].

I. Data Preprocessing

• Filter raw sequencing reads to remove low-quality sequences, generating a set of "clean reads" for assembly.

II. Genome Assembly

• Assemble the clean reads de novo using IDBA software [3]. This step reconstructs the mitochondrial genome without a reference.

III. Genome Annotation

• Annotate the assembled genome using the MITOS web server to identify protein-coding genes (PCGs) and rRNA gene fragments [3].
• Validate PCG annotations by performing a BLAST search against the GenBank database [3].
• Annotate rRNA gene fragments by comparing them to previously reported rRNA genes of related species (e.g., T. parva, T. orientalis) [3].
• Search for tRNA genes using tRNAscan-SE v.2.0 [3]. The result for T. velifera and other apicomplexans is typically the absence of tRNA genes.

IV. Data Analysis

• Calculate nucleotide composition and skewness using formulas:
  • AT-skew = (A - T) / (A + T)
  • GC-skew = (G - C) / (G + C) [3].
• Analyze relative synonymous codon usage (RSCU) using software like CodonW 1.4.2 [3].

Phylogenetic Analysis Protocol

This protocol details the procedure for constructing a phylogenetic tree to determine the evolutionary relationships of T. velifera [3].

I. Data Retrieval and Selection

• Retrieve nucleotide sequences of the cox1 and cob genes from a set of Apicomplexan parasites from GenBank.
• Omit the cox3 gene due to its high variability, which can confound deep-level phylogenetic analysis [3].

II. Sequence Alignment

• Concatenate the cox1 and cob gene sequences for each taxon.
• Align the concatenated nucleotide sequences using Clustal W as implemented in MEGA 11.0 software.
• Translate the aligned nucleotide sequences into corresponding amino acid sequences for analysis.

III. Tree Construction

• Use the maximum likelihood (ML) method in MEGA 11.0 to construct the phylogeny.
• Perform ML analysis with 1,000 bootstrap replicates to assess branch support, using models such as JTT with empirical base frequencies (Freqs) [3].
• Visualize the final tree using a tool like iTOL (Interactive Tree Of Life) [3].

The Scientist's Toolkit: Research Reagent Solutions

Table 4: Essential Reagents and Tools for Mitochondrial Genome Analysis

Reagent/Software	Function/Application
TIANamp Genomic DNA Kit	Extraction of high-quality total genomic DNA from blood samples [3].
Illumina NovaSeq 6000 Platform	High-throughput sequencing to generate raw genomic reads [3].
IDBA Software	De novo assembly of clean reads into contiguous sequences (contigs) [3].
MITOS Web Server	Automated annotation of the mitochondrial genome [3].
tRNAscan-SE	Identification of transfer RNA genes in genomic sequences [3].
MEGA 11.0 Software	Integrated suite for molecular evolutionary genetics analysis, including alignment and phylogenetic tree construction [3].
CodonW	Calculation of relative synonymous codon usage (RSCU) statistics [3].
A2A receptor antagonist 3	A2A Receptor Antagonist 3|Research Grade
2'-Deoxytubercidin 5'-triphosphate	2'-Deoxytubercidin 5'-triphosphate, MF:C11H17N4O12P3, MW:490.19 g/mol

Methodological Workflows: From Sample to Annotated Mitochondrial Genome

Best Practices for Parasite Sample Collection and High-Quality DNA Extraction

High-quality DNA is a prerequisite for successful downstream applications such as mitochondrial genome assembly, a cornerstone of modern parasite taxonomy and phylogenetic research [3]. The integrity of this genetic material is profoundly influenced by initial sample collection strategies and the DNA extraction methods employed. Suboptimal procedures can lead to degraded DNA, contaminant carryover, and inhibitor presence, ultimately compromising the reliability of genomic data [29]. This document outlines standardized, effective protocols for the collection of parasite samples and the subsequent extraction of high-quality DNA, contextualized within the framework of mitochondrial genomics.

Best Practices for Parasite Sample Collection

Proper sample collection and preservation are critical first steps to ensure the integrity of parasite DNA for mitochondrial genome sequencing.

• Sample Types: Common samples include whole ticks, blood from infected hosts, or isolated parasite stages [29] [3].
• Preservation Method: Immediate preservation is essential to prevent DNA degradation. Samples should be stored in absolute ethanol at 4°C or -20°C [29]. It is critical to ensure samples remain fully submerged to avoid evaporation, which can compromise DNA quality over long-term storage.
• Documentation: Maintain meticulous records of sample origin, host species, date, and geographical location for accurate taxonomic classification.

DNA Extraction Methodologies

The selection of a DNA extraction method balances cost, simplicity, and the requirements for downstream analytical success. The following table summarizes the performance of different methods evaluated on challenging, sub-optimally stored Ixodes ricinus ticks [29].

Table 1: Comparison of DNA Extraction Methods for Parasite Samples

Method	Key Steps	Avg. A260/280 Purity	Median DNA Yield (ng)	Inhibition (qPCR)	Relative Cost
Ammonium Hydroxide (Intact Tick) [29]	Incubation of intact tick in NH₄OH at 99°C [29].	~1.44 [29]	151 [29]	None detected [29]	Very Low
Ammonium Hydroxide (Crushed Tick) [29]	Homogenization with beads, then NHâ‚„OH hydrolysis [29].	~1.44 [29]	151 [29]	9/50 samples inhibitory [29]	Very Low
QIAGEN Blood & Tissue Kit [29]	Bead-beating homogenization, extended enzymatic lysis, silica-membrane purification [29].	1.63 [29]	151 [29]	None detected [29]	High
QIAGEN Mini Kit [29]	Bead-beating homogenization, silica-membrane purification [29].	~1.44 [29]	151 [29]	None detected [29]	Medium

Detailed Protocols

Protocol A: Ammonium Hydroxide Hydrolysis (for Intact Ticks)

This cheap and simple method is sufficient for qPCR-based pathogen detection and is as effective as commercial kits for this purpose [29].

• Placement: Transfer a single, intact tick into a 1.5 mL microcentrifuge tube.
• Lysis: Add 100 µL of a 0.7 M ammonium hydroxide solution to the tube.
• Incubation: Incubate the tube at 99°C for 30 minutes.
• Evaporation: Place the tube in a heating block at 99°C with the lid open for 10-15 minutes to evaporate the ammonia.
• Resuspension: Resuspend the resulting lysate in 100 µL of nuclease-free water or TE buffer.
• Storage: Store the extracted DNA at -20°C.

Protocol B: Silica-Membrane Kit Protocol (with Homogenization)

This method is recommended for applications requiring higher purity DNA, such as next-generation sequencing for mitochondrial genome assembly [29] [3].

• Homogenization: Place the individual tick in a 2 mL tube containing 180 µL of the provided lysis buffer (e.g., RLT from QIAGEN) and three 2.5 mm stainless steel beads.
• Disruption: Agitate the tube for 5 minutes at 10,000 Hz using a bead-mill homogenizer.
• Lysis: Following homogenization, extend the enzymatic lysis step to 16 hours (overnight) to ensure complete disruption of tough parasite structures.
• Purification: Follow the manufacturer's instructions for binding, washing, and eluting the DNA on the silica membrane.
• Elution: Elute the DNA in two successive 50 µL volumes of elution buffer to maximize final yield, resulting in a total volume of 100 µL [29].
• Storage: Store the purified DNA at -20°C or -80°C for long-term preservation.

Workflow for Mitochondrial Genome Assembly

The process from sample to assembled mitochondrial genome involves a series of interconnected steps, visualized in the following workflow.

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Reagents and Kits for Parasite DNA Research

Item	Function/Application
TIANamp Genomic DNA Kit [3]	Silica-membrane-based purification of total genomic DNA from various sample types, including blood.
QIAGEN DNeasy Blood & Tissue Kit [29]	Standardized protocol for efficient purification of DNA from tough-to-lyse samples, including ticks.
Ammonium Hydroxide (0.7 M) [29]	Single-reagent, low-cost hydrolysis method for rapid DNA release, suitable for PCR-based assays.
Stainless Steel Beads (2.5 mm) [29]	Used with a homogenizer for the mechanical disruption of tough parasite exoskeletons and cells.
Absolute Ethanol [29]	A key preservative for field-collected samples; also used in wash steps of many silica-membrane kits.
Illumina Novoseq 6000 Platform [3]	High-throughput sequencing platform for generating the short-read data used in genome assembly.
Oxford Nanopore Technologies (ONT) MinION [10]	Portable sequencer capable of producing long reads, useful for assembling through repetitive regions.
M8 metabolite of Carvedilol-d5	M8 metabolite of Carvedilol-d5, MF:C12H19NO4, MW:246.31 g/mol
DMT-dC(bz) Phosphoramidite-13C9,15N3	DMT-dC(bz) Phosphoramidite-13C9,15N3, MF:C46H52N5O8P, MW:845.8 g/mol

Quality Assessment and Downstream Application

Evaluating the success of DNA extraction is crucial before proceeding to costly sequencing.

• Spectrophotometry (e.g., NanoDrop): Provides a rapid assessment of DNA concentration and purity via A260/280 ratio. However, this method can be inaccurate for crude lysates, with ratios often below the ideal of 1.8-2.0 [29].
• Fluorometry (e.g., Qubit): A more DNA-specific quantification method that is less affected by contaminants like RNA or protein [29].
• Gel Electrophoresis: Assesses DNA integrity by visualizing high-molecular-weight fragments, indicating minimal degradation [29].
• qPCR: The most functional assessment, confirming that the DNA is amplifiable and free of inhibitors. This can be coupled with pathogen-specific assays to confirm infection status in the sample [29].

For mitochondrial genome assembly, the extracted DNA is used to generate sequencing libraries. Assembly can be performed de novo using tools like IDBA [3], and the resulting genome is annotated using specialized servers like MITOS [3]. The linear monomer structure of apicomplexan mitochondrial genomes, typically ~6,000 bp and encoding three core protein-coding genes (cox1, cob, cox3), serves as a valuable molecular marker for constructing robust phylogenies and clarifying taxonomic relationships [3].

Within parasite taxonomy research, the assembly of mitochondrial (mt) genomes is a fundamental tool for resolving phylogenetic relationships and understanding evolutionary histories [30]. The selection of an appropriate sequencing platform is critical, as it directly impacts the accuracy and completeness of the genomic data upon which these taxonomic conclusions are drawn. The two predominant technologies, Illumina (short-read) and Oxford Nanopore Technologies (ONT; long-read), offer distinct advantages and limitations. This application note provides a structured comparison of these platforms, focusing on coverage and accuracy within the specific context of mt-genome assembly for parasitic organisms. The guidance herein is designed to assist researchers in making an informed selection based on their specific project goals, whether for high-accuracy variant detection or for resolving complex genomic structures.

Platform Comparison: Technical Specifications and Performance

The fundamental differences in chemistry and data output between Illumina and Nanopore sequencing directly influence their performance in genomic applications. The following section quantifies these differences and summarizes their implications for mt-genome projects.

Table 1: Sequencing Platform Performance Metrics

Metric	Illumina	Oxford Nanopore
Read Length	Short-reads (up to 2x300 bp for MiSeq/NovaSeq X) [31]	Long-reads (full-length 16S rRNA ~1,500 bp; ultra-long reads >100 kb) [31] [32]
Raw Read Accuracy	Very High (>99.9%, Q30) [33]	High (Q20+ chemistry: >99%, up to 99.9% with latest basecalling) [32] [34]
Primary Error Mode	Substitution errors	Higher single-read error rate (5–15% historically), indel errors more common [31] [35]
Variant Calling (SNV) Accuracy	Exceptionally high; 6x fewer SNV errors than Ultima UG 100 in one benchmark [36]	Comparable to short-reads for SNPs with latest chemistry; F1 score is key metric [32]
Variant Calling (Indel) Accuracy	Excellent; 22x fewer indel errors than UG 100 platform in one benchmark [36]	Lower in homopolymers; indel accuracy decreases in homopolymers >10 bp [36]
Coverage of Challenging Regions	May struggle with repetitive regions, homopolymers [32]	Excels in complex regions; reduces "dark" areas of genome by 81% [32]
Consensus / Assembly Accuracy	High per-base accuracy ideal for mapping assemblies	Very high consensus accuracy (e.g., Q50 for bacterial assembly); long reads resolve repeats [32]
Epigenetic Modification Detection	Requires bisulfite conversion	Direct, real-time detection of base modifications (e.g., 5mC accuracy 99.5%) [32]

Interpreting the Data for Mitochondrial Genome Assembly

• Accuracy Definitions: For mt-genome assembly, it is crucial to distinguish between raw read accuracy and consensus accuracy. While Illumina provides highly accurate single reads, ONT's longer reads can be assembled into a consensus sequence that achieves very high accuracy (e.g., Q50 or higher), which is often the primary goal in a taxonomic study [32] [35].
• Coverage Comprehensiveness: A significant advantage of long-read technology is its ability to span repetitive and GC-rich regions that are problematic for short reads. This is critical for obtaining a complete, rather than a fragmented, mt-genome assembly [32]. Illumina's coverage can drop in mid-to-high GC-rich regions, potentially missing biologically relevant genes [36].
• Application-Specific Performance: A study on Clostridioides difficile found that while Illumina was superior for high-resolution epidemiological surveillance due to lower error rates, Nanopore performed satisfactorily for virulence gene detection and sequence type assignment, offering a faster, cheaper alternative for less detailed analyses [35]. This trade-off is directly applicable to parasite taxonomy, where the research question (e.g., species identification vs. strain-level variation) should guide platform choice.

Experimental Protocols for Mitochondrial Genome Sequencing

The following protocols outline proven methods for generating mt-genome data using Illumina and Nanopore platforms, as applied in recent taxonomic research.

Protocol A: Illumina-Based Mt-Genome Assembly via NGS

This protocol is adapted from the methodology used to sequence the mitochondrial genome of Dugesia cantonensis [26].

1. DNA Extraction:

• Material: Use the E.Z.N.A. Mollusc DNA Isolation Kit or a similar kit appropriate for the parasite sample.
• Procedure: Extract total genomic DNA from ~10 mg of tissue or 10 whole animals. Follow the manufacturer's protocol, with an optional extended lysis step (30-60 minutes) to ensure high molecular weight DNA.

2. Library Preparation and Sequencing:

• Technology: Illumina NovaSeq 6000 platform.
• Process: The extracted DNA is fragmented, and a library is constructed without prior enrichment for mitochondrial DNA. This whole-genome sequencing approach allows for the simultaneous recovery of nuclear and mitochondrial data.
• Coverage: Sequence to a high depth of coverage (e.g., 30x for the nuclear genome, which will result in ultra-high coverage for the mt-genome).

3. Data Processing and Assembly:

• Quality Control: Assess raw read quality using FastQC. Trim adapters and low-quality bases using Cutadapt [26].
• De Novo Assembly: Assemble the mitochondrial genome from the trimmed reads using a dedicated organelle assembler like NOVOPlasty, which is designed for assembling circular genomes [26].
• Annotation: Annotate the assembled genome using the MITOS2 online server or a similar tool to identify protein-coding genes, rRNAs, and tRNAs [26] [30].

Protocol B: Nanopore-Based Mt-Genome Assembly

This protocol is derived from methods used for high-accuracy assembly and can be applied to parasite samples.

1. DNA Extraction for Long Reads:

• Material: Use kits designed for high-molecular-weight (HMW) DNA extraction, such as the DNeasy PowerSoil Pro Kit, which includes mechanical lysis via bead beating [35].
• Critical Note: The integrity of the DNA is paramount. Avoid excessive vortexing and use wide-bore tips to prevent shearing.

2. Library Preparation and Sequencing:

• Technology: Oxford Nanopore MinION or PromethION sequencers.
• Library Kit: Ligation Sequencing Kit (e.g., SQK-LSK114) is commonly used. For the highest accuracy assemblies requiring ultra-long reads, the Ultra-long Sequencing Kit (ULK) is available [32].
• Basecalling: Perform basecalling in real-time or post-run using the Dorado basecaller with a Super Accuracy (SUP) model, which is the most computationally intense but recommended for de novo assembly [32].

3. Data Processing and Assembly:

• Quality Filtering: Filter raw reads for quality and length.
• Assembly: Perform de novo assembly using long-read assemblers such as Flye or Verkko. For the highest quality "telomere-to-telomere" assemblies, a combination of ultra-long reads, Pore-C, and assembly polishing data can be used [32].
• Polishing: While long reads alone can produce high-quality assemblies, the consensus accuracy can be further improved by polishing the assembly with the same set of reads using tools like Medaka [32].

Diagram 1: Mt-genome Sequencing Workflow

The Scientist's Toolkit: Essential Research Reagents

Table 2: Key Reagents and Kits for Mitochondrial Genome Sequencing

Item	Function	Application Note
E.Z.N.A. Mollusc DNA Kit	Genomic DNA extraction from difficult samples.	Optimal for diverse parasite tissues; used successfully for planarian DNA prep [26].
DNeasy PowerSoil Pro Kit	HMW DNA extraction with mechanical lysis.	Bead beating step is effective for tough parasite cell walls; used for bacterial DNA [35].
Illumina Nextera XT Kit	Library preparation for Illumina sequencing.	Facilitates rapid library construction for fragmented DNA; standard for WGS [35].
ONT Ligation Seq Kit	Prepares DNA libraries for nanopore sequencing.	Standard kit for generating sequencing-ready libraries from HMW DNA [32].
Dorado Basecaller	Converts raw nanopore signals to nucleotide sequences.	Use SUP model for highest accuracy in de novo assembly projects [32].
NOVOPlasty	De Novo assembler for organelle genomes.	Specifically designed for assembling circular mt-genomes from NGS reads [26].
MITOS2	Automated annotation of metazoan mt-genomes.	Critical for identifying and annotating genes in the newly assembled sequence [26] [30].
Sumatriptan hydrochloride	Sumatriptan Hydrochloride	Sumatriptan hydrochloride, a selective 5-HT1B/1D receptor agonist. For research into migraine mechanisms. This product is for Research Use Only.
PROTAC BRD9 Degrader-3	PROTAC BRD9 Degrader-3, MF:C41H47ClN6O6, MW:755.3 g/mol	Chemical Reagent

The choice between Illumina and Nanopore technologies for mitochondrial genome assembly in parasite taxonomy is not a matter of which platform is universally superior, but which is most fit-for-purpose.

• Select Illumina sequencing when the research priority is maximum per-base accuracy for detecting single-nucleotide variants (SNVs) and small indels within well-characterized genomic regions. It is the established choice for high-resolution phylogenetic analysis based on point mutations and for projects requiring high throughput at a lower cost per sample [31] [35].
• Select Nanopore sequencing when the research involves structurally complex genomes, requires the resolution of repetitive regions, or aims to produce a complete, gap-free mitochondrial genome. Its long reads are invaluable for de novo assembly, and its ability to natively detect epigenetic modifications can provide an additional layer of taxonomic information [32] [30]. Furthermore, its portability and rapid turnaround time make it ideal for fieldwork or diagnostics.

For the most robust outcomes, a hybrid approach is increasingly favored. This strategy leverages the high accuracy of Illumina short reads to polish consensus sequences generated from Nanopore long reads, thereby combining the strengths of both technologies to produce a highly accurate and complete mitochondrial genome assembly [31].

Mitochondrial genome assembly represents a critical frontier in parasite taxonomy and drug development research, providing insights into evolutionary history, species identification, and potential therapeutic targets. Unlike nuclear genomes, mitochondrial genomes present unique assembly challenges due to their dynamic structure, including the presence of extensive repetitive regions and frequent homologous recombination events [37]. The selection of appropriate k-mer sizes—short DNA subsequences of length k used for genome reconstruction—plays a pivotal role in determining assembly success, particularly in navigating these complex repetitive landscapes [38]. For researchers investigating parasitic taxa, where mitochondrial genomes may exhibit unconventional structures and high mutation rates, optimized de novo assembly strategies are indispensable for generating accurate genomic resources that support taxonomic classification and reveal essential biological functions. This protocol details comprehensive strategies for addressing k-mer selection and repeat resolution, specifically contextualized within parasite mitochondrial genome research.

Comparative Performance of Assembly Tools and Strategies

Current methodologies for plant mitochondrial genome assembly can be broadly categorized into three algorithmic approaches: reference-based assembly, de novo assembly, and iterative mapping and extension (IME) [39]. Reference-based methods align sequencing reads to a known reference genome, which can be ineffective for non-model parasites with distant relatives. De novo assembly reconstructs genomes without prior knowledge, making it ideal for novel parasitic organisms, while IME methods iteratively refine assemblies through repeated mapping. Among 416 analyzed articles on plant mitochondrial genomes, 333 utilized de novo assembly, establishing it as the predominant strategy for organelle genome reconstruction [39].

Tool Performance Evaluation

A comprehensive evaluation of 11 frequently used assembly tools over the past five years, along with two newly developed tools (TIPPo and Oatk), revealed significant performance variations [39]. The assessment considered completeness, contiguity, and correctness of assembled mitochondrial genomes. SMARTdenovo, NextDenovo, and Oatk demonstrated superior performance in terms of contiguity and completeness, generating longer contiguous sequences (contigs) with fewer gaps. Meanwhile, GetOrganelle and Flye excelled in correctness, producing assemblies with fewer misassemblies and errors [39]. Tools specifically designed for mitochondrial assembly, such as PMAT and MitoHiFi, leverage long-read data to better resolve complex repetitive structures [39].

Table 1: Performance Evaluation of Mitochondrial Genome Assembly Tools

Assembly Tool	Optimal Sequencing Data	Strengths	Key Applications in Literature
SMARTdenovo	Long-read (PacBio, Nanopore)	Superior contiguity	General plant mitochondrial assembly [39]
NextDenovo	Long-read	High completeness	General plant mitochondrial assembly [39]
Oatk	Long-read	Excellent contiguity	General plant mitochondrial assembly [39]
GetOrganelle	Short-read (Illumina)	High correctness	General plant mitochondrial assembly [39]
Flye	Long-read	High correctness, handles repeats	Oak mitogenome assembly [19]
NOVOPlasty	Short-read	Seed-and-extend algorithm	Aria alnifolia mitogenome [37]
Norgal	Short-read WGS	k-mer frequency analysis, no reference needed	Panda, brown algae, butterfly mitogenomes [40]
MitoHiFi	Long-read	Automated annotation	Aria alnifolia mitogenome [37]
PMAT	Long-read	Optimized for plant mitogenomes	Strobilanthes sarcorrhiza, Pueraria montana [39] [41] [42]

k-mer Selection Strategies for Mitochondrial Reads Extraction

Fundamentals of k-mer Analysis

k-mers serve as fundamental units in de novo assembly, providing a computationally efficient method for processing large sequencing datasets [38]. The frequency distribution of k-mers in sequencing data directly correlates with genomic depth, enabling discrimination between nuclear and organellar DNA based on copy number variation [40]. Mitochondrial genomes typically exhibit 10-100 times higher copy numbers than nuclear genomes, resulting in proportionally higher k-mer frequencies that facilitate their computational separation from nuclear reads [40].

Implementation in Assembly Pipelines

The Norgal pipeline exemplifies a specialized approach for extracting mitochondrial DNA from whole-genome sequencing data using k-mer frequency analysis [40]. Its methodology capitalizes on the differential abundance of mitochondrial DNA without requiring reference sequences, making it particularly valuable for non-model parasites.

Table 2: k-mer Selection Parameters in Assembly Tools

Tool	Default k-mer Size	Selection Strategy	Impact on Assembly
Norgal	31	Frequency threshold based on nuclear depth	Higher thresholds reduce nuclear contamination but may exclude low-coverage mitochondrial regions [40]
MEGAHIT (within Norgal)	21, 49, 77, 105	Multiple k-mer assembly	Larger k-mers resolve repeats but require higher coverage [40]
idba_ud (within Norgal)	20, 40, 60, 80, 100	Iterative assembly with increasing k	Progressive increase improves repeat resolution [40]
NOVOPlasty	Adaptive	Seed-based extension	k-mer size adapts to local sequence complexity [43]

Protocol: k-mer-Based Mitochondrial Reads Extraction Using Norgal

• Input Preparation: Obtain whole-genome sequencing (WGS) data from parasitic organisms. Preprocess reads by removing adapters and trimming low-quality bases using AdapterRemoval with parameters --minlength 30 [40].

• Nuclear Depth Threshold Estimation:

  • Perform an initial de novo assembly of all WGS reads using MEGAHIT with default parameters and k-mer range: 21, 49, 77, and 105.
  • Identify the longest assembled contig (assumed nuclear origin).
  • Map reads back to the longest contig using bwa mem.
  • Calculate nuclear depth threshold (ND threshold) using the formula:
    where d_i represents read depth at position i, and n is the number of non-zero depths in the percentile range [40].

• k-mer Counting and Read Binning:

  • Count 31-mer occurrences in all reads using BBTools.
  • Extract reads containing at least one 31-mer with frequency exceeding the ND threshold.
  • This step effectively bins mitochondrial-rich reads based on copy number differentiation.

• Mitochondrial Assembly:

  • Perform de novo assembly of the binned reads using idba_ud with multiple k-mer sizes (20, 40, 60, 80, 100).
  • Extract the longest contig as the putative mitochondrial genome candidate.

• Validation:

  • Test contig circularity by identifying overlapping ends.
  • Annotate potential mitochondrial genes using BLAST against reference databases.
  • Validate assembly by mapping reads back to the extracted mitochondrial sequence [40].

Figure 1: k-mer-Based Mitochondrial Read Extraction Workflow. This diagram illustrates the sequential process for extracting mitochondrial DNA from whole-genome sequencing data using k-mer frequency analysis, as implemented in the Norgal pipeline [40].

Strategies for Resolving Complex Repeat Regions

Nature of Repetitive Sequences in Mitochondrial Genomes

Plant mitochondrial genomes contain abundant repetitive sequences that facilitate frequent homologous recombination, leading to alternative genomic conformations including circular, linear, and branched molecules [37]. These dynamic structures present significant assembly challenges, particularly when repetitive regions exceed sequencing read lengths, causing misassemblies and incorrect genome size estimation [39]. In parasitic taxa, these challenges may be exacerbated by limited genomic resources and unusual architectures.

Repeat Identification and Classification

Comprehensive repeat analysis involves identifying various repeat types:

• Simple Sequence Repeats (SSRs): Short, tandemly repeated sequences of 1-6 base pairs.
• Dispersed Repeats: Longer repetitive sequences (≥30 bp) distributed throughout the genome, including forward, reverse, palindrome, and complementary repeats [41].
• Tandem Repeats: Direct repeats occurring adjacent to each other.

Protocol: Comprehensive Repeat Analysis in Assembled Mitogenomes

• SSR Identification:

  • Utilize MISA with minimum repeat thresholds:
    • Mononucleotide: 10 repeats
    • Dinucleotide: 5 repeats
    • Trinucleotide: 4 repeats
    • Tetranucleotide: 3 repeats
    • Pentanucleotide: 3 repeats
    • Hexanucleotide: 3 repeats [41]

• Dispersed Repeat Detection:

  • Employ REPuter with parameters:
    • Minimum repeat size: 30 bp
    • Hamming distance: 3
    • Maximum computed repeats: 5,000
  • Classify repeats into four types: forward (F), reverse (R), palindrome (P), and complementary (C) [41]

• Repeat-Mediated Recombination Analysis:

  • Visualize assembly graphs using Bandage to identify double-bifurcating structures (DBSs) indicating alternative genomic arrangements [37].
  • Map long reads across repetitive regions to validate correct assembly paths.
  • For each DBS, extract sequences and map HiFi reads completely encompassing these regions to confirm proper connections [37].

Advanced Strategies for Repeat Resolution

Long-read sequencing technologies (PacBio HiFi, Nanopore) substantially improve repeat resolution through several mechanisms:

• Read Length Advantage: HiFi reads often exceed 15-20 kb, spanning most repetitive elements to provide unique flanking sequences for unambiguous placement [19] [44].

• Graph-Based Assembly:

  • Implement Flye assembler with default parameters to generate graphical fragment assemblies (GFA).
  • Identify mitochondrial contigs based on coverage depth and sequence similarity to known mitochondrial genes.
  • Manually resolve complex regions by examining supporting read evidence [44].

• Multi-Platform Integration:

  • Combine Illumina short reads with PacBio HiFi data for error correction.
  • Use short reads for high-base accuracy and long reads for scaffold continuity.
  • Employ hybrid assembly pipelines like MitoHiFi for integrated analysis [37].

Figure 2: Comprehensive Repeat Resolution Workflow. This diagram outlines the multi-faceted approach for identifying and resolving complex repetitive regions in mitochondrial genomes, incorporating both computational tools and long-read sequencing validation [41] [37].

Case Studies in Parasite-Related Taxa

Strobilanthes sarcorrhiza Mitogenome Assembly

The mitochondrial genome of the medicinal plant Strobilanthes sarcorrhiza was assembled using PMAT v1.5.3 with PacBio HiFi data, employing the "autoMito" model parameterized with "-st hifi -g 820m -m" [41]. Despite a relatively large genome size (617,134 bp) with linear structure, the assembly successfully resolved 1,482 pairs of dispersed repeats accounting for 17.58% of the entire mitogenome [41]. This case exemplifies the challenges of assembling large mitochondrial genomes with abundant repeats and demonstrates the efficacy of long-read strategies in medicinal plants with potential parasitic relatives.

Camellia oleifera and C. lanceoleosa Multipartite Structures

A comparative analysis of two closely related Camellia species revealed extensive genome rearrangements and multipartite structures [44]. Researchers used Flye assembler with HiFi long reads, followed by meticulous source identification of contigs based on coverage depth and sequence similarity. The resulting assemblies (1,039,838 bp for C. oleifera and 934,155 bp for C. lanceoleosa) confirmed multiple-branched configurations rather than conventional circular molecules [44]. This study highlights the importance of graph-based assembly approaches and coverage-based binning for resolving complex mitochondrial architectures.

Aria alnifolia Recombination Analysis

In Aria alnifolia, researchers identified 12 double-bifurcating structures within the mitochondrial assembly graph, indicating potential sites for repeat-mediated homologous recombination [37]. By mapping both Illumina short reads and PacBio HiFi reads to these regions, they confirmed the presence of alternative genomic conformations and reconstructed the master circle configuration (455,361 bp) [37]. This approach demonstrates how integration of multiple sequencing technologies enables resolution of dynamic mitochondrial structures.

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Essential Research Reagents and Computational Tools for Mitochondrial Genome Assembly

Category	Item/Reagent	Specification/Function	Application Example
Sequencing Kits	PacBio SMRTbell prep kit 3.0	Library preparation for HiFi sequencing	Oak mitogenome assembly [19]
	DNeasy Plant Mini Kit	DNA extraction from plant tissues	Camellia mitogenome assembly [44]
DNA Quality Assessment	Agilent 4200 Bioanalyzer	DNA integrity assessment	Camellia mitogenome assembly [44]
	Qubit Fluorometer	DNA concentration measurement	Strobilanthes sarcorrhiza [41]
Computational Tools	BBTools Suite	k-mer counting and read binning	Norgal pipeline [40]
	Bandage v0.8.1	Visualization of assembly graphs	Aria alnifolia, Camellia [41] [44]
	MISA	Microsatellite identification	Strobilanthes sarcorrhiza [41]
	REPuter	Dispersed repeat detection	Strobilanthes sarcorrhiza [41]
Annotation Resources	GeSeq	Organelle genome annotation	Aria alnifolia, Strobilanthes sarcorrhiza [41] [37]
	tRNAscan-SE	tRNA gene identification	Strobilanthes sarcorrhiza [41]
Reference Databases	RefSeq Mitochondrial	Reference sequences for annotation	Norgal validation [40]
Ac-Gly-Ala-Val-Ile-Leu-Arg-Arg-NH2	Ac-Gly-Ala-Val-Ile-Leu-Arg-Arg-NH2, MF:C36H68N14O8, MW:825.0 g/mol	Chemical Reagent	Bench Chemicals
S-(N-Methylcarbamoyl)glutathione-d3	S-(N-Methylcarbamoyl)glutathione-d3, MF:C12H20N4O7S, MW:367.40 g/mol	Chemical Reagent	Bench Chemicals

Effective de novo assembly of mitochondrial genomes for parasite taxonomy research requires careful consideration of k-mer selection strategies and implementation of comprehensive repeat resolution protocols. The integration of long-read sequencing technologies with sophisticated computational approaches enables researchers to overcome historical challenges associated with complex mitochondrial architectures. By applying the detailed protocols and comparative analyses presented herein, research scientists and drug development professionals can generate high-quality mitochondrial genomic resources that support accurate taxonomic classification and provide insights into fundamental biological processes of parasitic organisms. As sequencing technologies continue to evolve, further refinements to these strategies will undoubtedly enhance our ability to decipher even the most complex mitochondrial genomes across diverse taxonomic groups.

Mitochondrial genome assembly and annotation serve as a cornerstone for modern parasite taxonomy and phylogenetics, providing critical insights into evolutionary relationships and molecular evolution [3]. For apicomplexan parasites, such as Theileria and Babesia species, the mitochondrial genome represents an essential molecular marker due to its higher evolutionary rate compared to nuclear DNA and greater reliability for discriminating between closely related species [3]. The annotation of these genomes, however, presents significant challenges, including the accurate identification of protein-coding genes (PCGs), structural RNAs, and the determination of gene boundaries amidst unusual genetic codes and compositional biases [45]. This application note provides detailed protocols for utilizing two complementary tools—MITOS and BLAST—to address these challenges and produce high-quality mitochondrial genome annotations specifically tailored for parasite research.

The selection of appropriate annotation tools is critical for generating accurate and biologically meaningful mitochondrial genome annotations. MITOS (MITOchondrial genome annotation Server) is an automated pipeline designed for de novo annotation of metazoan mitochondrial genomes, leveraging curated covariance models for structural RNAs and sophisticated similarity searches for protein-coding genes [46]. In parallel, BLAST (Basic Local Alignment Search Tool) finds regions of local similarity between sequences, enabling researchers to infer functional and evolutionary relationships and identify members of gene families through comparison with sequence databases [47] [48].

Table 1: Comparison of Mitochondrial Genome Annotation Tools

Tool	Primary Approach	Strengths	Limitations	Ideal Use Cases
MITOS	De novo annotation using profile HMMs and covariance models [46]	Consistent annotation strategy; automated tRNA and rRNA identification; no close relative requirement [46]	Optimized for metazoans; limited capability for intron-rich genes [45]	Initial structural annotation of novel parasite mitogenomes
BLAST	Similarity-based search using local alignment algorithms [48]	Extremely versatile; can identify novel homologs; confirms gene identity and function [48]	Requires existing reference sequences; prone to propagating historical errors [45]	Validation of MITOS predictions; identification of horizontal gene transfer events

Two fundamental annotation strategies exist: next-neighbour-guided annotation, which transfers annotations from closely related species using BLAST-like algorithms, and ab-initio inference, which uses probabilistic methods like profile Hidden Markov Models (HMMs) to identify evolutionarily conserved signatures without requiring close relatives [45]. MITOS primarily employs an ab-initio approach, making it particularly valuable for parasite taxa with no closely related annotated mitogenomes available.

Experimental Protocols

Protocol 1: De Novo Annotation Using the MITOS Web Server

This protocol details the steps for annotating a mitochondrial genome using the MITOS web interface, which is optimized for metazoan sequences and requires no prior programming knowledge.

Input Preparation and Submission

• Prepare Sequence File: Ensure your mitochondrial genome assembly is in FASTA format. The sequence header should contain relevant identifiers (e.g., species name, isolate code).
• Access MITOS: Navigate to the MITOS web server at http://mitos.bioinf.uni-leipzig.de/.
• Submit Job:
  • Upload your FASTA file or paste the sequence directly into the input field.
  • Select the appropriate Genetic Code. For most apicomplexan parasites, the "Invertebrate Mitochondrial" or "Protozoan Mitochondrial" code is applicable.
  • Retain the default parameters for the initial analysis. The "Reference" option can be set to "RefSeq 63" for a comprehensive gene set.
  • Initiate the annotation by clicking the "Submit" button. Processing time varies from minutes for small genomes to a few hours for larger or more complex assemblies.

Interpretation of Results

• Summary Page: Upon completion, MITOS presents a summary table listing all predicted genes, their positions, strands, and other features.
• Gene Annotations:
  • Protein-Coding Genes (PCGs): Verify the predicted PCGs (e.g., cox1, cob, cox3 for apicomplexans [3]) by examining their length and comparing them to known gene lengths in related species to avoid truncated or extended predictions.
  • Structural RNA Genes: MITOS uses covariance models to annotate rRNAs and tRNAs with high accuracy [46]. Note that some apicomplexan mitogenomes, like those of Theileria, lack tRNA genes and contain fragmented rRNA genes [3].
  • Genetic Code Verification: Confirm that start and stop codons conform to the expected patterns for your parasite clade. For example, start codons in Theileria and Babesia spp. predominantly include ATN, GTN, and TTN, while end codons are mainly TAA, TAG, and TGA [3].
• Output Files: Download the complete result set, typically including a annotated GenBank file, a GFF3 file for genomic viewers, and detailed reports for tRNAs and rRNAs.

Protocol 2: Gene Validation and Analysis Using BLAST

This protocol describes how to use BLAST to validate gene predictions from MITOS and investigate specific gene characteristics, such as potential horizontal gene transfer (HGT) events, which are common in parasitic plants [49].

Single Gene Validation via blastp/tblastn

• Extract Gene Sequence: From the MITOS output, extract the nucleotide or amino acid sequence of the gene requiring validation (e.g., cox1).
• Access BLAST: Use the NCBI BLAST web service (https://blast.ncbi.nlm.nih.gov/Blast.cgi) or an integrated platform like Geneious.
• Configure Search:
  • For Amino Acid Sequences (blastp): Use the "blastp" algorithm against the "Non-redundant protein sequences (nr)" database to find homologous proteins.
  • For Nucleotide Sequences (tblastn): If the gene is poorly conserved at the nucleotide level, use "tblastn" to query your protein sequence against a translated nucleotide database. This is often more sensitive for divergent sequences.
  • Set Search Parameters: Adjust the "Expect threshold" (E-value) to 0.05 or lower to restrict results to statistically significant matches [48]. Use an "Entrez Query" to limit searches to specific taxonomic groups (e.g., Apicomplexa).
• Analyze Results:
  • Assess the top hits based on E-value, Percent Identity, and Query Coverage [48]. A significant E-value (e.g., < 1e-10) and high query coverage strongly support the annotation.
  • Inspect the alignment viewer to check for correct gene boundaries, the absence of frameshifts, and conserved functional domains.

Advanced Analysis: Screening for HGT

• Perform Broad BLAST Search: Run a blastp or tblastn search without taxonomic restrictions using a mitochondrial gene of interest.
• Identify Anomalous Hits: Construct a phylogenetic tree or carefully examine the taxonomy of the top hits. Genes acquired via HGT will show unexpectedly high similarity to genes from distantly related taxa or potential host species [49].
• Confirm with Additional Genes: HGT events are more robustly supported when evidence is found for multiple genes from the same donor lineage.

Table 2: BLAST Algorithms and Their Applications in Mitogenome Annotation

Algorithm	Query Type	Database Type	Primary Application in Annotation
blastn	Nucleotide	Nucleotide	Initial identification of conserved genes; validation of rRNA fragments
Megablast	Nucleotide	Nucleotide	Fast, high-similarity searches for well-conserved genes within a clade
blastp	Protein	Protein	Standard for validating predicted PCGs; functional inference
tblastn	Protein	Nucleotide (translated)	Sensitive identification of divergent or novel PCGs; finding potential HGT events
blastx	Nucleotide (translated)	Protein	Identifying protein-coding regions in unannotated genomic sequence

Table 3: Key Research Reagents and Computational Resources for Mitochondrial Genome Annotation

Item/Resource	Function/Description	Example/Application in Parasite Taxonomy
MITOS Web Server	De novo annotation pipeline for metazoan mitochondrial genomes [46].	Primary structural annotation of PCGs, rRNAs, and tRNAs in a novel parasite mitogenome.
NCBI BLAST Suite	Toolsuite for finding sequence similarities and inferring homology [47] [48].	Validating MITOS-predicted cox1 gene; screening for horizontally acquired genes.
Genetic Code Table (Protozoan/Invertebrate)	Specifies codon assignments for translation.	Correctly interpreting start/stop codons and coding sequences in apicomplexan mitochondrial genes.
Reference Mitogenomes (RefSeq)	Curated, high-quality genomic sequences from databases like NCBI RefSeq.	Used as references for BLAST searches and for comparative analysis of gene order and content.
Sequence Assembly Software (e.g., IDBA)	Assembles sequencing reads into contiguous sequences (contigs) [3].	Generating a complete mitochondrial genome assembly from Illumina or Nanopore reads.
Multiple Sequence Alignment Tool (e.g., ClustalW)	Aligns three or more biological sequences to identify regions of similarity [3].	Preparing aligned sequences of cox1 and cob genes for phylogenetic analysis.

Application in Parasite Taxonomy: A Case Study onTheileria velifera

The utility of this integrated MITOS-BLAST pipeline is exemplified by its application in the characterization of the mitochondrial genome of Theileria velifera, an apicomplexan parasite [3]. The study followed a streamlined workflow: after sequencing and assembly, the mitochondrial genome was annotated using the MITOS web server [3]. The MITOS output provided the initial gene models, which were subsequently verified and refined by BLAST searches against the GenBank database to confirm the identity of homologous proteins [3].

This combined approach successfully identified the three characteristic apicomplexan PCGs—cox1, cox3, and cob—and the five fragmented large subunit (LSU) rRNA genes, while confirming the absence of tRNA genes [3]. The resulting annotation, 6,125 bp in length, was deposited in GenBank (ON684327) and served as the foundation for downstream comparative and phylogenetic analyses. These analyses involved aligning the cox1 and cob gene sequences with those from other apicomplexans using ClustalW in MEGA software and constructing a maximum likelihood phylogenetic tree, which clearly resolved the evolutionary position of T. velifera relative to other Theileria species [3].

The synergistic use of the MITOS pipeline for de novo gene prediction and BLAST for homology-based validation provides a robust, efficient, and accessible framework for annotating mitochondrial genomes. This integrated approach is particularly powerful in the field of parasite taxonomy, where it enables the accurate resolution of evolutionary relationships, even in non-model organisms with limited prior genomic information. As demonstrated in the Theileria velifera case study, the annotations generated through these protocols form a critical foundation for subsequent comparative genomics, population genetics, and phylogenetic studies, ultimately advancing our understanding of parasite evolution and aiding in the development of targeted control strategies.

Integrative taxonomy combines multiple lines of evidence to delineate species boundaries and establish robust taxonomic classifications. For parasite research, this approach is particularly valuable, as many parasites exhibit conservative morphology with cryptic genetic diversity. The mitochondrial genome serves as a cornerstone in these investigations due to its maternal inheritance, lack of recombination, and predictable mutation rate, which provide reliable phylogenetic signals. This protocol details the methodology for combining mitochondrial genomic data with morphological and histological analyses, creating a comprehensive framework for parasite taxonomy and drug target identification.

The application of integrative taxonomy has become increasingly important in haemosporidian parasite studies (phylum Apicomplexa, order Haemosporida), which include the agents of malaria. These parasites infect various reptiles, mammals, and birds worldwide, and their accurate identification is crucial for understanding disease dynamics and developing control strategies. Mitochondrial genes, particularly the cytochrome b gene (cytb), have become a de facto DNA barcode for these parasites, but full mitochondrial genome sequencing provides substantially more phylogenetic information for resolving complex taxonomic relationships [7].

Experimental Protocols and Workflows

Mitochondrial Genome Assembly Protocol

Sample Preparation and DNA Extraction

• Source Material: Begin with 10-20 mg of parasite-rich tissue or purified parasites from host blood samples.
• DNA Extraction: Use the E.Z.N.A. Mollusc DNA Isolation Kit (or similar specialized kit for parasite DNA) according to manufacturer protocols [26]. Include negative controls to monitor contamination.
• DNA Quantification: Assess DNA quality and quantity using spectrophotometry (NanoDrop) and fluorometry (Qubit). High-molecular-weight DNA with A260/A280 ratios of 1.8-2.0 is ideal for long-read sequencing.

Mitochondrial Genome Amplification and Sequencing Two complementary approaches can be employed based on available technology and sample quality:

• Long-Range PCR and PacBio HiFi Sequencing (Ideal for detecting mixed infections) [7]:

  • Primer Design: Use slightly modified AE170 (5′-GAT TCT CTC CAC ACT TCA ATT CGT ACT TC-3′) and AE171 (5′-GAA AAT WAT AGA CCG AAC CTT GGA CTC-3′) primers with 5′ barcode sequences for multiplexing.
  • PCR Conditions: Denaturation at 94°C for 1 min, followed by 30 cycles of 98°C for 10 s, annealing at 55°C for 15 s, and extension at 68°C for 5 min.
  • Library Preparation: Multiplex up to 192 specimens in an SMRTbell library preparation for PacBio HiFi sequencing, which generates long reads (10-25 kb) with approximately 99.5% accuracy.

• Shotgun Sequencing and Assembly (Alternative approach) [26]:

  • Library Preparation: Use Illumina Novaseq 6000 for next-generation sequencing with 150-250 bp insert sizes.
  • Genome Assembly: Perform de novo assembly using NOVOPlasty software version 4.3.1 or similar specialized mitogenome assemblers.
  • Gap Filling: Design primers at both ends of linear sequences to fill gaps and validate assembly through PCR and Sanger sequencing.

Mitochondrial Genome Annotation

• Automated Annotation: Use Mitos2 for initial gene annotation [26].
• tRNA Validation: Predict secondary structures of mitochondrial transfer RNA genes using tRNAscan-SE.
• Manual Curation: Verify gene boundaries, identify non-coding regions, and confirm start/stop codons through comparison with closely related species.
• Visualization: Generate circular mitochondrial maps using Organellar Genome DRAW (OGDRAW) version 1.3.1 [26].

Morphological and Histological Analysis Protocol

Specimen Collection and Preservation

• Collect parasites from natural infections or experimental hosts using appropriate methods (blood draws, tissue biopsies).
• Fix specimens simultaneously for molecular (96% ethanol) and morphological (4% formaldehyde, 2.5% glutaraldehyde) analyses.
• Document collection location, date, host species, and tissue type for each specimen.

Morphological Characterization

• Light Microscopy: Examine fresh or fixed specimens under compound and stereomicroscopes.
• Measurement Protocol: Capture digital images and use ImageJ software version 1.53 or similar to measure key morphological characters [26].
• Documentation: Record quantitative morphological data including length, width, organ proportions, and distinctive features.

Histological Processing and Analysis

• Tissue Processing: Dehydrate specimens through ethanol series, clear with xylene, and embed in paraffin wax.
• Sectioning: Cut 5-7 μm thick sections using a rotary microtome.
• Staining: Employ standard staining protocols (Haematoxylin and Eosin, Giemsa) and specialized stains for specific structures.
• Microscopy: Examine sections under light microscope, capturing digital images of key anatomical features.
• Reproductive Analysis: Carefully examine histological sections for presence/absence of reproductive organs to confirm asexual lineages [26].

Integrative Data Analysis Workflow

Figure 1. Integrative taxonomy workflow for parasite classification combining mitochondrial genomics with morphological and histological data.

Data Presentation and Analysis

Quantitative Data Analysis in Integrative Taxonomy

Table 1. Key Quantitative Measurements for Integrative Taxonomy of Parasites

Data Category	Specific Metrics	Measurement Method	Application in Taxonomy
Mitogenomic Features	Genome length, GC content, gene order, AT/GC skew	Sequencing assembly, Mitos2 annotation	Phylogenetic placement, evolutionary relationships
Gene Sequences	COI, cytb, 18S rDNA, 28S rDNA sequences	PCR amplification, Sanger/NGS sequencing	DNA barcoding, species delimitation
Morphometrics	Body length, width, organ proportions, cell counts	Digital imaging, ImageJ analysis	Species characterization, diagnostic features
Karyological Data	Chromosome number, ploidy, centromere position	Karyotyping, flow cytometry	Evolutionary history, kinship relationships
Histological Features	Tissue organization, reproductive structures, gland patterns	Histological staining, microscopic examination	Confirmation of reproductive mode, structural analysis

Table 2. Mitochondrial Genome Characteristics of Representative Parasite Taxa

Parasite Group	Genome Size (bp)	GC Content (%)	Gene Content	Unique Features	References
Haemosporida	~6,000	20-30%	3 COX genes, cytb, rRNAs	Linear conformation, multicopy	[7]
Plasmodium spp.	5,966-6,009	24.8-30.4%	cox1, cox3, cytb, rRNAs	Conserved gene order	[7]
Dugesia cantonensis	18,125	-	36 genes, lacks atp8	Circular conformation	[26]
Leucocytozoon spp.	~6,000	~25%	cox1, cox3, cytb, rRNAs	Distinct from Plasmodium	[7]

Phylogenetic Analysis and Species Delimitation

Molecular Phylogenetics

• Sequence Alignment: Use PhyloSuite v.1.2.3 or similar for sequence extraction and concatenation [26].
• Tree Construction: Perform maximum likelihood analysis with 1000 bootstrap replicates under the GTR + I + G model using MEGA version 6 or similar software [26].
• Multi-locus Analysis: Combine mitochondrial markers (COI, cytb, complete mitogenome) with nuclear markers (18S rDNA, 28S rDNA) for robust phylogenies.

Species Delimitation Methods

• Genetic Distance: Calculate pairwise distances for COI and other markers to establish divergence thresholds.
• Phylogenetic Species Concept: Identify monophyletic clades with strong bootstrap support.
• Integrative Decision Framework: Combine genetic data with morphological gaps and biological characteristics.

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3. Research Reagent Solutions for Integrative Taxonomy of Parasites

Reagent/Material	Specification	Application	Key Considerations
DNA Extraction Kit	E.Z.N.A. Mollusc DNA Kit, DNeasy Blood & Tissue	High-quality DNA from various sample types	Optimize for parasite vs. host DNA recovery
PCR Primers	AE170/AE171, COI primers, 18S/28S rDNA primers	Mitochondrial genome amplification, barcoding	Test specificity for target parasite group
Sequencing Technology	PacBio HiFi, Illumina Novaseq, Sanger sequencing	Genome assembly, variant detection	Choose based on required resolution and budget
Histological Stains	Haematoxylin & Eosin, Giemsa, specialized stains	Tissue morphology, cellular structure	Standardize staining protocols for consistency
Image Analysis Software	ImageJ, commercial morphometrics packages	Quantitative morphological measurements	Calibrate measurements across users
Phylogenetic Software	MEGA, PhyloSuite, MrBayes, RAxML	Evolutionary analysis, tree construction	Use appropriate models of sequence evolution
3-Hydroxy-3-methylvaleric acid-d5	3-Hydroxy-3-methylvaleric acid-d5, MF:C6H12O3, MW:137.19 g/mol	Chemical Reagent	Bench Chemicals

Advanced Applications in Parasite Research

Mitochondrial Genome Analysis for Mixed Infections

The PacBio HiFi protocol enables detection and characterization of mixed infections and co-infections, which are common in wildlife parasites [7]. This approach uses a machine-learning pipeline (HmtG-PacBio Pipeline) that integrates multiple sequence alignments with modified variational autoencoders and clustering methods to identify mitochondrial haplotypes and species in a sample.

Key Parameters for Mixed Infection Analysis:

• Minimum Coverage: 30X per haplotype based on detected error rates [7]
• Error Rate: Averages 0.2% per read (range 0.03%-0.46%) [7]
• Detection Sensitivity: Capable of identifying parasite lineages present in very low parasitemia

Data Integration and Visualization Framework

Figure 2. Data integration framework for taxonomic decision-making in parasite systematics.

Integrative taxonomy combining mitochondrial genomics with morphology and histology provides a robust framework for parasite classification that resolves limitations of single-method approaches. The protocols outlined here enable researchers to generate comprehensive datasets that capture both genetic and phenotypic dimensions of parasite diversity. As sequencing technologies continue to advance, particularly long-read methods like PacBio HiFi, the capacity to resolve complex taxonomic questions in parasitology will further improve.

This integrated approach has significant implications for understanding parasite evolution, host-parasite interactions, and ultimately for developing targeted interventions against parasitic diseases. The methodological framework presented can be adapted to various parasite taxa, providing a standardized approach for taxonomic studies that facilitates comparison across groups and geographic regions.

Overcoming Assembly Challenges: Repeats, NUMTs, and Genome Instability

Addressing Repeat-Induced Assembly Collapse in Complex Mitochondrial Genomes

The assembly of complex mitochondrial genomes, particularly for parasite taxonomy research, is significantly hampered by the presence of abundant repetitive sequences. These repeats induce assembly collapse through misassembly during the merging of sequencing reads, leading to incorrect genome structures and size estimations [39]. Plant mitochondrial genomes exhibit remarkable structural diversity, including circular, linear, and branched configurations, with sizes ranging from 66 kb in Viscum scurruloideum to 18.99 Mb in Cathaya argyrophylla [39] [50]. This extensive size variation is largely driven by repetitive elements that facilitate frequent recombination events, creating substantial challenges for assembly algorithms not specifically designed to handle such complexity [50].

In parasite genomics, these challenges are compounded by the need for accurate taxonomic classification, where misassembled mitochondrial genomes can lead to incorrect phylogenetic placements and misunderstood evolutionary relationships. The presence of nuclear mitochondrial DNA (NUMT) and mitochondrial plastid DNA (MTPT) contamination further complicates the assembly process, as these sequences can be mistakenly incorporated into mitochondrial assemblies [39]. Overcoming repeat-induced assembly collapse is therefore critical for generating reliable mitochondrial references that can advance parasite taxonomy and drug development research.

Computational Tools and Performance Evaluation

Assembly Algorithm Classification and Selection

Mitochondrial genome assembly approaches generally fall into three algorithmic categories, each with distinct strengths for handling repetitive elements. Reference-based assembly utilizes closely related mitochondrial sequences as templates but is limited by the low interspecific sequence conservation in mitochondrial genomes [39]. De novo assembly has been the predominant approach, used in 333 of 387 analyzed studies, and reconstructs genomes without reference bias, making it suitable for novel parasite lineages [39]. Iterative mapping and extension methods, employed in 48 studies, build consensus sequences through repeated read mapping steps [39].

The performance of assembly tools varies significantly in handling repetitive regions. Tools specifically designed for mitochondrial genomes, such as GetOrganelle and Flye, excel in assembly correctness, while SMARTdenovo, NextDenovo, and Oatk demonstrate superior contiguity and completeness [39]. The recently developed PMAT (Plastid and Mitochondrial Assembly Tool) utilizes copy number differences among organellar genomes without pre-assembly filtering of mitochondrial reads, reducing susceptibility to NUMT and MTPT interference [50]. For parasite researchers, selection criteria should prioritize tools with demonstrated efficacy on genomes with high repeat content and those capable of resolving complex structural variants prevalent in mitochondrial genomes.

Table 1: Performance Evaluation of Mitochondrial Genome Assembly Tools

Tool	Algorithm Type	Strengths	Limitations	Optimal Use Cases
GetOrganelle	De novo	High correctness	Requires k-mer based read enrichment	Well-characterized parasites with reference data
Flye	De novo	Excellent correctness, handles long repeats	Computationally intensive	Novel parasite genomes with complex repeats
SMARTdenovo	De novo	Superior contiguity	May require additional polishing	Taxonomic studies requiring complete gene sets
NextDenovo	De novo	High completeness, handles long reads	Limited for short-read data	Large mitochondrial genomes with extensive repeats
Oatk	De novo	Excellent contiguity	Newer tool with less validation	Challenging repeat structures in uncharacterized parasites
PMAT	De novo	Direct assembly without read filtering; reduced NUMT/MTPT interference	Designed specifically for long-read data	Parasite taxa with potential plastid integrations

Advanced Tools for Complex Repeats

Specialized tools have emerged to address specific repeat-induced challenges in mitochondrial genome assembly. TIPPo and Oatk incorporate novel algorithms for resolving intermediate-sized repeats (50 bp-1 kb) that frequently cause assembly collapse in conventional pipelines [39]. For parasite taxa with potential plastid integrations, PMAT's approach of leveraging copy number differences without pre-filtering mitochondrial reads provides significant advantages in avoiding misclassification and structural loss [50] [51].

In application to parasite taxonomy, these tools enable the resolution of cryptic species complexes through complete mitochondrial assembly. For example, in avian haemosporidian parasites, Nanopore sequencing coupled with advanced assembly tools successfully resolved co-infections of Haemoproteus and Plasmodium lineages that would remain undetected with Sanger sequencing [52]. The implementation of these tools in taxonomic workflows significantly improves the detection of structural variants that define parasite lineages.

Experimental Design and Workflow Visualization

Integrated Assembly Strategy

A robust mitochondrial assembly protocol incorporates complementary technologies and algorithms to address repetitive elements. The workflow begins with DNA extraction using high-quality kits such as the Hi-DNAsecure Plant Kit or E.Z.N.A. Mollusc DNA Isolation Kit, optimized for mitochondrial preservation [26] [50]. Sequencing technology selection is critical, with third-generation long-read platforms (PacBio HiFi, Nanopore) providing the read lengths necessary to span repetitive regions that collapse short-read assemblies [39] [50].

A hybrid assembly approach implementing multiple algorithms on the same dataset significantly improves assembly completeness and correctness. For the oak mitochondrial genome project, researchers applied four complementary assembly strategies to reconstruct six distinct structural variants ranging from 339 to 622 kb [53]. This multi-algorithm approach enables the identification of assembly artifacts versus genuine biological structures, particularly in repetitive regions susceptible to misassembly.

Table 2: Research Reagent Solutions for Mitochondrial Genome Assembly

Reagent/Category	Specific Examples	Function in Workflow	Parasitology Applications
DNA Extraction Kits	Hi-DNAsecure Plant Kit, E.Z.N.A. Mollusc DNA Isolation Kit, Plant DNAzol Reagent	High-quality, high-molecular-weight DNA preservation	Parasite isolation from host tissues, mixed DNA populations
Long-read Sequencing	PacBio Revio (HiFi), Nanopore PromethION	Span repetitive elements, resolve structural variants	Cryptic parasite detection, co-infection resolution
Assembly Software	PMAT, Oatk, TIPPo, GetOrganelle	Specialized mitochondrial assembly, repeat resolution	Taxon-specific repeat profiles, NUMT exclusion
Validation Tools	Bandage, hifisr variant frequency estimation	Assembly graph visualization, structural validation	Verification of parasite-specific structural variants
Bait Enrichment	Custom mitochondrial bait panels	Mitochondrial read enrichment from total DNA	Host DNA depletion in host-parasite systems

Visualizing the Assembly Workflow

The following diagram illustrates a comprehensive workflow for addressing repeat-induced assembly collapse in complex mitochondrial genomes, integrating both experimental and computational components:

Diagram 1: Comprehensive Mitochondrial Assembly Workflow. This workflow integrates experimental and computational steps with specialized repeat handling to overcome assembly collapse. The color-coded nodes represent process categories: yellow for wet lab, blue for computation, green for analysis, and red for validation.

Detailed Methodological Protocols

Protocol 1: Hybrid Sequencing and Assembly for Repetitive Regions

Principle: This protocol combines long-read and short-read technologies to leverage their complementary strengths for resolving repetitive regions while maintaining base-level accuracy [50] [54] [18].

Step-by-Step Procedure:

• DNA Extraction and Quality Control
  • Use the Hi-DNAsecure Plant Kit or similar optimized for mitochondrial DNA preservation
  • Assess DNA quality via agarose gel electrophoresis (1.0%) and Nanodrop spectrophotometry
  • Require A260/A280 ratio of 1.8-2.0 and A260/A230 ratio >2.0
  • Verify high molecular weight DNA (>20 kb) using pulsed-field gel electrophoresis

• Library Preparation and Sequencing

  • Prepare 15-kb SMRTbell libraries using SMRTbell Express Template Prep Kit 2.0 (PacBio)
  • Construct Illumina libraries with 350-500 bp insert sizes using TruSeq DNA PCR-Free kit
  • Sequence on PacBio Revio platform for HiFi reads (minimum 7 Gb data, >7 kb average length)
  • Sequence on Illumina NovaSeq 6000 for 150 bp paired-end reads (minimum 30x coverage)

• Data Preprocessing

  • Process PacBio raw data with CCS algorithm (minPasses=3, minPredictedAccuracy=0.99)
  • Filter Illumina reads using Fastp v0.23.4 (remove adapters, quality cutoff Q5)
  • Filter Nanopore data using FilLong v0.2.1 (--minlength 1000, --minmean_q 7)

• Multi-Tool Assembly Implementation

  • Execute PMAT v1.5.3 with parameters: -g 3.8G, -t hifi, -m, -T 50
  • Run Flye v2.9 with parameters: --nano-hq for PacBio data or --nano-raw for Nanopore
  • Perform iterative assembly with NOVOPlasty v4.3.1 using mitochondrial core genes as seed
  • Visualize and compare assembly graphs using Bandage v0.8.1

• Repeat-Specific Resolution

  • Identify repetitive elements using MISA v1.0 and Tandem Repeats Finder
  • Resolve repetitive regions through manual inspection of read alignments
  • Verify repeat boundaries through PCR amplification and Sanger sequencing
  • Validate assembly continuity across repeats using read depth analysis

Validation and Quality Assessment:

• Estimate variant frequencies using hifisr pipeline to identify potential misassemblies [55]
• Calculate assembly completeness using BUSCO with viridiplantae_odb10 dataset
• Verify structural accuracy through PCR amplicon sequencing across repeat junctions
• Assess base-level accuracy by mapping Illumina reads to the final assembly

Protocol 2: NUMT and MTPT Exclusion Strategy

Principle: This protocol specifically addresses the challenge of nuclear mitochondrial DNA (NUMT) and mitochondrial plastid DNA (MTPT) contamination that disproportionately affects parasite mitochondrial assemblies [39] [50].

Step-by-Step Procedure:

• Mitochondrial Read Enrichment
  • Align sequencing reads to plant mitochondrial core gene database using Minimap2
  • Extract reads with alignment length >50 bp and multiple core gene matches as candidate mitochondrial reads
  • Perform k-mer frequency analysis to distinguish mitochondrial from nuclear and plastid reads
  • Apply copy-number based filtering using KMCP or similar tools

• Assembly with NUMT-Aware Tools

  • Utilize PMAT which employs copy number differences without pre-filtering to reduce NUMT incorporation
  • Implement MitoHiFi for specialized mitochondrial assembly with integrated contamination screening
  • Execute Oatk with k-mer based mitochondrial read enrichment followed by repeat resolution

• Contamination Identification and Removal

  • Align assembly contigs to reference nuclear and plastid genomes using BLASTn
  • Identify potential NUMTs/NUPTs through sequence similarity and read depth discrepancies
  • Verify mitochondrial origin through codon usage and substitution pattern analysis
  • Perform PCR validation of ambiguous regions using mitochondrial-specific primers

• Consensus Assembly Generation

  • Compare assemblies from multiple tools using QUAST-LG or similar comparison software
  • Generate consensus structures supported by at least two independent assembly methods
  • Resolve discrepancies through manual inspection of read alignments and coverage patterns
  • Verify assembly continuity through optical mapping or Hi-C data where available

Application Notes: This protocol is particularly valuable for parasite taxonomy studies where NUMT contamination can lead to incorrect phylogenetic inferences. The implementation requires balancing sensitivity (retaining genuine mitochondrial sequences) with specificity (excluding contaminants), which may require taxon-specific parameter optimization.

Case Studies and Applications in Parasite Taxonomy

Avian Haemosporidian Parasite Discrimination

The application of advanced mitochondrial assembly protocols has revolutionized detection and discrimination of avian haemosporidian parasites. In Lophura swinhoii, Nanopore sequencing enabled resolution of co-infections with multiple parasite lineages that remained undetected with conventional Sanger sequencing [52]. The implementation of long-read sequencing followed by careful assembly using repeat-aware algorithms identified two novel Haemoproteus lineages (hLOPSWI01 and hLOPSWI02) and one Plasmodium lineage (pNILSUN01) within the same host specimen [52].

The mitochondrial genome assembly provided sufficient phylogenetic signal to resolve the Haemoproteus lineages within the Parahaemoproteus clade and the Plasmodium lineage within the Giovannolaia-Haemamoeba clade [52]. This discrimination at the subgeneric level demonstrates the taxonomic precision enabled by complete mitochondrial genomes assembled with protocols that address repeat-induced collapse. For drug development professionals, this resolution enables targeted development of interventions against specific parasite lineages with different host specificities and pathogenic potentials.

Structural Variants as Taxonomic Markers

In oak species (Quercus), comprehensive mitochondrial assembly revealed six distinct structural variants ranging from 339 to 622 kb, with dispersed repeats identified as the primary drivers of mitochondrial genome expansion and structural dynamism [53]. The assembly of complete mitochondrial genomes from 15 phylogenetically representative species revealed remarkable intragenus variation in organellar gene inventories, encoding 34-41 genes, 20-28 tRNAs, and 2-5 rRNAs [53].

Phylogenomic analysis of 39 mitochondrial genes resolved deep evolutionary relationships in Quercus with clear cytonuclear discordance compared to nuclear phylogenies [53]. This demonstrates the value of mitochondrial structural variants as complementary taxonomic markers that can reveal distinct evolutionary histories and potential divergent selection pressures on cytoplasmic versus nuclear genomes. For parasite taxonomy, similar approaches could resolve cryptic species complexes where morphological differentiation is minimal but biological differences significant.

Addressing repeat-induced assembly collapse in complex mitochondrial genomes requires integrated experimental and computational approaches specifically designed to handle repetitive elements and structural complexity. The protocols outlined here, incorporating long-read sequencing technologies, multi-algorithm assembly strategies, and specialized repeat-resolution methods, provide robust solutions for generating complete and accurate mitochondrial genomes essential for parasite taxonomy research.

Future developments in this field will likely include more sophisticated algorithms capable of resolving heteroplasmy and structural variants at the population level, improved tools for distinguishing genuine mitochondrial sequences from NUMTs in diverse parasite taxa, and standardized validation metrics for assessing assembly completeness and accuracy. For researchers in parasite taxonomy and drug development, implementing these advanced assembly protocols will enable more precise species delineation, deeper understanding of evolutionary relationships, and identification of potential molecular targets for intervention.

Mitigating Nuclear Mitochondrial DNA Sequences (NUMTs) Contamination

Nuclear Mitochondrial DNA segments (NUMTs) are fragments of the mitochondrial genome (mtDNA) that have been inserted into the nuclear genome [56]. These sequences pose a significant challenge for mitochondrial genome assembly in parasite taxonomy research, as they can be co-amplified during PCR, leading to the misassembly of nuclear pseudogenes as mitochondrial sequences [57] [58]. This contamination risks incorrect phylogenetic tree inference and false identification of heteroplasmic variants, ultimately compromising taxonomic classification [56].

The transfer of mtDNA into the nuclear genome is an ongoing process, facilitated by mechanisms such as the repair of double-stranded DNA breaks via non-homologous end joining [56]. NUMTs can range in size from 24 base pairs to nearly the entire mitochondrial genome, and their sequence similarity to authentic mtDNA makes them particularly problematic for molecular studies [56]. For research focused on the mitochondrial genome assembly of parasites, implementing robust strategies to mitigate NUMT contamination is therefore a critical prerequisite for ensuring data integrity.

Experimental Protocols for NUMTs Prevention

Proactive, bench-level methods are essential for preventing NUMT contamination. The following protocols have been demonstrated to significantly reduce the co-amplification of nuclear pseudogenes.

Pre-PCR Dilution of DNA Template

This method leverages the natural abundance of mtDNA copies relative to nuclear NUMTs within a cell [57] [58].

• Principle: A single cell contains hundreds to thousands of mtDNA copies, while any specific NUMT locus is present in one or two copies per diploid nucleus. Extreme dilution of the DNA template can reduce the probability of amplifying low-copy-number NUMTs below the PCR detection threshold, while still allowing amplification of the high-copy-number authentic mtDNA [58].
• Protocol:
  • Quantify the genomic DNA extract using a fluorometric method.
  • Perform a series of template dilutions. Testing a range from 1-10 ng/μL down to 0.01-0.1 ng/μL is recommended [58].
  • Perform PCR amplification using the diluted templates.
  • The optimal dilution is identified as the highest dilution that still yields a robust PCR product for the target mtDNA fragment.
• Advantages: A simple, low-cost, and robust method that does not require specialized reagents [57] [58].

Long-Range PCR Amplification

This approach exploits the structural difference between the intact, circular mitochondrial genome and the typically short, fragmented NUMTs [58].

• Principle: Native mtDNA is a full-length, circular molecule (e.g., ~6 kb in apicomplexan parasites [3], ~16 kbp in metazoans), whereas most NUMTs are short fragments (mean length of 100-300 bp) [58] [56]. Designing primers to amplify a large, multi-kb fragment of the mtDNA selectively targets the genuine, full-length genome.
• Protocol:
  • Design primers targeting distant sites on the mitochondrial genome to generate amplicons of >3 kb.
  • Use a high-fidelity DNA polymerase blend optimized for long-range PCR.
  • A semi-nested or nested PCR approach can subsequently be performed using the long-range product as a template to increase specificity and yield for sequencing [58].
• Advantages: Directly targets the physical structure of the authentic mtDNA and can be combined with other methods.

Mitochondrial DNA Enrichment via Organelle Isolation

This method physically separates mitochondria from nuclei before DNA extraction, thereby removing the source of NUMTs [56].

• Principle: By isolating intact mitochondria through differential centrifugation, the nuclear genomic DNA (including NUMTs) is excluded. Subsequent DNA extraction yields a sample highly enriched for authentic mtDNA [56].
• Protocol:
  • Homogenize cell or tissue samples in an isotonic buffer.
  • Perform a series of differential centrifugations (e.g., low speed to remove nuclei and debris, followed by high speed to pellet mitochondria).
  • Use a commercial mitochondrial DNA isolation kit or alkaline lysis to purify mtDNA from the organelle pellet.
  • The enriched mtDNA can then be used for downstream applications like PCR or next-generation sequencing [56].
• Advantages: Provides a direct physical separation of mtDNA from nuclear DNA.
• Disadvantages: More labor-intensive, requires fresh or specially preserved tissue, and cross-contamination with nuclear DNA can occur if isolation is not optimal [57] [58].

cDNA Synthesis from RNA Template

This method targets the transcribed mitochondrial genome, as NUMTs are generally not transcribed [58].

• Principle: Functional mitochondrial genes are actively transcribed, while NUMTs integrated into the nuclear genome are typically silent. Reverse transcribing RNA into cDNA provides a template for PCR that is devoid of NUMTs [58].
• Protocol:
  • Extract total RNA from the sample, ensuring to remove any contaminating genomic DNA with DNase I treatment.
  • Use a reverse transcriptase enzyme and random hexamers or gene-specific primers to synthesize cDNA.
  • Perform standard PCR amplification using the cDNA as template.
• Advantages: Effectively avoids amplification of non-transcribed NUMTs.
• Disadvantages: Limited to profiling the expressed, coding regions of the mitochondrial genome; requires high-quality RNA [58].

The table below provides a comparative summary of these key experimental methods.

Table 1: Comparison of Experimental Methods for Preventing NUMT Contamination

Method	Principle	Key Advantage	Key Limitation
Pre-PCR Dilution [57] [58]	Dilutes single-copy NUMTs below PCR threshold	Simple, robust, and cost-effective	Requires empirical optimization of dilution factor
Long-Range PCR [58]	Amplifies long fragments absent in fragmented NUMTs	Targets structural integrity of mtDNA	Requires intact, high-quality DNA template
mtDNA Enrichment [56]	Physical separation of mitochondria from nuclei	Directly removes nuclear DNA (NUMTs source)	Labor-intensive; risk of nuclear contamination
cDNA Amplification [58]	Amplifies from transcribed mtDNA, not genomic NUMTs	Effectively avoids non-transcribed NUMTs	Limited to coding regions; requires high-quality RNA

Figure 1: An integrated experimental and computational workflow for mitigating NUMT contamination in mitochondrial genome studies.

Computational Identification of NUMTs

Despite preventive wet-lab measures, computational checks remain essential for identifying any residual NUMT contamination in sequencing data, especially with the sensitivity of next-generation sequencing [56].

Post-Sequencing Quality Control

The following analyses can be performed on putative mtDNA sequences to identify anomalies typical of NUMTs.

• Analysis of Coding Sequences: Examine protein-coding genes for the presence of indels that cause frameshifts or the appearance of premature stop codons. As NUMTs are non-functional, they accumulate such disruptive mutations, which are purged from functional mtDNA [58].
• Codon Position Substitution Bias: Calculate the distribution of nucleotide substitutions across codon positions. Authentic mtDNA concentrates synonymous changes in the 3rd codon position due to functional constraints. NUMTs, under no such constraint, accumulate mutations randomly across all positions [58].
• Phylogenetic Incongruence: Reconstruct a phylogenetic tree with the putative mtDNA sequences and compare it to a trusted tree based on morphology or nuclear markers. NUMTs often appear as outliers or form aberrant sister groups due to their different evolutionary history [58].

Bioinformatic Filtering Pipelines

For large-scale sequencing projects, systematic bioinformatic pipelines are required.

• Alignment-Based Filtering: Map sequencing reads to both a reference nuclear genome and a reference mitochondrial genome. Reads that map equally or better to the nuclear genome (especially to known NUMT regions) should be flagged or removed [56].
• Variant Filtering: In heteroplasmy studies, false positive variants from NUMTs can be identified by their characteristics. They often appear at a low variant allele frequency (VAF) and may have a lower sequencing quality score. Filtering out variants with VAF below a certain threshold (e.g., 1-2%) can help, though this may also remove true low-level heteroplasmies [56].
• K-mer-Based Detection: Advanced methods use k-mer profiles to distinguish the subtle sequence composition differences between true mtDNA and recently integrated NUMTs without relying solely on alignment [56].

Table 2: Computational Strategies for Identifying NUMT-Derived Sequences

Method	Key Indicator of NUMT	Applicable Scenario
Coding Sequence Analysis [58]	Frameshifts, premature stop codons	Analysis of assembled protein-coding genes
Codon Bias Analysis [58]	Elevated non-synonymous substitutions in 1st/2nd codon positions	Population genetics and evolutionary studies
Phylogenetic Incongruence [58]	Aberrant placement in mtDNA tree	Taxonomic and phylogenetic research
Read Alignment & Filtering [56]	Reads mapping to nuclear NUMT loci	All NGS-based mtDNA studies
Variant Allele Frequency (VAF) Filtering [56]	Low VAF variants; inconsistent with mtDNA copy number	Heteroplasmy and disease association studies

The Scientist's Toolkit: Essential Research Reagents & Materials

Successful mitigation of NUMTs requires a combination of specific laboratory reagents and bioinformatic tools. The following table details key solutions for designing a robust NUMT-contamination-free workflow.

Table 3: Research Reagent and Tool Solutions for NUMT Mitigation

Reagent / Tool	Function / Description	Application in NUMT Mitigation
High-Fidelity DNA Polymerase	Enzyme with proofreading activity for accurate long-range amplification.	Essential for Long-Range PCR protocol to generate full-length mtDNA amplicons without errors [58].
DNase I (RNase-free)	Enzyme that degrades single- and double-stranded DNA.	Critical for cDNA Synthesis protocol to remove contaminating genomic DNA from RNA samples prior to reverse transcription [58].
Mitochondrial Isolation Kit	Commercial reagent kit for purifying intact mitochondria via differential centrifugation.	Enables Mitochondrial Enrichment protocol, providing a physical barrier against NUMTs [56].
PacBio HiFi Reads	Long-read sequencing technology with high accuracy (>99.5%).	Ideal for sequencing long-range mtDNA amplicons, allowing for full-length haplotype reconstruction and easier detection of structural inconsistencies caused by NUMTs [7].
BLAST+ Suite	Basic Local Alignment Search Tool command-line version.	Foundational for Computational Identification; used to align putative mtDNA sequences or reads against nuclear genome databases to identify NUMT-homologous regions [58] [56].
MITOS2 / GeSeq	Web-based platforms for automated annotation of mitochondrial genomes.	Helps in the initial quality control by identifying and annotating genes, making it easier to spot frameshifts and stop codons indicative of NUMTs [30] [3].

Strategies for Resolving Regions with Low Coverage or Missing Genes (e.g., atp8)

Within the specific context of parasite taxonomy research, the assembly of complete and accurate mitochondrial genomes is a critical endeavor. A persistent challenge in this field is the accurate annotation of genes that are frequently missing from automated annotations or reside in regions of low sequencing coverage, with the ATP synthase F0 subunit 8 (atp8) gene being a prime example [59] [60] [61]. The inability to correctly identify such genes can hamper subsequent phylogenetic analyses and obscure the true genetic capabilities of the organism under study. This Application Note provides a consolidated strategic and technical guide for resolving these problematic genomic regions, enabling more reliable mitochondrial genome assembly for precise parasite classification and evolutionary studies.

The Challenge of Missing Genes and Low Coverage

The atp8 gene is notoriously difficult to annotate in mitochondrial genomes due to its small size, highly variable sequence, and divergent nature [59] [60]. Consequently, its absence from many published mitogenomes may often be an annotation artifact rather than a true biological deletion. For instance, a comprehensive analysis of Mytilidae mussels demonstrated that the atp8 gene, previously thought to be missing, could be identified through manual re-annotation, revealing characteristic transmembrane domains and hydropathy profiles [59] [60]. Similarly, the mitogenome of the acanthocephalan parasite Longicollum pagrosomi was reported to lack the atp8 gene [61]. These annotation gaps can lead to incomplete genetic characterizations, potentially misinforming taxonomic and phylogenetic inferences.

Low coverage regions in mitochondrial assemblies typically arise from technical limitations related to the sequencing technology employed. Homopolymer regions are particularly problematic for certain long-read sequencing technologies, leading to indels that compromise assembly accuracy [62]. Furthermore, complex repetitive sequences and structural variations can cause assembly algorithms to break, resulting in incomplete drafts [63]. In parasite research, these challenges are compounded by the difficulty of isolating pure mitochondrial DNA away from host tissue contamination [8].

Strategic Approaches and Experimental Protocols

A multi-faceted approach, combining advanced sequencing technologies with specialized bioinformatic tools and manual curation, is essential for producing high-quality mitochondrial assemblies.

Sequencing Technology and Assembly Selection

The choice of sequencing technology and assembler fundamentally influences the ability to resolve difficult regions. The table below benchmarks the capabilities of different technologies.

Table 1: Comparison of Sequencing and Assembly Strategies for Mitogenome Completion

Strategy	Typical Use Case	Advantages	Limitations	Key Tools/Examples
Illumina Short-Reads	Gold-standard reference assembly [62]	High base-level accuracy	Struggles with long repeats; time-consuming workflow [62]	GetOrganelle [62], NOVOPlasty [60]
Nanopore Long-Reads	Rapid in-situ assembly; resolving co-infections [52]	Long reads span repetitive regions; fast turnaround	Higher initial error rate, especially in homopolymers [62]	Customized de-novo & reference-based workflows [62]
Hybrid Assembly	Complex plant mitogenomes [63]	Leverages accuracy of short reads and continuity of long reads	Complex workflow; requires multiple data types	SPAdes (with --meta) [8]
Specialized Assemblers	Standardized, accurate mitogenome assembly	Optimized for organellar genomes; improves completeness	May not be suitable for all organism types	PMAT (for plants) [63], MITOS2 [62]

Experimental Protocol 1: Mitochondrial Genome Assembly using Low-Coverage Nanopore Sequencing This protocol is adapted from a study on the silky shark [62] and haemosporidian parasites [52], emphasizing its utility for rapid characterization.

• DNA Extraction: Use a high-quality tissue DNA extraction kit, preferably from purified parasites to minimize host contamination.
• Library Preparation & Sequencing: Perform low-pass whole genome sequencing using a portable MinION ONT platform. As little as 1x coverage of the genome can be sufficient for mitogenome assembly with specialized tools [63].
• In-situ De Novo Assembly: Assemble reads using a long-read assembler (e.g., SMARTdenovo, Canu). Subsequently, circularize the contig to form the complete mitochondrial chromosome.
• Benchmarking (Optional but Recommended): Compare the long-read assembly to a "gold-standard" mitogenome assembled from Illumina short-reads from the same sample to quantify accuracy, focusing on indel errors in homopolymer regions [62].

Manual Annotation and Verification of Missing Genes

When automated pipelines fail to identify genes like atp8, manual re-annotation is required.

Experimental Protocol 2: Manual Annotation of the atp8 Gene This protocol is based on the successful strategy employed for Mytilidae mussels [59] [60].

• Identify Candidate Regions: Scan intergenic regions of the assembled mitogenome using ORFfinder to locate all possible Open Reading Frames (ORFs), regardless of their length.
• Validate Start Codons: Manually inspect and correct the start codons of candidate ORFs by comparing them to atp8 sequences from closely related species.
• Predict Transmembrane Domains: Analyze the hydropathy profile of the predicted protein sequence using the PROTSCALE tool on ExPASy. Submit the sequence to TMHMM Server v.2.0 to identify the presence of characteristic transmembrane helices. A true atp8 should typically contain at least one transmembrane domain [60].
• Check for Charged C-Termini: Verify that the C-terminal region of the predicted protein is enriched with positively charged amino acids, a conserved feature of atp8 [60].
• Advanced Profile Search (Optional): For greater confidence, perform a Hidden Markov Model (HMM)-based search using HHblits. Construct an HMM for the candidate ORF and align it against a database of known atp8 HMMs [60].

Table 2: Key Tools for Manual Verification of Difficult Genes like atp8

Tool	Function	Application in Protocol
ORFfinder	Finds all possible Open Reading Frames	Identifies candidate atp8 sequences in intergenic regions [60]
TMHMM Server v.2.0	Predicts transmembrane helices	Confirms presence of a transmembrane domain, supporting atp8 identity [60]
PROTSCALE (ExPASy)	Calculates hydrophobicity profiles	Validates hydropathy profile consistent with atp8 [60]
HHblits	Performs HMM-HMM alignment	Provides robust, homology-based evidence for atp8 annotation [60]
MITOS2 Web Server	Automated mitogenome annotation	Provides a baseline annotation to be verified and corrected manually [8]

The following diagram illustrates the logical workflow for resolving missing genes and low-coverage regions, integrating both bioinformatic and experimental strategies.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Reagents and Kits for Mitochondrial Genome Assembly

Item	Function/Application	Example/Reference
TIANGEN Marine Animal Tissue DNA Kit	DNA extraction from parasitic or marine samples, critical for obtaining high-quality input DNA.	Used for Didymozoidae parasites from yellowfin tuna [8].
E.Z.N.A. Mollusc DNA Kit	DNA extraction optimized for molluscs and other challenging invertebrate taxa.	Used for genomic DNA extraction from planarians [26].
NEBNext Ultra DNA Library Prep Kit	Preparation of sequencing libraries for the Illumina platform for high-accuracy "gold-standard" assemblies.	Used for constructing libraries for Didymozoidae parasite sequencing [8].
Illumina NovaSeq 6000 Platform	Generating high-coverage, accurate short reads for benchmarking or hybrid assembly.	Used for sequencing the Didymozoidae parasite mitogenome [8].
MinION ONT Pocket-Sized Sequencer	Portable, real-time long-read sequencing for rapid in-situ mitogenome assembly.	Used for assembling the silky shark and haemosporidian parasite mitogenomes [62] [52].
Trimmomatic	Quality control tool for trimming adapter sequences and low-quality bases from raw sequencing reads.	Used for filtering raw data in Mytilidae mussel mitogenome studies [60].

Concluding Remarks

The path to a complete and accurate mitochondrial genome requires more than just automated pipelines. As demonstrated in studies across diverse taxa, from marine mussels to flatworm parasites, the strategic integration of long-read sequencing technologies, specialized assembly toolkits, and, crucially, meticulous manual annotation is paramount for resolving problematic regions and uncovering missing genes like atp8 [59] [60] [61]. For researchers in parasite taxonomy, adopting these comprehensive strategies will yield more reliable mitogenomes, thereby strengthening the foundation for phylogenetic analysis, species identification, and understanding of evolutionary adaptation.

The accurate assembly of mitochondrial genomes is a cornerstone of phylogenetic studies and parasite taxonomy. Traditional sequencing methods, particularly those relying on short-read technologies, face significant challenges in resolving structural variants, determining the phase of mutations, and distinguishing between genuine mitochondrial DNA and nuclear mitochondrial sequences (NUMTs) [64] [65]. These limitations can obscure the true genetic relationships between parasitic species and hinder the resolution of complex taxonomic classifications. The advent of CRISPR-Cas9-based enrichment coupled with long-read sequencing technologies now enables the acquisition of complete, phased mitochondrial genomes from native DNA molecules. This technical advance provides a powerful tool for researchers in parasitology, allowing for the precise delineation of species and strains based on full-length, haplotype-resolved mitochondrial data, free from the biases of amplification [66] [67].

Core Principle: Cas9-Guided Enrichment for Full-Length mtDNA Sequencing

The fundamental innovation of this method is the use of the RNA-guided Cas9 nuclease to selectively cleave the target mitochondrial genome, which simultaneously enriches it from the total genomic DNA and defines the start and end points for sequencing reads. This amplification-free approach is crucial for preserving the native state of the DNA and avoiding the recombination artifacts often introduced by PCR [64] [66].

The process capitalizes on the topology of the mitochondrial genome. By dephosphorylating all free 5' ends in the genomic DNA sample, ligation of sequencing adapters is blocked for the vast majority of nuclear DNA. A sequence-specific Cas9 cleavage is then performed on the circular mtDNA, creating a double-strand break with a defined 5' phosphate group. This allows for selective ligation of sequencing adapters only to the ends created by Cas9 on the mtDNA. Consequently, during a long-read sequencing run, the instrument is primed to sequence reads that initiate from the Cas9 cut site, continuing through the entire mitochondrial genome until the read concludes at the same cut site, thereby generating full-length, single-molecule sequences [65] [66]. This "cut-site as barcode" strategy can also be leveraged to multiplex samples from different parasites, each cleaved with a unique guide RNA, within a single sequencing run [66].

Experimental Protocol: A Step-by-Step Guide

Laboratory Workflow for Cas9-mtDNA Enrichment and Sequencing

The following protocol, adapted from current methodologies, details the steps for preparing a sequencing library for full-length mtDNA from parasite samples [64] [66].

Step 1: DNA Extraction and Quality Control

• Extract high-molecular-weight genomic DNA from the parasite sample using a method that preserves long DNA fragments (e.g., phenol-chloroform extraction or a commercial kit designed for long-read sequencing).
• Assess DNA integrity using pulsed-field gel electrophoresis or a Fragment Analyzer. For high-integrity DNA, an optional pre-enrichment step using Exonuclease V can be included to digest linear nuclear DNA, thereby enriching circular mtDNA and reducing NUMT interference [66].

Step 2: Dephosphorylation of Nuclear DNA

• Treat 3 µg of input gDNA with a quick dephosphorylation enzyme (e.g., Quick CIP from New England Biolabs) to remove 5' phosphate groups from all linear DNA fragments.
• Inactivate the enzyme according to the manufacturer's instructions [64]. This step is critical to prevent adapter ligation to nuclear DNA.

Step 3: Sequence-Specific Cleavage with Cas9

• Design guide RNA (gRNA): Design a crRNA targeting a unique ~20 nt sequence in the mitochondrial genome of the target parasite. The target site should be chosen to minimize off-target effects. For multiplexing, design a unique gRNA for each sample.
• Form Ribonucleoprotein (RNP) Complex: Mix the crRNA and tracrRNA (each at 10 µM) to form a duplex. Then combine with a high-fidelity Cas9 nuclease (e.g., HiFi Cas9) in a suitable buffer (e.g., 1X CutSmart Buffer).
• Cleave mtDNA: Add 10 µL of the 333 nM RNP complex to the dephosphorylated DNA. Incubate at 37°C for 30-60 minutes to allow for targeted cleavage [64] [66].

Step 4: Ligation of Sequencing Adapters

• Terminate the Cas9 reaction by adding Proteinase K to digest the nuclease.
• Purify the DNA and ligate long-read sequencing adapters (e.g., Oxford Nanopore Ligation Sequencing Kit, LSK-109) directly to the phosphorylated ends created by Cas9 cleavage [64].
• If multiplexing, pool the samples cleaved with different gRNAs at this stage.

Step 5: Sequencing and Data Collection

• Load the final library onto a long-read sequencer (e.g., Oxford Nanopore GridION or PromethION).
• Sequence for an appropriate number of hours to achieve sufficient coverage (e.g., >500x per mtDNA genome). Base-calling can be performed in real-time or post-run using the sequencer's software (e.g., Guppy) [64].

Bioinformatic Analysis Workflow

The raw sequencing data requires a specialized pipeline to demultiplex samples and call variants accurately [66].

• Demultiplexing: Align reads to a whole-genome reference using minimap2. Reads are assigned to a sample based on the proximity of their start and/or end coordinates to the expected Cas9 cut-site for each gRNA (within a 100 bp window). The most stringent demultiplexing strategy ("Both") selects only full-length reads that start and end near the same cut-site [66].
• Variant Calling and Phasing: Use a custom-developed tool like baldur [66] or a pipeline that incorporates long-read variant callers (e.g., Clair3 [66]) to identify single nucleotide variants (SNVs) and small indels. The long reads physically link variants, allowing for the determination of phase (haplotype) directly from the data.
• Structural Variant Detection: The same long reads are scanned for large deletions or other structural variants. The ability to span these events entirely allows for precise mapping of their breakpoints [64] [65].
• Annotation and Feature Table Generation: For taxonomic purposes, generate an annotated mitochondrial genome. The tool aln2tbl.py can create a feature table from a manually curated alignment, which is then used with tbl2asn to create a submission-ready file for public databases [68].

Diagram 1: Integrated laboratory and computational workflow for Cas9-based full-length mtDNA sequencing.

The Scientist's Toolkit: Essential Research Reagents and Materials

Successful implementation of this technique relies on a set of key reagents and tools, summarized in the table below.

Table 1: Key Research Reagent Solutions for Cas9-based mtDNA Sequencing

Item	Function/Description	Example Products/Suppliers
High-Fidelity Cas9	RNA-guided endonuclease for precise mtDNA cleavage.	HiFi Cas9 V3 (IDT) [64]
crRNA & tracrRNA	Custom guide RNA components that direct Cas9 to the target mtDNA site.	Integrated DNA Technologies (IDT) [64]
Dephosphorylation Enzyme	Removes 5' phosphates from linear DNA to prevent non-specific adapter ligation.	Quick CIP (New England Biolabs) [64]
Long-Read Sequencing Kit	Provides enzymes and adapters for library preparation.	Oxford Nanopore Ligation Sequencing Kit (LSK-109) [64]
Exonuclease V	Optional enzyme for pre-enrichment of circular mtDNA by degrading linear gDNA.	Available from molecular biology suppliers [66]
Bioinformatics Tools	Software for demultiplexing, variant calling, and genome annotation.	minimap2, baldur, Clair3, aln2tbl.py [64] [68] [66]

Application and Performance: Data from Model Studies

This method has been rigorously tested and shown to overcome the key limitations of short-read sequencing for mitochondrial genomics.

Resolving Heteroplasmy and Phasing Mutations

In a proof-of-concept study using blood from a patient with MELAS syndrome, long reads generated by Cas9-enrichment successfully phased two putative pathogenic mutations (m.1642A and m.5007A), revealing they existed on separate mtDNA molecules rather than co-segregating on a single, highly deleterious haplotype [64] [65]. This ability to determine phase is impossible with standard short-read sequencing and is critical for accurate genotype-phenotype correlation in parasitic organisms, where mixed infections or heteroplasmy are common.

Mapping Structural Variants

When applied to aged human muscle tissue, a known site of mtDNA deletion accumulation, the method readily identified and mapped the breakpoints of large deletions. This demonstrates its utility for discovering and characterizing structural variants in parasite genomes, which may be associated with drug resistance or adaptive evolution [64] [65].

The performance metrics of the optimized workflow, particularly when using the latest sequencing chemistries (e.g., Q20+), are superior to conventional methods.

Table 2: Quantitative Performance of Cas9-based mtDNA Sequencing

Performance Metric	Capability of Cas9-Based Method	Limitation of Short-Read Methods
Variant Phasing	Direct, physical phasing of variants across the entire genome [64] [66]	Statistical phasing with limited accuracy, impossible for distant variants [64]
Heteroplasmy Detection	Sensitivity down to <1% for single nucleotide variants [66]	Requires extremely high coverage and is confounded by NUMTs [64] [65]
Structural Variant Detection	Precise mapping of large deletion breakpoints and complex rearrangements [64] [66]	Ineffective for large deletions and complex structural variants [65]
Coverage Uniformity	Even coverage across the genome with "Both" demux strategy [66]	Significant coverage bias due to GC-content and amplification [64]
Multiplexing	Yes, using gRNA cut-sites as barcodes [66]	Requires separate index ligation steps

Within the field of parasite taxonomy research, a high-quality mitochondrial (mt) genome assembly is a prerequisite for reliable phylogenetic analysis and species identification. The compact nature, lack of introns, and general absence of recombinant events make the mitochondrial genome an ideal genetic marker for resolving evolutionary relationships [26]. However, the assembly process, often relying on tools designed for larger nuclear genomes, can introduce errors that compromise downstream biological interpretations. Therefore, rigorous quality control (QC) is not a final step but an integral part of the mitochondrial genome assembly pipeline. This protocol details the essential quantitative metrics and experimental methodologies for validating both the completeness and the gene content of mitochondrial genome assemblies, with a specific focus on applications in parasitology.

Key Quantitative Metrics for Assembly Quality Assessment

A multi-faceted approach to quality assessment is crucial. The metrics below provide a comprehensive picture of assembly integrity, from global architecture to base-level accuracy.

Table 1: Core Metrics for Genome Assembly Quality Assessment

Metric Category	Specific Metric	Definition and Interpretation	Optimal Value for mt-Genomes
Contiguity	Number of Contigs	Total sequences in the assembly.	1 (indicating a single, circular genome) [30] [8]
	N50/L50	N50: length of the shortest contig such that 50% of the total assembly length is contained in contigs of at least this size. L50: the number of contigs at the N50 length.	N50 should be equal or close to the full mt-genome length (~16-18 kb for many parasites) [69].
Completeness	Total Assembly Length	Total base pairs in the assembly.	Should match the expected size for the clade (e.g., ~16.5 kb for Didymozoidae [8], ~17.5 kb for Chaunocephalus ferox [30]).
	BUSCO Score	Percentage of universal single-copy orthologs found in the assembly [70]. Measures gene space completeness.	High C (Complete), low D (Duplicated), F (Fragmented), and M (Missing). A near-complete set of 12-14 PCGs is expected.
Gene Content	Gene Count	Number of protein-coding genes (PCGs), tRNAs, and rRNAs identified.	Typically 12-13 PCGs, 2 rRNAs, and a variable number of tRNAs (often 22-26) [30] [8].
	Missing Genes	Identification of commonly absent genes.	atp8 is frequently missing in many flatworm mt-genomes [26] [30] [8].
Base-Level Accuracy	QC Value / Read Mapping	A k-mer based measure of consensus quality [70]. Visual inspection of read coverage and variants.	QV > 50 is considered high-quality. Uniform read coverage and few variants suggest accurate assembly.

Interpreting Key Metric Outcomes

• Contiguity and Completeness: A perfect mitochondrial assembly is a single, circular contig with telomere-to-telomere (T2T) resolution, as achieved for several chromosomes in the Magnaporthe oryzae nuclear genome [71]. The presence of multiple contigs may indicate a fragmented assembly or, in some cases, nuclear mitochondrial sequences (NUMTs) misassembled as the genuine mt-DNA.
• BUSCO for Mitochondria: While BUSCO is often used for nuclear genomes, the principle can be applied to mitochondrial genomes using specific lineage datasets. It assesses whether the expected set of core mitochondrial genes is present and intact [70].
• Gene Content Anomalies: The consistent absence of the atp8 gene in trematodes like Chaunocephalus ferox and Didymozoidae parasites is a known feature and not an assembly error [30] [8]. However, the absence of other core genes like cox1 or cob typically indicates a problem with assembly completeness.

Experimental Protocols for Validation

Protocol 1: Assembly Completeness Assessment with BUSCO

This protocol evaluates the completeness of a genome assembly by searching for universal single-copy orthologs.

1. Software and Data Setup

• Tool: compleasm or BUSCO.
• Prerequisites: Assembled mt-genome in FASTA format.
• Database: A appropriate lineage dataset (e.g., eukaryota_odb12).

2. Step-by-Step Procedure

• Environment Setup: Load necessary modules and activate the software environment.

• Execution: Run the analysis specifying the assembly, database, and computational resources.
  Parameters: -t 8 uses 8 threads; -l and -L specify the lineage database; -a is the input assembly; -o is the output directory [70].

• Interpretation: Analyze the output short_summary.json file. Key results are:
  • Complete (C): The percentage of BUSCO genes found as single-copies. Expected to be very high for a complete mt-genome.
  • Fragmented (F): The percentage found as partial sequences.
  • Missing (M): The percentage not found in the assembly.

Protocol 2: K-mer Based Quality Assessment with Merqury

This protocol evaluates base-level accuracy and assembly completeness by comparing k-mer frequencies between the raw sequencing reads and the final assembly.

1. Software and Data Setup

• Tool: Merqury.
• Prerequisites: Assembled mt-genome in FASTA format; original long-read (HiFi) or short-read data used for assembly.
• K-mer Size: Determine the optimal k-mer size based on genome size.

2. Step-by-Step Procedure

• Optimal K-mer Size: Calculate the best k-value for your data.

• K-mer Counting: Count k-mers present in the original sequencing reads.

• Merqury Evaluation: Run Merqury to compare the assembly to the read k-mers.

• Interpretation: Merqury generates several outputs, including:
  • QV (Quality Value): A log-scaled probability of an error per base. A higher QV (e.g., >50) indicates higher accuracy.
  • Completeness: The percentage of k-mers from the reads that are found in the assembly, indicating assembly completeness [70].

Workflow Visualization

The following diagram illustrates the logical workflow for validating mitochondrial genome assembly quality, integrating the metrics and protocols described above.

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Research Reagents and Tools for mt-Genome Assembly & QC

Item / Resource	Function / Purpose	Example / Notes
DNA Extraction Kit	High-quality, high-molecular-weight DNA extraction from parasite tissue.	TIANGEN Marine Animal Tissue DNA Extraction Kit [8].
Long-Read Sequencing	Generates long reads capable of spanning repetitive regions.	PacBio HiFi sequencing provides high accuracy [71] [69].
Assembly Software	De novo assembly of long reads into contiguous sequences (contigs).	Hifiasm [70], NextDenovo, NECAT, Flye (balanced performance) [72].
Annotation Pipeline	Automated identification and annotation of genes in the assembled genome.	MITOS2 Web Server [30] [8].
Quality Assessment Tools	Evaluation of assembly contiguity, completeness, and accuracy.	BUSCO/compleasm [70], Merqury [70], QUAST.
Visualization Software	Graphical representation of the circular mitochondrial genome.	OGDRAW [26], CGView [8].

The reliability of mitochondrial genome assemblies for parasite taxonomy is non-negotiable. By systematically applying the quantitative metrics and experimental protocols outlined in this document—focusing on contiguity, completeness, gene content, and base-level accuracy—researchers can confidently validate their assemblies. This rigorous QC framework ensures that subsequent phylogenetic analyses and taxonomic conclusions are built upon a foundation of high-quality genomic data, thereby advancing the field of molecular parasitology.

Validation, Comparative Genomics, and Applications in Biomedical Research

Within the framework of mitochondrial genome assembly for parasite taxonomy research, the selection of appropriate molecular markers is fundamental for resolving evolutionary relationships. Mitochondrial protein-coding genes (PCGs), particularly cytochrome c oxidase subunit 1 (cox1) and cytochrome b (cob), have emerged as powerful tools for constructing robust phylogenetic trees [73]. Their utility stems from a balance of evolutionary rate characteristics and functional conservation. The maternal inheritance and general lack of recombination in mitochondrial DNA provide a clear lineage tracing path, while a faster mutation rate than nuclear DNA offers sufficient variation for distinguishing closely related species [73]. Constructing phylogenetic trees using concatenated sequences of these genes significantly enhances resolution and support for evolutionary branches, providing a reliable molecular basis for studies in population genetics, species identification, and taxonomic classification [3] [73].

The Scientific Rationale for Gene Selection

The selection of cox1 and cob is not arbitrary; it is grounded in their distinct biological properties and empirical performance in phylogenetic studies.

Functional Conservation and Evolutionary Rate

The cox1 and cob genes encode critical subunits of the mitochondrial electron transport chain. This essential function imposes selective constraints, ensuring that the genes remain conserved enough for alignment across diverse taxa while accumulating neutral substitutions useful for phylogenetics. Analysis of cox1 and cob genes in Apicomplexan parasites reveals they are generally subject to purifying selection, with Ka/Ks (non-synonymous to synonymous substitution rate) ratios less than 1, which preserves their protein function while allowing for the accumulation of phylogenetically informative synonymous substitutions [3] [74].

Superiority in Phylogenetic Resolution

Single-gene analyses can sometimes yield unresolved or weakly supported phylogenies. The concatenation of cox1 and cob into a single super-alignment provides a larger number of informative sites for phylogenetic analysis. This approach was successfully used to resolve the relationships among five critical Bipolaris phytopathogen species, resulting in a well-supported phylogenetic tree [74]. Similarly, in studies of Theileria parasites, a combined dataset of cox1 and cob was employed to clarify evolutionary relationships when a third gene, cox3, was excluded due to excessive sequence variation [3].

Table 1: Characteristics of cox1 and cob genes in various parasite taxa.

Taxonomic Group	Gene Length (bp)	Evolutionary Characteristic	Primary Use in Phylogenetics
Apicomplexan Parasites [3]	cox1: 1,428; cob: ~1,200	Highly conserved; under purifying selection	Resolving inter- and intra-species relationships
Bipolaris Fungi [74]	~1,500-1,600 each	Low genetic distance (conserved)	Species delineation and genus-level phylogeny
Plant-Parasitic Nematodes [73]	~1,000-1,500 each	Rapidly evolving; high sequence variation	Distinguishing cryptic species complexes

Wet-Lab Protocol: From Sample to Sequence

This section details a standard workflow for obtaining cox1 and cob sequences from parasite samples, incorporating methods from several studies.

Sample Collection and DNA Extraction

• Sample Source: The protocol can be applied to various samples, including parasitized blood [3], infected plant tissue [74], or environmental samples like soil containing nematodes [73].
• DNA Extraction: Use commercial kits, such as the TIANamp Genomic DNA Kit, following the manufacturer's instructions [3]. The quality and quantity of the extracted DNA should be assessed via spectrophotometry or fluorometry. For complex samples, a nematode suspension step is recommended to enrich the target organism before DNA extraction [73].

Mitochondrial Gene Amplification and Sequencing

Two primary strategies are employed:

• Polymerase Chain Reaction (PCR): Use dinoflagellate-specific PCR primers (or other taxon-specific primers) to isolate the cox1 and cob genes and produce corresponding cDNAs [75]. This is ideal for targeted sequencing of specific genes.
• High-Throughput Sequencing (HTS): For a more comprehensive approach, subject the extracted DNA to whole-genome sequencing on platforms like Illumina Novoseq 6000 [3] or use long-read technologies such as Oxford Nanopore [73]. This method is particularly powerful for metagenomic samples or for assembling complete mitogenomes from which cox1 and cob can be extracted bioinformatically.

In Silico Protocol: Phylogenetic Tree Construction

The following workflow outlines the computational steps from raw sequence data to a finalized phylogenetic tree.

Data Preparation and Multiple Sequence Alignment

• Genome Assembly: Assemble clean reads from HTS into contigs using assemblers like IDBA [3]. Annotate the resulting mitogenome using web servers such as MITOS to identify and extract the cox1 and cob gene sequences [3].
• Multiple Sequence Alignment: Perform alignment of the extracted cox1 and cob nucleotide or amino acid sequences using algorithms such as ClustalW, which is implemented in software like MEGA [3]. Manually inspect and refine alignments to ensure accuracy.

Concatenation and Evolutionary Model Selection

• Gene Concatenation: Combine the aligned cox1 and cob sequences into a single, partitioned dataset. The order is typically cox1 followed by cob.
• Model Selection: Determine the best-fit model of nucleotide or amino acid substitution for the concatenated dataset using programs like ProtTest for protein sequences [75]. Common models for mitochondrial proteins include the JTT (Jones-Taylor-Thornton) model with a heterogeneity parameter (e.g., "+G") [3].

Tree Building and Visualization

• Phylogenetic Inference: Construct the tree using robust statistical methods.
  • Maximum Likelihood (ML): Perform ML analysis with software like MEGA or RAxML, assessing branch support with 1,000 bootstrap replicates [3].
  • Bayesian Inference (BI): Use MrBayes to generate a tree with posterior probabilities, offering an alternative measure of clade credibility.
• Tree Visualization and Annotation: Visualize the final tree using tools like iTOL (Interactive Tree Of Life) [3] or R packages like phytools and ape [76] [77]. Use color-coding for branches or clades to represent different taxonomic groups or reconstructed ancestral states.

The Scientist's Toolkit

Table 2: Essential research reagents and software for phylogenetic analysis of mitochondrial PCGs.

Item Name	Function/Application	Specific Example/Use Case
TIANamp Genomic DNA Kit [3]	Extraction of high-quality total genomic DNA from parasite samples.	Used for DNA isolation from bovine blood infected with Theileria velifera.
Illumina Novoseq 6000 [3]	High-throughput sequencing platform for generating raw sequence reads.	Sequencing the complete mitochondrial genome of Theileria velifera.
MITOS Web Server [3]	Automated annotation of mitochondrial genomes.	Annotating PCGs, rRNAs, and identifying open reading frames.
ClustalW [3]	Multiple sequence alignment of nucleotide or amino acid sequences.	Aligning cox1 and cob sequences from multiple Apicomplexan parasites.
MEGA 11.0 Software [3]	Integrated tool for sequence alignment, model selection, and phylogenetic tree construction.	Performing Maximum Likelihood analysis with bootstrap validation.
JTT Model [3]	An empirical model of amino acid substitution.	Selected as the best-fit model for phylogenetic analysis of concatenated COX1 and COB proteins.
iTOL (Interactive Tree Of Life) [3]	Online tool for the display, annotation, and management of phylogenetic trees.	Visualizing and annotating the final phylogenetic tree for publication.

Troubleshooting and Data Interpretation

Even with a robust protocol, challenges in phylogenetic analysis are common. Below are key considerations for data interpretation and troubleshooting.

• Weak Branch Support: If bootstrap values or posterior probabilities are low, consider increasing the number of bootstrap replicates, re-checking the multiple sequence alignment for errors, or testing if the evolutionary model is still appropriate. In some cases, adding more taxa can help break up long branches and improve resolution.
• Handling Missing Data and Sequence Variation: The cox3 gene is sometimes excluded from phylogenetic analyses in Theileria and Babesia due to "significant variation" which can complicate alignment and model fitting [3]. A similar approach can be considered for other variable genes or regions that disrupt the phylogenetic signal.
• Color Mapping and Visualization: When coloring trees by traits or ancestral state reconstructions, ensure the color vector correctly maps to all possible states. If a state is absent from the tips but present in ancestors, explicitly define the state sequence (e.g., seq = LETTERS[1:5]) in R to prevent color misassignment [77].
• Introns and Genome Structure: Be aware that introns in fungal mitogenomes, particularly in the cox1 gene, can experience frequent gain/loss events and contribute significantly to size variation [74]. For PCR-based approaches, this can lead to amplification failures, making HTS a more reliable option for such taxa.

Application Note & Protocol

Framed within a thesis on mitochondrial genome assembly for parasite taxonomy research

Experimental Workflow for Comparative Mitogenomics

The following diagram illustrates the integrated bioinformatic pipeline for comparative mitogenomic analysis, from initial assembly to evolutionary insights, specifically contextualized for parasite research.

Quantitative Analysis of Nucleotide Composition

Protocol for Nucleotide Composition and Skew Analysis

Principle: Nucleotide composition skews are calculated to assess strand asymmetry and provide insights into mutational biases and replication mechanisms [78]. These analyses are particularly valuable for understanding evolutionary pressures on parasite mitochondrial genomes.

Procedure:

• Data Preparation: Extract complete mitochondrial genome sequences in FASTA format.
• Composition Calculation: Use bioinformatics tools (e.g., DNASTAR Lasergene, MEGA, or custom scripts) to calculate the percentage of each nucleotide (A, T, G, C) across the entire genome and for specific genetic elements [78] [79].
• Skew Analysis: Apply the following standard formulas to determine strand asymmetry [78]:
  • AT-skew = (A - T) / (A + T)
  • GC-skew = (G - C) / (G + C)
• Comparative Assessment: Calculate these metrics for multiple species and compare values across taxonomic groups to identify lineage-specific patterns.

Technical Notes: Skew values typically range from -1 to +1. Positive AT-skew indicates an excess of adenine over thymine, while negative GC-skew indicates a deficit of guanine relative to cytosine [79]. In Tortricidae mitogenomes, for example, consistently negligible AT-skews and negative GC-skews have been observed [79].

Nucleotide Composition Across Taxonomic Groups

Table 1: Comparative nucleotide composition and skew analysis across diverse taxonomic groups

Taxonomic Group	A+T Content (%)	AT-skew	GC-skew	Reference
Tortricidae (Lepidoptera)	80.7	~0.004	Negative	[79]
Theileria velifera (Apicomplexan parasite)	Not specified	Not specified	Not specified	[3]
Camallanus cotti (Nematode parasite)	Not specified	Not specified	Not specified	[80]
Archipini (Tortricidae tribe)	80.8	Not specified	Not specified	[79]
Olethreutini (Tortricidae tribe)	80.2	Not specified	Not specified	[79]

Interpretation: The high A+T content observed in Tortricidae (80.7%) is typical for insect mitogenomes and reflects mutational biases common in arthropods [79]. Similar compositional biases are observed in parasite mitogenomes, though specific values vary by taxonomic group [3] [80].

Protocol for Synteny Analysis

Synteny Analysis Workflow

The following diagram details the procedural workflow for conducting synteny analysis to identify genomic rearrangements, a method particularly relevant for studying genome evolution in parasites.

Procedure:

• Data Preparation: Collect complete, annotated mitochondrial genomes from multiple species in FASTA format.
• Multiple Genome Alignment: Perform whole-genome alignment using Mauve v2.4.0 or similar software (e.g., progressiveMauve for larger datasets) [78].
• Visualization: Examine the output visualization to identify conserved genomic blocks (locally collinear blocks) and rearrangement breakpoints.
• Rearrangement Mapping: Note the position and orientation of rearranged blocks, including inversions, translocations, and duplications.
• Evolutionary Analysis: Map rearrangement events onto phylogenetic trees to understand their evolutionary timing and significance.

Application in Parasite Research: In nematode parasites like Camallanus cotti, synteny analysis has revealed exceptionally high rates of gene rearrangement, including duplicated protein-coding genes and tRNAs, providing insights into the mechanisms of mitochondrial genome evolution in parasitic lineages [80].

Gene Order Conservation Metrics

Table 2: Synteny and gene order conservation across taxonomic groups

Taxonomic Group	Gene Order Conservation	Notable Rearrangements	Phylogenetic Utility
Tortricidae (Lepidoptera)	Highly conserved, typical Lepidoptera pattern	Minimal rearrangements	High for deep phylogeny
Camallanus cotti (Nematoda)	Extremely derived	Multiple gene duplications (6 PCGs, 6 tRNAs)	Lineage-specific signatures
Ganoderma (Fungi)	Highly conserved gene order	Collinearity: 82.93-92.02% similarity	Reliable for species delimitation
Theileria (Apicomplexan parasites)	Linear genome structure with terminal inverted repeats	Unique among eukaryotes	Distinguishes closely related species

Technical Notes: The percentage similarity in collinearity analysis for Ganoderma was calculated by comparing newly assembled mitogenomes with reference mitogenomes at the nucleotide level [81]. Gene order conservation is typically high within animal phyla but shows significant divergence in certain parasitic groups like nematodes [80].

The Scientist's Toolkit: Essential Research Reagents & Software

Table 3: Key research reagents and bioinformatic tools for comparative mitogenomics

Tool/Resource	Type	Primary Function	Application in Parasite Taxonomy
MEGA 6.0/11.0	Software	Evolutionary genetics analysis	Genetic distance calculation, phylogenetic tree construction [78] [3]
Mauve v2.4.0	Software	Multiple genome alignment	Synteny analysis and rearrangement detection [78]
MITOS/MITOS2	Web server	Automated mitogenome annotation	Gene boundary prediction, tRNA identification [3] [68]
SPAdes	Software	De novo genome assembly	Mitogenome assembly from NGS data [82]
Geneious Prime	Software	Sequence analysis platform	Reference-based assembly, annotation, visualization [82]
Aliview/Seaview	Software	Sequence alignment editor	Manual curation of gene annotations [68]
tbl2asn	Command-line tool	GenBank submission	Preparation of annotated genomes for submission [68]
Illumina NovaSeq	Sequencing platform	High-throughput sequencing	Mitogenome sequencing from total DNA [3]

Case Study: Apicomplexan Parasite Mitogenomics

Background: The mitochondrial genomes of Apicomplexan parasites like Theileria velifera exhibit highly derived characteristics that provide valuable taxonomic markers [3].

Methodology Implementation:

• Sequencing: Illumina NovaSeq 6000 platform for high-coverage sequencing of total genomic DNA from infected bovine blood samples [3].
• Assembly: IDBA software for de novo assembly of the mitochondrial genome from clean reads [3].
• Annotation: MITOS web server for initial annotation, followed by BLAST searches against GenBank for homologous sequence identification [3].
• Composition Analysis: DNAStar V7.1 for nucleotide composition analysis and CodonW 1.4.2 for relative synonymous codon usage (RSCU) calculation [3].

Key Findings: The mitochondrial genome of T. velifera is a linear molecule of 6,125 bp containing 3 protein-coding genes (cox1, cob, and cox3), 5 large subunit rRNA gene fragments, and terminal inverted repeats at both ends [3]. This structure is typical for Apicomplexan parasites and distinguishes them from the circular mitogenomes of most other eukaryotes.

Taxonomic Application: Phylogenetic analysis using concatenated amino acid sequences of cox1 and cob successfully resolved evolutionary relationships among Theileria species, providing a molecular framework for parasite classification and identification [3].

Advanced Protocol: Integrating Composition and Synteny Data for Phylogenetic Inference

Procedure:

• Data Matrix Preparation: Combine nucleotide composition metrics (A+T content, skew values) with synteny data (presence/absence of specific gene arrangements) into a phylogenetic character matrix.
• Phylogenetic Analysis: Perform maximum likelihood or Bayesian analysis using concatenated protein-coding gene sequences (e.g., 15 core PCGs) with appropriate substitution models [81].
• Character Mapping: Map compositional and syntenic characters onto the resulting phylogenetic tree to identify evolutionarily significant patterns.
• Statistical Testing: Assess the congruence between different data partitions (composition vs. gene order vs. sequence data) using likelihood-based tests.

Application in Parasite Taxonomy: This integrated approach has revealed that NUMTs (nuclear mitochondrial segments) exhibit non-random origins from mtDNA and are preferentially located in transposon-rich regions, providing insights into the evolutionary dynamics of mitochondrial sequence transfer to the nucleus in mammalian parasites [83].

Linking Mitochondrial Lineages to Clinical Outcomes and Drug Resistance Phenotypes

Mitochondrial genomics has emerged as a pivotal field for understanding the complex interplay between genetic lineages, clinical disease manifestations, and therapeutic resistance. This relationship is particularly pronounced in parasitology, where mitochondrial genome variations serve as crucial biomarkers for species identification, virulence assessment, and treatment response prediction. The mitochondrial genome offers distinct advantages for these studies, including higher copy number per cell than nuclear DNA, minimal recombination, and rapid evolution rates that provide robust phylogenetic signals.

This Application Note provides a comprehensive framework for utilizing mitochondrial genome assembly and analysis to link genetic lineages to clinical outcomes and drug resistance phenotypes, with specific application to parasite taxonomy research. We present standardized protocols for generating high-quality mitochondrial genomic data, analytical methods for correlating haplotypes with phenotypic traits, and visualization tools for interpreting complex relationships in mitochondrial biology.

Applications in Clinical and Drug Resistance Research

Mitochondrial lineage analysis provides critical insights for clinical management and therapeutic development. Several studies have demonstrated the utility of this approach across multiple disease contexts, with particular relevance for parasitic infections and cancer research.

Table 1: Applications of Mitochondrial Lineage Analysis in Disease Research

Application Area	Specific Utility	Research Context
Parasite Taxonomy & Biodiversity	Species identification, discovery of cryptic species, and understanding evolutionary relationships [84] [85]	Haemosporidian parasites (Plasmodium, Haemoproteus, Leucocytozoon) and nematodes (Heterakis)
Chemotherapy Resistance	Prognostic biomarker identification and resistance mechanism characterization [86] [87]	Laryngeal squamous cell carcinoma and various cancer cell lines
Infection Tracking	Monitoring emerging outbreaks and zoonotic transmission dynamics [88]	Mpox virus Clade Ib and H5N1 influenza
Drug Discovery	Identifying novel therapeutic targets in mitochondrial processes [89] [87]	Mitochondrial dynamics proteins (MFN1/2, DRP1, OPA1) and mitophagy pathways

The molecular mechanisms linking mitochondrial dynamics to drug resistance involve complex regulatory pathways. Mitochondrial fusion, fission, and mitophagy processes have been identified as significant contributors to treatment resistance across various cancer types [87]. For instance, dysregulation of mitochondrial dynamin-related proteins (MFN1, MFN2, and DRP1) correlates with proliferation and chemoresistance in multiple tumors. Similarly, mitophagy induction enables tumor cell survival under therapeutic pressure by maintaining mitochondrial homeostasis [87].

In parasitic diseases, mitochondrial genome analysis enables precise tracking of resistant strains. A recent study on haemosporidian parasites demonstrated that complete mitochondrial sequencing could resolve mixed infections and co-infections that would be mischaracterized using standard barcoding approaches [84]. This capability is crucial for understanding treatment failures and emerging resistance in endemic regions.

Protocol: Mitochondrial Genome Assembly for Parasite Taxonomy

This protocol describes a comprehensive method for obtaining complete mitochondrial genomes from parasite samples using long-read sequencing technology, optimized for haplotyping and lineage analysis.

Sample Preparation and DNA Extraction

• Sample Collection: Obtain biological samples containing parasite material (blood, tissue, or cultured isolates). Preserve immediately in appropriate storage buffers (e.g., 80% ethanol for nematodes) [85] or process directly for DNA extraction.
• DNA Extraction: Use column-based genomic DNA isolation kits (e.g., DNeasy Blood & Tissue Kit, Qiagen) following manufacturer's protocols with minor modifications:
  • Incubate samples with proteinase K for 3-6 hours to ensure complete lysis
  • Include RNase A treatment (10-20 μg/mL) for 10 minutes at room temperature
  • Elute DNA in nuclease-free water or low-EDTA TE buffer
  • Assess DNA quality and quantity using fluorometry (e.g., Qubit) and fragment analysis (e.g., TapeStation)

Mitochondrial Genome Amplification

• Primer Design: Utilize previously validated primers AE170 (5′-GAT TCT CTC CAC ACT TCA ATT CGT ACT TC-3′) and AE171 (5′-GAA AAT WAT AGA CCG AAC CTT GGA CTC-3′) that target approximately 6 kb of the haemosporidian mitochondrial genome [84]. For other parasite taxa, design primers based on conserved mitochondrial regions.
• PCR Amplification: Perform reactions in 50 μL volumes using high-fidelity polymerase (e.g., TaKaRa LA Taq Polymerase) with the following cycling conditions [84]:
  • Initial denaturation: 94°C for 1 minute
  • 30 cycles of:
    • Denaturation: 94°C for 30 seconds
    • Annealing/Extension: 68°C for 7 minutes
  • Final extension: 72°C for 10 minutes
• Product Purification: Visualize PCR products on 1% agarose gels, excise bands of correct size (~6 kb), and purify using gel extraction kits (e.g., QIAquick).

Library Preparation and Sequencing

• Library Construction: Prepare SMRTbell libraries using the PacBio HiFi library preparation kit with the following specifications:
  • Damage repair and end-repair of amplified PCR products
  • Ligation of universal hairpin adapters
  • Size selection using BluePippin or similar systems (select 4-8 kb fragments)
  • Final library quality assessment using Fragment Analyzer or Bioanalyzer
• Sequencing: Perform PacBio HiFi sequencing on Sequel IIe or Revio systems with [84]:
  • 30-hour movie times
  • On-instrument circular consensus sequencing (CCS) generation
  • Minimum 30X coverage per haplotype for accurate variant calling

Bioinformatic Analysis

• Read Processing: Use the HmtG-PacBio Pipeline [84] which incorporates:
  • Demultiplexing of barcoded samples
  • CCS read generation with minimum QV ≥ 30
  • Read quality filtering and length selection
• Haplotype Reconstruction: Employ machine learning algorithms (modified variational autoencoders) and clustering methods to identify mitochondrial haplotypes/species in each sample, specifically designed to detect mixed infections [84].
• Genome Assembly and Annotation:
  • Assemble mitochondrial genomes using GetOrganelle v1.7.7.0 [85]
  • Annotate protein-coding genes, rRNAs, and tRNAs using MITOS2 web server and MitoZ v3.6 [85]
  • Manually verify open reading frames using ORF Finder
  • Identify tRNA genes using BLAST against nematode tRNA databases [85]
• Phylogenetic Analysis:
  • Perform multiple sequence alignment of mitochondrial protein-coding genes (e.g., MUSCLE, Kalign)
  • Construct phylogenetic trees using Maximum Likelihood (IQ-TREE) and Bayesian Inference (MrBayes) methods [85]
  • Assess node support with bootstrap analysis (1000 replicates) and posterior probabilities

Experimental workflow for mitochondrial genome assembly and analysis

Data Analysis and Integration Methods

Mitochondrial Lineage Classification

Effective classification of mitochondrial lineages forms the foundation for correlating genetic variation with clinical phenotypes. The following analytical approaches are recommended:

• Barcode Index Number (BIN) System: Assign unique identifiers to mitochondrial lineages based on COI sequence clusters using the BOLD Systems platform [90]. Generate BIN Discordance Reports to identify taxonomic conflicts that may indicate misidentification or cryptic diversity.
• Species Delimitation: Apply the Assembled Species by Automatic Partitioning (ASAP) method using the Kimura (k80) ts/tv 2.0 model to objectively delineate species boundaries based on genetic distances [85]. Use multiple genetic markers (ITS, cox1, cox2, 12S) for robust delimitation.
• Distance Analyses: Calculate intra-specific and inter-specific genetic distances using appropriate substitution models (e.g., K2P, p-distance). Generate Barcode Gap Summaries to visualize the distribution of genetic variation within and between species [90].

Correlation with Clinical and Resistance Phenotypes

• Morphological Analysis: For parasitic nematodes, combine molecular data with detailed morphological characterization using light and scanning electron microscopy [85]. Document key diagnostic features and deposit voucher specimens in accessible collections.
• Drug Sensitivity Profiling: For cancer applications, correlate mitochondrial genomic features with drug response data using the "oncoPredict" R package to calculate half-maximal inhibitory concentration (IC50) values [86]. Compare sensitivity patterns between mitochondrial lineage groups.
• Immune Microenvironment Assessment: Analyze differences in immune cell infiltration and function between mitochondrial lineage groups using CIBERSORT and single-sample Gene Set Enrichment Analysis (ssGSEA) [86]. Calculate Tumor Immune Dysfunction and Exclusion (TIDE) scores to predict immunotherapy response.

Table 2: Key Analytical Tools for Mitochondrial Data Integration

Tool/Platform	Primary Function	Application Context
BOLD Systems [91] [90]	DNA barcode data management, BIN assignment, and distance analysis	Species identification and biodiversity studies
HmtG-PacBio Pipeline [84]	Haplotype reconstruction from long-read data using machine learning	Detecting mixed infections and co-infections
MITOS2 & MitoZ [85]	Mitochondrial genome annotation and visualization	Gene identification and genome feature mapping
oncoPredict R package [86]	Drug sensitivity prediction from genomic data	Correlating mitochondrial variants with chemoresistance
CIBERSORT & ssGSEA [86]	Immune cell infiltration quantification	Tumor microenvironment analysis

Mitochondrial Dynamics in Drug Resistance: Mechanisms and Visualization

The role of mitochondrial dynamics in drug resistance represents a crucial connection between organelle biology and therapeutic outcomes. Understanding these mechanisms provides insights for overcoming treatment resistance.

Mitochondria undergo constant fusion and fission processes regulated by specific protein complexes. Fusion is mediated by mitofusins (MFN1, MFN2) in the outer membrane and OPA1 in the inner membrane, while fission is primarily driven by DRP1 [87]. These dynamic processes maintain mitochondrial health and function, but their dysregulation contributes significantly to drug resistance.

In chemotherapy resistance, mitochondrial fusion enables content complementation between damaged and healthy mitochondria, allowing cancer cells to survive therapeutic insult. Conversely, fission facilitates the segregation of damaged components for selective removal through mitophagy [87]. Both processes have been implicated in various resistance mechanisms:

• Metabolic Adaptation: Elongated mitochondrial networks promote efficient ATP production under nutrient stress, enhancing tumor cell survival during treatment [87].
• Apoptosis Evasion: Balanced mitochondrial dynamics maintain membrane potential and prevent cytochrome c release, circumventing drug-induced apoptosis [92] [87].
• Mitophagy-mediated Quality Control: Selective removal of drug-damaged mitochondria through PINK1-Parkin mediated mitophagy enables cellular recovery post-treatment [92] [87].

Mitochondrial dynamics in drug resistance mechanisms

Research Reagent Solutions

Table 3: Essential Research Reagents for Mitochondrial Lineage Studies

Reagent/Category	Specific Examples	Application and Function
Primers for Mitochondrial Amplification	AE170/AE171 [84], genus-specific primers [85]	Amplification of mitochondrial targets from various parasite species
DNA Extraction Kits	DNeasy Blood & Tissue Kit (Qiagen) [84] [85]	High-quality genomic DNA extraction from diverse sample types
PCR Reagents	TaKaRa LA Taq Polymerase [84]	Long-range, high-fidelity amplification of mitochondrial genomes
Library Prep Kits	SMRTbell Express Template Prep Kit [84]	Preparation of libraries for PacBio HiFi sequencing
Mitochondrial Stains	MitoTracker dyes (Red, Green) [92], TOM20 antibodies [92] [89]	Visualization of mitochondrial morphology and network structure
Image Analysis Tools	MoDL (deep learning algorithm) [89]	Automated mitochondrial segmentation and function prediction
Bioinformatics Tools	HmtG-PacBio Pipeline [84], MITOS2 [85], MitoZ [85]	Mitochondrial genome assembly, annotation, and haplotype analysis

The integration of mitochondrial genome assembly with clinical outcome data provides a powerful framework for understanding treatment resistance and disease pathogenesis. The protocols and分析方法 outlined in this Application Note enable researchers to establish robust correlations between mitochondrial lineages and phenotypic traits, particularly in the context of parasite taxonomy and evolution.

Future directions in this field will likely include the development of standardized mitochondrial variant reporting frameworks, multi-omics integration approaches, and machine learning models for predicting clinical outcomes based on mitochondrial genomic features. The growing accessibility of long-read sequencing technologies will further enhance our ability to resolve complex mitochondrial haplotypes and their association with drug resistance across diverse pathological contexts.

Identifying Novel Drug Targets Through Essential Mitochondrial Pathways

Mitochondria are essential organelles that perform crucial functions beyond energy production, including calcium homeostasis, regulation of apoptosis, and biosynthesis of key metabolites [93] [94]. Over the past decade, mitochondrial dysfunction has been implicated in a wide spectrum of human diseases, making mitochondrial proteins (MPs) increasingly appealing targets for therapeutic intervention [95]. Current research indicates that approximately 20% of the mitochondrial proteome (312 out of an estimated 1,500 MPs) has known interactions with small molecules, suggesting MPs are highly targetable [95]. The unique structural and functional characteristics of mitochondria, particularly the electrochemical gradient across the inner mitochondrial membrane (IMM), enable selective targeting of drugs designed to modulate the function of this organelle for therapeutic gain [96]. Mitochondrial drug-targeting strategies open new avenues for manipulating mitochondrial functions, allowing for selective protection or eradication of cells in various diseases, including cancer, neurodegenerative disorders, cardiovascular conditions, and metabolic syndromes [96] [93] [94].

Table 1: Mitochondrial Functions and Associated Therapeutic Opportunities

Mitochondrial Function	Biological Process	Therapeutic Opportunity
Energy Production	Oxidative Phosphorylation	Enhance ATP synthesis in metabolic diseases
Calcium Homeostasis	Calcium Signaling	Modulate calcium-dependent cell signaling
Apoptosis Regulation	Permeability Transition	Induce apoptosis in cancer; inhibit in neurodegeneration
ROS Production	Redox Signaling	Antioxidant delivery for oxidative stress-related conditions
Metabolic Integration	TCA Cycle, β-oxidation	Regulate substrate utilization in metabolic syndrome

Mitochondrial Genome Insights from Parasite Taxonomy Research

The study of mitochondrial genomes in parasitic organisms provides critical insights into evolutionary relationships and identifies conserved, essential pathways that represent promising drug targets [30] [8] [3]. Recent advances in sequencing technologies have enabled the complete assembly and annotation of mitochondrial genomes from various parasites, revealing both conserved and unique features that can be exploited for therapeutic development.

In trematode parasites such as Chaunocephalus ferox, the complete mitochondrial genome spans 17,482 bp and encodes 12 protein-coding genes, 22 tRNAs, and 2 rRNAs, with notable absence of the atp8 gene [30]. Similarly, studies on Didymozoidae parasites from yellowfin tuna reveal mitochondrial genomes of 16,468 bp with 12 protein-coding genes and 19 tRNA genes, also lacking the atp8 gene [8]. This consistent absence of specific genes across parasite lineages highlights potentially divergent metabolic requirements and identifies lineage-specific adaptations that could serve as selective drug targets.

Apicomplexan parasites, including Theileria velifera, exhibit particularly streamlined mitochondrial genomes, with T. velifera possessing a linear monomer mitochondrial genome spanning 6,125 bp that encodes only 3 protein-coding genes (cox1, cob, and cox3) and contains 5 large subunit rRNA gene fragments [3]. The significant reduction in mitochondrial gene content in these parasites suggests increased reliance on host metabolic processes, offering potential targets for disrupting this parasitic dependency.

Figure 1: Mitochondrial Genome Analysis for Drug Target Identification. Parasite mitochondrial genomes reveal structural variations, gene content patterns, and unique features that inform conserved and lineage-specific drug targets.

Key Mitochondrial Drug Targeting Strategies

Delocalized Lipophilic Cations for Mitochondrial Accumulation

Delocalized lipophilic cations (DLCs) represent a primary strategy for targeting bioactive molecules to mitochondria [96]. These compounds preferentially accumulate in the mitochondrial matrix driven by the high mitochondrial membrane potential (Δψm), typically -150 to -180 mV [96] [94]. The accumulation ratio can reach several hundred-fold compared to the extracellular concentration, making DLCs exceptionally efficient for mitochondrial targeting [96]. This approach is particularly effective in cancer cells, which often exhibit higher Δψm compared to normal cells, enabling selective targeting [96].

Notable examples of DLC-based therapeutics include:

• MitoQ: A triphenylphosphonium (TPP+) conjugate of ubiquinone that functions as a potent mitochondrial-targeted antioxidant [96] [94]. MitoQ has demonstrated efficacy in models of ischemia-reperfusion injury and is currently in clinical trials for Parkinson's disease [96].
• Rhodamine 123 and MKT-077: These DLCs selectively localize to carcinoma cells and can be used to deliver pro-apoptotic compounds specifically to cancer cell mitochondria [96].
• SkQ1: A plastoquinone conjugated to TPP+ that exhibits potent antioxidant activity and has shown beneficial effects in various animal models of age-related diseases [94].

Targeting Mitochondrial Electron Transport Chain (ETC)

The electron transport chain represents a key target for modulating mitochondrial function, with multiple complexes offering distinct therapeutic opportunities [96] [94]. Complex I and III are primary sites of reactive oxygen species (ROS) production, making them targets for antioxidant strategies [94]. Additionally, inhibition of specific ETC complexes can trigger apoptosis in cancer cells, while partial uncoupling may reduce ROS production in neurodegenerative conditions [96].

Modulation of Mitochondrial Permeability Transition (MPT)

The mitochondrial permeability transition pore (mPTP) plays a critical role in cell death pathways and represents a promising target for cytoprotective therapies [96] [94]. Cyclosporin A (CsA) and its analogues inhibit MPT by binding to cyclophilin D, providing protection against ischemia-reperfusion injury in heart and brain [96]. Other compounds including sangliferhin and Ro 68-3400 also target MPT, demonstrating the therapeutic potential of this pathway [96].

Figure 2: Strategic Approaches to Mitochondrial Drug Targeting. Multiple strategies including delocalized lipophilic cations, electron transport chain modulation, permeability transition regulation, and apoptotic pathway targeting enable precise therapeutic interventions.

Experimental Protocols for Mitochondrial Drug Discovery

Protocol 1: Mitochondrial Genome Assembly and Analysis for Target Identification

Purpose: To identify essential and conserved mitochondrial pathways in parasites that may represent novel drug targets.

Materials and Reagents:

• Parasite samples or infected host tissues
• DNA extraction kit (e.g., TIANGEN Marine Animal Tissue DNA Extraction Kit)
• Illumina NovaSeq 6000 platform or Oxford Nanopore Technologies for sequencing
• MITOS2 or GeSeq annotation platforms
• MEGA software for phylogenetic analysis

Procedure:

• Sample Preparation and DNA Extraction: Isolate parasites from host tissues under sterile conditions. Extract genomic DNA using specialized kits, ensuring high molecular weight and purity suitable for long-read sequencing [8] [3].
• Library Preparation and Sequencing: Prepare sequencing libraries using appropriate kits (e.g., NEBNext Ultra DNA Library Prep Kit for Illumina). Sequence using either short-read (Illumina) or long-read (Oxford Nanopore, PacBio) platforms. For complex samples or co-infections, long-read technologies are preferred for resolving individual haplotypes [97].
• Genome Assembly and Annotation: Assemble clean reads into mitochondrial genomes using specialized assemblers (e.g., SPAdes, GetOrganelle). Annotate the assembled genome using MITOS2 web server or similar platforms to identify protein-coding genes, tRNAs, and rRNAs [8] [98] [85].
• Comparative Analysis and Target Identification: Compare gene content, arrangement, and sequences across related species. Identify conserved essential genes and lineage-specific adaptations. Perform phylogenetic analysis using concatenated protein sequences to resolve evolutionary relationships and identify clade-specific features [30] [3].

Protocol 2: High-Throughput Screening of Mitochondria-Targeted Compounds

Purpose: To identify and validate compounds that selectively modulate mitochondrial functions.

Materials and Reagents:

• Cell culture models (primary cells or established cell lines)
• Mitochondria-specific fluorescent dyes (TMRE, JC-1, MitoTracker)
• Oxygen consumption measurement systems (Seahorse XF Analyzer)
• Mitochondria-targeted compound libraries
• MitoQ and other reference compounds as controls

Procedure:

• Cell Culture and Treatment: Culture appropriate cell models (cancer cells for pro-apoptotic compounds, neuronal cells for protective compounds). Treat with candidate compounds across a concentration range and time course [96] [94].
• Assessment of Mitochondrial Membrane Potential: Stain cells with potential-sensitive dyes (TMRE, JC-1) after compound treatment. Quantify fluorescence shifts using flow cytometry or fluorescence microscopy to detect depolarization or hyperpolarization [96].
• Metabolic Profiling: Measure oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) using Seahorse XF Analyzer or similar systems. Calculate basal respiration, ATP-linked respiration, proton leak, and spare respiratory capacity [93] [94].
• Viability and Apoptosis Assessment: Perform MTT or Alamar Blue assays to assess cell viability. Use Annexin V/propidium iodide staining to quantify apoptosis induction. For antioxidant compounds, measure ROS production using H2DCFDA or MitoSOX Red [96] [94].
• Mechanistic Studies: Investigate compound effects on specific pathways including ETC complex activities, calcium homeostasis, MPT opening, and apoptotic protein release [96] [94].

Table 2: Key Research Reagent Solutions for Mitochondrial Drug Discovery

Reagent/Category	Specific Examples	Function/Application
Sequencing Platforms	Illumina NovaSeq 6000, Oxford Nanopore	Mitochondrial genome assembly and variation analysis
Annotation Tools	MITOS2, GeSeq, MitoZ	Gene identification and genome annotation
Targeting Moieties	Triphenylphosphonium (TPP+), Rhodamine	Mitochondrial accumulation of conjugated compounds
Reference Compounds	MitoQ, SkQ1, Cyclosporin A	Benchmark molecules for specific mitochondrial targets
Assessment Dyes	TMRE, JC-1, MitoTracker, MitoSOX	Measurement of membrane potential, mass, and ROS

Mitochondrial Signaling Pathways as Therapeutic Targets

Apoptotic Regulation Pathways

Mitochondria play a central role in the intrinsic apoptosis pathway, making them prime targets for therapeutic modulation in cancer and degenerative diseases [96] [94]. Key targets include Bcl-2 family proteins, which regulate outer mitochondrial membrane permeability, and the permeability transition pore, which controls the release of cytochrome c and other pro-apoptotic factors [96]. Experimental approaches targeting these pathways include:

• BH3 mimetics such as SAHB (stabilized alpha-helix of BCL-2 domains) that directly activate Bax to trigger apoptosis [96].
• Compounds targeting the benzodiazepine receptor including PK11195 that modulate the interaction between Bcl-2 family proteins and the permeability transition pore [96].
• SMAC-DIABLO mimetics that inhibit XIAP (inhibitor of apoptosis protein) to promote cell death in cancer cells [96].

ROS Signaling and Antioxidant Defense Pathways

Mitochondria are both producers and targets of reactive oxygen species, creating a signaling network that can be therapeutically modulated [96] [94]. While excessive ROS production causes oxidative damage, controlled ROS generation serves as important signaling molecules in processes including apoptosis, inflammasome activation, and calcium signaling [94]. Mitochondria-targeted antioxidants including MitoQ and MitoPBN directly quench mitochondrial ROS, protecting against oxidative damage in neurodegenerative diseases, ischemia-reperfusion injury, and diabetes [96]. These compounds typically consist of an antioxidant moiety (ubiquinone, PBN) conjugated to a lipophilic cation that facilitates mitochondrial accumulation [96].

Metabolic Integration Pathways

Mitochondria serve as central hubs for cellular metabolism, integrating carbohydrate, lipid, and amino acid metabolism through the TCA cycle and oxidative phosphorylation [93] [94]. In metabolic diseases including type 2 diabetes and obesity, mitochondrial oxidative capacity is often impaired, leading to accumulation of lipid intermediates and insulin resistance [94]. Therapeutic strategies targeting mitochondrial metabolism include:

• Uncoupling protein activators that reduce mitochondrial membrane potential, thereby decreasing ROS production while increasing energy expenditure [96].
• Substrate utilization modulators that shift fuel preference between glucose, fatty acids, and amino acids [93].
• Enzyme replacement approaches for mitochondrial metabolic disorders [96].

Figure 3: Mitochondrial Signaling Pathways as Therapeutic Targets. Key pathways including apoptotic regulation, ROS signaling, metabolic integration, and calcium homeostasis offer multiple entry points for pharmacological intervention.

Advanced Therapeutic Approaches

Mitochondrial Transplantation and Component-Based Therapies

Recent advances have demonstrated that functionally intact mitochondria can be transferred between cells, opening new therapeutic possibilities [93]. Mitochondrial transplantation has shown promise in preclinical models of various diseases, including ischemia-reperfusion injury and metabolic disorders [93]. Additionally, mitochondrial components including mtDNA, mitochondria-located microRNA, and specific proteins can function as therapeutic agents to augment mitochondrial function in immunometabolic diseases and tissue injuries [93].

Gene Therapy Approaches for Mitochondrial Diseases

Given that mitochondrial dysfunction often has genetic causes, gene therapy represents a promising approach for primary mitochondrial diseases [93]. Strategies include:

• Allotopic expression of mutant mitochondrial genes from the nuclear genome, with subsequent import of the functional protein into mitochondria [93].
• Mitochondrial genome editing using emerging technologies to correct pathogenic mutations [93].
• Modulation of mitochondrial biogenesis through targeting master regulators such as PGC-1α [93].

The strategic targeting of essential mitochondrial pathways represents a promising frontier in drug discovery, with potential applications across a broad spectrum of human diseases. Insights from parasite mitochondrial genome research reveal both conserved essential pathways and lineage-specific adaptations that can be exploited for selective therapeutic intervention. Current mitochondria-targeted compounds, including MitoQ and SkQ1, have demonstrated efficacy in preclinical models and are advancing in clinical trials, validating the overall approach [96] [94]. Future directions in mitochondrial drug discovery will likely include more sophisticated targeting strategies, personalized approaches based on individual mitochondrial phenotypes, and combination therapies that simultaneously modulate multiple mitochondrial processes. As our understanding of mitochondrial biology continues to advance, particularly through comparative genomics of diverse organisms including parasites, new therapeutic opportunities will undoubtedly emerge, making mitochondria an increasingly important target for pharmacological intervention in human disease.

Didymozoidae trematodes are significant parasitic pathogens impacting global aquaculture, particularly infecting high-value marine fishes such as tunas and groupers. The morphological similarity among didymozoid species complicates accurate identification, leading to challenges in disease diagnosis and management within aquaculture settings. This application note outlines a mitogenomic framework for precise species differentiation of these parasites. We present a validated experimental protocol for assembling complete mitochondrial genomes from parasite samples, enabling high-resolution phylogenetic analysis. This approach supports the development of targeted control strategies, thereby mitigating economic losses in aquaculture operations.

The protocol leverages third-generation sequencing technologies (PacBio HiFi) to overcome challenges associated with complex mitochondrial genome structures. By providing comprehensive genetic data, this method addresses the critical need for molecular tools in parasite identification, previously hindered by scarce genetic information for didymozoids [99] [100]. Implementation of this mitogenomic approach will enhance diagnostic accuracy, inform treatment decisions, and contribute to sustainable aquaculture health management.

The family Didymozoidae (Digenea: Hemiuroidea) comprises over 270 species of trematodes that parasitize marine fishes, especially members of the Scombridae family (e.g., tunas) [99]. These parasites form characteristic cysts or capsules in host tissues including gills, skin, and internal organs, potentially causing reduced growth, tissue damage, and increased susceptibility to secondary infections in farmed species. The taxonomic classification of didymozoids has historically relied on morphological characteristics such as body shape, reproductive organ arrangement, and cyst structure. However, these features often show considerable intraspecific variation and overlap between taxa, leading to misidentification and taxonomic confusion [99] [101].

Integrative taxonomy, which combines morphological, genetic, and ecological data, has emerged as the standard for reliable species delimitation in parasitic flatworms. Current molecular approaches for didymozoid identification primarily utilize the nuclear 28S ribosomal DNA (rDNA) and Internal Transcribed Spacer 2 (ITS2) regions [99] [100]. While these markers provide valuable phylogenetic information, they may lack sufficient resolution for distinguishing closely related species or understanding population-level dynamics. Mitochondrial genomes offer a powerful alternative with several advantages for species differentiation, including higher mutation rates, absence of introns, and haploid inheritance without recombination.

The mitogenomic approach detailed in this study addresses a significant gap in didymozoid research. Previous studies have highlighted the scarcity of genetic data for morphologically characterized didymozoids, with many clades represented by only one or two sequences in public databases [99]. This protocol enables researchers to generate complete mitochondrial genomes for didymozoid species, providing a robust foundation for phylogenetic analysis, species delimitation, and the development of specific diagnostic markers for aquaculture health monitoring.

Experimental Design and Workflow

The experimental design incorporates a holistic approach to parasite characterization, integrating both morphological and molecular methodologies to ensure comprehensive species identification. The workflow proceeds through four critical phases: (1) sample collection and morphological examination; (2) high-quality DNA extraction; (3) mitochondrial genome sequencing and assembly; and (4) comparative genomic and phylogenetic analysis.

Sample Collection and Preservation

Parasite specimens should be collected from freshly sacrificed host fishes during routine health inspections or disease outbreaks. Didymozoids typically form visible yellow capsules on host tissues, particularly under the skin, near the caudal fin, and in the branchial cavity [99]. For each parasite specimen, implement a split-sample protocol: one portion should be preserved for morphological analysis while the other is dedicated to genetic studies.

• Morphological Preservation: Fix parasites in Alcohol-Formalin-Acetic acid (AFA) solution, with or without compression, then transfer to 70-80% ethanol for long-term storage. Subsequently, stain specimens with Langeron's alcoholic-acid carmine, dehydrate through a graded ethanol series, clear in clove oil, and mount permanently in Canada balsam for microscopic examination and morphological characterization [99].
• Genetic Preservation: Preserve tissue samples (minimum 10-20 mg) in 95-100% ethanol, ensuring complete tissue immersion. Store at -20°C or -80°C until DNA extraction. This preservation method inhibits DNase activity and maintains DNA integrity for subsequent sequencing applications.

Integrative Taxonomic Approach

The integrated methodology combines morphological examination with genetic analysis to establish robust species identifications. Morphological characterization should focus on diagnostic features including body shape, sucker morphology, reproductive structures (testes and ovary arrangement), and cyst formation [99]. Genetic analysis then complements these morphological data through mitochondrial genome sequencing, enabling phylogenetic placement and validation of morphological identifications.

Detailed Protocols

High Molecular Weight DNA Extraction

High-quality DNA with minimal fragmentation is essential for successful mitochondrial genome assembly, particularly for leveraging long-read sequencing technologies.

• Reagents and Equipment:

  • TRIzol Reagent (Thermo Fisher Scientific) or CTAB extraction buffer
  • Wizard SV Genomic DNA Purification System (Promega)
  • Proteinase K (10 mg/mL)
  • DNeasy Plant MiniKit (QIAGEN)
  • NanoDrop spectrophotometer or Qubit fluorometer
  • Pulsed-field gel electrophoresis system (for quality assessment)

• Procedure:

  • Grind 20-30 mg of parasite tissue in liquid nitrogen using a sterile mortar and pestle.
  • Transfer powdered tissue to a microcentrifuge tube and add 1 mL of TRIzol reagent or CTAB buffer.
  • Incubate with 10 mg/mL Proteinase K at 56°C for 12-24 hours to ensure complete tissue lysis [99].
  • Extract genomic DNA using the Wizard SV Genomic DNA Purification System according to the manufacturer's instructions.
  • Assess DNA purity spectrophotometrically (ideal A260/A280 ratio: 1.8-2.0; A260/A230 ratio: >2.0).
  • Evaluate DNA integrity by pulsed-field gel electrophoresis or Agilent Femto Pulse System, confirming fragment sizes >50 kb for long-read sequencing [19].
  • Quantify DNA using Qubit fluorometer with dsDNA HS Assay kit.

Mitochondrial Genome Sequencing

This protocol utilizes a hybrid sequencing approach, combining long-read PacBio HiFi technology for structural resolution with Illumina short-reads for accuracy validation.

• Library Preparation and Sequencing:
  • For PacBio HiFi sequencing: Prepare 15-20 kb SMRTbell libraries using the SMRTbell Express Template Prep Kit 2.0 with 7 μg of high molecular weight DNA [51] [19].
  • Size-select libraries using the BluePippin System (15-20 kb cutoff).
  • Sequence on PacBio Revio platform with 30-hour movie times to generate HiFi reads with >99.9% accuracy.
  • For Illumina sequencing: Prepare 350-550 bp insert libraries using Nextera XT DNA Library Prep Kit.
  • Sequence on Illumina MiSeq or NovaSeq platform (2×150 bp or 2×250 bp).

Table 1: Sequencing Platform Comparison for Mitochondrial Genome Assembly

Platform	Read Length	Accuracy	Advantages	Limitations	Cost
PacBio HiFi	15-20 kb	>99.9%	Resolves complex repeats, structural variants	Higher input DNA requirement, cost	$$$
Illumina	150-300 bp	>99.9%	Cost-effective, high coverage	Short reads struggle with repeats	$
Oxford Nanopore	10 kb - 2 Mb	~95-97%	Very long reads, direct epigenetics	Higher error rate requires correction	$$

Mitochondrial Genome Assembly

The assembly process employs specialized tools designed to handle the complex repetitive structures and variable configurations characteristic of mitochondrial genomes.

• Software Requirements:

  • PMAT (Plastid and Mitochondrial Assembly Tool)
  • Flye v2.9 or newer
  • GetOrganelle v1.7.5 or newer
  • Bandage v0.8.1 or newer
  • BLAST+ v2.12 or newer
  • MITOS2 WebServer or GeSeq for annotation

• Assembly Procedure:

  • Quality Control: Filter and trim raw reads using Fastp v0.23.2 with default parameters.
  • Initial Assembly:
    • For PacBio HiFi data: Use PMAT with parameters 'autoMito -t hifi -m -T 50' for mitochondrial-specific assembly [51] [39].
    • Alternatively, use Flye with '--meta' and '--plasmids' flags for metagenomic-style assembly.
  • Graph-based Resolution: Visualize assembly graphs using Bandage to identify circular contigs and resolve complex structures [51].
  • Hybrid Assembly: For data from multiple platforms, perform hybrid assembly using Unicycler v0.5.0 with '--mode normal' for optimal integration of long and short reads.
  • Contig Validation: Validate mitochondrial contigs through BLASTn search against the NCBI nt database and check for mitochondrial gene content.
  • Circularization Verification: Confirm circular genome structure by identifying overlapping termini or through PCR validation across putative junction regions.

Table 2: Performance Comparison of Mitochondrial Genome Assembly Tools

Tool	Algorithm Type	Read Type	Strengths	Limitations
PMAT	De novo	Long reads	Specifically designed for plant mitogenomes, handles complex structures	Limited to long-read data
Flye	De novo	Long reads	Excellent for complex repeats, produces high-contiguity assemblies	May require high coverage
GetOrganelle	IME	Short/Long reads	Effective organelle genome assembly, minimizes nuclear DNA contamination	May produce fragmented assemblies for complex mitogenomes
SMARTdenovo	De novo	Long reads	Fast assembly, good contiguity	Less accurate for repetitive regions
NextDenovo	De novo	Long reads	High accuracy, good for complex genomes	Computationally intensive

Genome Annotation and Phylogenetic Analysis

Comprehensive annotation and analysis transform assembled sequences into biologically meaningful information for species differentiation.

• Genome Annotation:

  • Protein-Coding Genes: Annotate using the online platform PMGA (Plant Mitochondrial Genome Annotation) or MITOS2 WebServer with metazoan genetic code.
  • RNA Genes: Identify tRNA genes using tRNAscan-SE v2.0 with organellar genetic code settings. Locate rRNA genes through BLASTn against reference mitochondrial rRNAs.
  • Manual Curation: Verify all gene models through alignment with known didymozoid mitochondrial genes and manual inspection in software such as MacVector or Geneious.
  • Visualization: Generate circular genome maps using PMGmap or Circos v0.69-9.

• Phylogenetic Analysis:

  • Gene Selection: Extract and align conserved mitochondrial protein-coding genes (e.g., cox1, cytb, nad1, nad4, nad5) from target and reference species.
  • Sequence Alignment: Perform multiple sequence alignment using MAFFT v7.490 with L-INS-i algorithm.
  • Partitioning Scheme: Determine optimal partitioning scheme and substitution models using PartitionFinder v2.1.1 or ModelTest-NG.
  • Tree Reconstruction:
    • Maximum Likelihood: Use IQ-TREE v2.2.0 with 1000 ultrafast bootstrap replicates.
    • Bayesian Inference: Perform analysis in MrBayes v3.2.7 with two parallel runs of 10 million generations.

Expected Results and Interpretation

Successful implementation of this protocol will yield the complete mitochondrial genome of the target didymozoid parasite, typically ranging between 14-18 kb in length and encoding 37-42 genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs). The genome structure should be verified as circular, with potential observation of multiple conformations or isomeric forms due to repetitive elements.

Comparative Genomic Analysis

Comparative analysis of the assembled mitogenome against reference sequences enables identification of species-specific genetic markers and phylogenetic placement:

• Sequence Divergence: Calculate pairwise genetic distances (p-distance or Kimura 2-parameter) for protein-coding genes, particularly cox1, to establish divergence thresholds for species delimitation. Cox1 typically shows 10-15% divergence between congeneric didymozoid species.
• Synonymous vs. Non-synonymous Substitutions: Determine Ka/Ks ratios for protein-coding genes to identify evolutionary pressures. Most mitochondrial genes should show Ka/Ks < 1, indicating purifying selection and functional constraint.
• Gene Order Conservation: Compare gene arrangement with related species. While mitochondrial gene order is generally conserved within families, rearrangements can provide phylogenetic signals for deep divergences.
• Repetitive Element Analysis: Identify simple sequence repeats (SSRs) and dispersed repeats that may serve as potential molecular markers for population-level studies.

Table 3: Key Mitochondrial Genes for Didymozoid Phylogenetics and Diagnostics

Genetic Marker	Sequence Length	Evolutionary Rate	Application	Technical Considerations
cox1	~720 bp	Medium	Species barcoding, population genetics	Universal primers available
cytb	~1,100 bp	Medium	Species identification, phylogenetics	Informative at species level
nad1	~900 bp	Medium-Fast	Population studies, species delimitation	More variable than cox1
rrnL	~950 bp	Slow	Deep phylogeny, family-level relationships	Secondary structure important
Complete Mitogenome	14-18 kb	Variable	Comprehensive phylogeny, genome evolution	Requires advanced bioinformatics

Phylogenetic Resolution

The phylogenetic analysis should robustly resolve the taxonomic position of the target didymozoid with strong statistical support (bootstrap values >90%, posterior probabilities >0.95). The resulting phylogeny will:

• Establish monophyletic clustering with congeners or morphologically similar species
• Confirm or refute species hypotheses based on morphological data
• Reveal potential cryptic species complexes through distinct, well-supported clades
• Provide insights into host-parasite co-evolutionary patterns

Previous studies using mitochondrial markers have successfully resolved phylogenetic relationships within the Didymozoidae, such as distinguishing between Platodidymocystis yamagutii n. gen., n. sp. and related genera like Platocystis and Didymocystis [99] [100]. The complete mitogenome approach provides substantially greater phylogenetic signal through concatenated analysis of all protein-coding genes.

The Scientist's Toolkit

Table 4: Essential Research Reagents and Tools for Didymozoid Mitogenomics

Item	Specification	Application	Notes
DNA Extraction	Wizard SV Genomic DNA Purification System	High molecular weight DNA isolation	Preserve high molecular weight
Long-read Sequencing	PacBio SMRTbell Express Template Prep Kit 3.0	Library prep for HiFi sequencing	Enables 15-20 kb inserts
Short-read Sequencing	Illumina Nextera XT DNA Library Prep Kit	Complementary short-read data	350-550 bp inserts
Assembly Software	PMAT v1.5 or newer	Mitochondrial-specific assembly	Optimized for organelle genomes
Assembly Visualization	Bandage v0.8.1+	Assembly graph inspection	Critical for complex structures
Gene Annotation	MITOS2 WebServer	Automated gene finding	Uses metazoan genetic code
Sequence Alignment	MAFFT v7.490	Multiple sequence alignment	L-INS-i algorithm recommended
Phylogenetic Analysis	IQ-TREE v2.2.0+	Maximum likelihood tree building	Ultrafast bootstrap
Morphological Analysis	Alum carmine stain	Tissue staining for morphology	Follow standardized protocols

Applications in Aquaculture Health

Implementation of this mitogenomic protocol directly benefits aquaculture health management through:

• Accurate Parasite Identification: Enable precise species-level identification of didymozoid infestations in farmed fish populations, facilitating targeted treatment strategies.
• Disease Surveillance: Develop PCR-based diagnostic assays targeting species-specific mitochondrial markers for routine monitoring of parasite prevalence in aquaculture facilities.
• Treatment Efficacy Assessment: Monitor genetic diversity and potential emergence of treatment-resistant parasite populations through mitochondrial marker analysis.
• Biosecurity Protocols: Identify parasite sources and transmission pathways to implement effective biosecurity measures and prevent introduction into aquaculture systems.
• Epidemiological Tracking: Use mitochondrial sequence variability to trace outbreak origins and understand parasite dispersal patterns within and between aquaculture sites.

The mitogenomic approach outlined here represents a significant advancement over traditional morphological identification, providing aquaculture professionals with powerful molecular tools for proactive parasite management. This methodology supports the development of evidence-based health management strategies, ultimately contributing to reduced economic losses and improved sustainability of aquaculture operations.

Conclusion

Mitochondrial genome assembly has emerged as an indispensable tool for precise parasite taxonomy, resolving species complexes that morphological alone cannot distinguish. The integration of advanced sequencing technologies, robust bioinformatic pipelines, and comparative genomic frameworks provides unprecedented resolution for phylogenetic studies and population genetics. Future directions will leverage these detailed mitochondrial blueprints to identify essential, parasite-specific metabolic pathways, directly fueling target-based drug design and the development of novel therapeutics against neglected parasitic diseases. This field stands to make significant contributions to both evolutionary biology and translational clinical research.

References

• 1. Divergence of the mitochondrial structure in the... genome [https://pubmed.ncbi.nlm.nih.gov/20034997/]
• 2. Novel type of linear with dual flip-flop... mitochondrial genomes [https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-13-622]
• 3. Insights into the mitochondrial genome structure and ... [https://www.nature.com/articles/s41598-025-92939-0]
• 4. -encoded mRNAs are a distinct subpopulation... Circular mitochondrial [https://www.nature.com/articles/s41598-023-34255-z?error=cookies_not_supported&code=a8062fd2-7403-454c-9e95-2f48f887da9b]
• 5. 'Why genomes in pieces?' revisited: Sucking lice do their own ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3647531/]
• 6. Application of mitochondrial genomes to species ... [https://www.sciencedirect.com/science/article/pii/S0717345825000193]
• 7. A new long-read mitochondrial-genome protocol (PacBio HiFi ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11069185/]
• 8. Comprehensive analysis of the complete mitochondrial ... [https://ija.scholasticahq.com/article/142194]
• 9. Mitochondrial genome diversity across the subphylum ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10701893/]
• 10. Assessing de novo parasite genomes assembled using ... [https://www.sciencedirect.com/science/article/pii/S258900422401839X]
• 11. Mitochondrial genome diversity in parasites [https://www.sciencedirect.com/science/article/abs/pii/S0020751999001903]
• 12. Mitochondrial DNA [https://en.wikipedia.org/wiki/Mitochondrial_DNA]
• 13. The complete mitochondrial genomes of three parasitic ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3693896/]
• 14. Characterization of fragmented mitochondrial ribosomal ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC140509/]
• 15. Improved annotation of protein-coding genes boundaries ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6847864/]
• 16. The complete mitochondrial genome of a parasite at the ... [https://parasitesandvectors.biomedcentral.com/articles/10.1186/s13071-020-3926-5]
• 17. Comparative analysis of the mitochondrial genomes of four ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2024.1429545/full]
• 18. Assembly and comparative analysis of the complete ... [https://www.nature.com/articles/s41598-025-93762-3]
• 19. Mitochondrial Genome Assembly and Comparative ... [https://www.mdpi.com/2311-7524/11/10/1231]
• 20. Insights into the mitochondrial genome structure and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11950482/]
• 21. Red Algae Lose Key Mitochondrial Genes in Response to ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3014286/]
• 22. Transposable element abundance correlates with mode of ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7313128/]
• 23. Extensive intraspecific gene order and gene structural ... [https://www.nature.com/articles/s41467-019-10820-x]
• 24. Deciphering the dynamics of active autonomous terminal ... [https://www.sciencedirect.com/science/article/pii/S2662173825002164]
• 25. Evolutionary genetics of the mitochondrial genome: insights ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10324950/]
• 26. The Integrative Taxonomy and Mitochondrial Genome ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12192195/]
• 27. Extensive analysis of mitochondrial DNA quantity and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8730313/]
• 28. Sequence analyses at mitochondrial and nuclear loci reveal a ... [https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0225822]
• 29. Making the best of a bad sample: Comparison of DNA ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12121741/]
• 30. Complete mitochondrial genome and phylogenetic ... [https://www.sciencedirect.com/science/article/pii/S221322442500094X]
• 31. Comparative analysis of illumina and oxford nanopore ... [https://www.nature.com/articles/s41598-025-18768-3]
• 32. Nanopore sequencing accuracy [https://nanoporetech.com/platform/accuracy]
• 33. Sequencing Quality Scores [https://www.illumina.com/science/technology/next-generation-sequencing/plan-experiments/quality-scores.html]
• 34. What does accuracy mean to you? [https://nanoporetech.com/blog/what-does-accuracy-mean-to-you-]
• 35. Comparison of Illumina and Oxford Nanopore sequencing ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-025-11267-9]
• 36. Ultima Genomics UG 100 platform vs. Illumina NovaSeq X ... [https://www.illumina.com/science/genomics-research/articles/ultima-ug-100-vs-illumina-novaseq-x-wgs-data.html]
• 37. De Novo assembly and characterization of Aria alnifolia ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-025-11892-4]
• 38. A survey of k-mer methods and applications in bioinformatics [https://www.sciencedirect.com/science/article/pii/S2001037024001703]
• 39. Advance in the assembly of the plant mitochondrial ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12576441/]
• 40. Norgal: extraction and de novo assembly of mitochondrial ... [https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-017-1927-y]
• 41. De-novo assembly and comparative analysis of the complete ... [https://bmcplantbiol.biomedcentral.com/articles/10.1186/s12870-025-06731-3]
• 42. De Novo Assembly and Comparative Analysis ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11171644/]
• 43. Optimizing the procedure of mitochondrial genome assembly ... [https://cepdnaclk.github.io/e15-4yp-Optimizing-Mitochondria-Genome-Assembly-And-Annotation-With-Skim-Sequencing-Data/]
• 44. De novo assembly of the complete mitochondrial genomes ... [https://www.nature.com/articles/s41598-025-86411-2]
• 45. Mitochondrial genome annotation with MFannot: a critical ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10352661/]
• 46. MITOS: Improved de novo metazoan mitochondrial ... [https://www.sciencedirect.com/science/article/abs/pii/S1055790312003326]
• 47. BLAST: Basic Local Alignment Search Tool [https://blast.ncbi.nlm.nih.gov/Blast.cgi]
• 48. BLAST Searching [https://www.geneious.com/tutorials/blast-searching]
• 49. Miniaturized mitogenome of the parasitic plant Viscum ... - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC4500244/]
• 50. Assembly and comparative analysis of the mitochondrial ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2025.1704888/full]
• 51. Assembly and comparative analysis of the mitochondrial ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12597970/]
• 52. Characterization of co-infections of haemosporidian ... [https://www.sciencedirect.com/science/article/pii/S2667114X25000731]
• 53. Structural variations and phylogenetic implications of ... [https://www.sciencedirect.com/science/article/pii/S0926669025013639]
• 54. Mitochondrial genome assembly and comparative analysis ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2025.1556379/full]
• 55. The evolutionary dynamics of organellar pan-genomes in ... [https://genomebiology.biomedcentral.com/articles/10.1186/s13059-025-03717-0]
• 56. The Mighty NUMT: Mitochondrial DNA Flexing Its Code in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10216076/]
• 57. Preventing the pollution of mitochondrial datasets with ... [https://pubmed.ncbi.nlm.nih.gov/21047564/]
• 58. Preventing the pollution of mitochondrial datasets with ... [https://www.sciencedirect.com/science/article/abs/pii/S1567724910001881]
• 59. Mitochondrial genomic analyses provide new insights into ... [https://pubmed.ncbi.nlm.nih.gov/36324074/]
• 60. Mitochondrial genomic analyses provide new insights into the ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9628169/]
• 61. Characterization of the complete mitochondrial genome of ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12644934/]
• 62. Are we there yet? Benchmarking low-coverage nanopore long ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-022-08482-z]
• 63. PMAT: an efficient plant mitogenome assembly toolkit ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10925850/]
• 64. Long read mitochondrial genome sequencing using Cas9 ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9399971/]
• 65. Long read mitochondrial genome sequencing using Cas9 ... [https://www.sciencedirect.com/science/article/abs/pii/S1567724922000514]
• 66. A method for multiplexed full-length single-molecule ... [https://www.nature.com/articles/s41467-022-33530-3]
• 67. Novel method for multiplexed, full-length single-molecule ... [https://nanoporetech.com/resource-centre/lc22-novel-method-for-multiplexed-full-length-single-molecule-sequencing-of-human-mitochondrial-dna-using-cas9-mediated-enrichment]
• 68. Aln2tbl: building a mitochondrial features table from ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8386716/]
• 69. Chromosome-level genome assembly and annotation of ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12639136/]
• 70. Data Preparation and Quality Assessment in Genome Assembly [https://scinet.usda.gov/events/2025-04-29-data]
• 71. An annotated near-complete sequence assembly of the ... [https://www.nature.com/articles/s41597-025-05116-3]
• 72. Research Article Benchmarking long-read assembly tools ... [https://www.sciencedirect.com/science/article/pii/S2215017X2500058X]
• 73. High-Throughput Sequencing Enables Rapid Analyses of ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11944570/]
• 74. The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny | IMA Fungus | Full Text [https://imafungus.biomedcentral.com/articles/10.1186/s43008-024-00149-6]
• 75. Mitochondrial cob and cox1 genes and editing of the corresponding mRNAs in Dinophysis acuminata from Narragansett Bay, with special reference to the phylogenetic position of the genus Dinophysis - PubMed [https://pubmed.ncbi.nlm.nih.gov/18165361/]
• 76. Coloring the edges of a plotted tree by the marginal ... [https://blog.phytools.org/2023/04/coloring-edges-of-plotted-tree-by_18.html]
• 77. Specifying colours in phylogenetic ancestral state ... [https://stackoverflow.com/questions/44047964/specifying-colours-in-phylogenetic-ancestral-state-reconstruction]
• 78. 2.4. Mitogenomic of Six Ophiocordyceps Species Comparison [https://bio-protocol.org/exchange/minidetail?id=14847798&type=30]
• 79. The mitochondrial genomes of Tortricidae: nucleotide ... composition [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-021-08041-y]
• 80. The complete mitochondrial genome of parasitic nematode ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-017-4237-x]
• 81. Comparative Mitogenomics Provides Valuable Insights for ... [https://www.mdpi.com/2309-608X/10/11/769]
• 82. Annotated Bioinformatic Pipelines for Genome and... Assembly [https://bio-protocol.org/en/bpdetail?id=5231&type=0]
• 83. Genome Microsynteny Illuminates the Fast Evolution of... Comparative [https://pubmed.ncbi.nlm.nih.gov/38124445/]
• 84. A new long-read mitochondrial - genome (PacBio HiFi) for... protocol [https://malariajournal.biomedcentral.com/articles/10.1186/s12936-024-04961-8]
• 85. Frontiers | Morphology, complete mitochondrial and... genome [https://www.frontiersin.org/journals/veterinary-science/articles/10.3389/fvets.2025.1519092/full]
• 86. Mitochondrial-related drug resistance lncRNAs as ... [https://link.springer.com/article/10.1007/s12672-024-01690-x]
• 87. The molecular mechanisms of mitochondrial dynamics and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12492677/]
• 88. GISAID - gisaid.org [https://gisaid.org/]
• 89. Mitochondrial segmentation and function prediction in live ... [https://www.nature.com/articles/s41467-025-55825-x]
• 90. Extended Analysis Tools [https://v3.boldsystems.org/index.php/resources/handbook?chapter=5_analysis_extended.html]
• 91. Barcode ID [https://boldsystems.org/barcode-id/]
• 92. Best Image Award 2025 Winners Announced [https://wms-site.com/alert-on-mitochondria/1352-best-image-award-2025-winners-announced]
• 93. Mitochondrial dysfunction: mechanisms and advances in ... [https://www.nature.com/articles/s41392-024-01839-8]
• 94. Mitochondria-Targeted Drugs - PMC - PubMed Central [https://pmc.ncbi.nlm.nih.gov/articles/PMC6875871/]
• 95. Mitochondrial Targets for Pharmacological Intervention in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC4286170/]
• 96. Pharmacological targeting of mitochondria in disease - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC2189819/]
• 97. Characterization of co-infections of haemosporidian parasites ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12419108/]
• 98. Next-generation sequencing, assembly, annotation, and ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-025-12171-y]
• 99. Integrative Taxonomy of Didymozoids Parasitizing ... [https://www.mdpi.com/2076-0817/14/4/359]
• 100. Integrative Taxonomy of Didymozoids Parasitizing ... [https://pubmed.ncbi.nlm.nih.gov/40333176/]
• 101. Detection and molecular analysis of Philopinna higai ... [https://pubmed.ncbi.nlm.nih.gov/37562581/]