Skip to main content

Mitochondrial metabolism of the facultative parasite Chilodonella uncinata (Alveolata, Ciliophora)

Abstract

Background

Chilodonella uncinata is an aerobic ciliate capable of switching between being free-living and parasitic on fish fins and gills, causing tissue damage and host mortality. It is widely used as a model organism for genetic studies, but its mitochondrial metabolism has never been studied. Therefore, we aimed to describe the morphological features and metabolic characteristics of its mitochondria.

Methods

Fluorescence staining and transmission electron microscopy (TEM) were used to observe the morphology of mitochondria. Single-cell transcriptome data of C. uncinata were annotated by the Clusters of Orthologous Genes (COG) database. Meanwhile, the metabolic pathways were constructed based on the transcriptomes. The phylogenetic analysis was also made based on the sequenced cytochrome c oxidase subunit 1 (COX1) gene.

Results

Mitochondria were stained red using Mito-tracker Red staining and were stained slightly blue by DAPI dye. The cristae and double membrane structures of the mitochondria were observed by TEM. Besides, many lipid droplets were evenly distributed around the macronucleus. A total of 2594 unigenes were assigned to 23 functional classifications of COG. Mitochondrial metabolic pathways were depicted. The mitochondria contained enzymes for the complete tricarboxylic acid (TCA) cycle, fatty acid metabolism, amino acid metabolism, and cytochrome-based electron transport chain (ETC), but only partial enzymes involved in the iron-sulfur clusters (ISCs).

Conclusions

Our results showed that C. uncinata possess typical mitochondria. Stored lipid droplets inside mitochondria may be the energy storage of C. uncinata that helps its transmission from a free-living to a parasitic lifestyle. These findings also have improved our knowledge of the mitochondrial metabolism of C. uncinata and increased the volume of molecular data for future studies of this facultative parasite.

Graphical Abstract

Background

Mitochondria are thought to have arisen as a consequence of endosymbiosis of an ancestral α-proteobacterium into an archaeon [1, 2]. Mitochondria can be categorized into five classes based on energy metabolism: aerobic mitochondria, anaerobic mitochondria, H2-producing mitochondria, hydrogenosomes, and mitosomes [3]. For example, anaerobic flagellates such as trichomonads have hydrogenosomes that lack an cytochrome-based electron transport chain (ETC) and can produce H2 [4]. The armophorean ciliate like Nyctotherus ovalis has a hydrogen-producing mitochondrion [5]. The amoebozoan Entamoeba istolytica contains mitosomes, which have lost the ability to produce ATP [6]. The free-living oligohymenophorean aerobic ciliates belonging to the genera Tetrahymena and Paramecium have more typical aerobic mitochondria [7]. These mitochondria-related organelles (MROs) have evolved several times along with the adaptation of ciliates to the anoxic and oxygen-depleted habitats [8]. The MROs of different anaerobic ciliate lineages (class Armophorea, class Plagiopylea, class Muranotrichea, etc.) have common features and possess unique characteristic as well [8,9,10]. Compared with these anaerobic ciliates, the understanding of mitochondrial metabolism in aerobic parasitic ciliates is still not enough.

The aerobic ciliate Chilodonella uncinata is considered to be a facultative parasite, which can be infective to goldfish and cause tissue damage and even host mortality [11, 12]. Other Chilodonella species are either free-living or obligate parasitic [13, 14]. The ability of C. uncinata to adapt to both free-living and parasitic environment indicates that it can use different nutrients to produce energy. Mitochondria play a central role in ATP production via oxidative phosphorylation, so it is necessary to study the mitochondria of C. uncinata to explore their metabolic characteristics and lifestyle adaptability. Our previous work found that there was no difference in the mitochondrial proteins between the parasitic C. uncinata and the free-living C. uncinata [15]. However, the metabolic characteristics of C. uncinata's mitochondria are still unclear.

In the present study, multiple methods were combined to study the mitochondria of C. uncinata. First, fluorescence staining was used to observe the mitochondria in vivo. Then, the morphology of mitochondria was revealed by transmission electron microscopy (TEM). Single-cell transcriptome data of C. uncinata were used to identify the metabolic pathways and make the Clusters of Orthologous Groups (COG) annotation. We also sequenced the cytochrome c oxidase subunit 1 (COX1) gene and performed the phylogenetic analysis.

Methods

Specimen collection

Specimens used in this study were all collected in our laboratory since an in vitro culture and a successful infection model of C. uncinata had been formerly established [15]. Since the expression of mitochondrial proteins did not significantly differ between the parasitic and free-living C. uncinata [15], we used free-living specimens (from the culture medium) for Mito-tracker and DAPI staining, TEM, and DNA extraction for the convenience of sampling. Both free-living and parasitic C. uncinata were collected for the transcriptome sequencing.

Mito-tracker and DAPI staining

Mito-tracker and DAPI stainings were carried out to locate and visualize DNA and mitochondria. First, C. uncinata cells were washed in sterile water, and then the Mito-tracker Red Dye was dissolved in DMSO for a stock solution. Finally, the stock solution was diluted with PBS (phosphate-buffered saline) for a working solution. Chilodonella uncinata specimens were subsequently stained with a few drops of the working solution for 30 min at room temperature. After the Mito-tracker Red staining, 0.02% DAPI (4′,6-diamidino-2-phenyl-indole) was added to the cell suspension for 10 min. The cells were visually inspected and photographed using Axioplan 2 imaging and Axiophot 2 (Zeiss, Oberkochen, Germany).

Transmission electron microscopy (TEM)

For the transmission electron microscopy, C. uncinata cells were fixed in 2.5% glutaraldehyde. Then, the samples were postfixed in 1% (v/v) osmium tetroxide in PBS for 2 h at 4 ℃ followed by dehydration in a gradient acetone series and embedded in Araldite. Ultrathin sections were cut on a Leica Ultracut R ultramicrotome (Leica, Germany) and stained with uranyl acetate and lead citrate before being observed using a JEM-1230 Transmission Electron Microscope (JEOL, Japan).

Transcriptome analysis and mitochondrial protein prediction

Transcriptome analysis was conducted based on our former sequencing and assembly results, with reads deposited in GenBank under the BioSample number PRJNA842862 [15]. Transdecoder v5.5.0 (http:transdecoder.github.io) was used to predict open-reading frames (ORFs). After prediction, genes were annotated using the COG database. Putative mitochondrial proteins were detected by BLASTP [16] using homologs from other ciliates and well-described mitochondrial proteomes as queries: Homo sapiens [17], Saccharomyces cerevisiae [18], and Tetrahymena thermophila [19]. TargetP [20] and MitoFates [21] were also used to predict the mitochondrial-targeting signals. BlastKOALA (www.kegg.jp/blastkoala/) server [22] was then used to annotate these putative mitochondrial proteins. Mitochondrial metabolic pathways were visualized with Adobe illustrator. All other visualizations were completed in R v4.1.1 [23].

DNA extraction, amplification, and sequencing

About 500 C. uncinata were collected from the culture medium, suspended in lysis buffer (10 mM Tris–HCl, pH 8.0; 1 M EDTA, pH 8.0; 0.5% sodium dodecyl sulfate; 60 μg/ml proteinase K), and incubated at 55 °C for 12–20 h. DNA was extracted using the REDExtract-N-Amp Tissue PCR Kit (Sigma, St. Louis, MO, USA). The COX1 gene was amplified with two primers listed in Additional file 1: Table S1. Polymerase chain reaction (PCR) was conducted in a 20-µl reaction mixture containing 7.4 µl ddH2O, 10 µl 2 × PCR buffer (Mg2+, dNTP plus, Takara, Dalian, China), 0.6 µl of each primer, 0.4 µl recombinant Taq DNA polymerase (250 U/µl, Takara, Dalian, China), and 1 µl DNA template. PCR conditions were as follows: initial denaturation at 98 °C for 2 min, followed by 40 cycles at 98 °C for 10 s, 50 °C for 15 s, 68 °C for 1 min/kb, and a final extension at 68 °C for 10 min. PCR amplified products were sequenced on an ABI PRISM® 3730 DNA Sequencer (Applied Biosystems, USA).

Sequence comparisons and phylogenetic analysis

BLAST [16] was used to search the homology from the NCBI database. Amino acid sequences were deduced using ExPASy translate tool (http://web.expasy.org/translate/) [24]. The secondary structure of the COX1 in C. uncinata was determined by PSIPRED Protein Sequence Analysis Workbench (http://bioinf.cs.ucl.ac.uk/psipred) [25] to predict coil, strand, and helix chains of the COX1 protein. The transmembrane segment prediction was conducted and visualized by the EMBOSS packages tmap and pepwheel [26]. The 3D structure of COX1 protein was predicted by Swiss-Model online server (http://swissmodel.expasy.org/interactive) [27]. Swiss-PDB viewer [28] was used to check the stereochemical quality of the structure.

Aside from the COX1 gene sequence of C. uncinata obtained in this study, other sequences of COX1 were retrieved from the GenBank database using PhyloSuite v1.1.15 [29]. Sequences of the class Heterotrichea were selected as outgroups. Sequence alignments, best-fit model selection, and phylogenetic tree reconstruction were conducted using plug-in programs implemented in PhyloSuite: MAFFT v7.313 [30], ModelFinder [31], MrBayes v3.1.2 [32], and IQ-TREE v1.6.8 [33]. The GTR + I + F + G4 model was selected as the best model, which was used for both maximum likelihood (ML) and Bayesian inference (BI) methods. The ML analysis was performed by the IQ-TREE v1.6.8 [33] with 1000 bootstrap replicates. BI analysis was performed using MrBayes 3.2.6 [32] with 1,000,000 generations. Trees were sampled every 1000 generations with the initial 25% of sampled data discarded as burn-in. The ML and BI trees were visualized and edited using the Interactive Tree Of Life (iTOL v4) tool [34], with annotation files generated by PhyloSuite.

Results

Mito-tracker and DAPI staining

Different staining methods were used to observe mitochondria of C. uncinata in vivo. Living specimens were photographed from the dorsal side (Fig. 1A) and ventral side (Fig. 1B right) to show the gross morphology of C. uncinata. Cilia, cytostome, oral basket, macronucleus, and micronucleus are clearly visible. The macronucleus and micronucleus were stained blue by DAPI dye. Mitochondria were also slightly blue (Fig. 1C), which indicates that they contain traces of DNA. Mitochondria were stained red using the Mito-tracker Red staining, which revealed that they are evenly distributed within the cell (Fig. 1D).

Fig. 1
figure 1

Overall morphology of C. uncinata and Mito-tracker Red and DAPI staining of its mitochondria. A Differential interference contrast (DIC) image of a living specimen from the dorsal side. B The macronucleus and micronucleus of a living specimen from the dorsal side (left) and cytostome (arrow) from the ventral side (right). C DAPI staining, showing the macronucleus and micronucleus of C. uncinata; arrow shows mitochondria. D The mitochondria of C. uncinata stained red with Mito-tracker Red

Transmission electron microscopy (TEM)

The macronucleus, mitochondria, and other organelles are shown in Fig. 2A. There was a circle of lipid droplets around the macronucleus, which were round in outline and had low density areas (Fig. 2B). Mitochondria were unevenly distributed among lipid droplets and their morphologies were slightly different, ranging from oval to round (Fig. 2B). Some free ribosomes were present around the mitochondria (Fig. 2C). Further magnifications allowed us to observe mitochondrial cristae and the double membrane (Fig. 2D).

Fig. 2
figure 2

Transmission electron microscope pictures of C. uncinata. A The whole section of C. uncinata, where N is the nucleolus of the macronucleus, M is a mitochondrion, and the star sign denotes lipid droplets (Ldp). B Lipid droplets and mitochondria of various shapes are evenly distributed around the macronucleus. C A detail showing mitochondria and ribosomes (white arrow) around them. D Cristae and the double membrane of mitochondria

COG annotation

A total of 2594 unigenes were assigned to 23 functional classifications of COG (Fig. 3). Among them, 680 genes were assigned to the ‘information storage and processing’ category; 555 genes were assigned to the ‘cellular processes and signaling’ category; 1075 genes were assigned to the ‘metabolism’ category; 284 genes were assigned to the ‘poorly characterized’ category. Within the information storage and processing category, translation, ribosomal structure, and biogenesis were the most abundant sub-categories that accounted for 63.53% of unigenes (n = 432) (Fig. 3A). Within the cellular processes and signaling category, the most abundant sub-category was represented by posttranslational modification, protein turnover, and chaperones, which accounted for 38.20% of unigenes (n = 212). There were no assigned genes in extracellular structures and nuclear structure categories (Fig. 3B). Within the metabolism category, the top three identified functional sub-categories were energy production and conversion, amino acid transport and metabolism, and lipid transport and metabolism, which accounted for 58.60% (n = 630) of unigenes (Fig. 3C). As for the category named simply ‘poorly characterized,’ one of the most abundant sub-categories embraced was the general function prediction category, which accounted for 68.66% (n = 195) of unigenes (Fig. 3D). All these data are summarized in Additional file 1: Table S2.

Fig. 3
figure 3

Clusters of orthologous groups (COG) annotation of unigenes. A Information storage and processing; B cellular processes and signaling; C metabolism; D poorly characterized

Fig. 4
figure 4

Predicted mitochondrial metabolic pathway map of C. uncinata based on the molecular data. Circles in different colors represent different biological functions. The subunits that were not identified are outlined by a solid line. The cycles comprise: glycolysis (purple); fatty acid metabolism (grey); mitochondrial pyruvate carrier (pink); protein import and folding (brown); iron-sulfur cluster (ISC) assembly system (green); tricarboxylic acid (TCA) cycle (yellow); amino acid metabolism (light blue); superoxide catalytic reaction (plum-red); ETC complex (dark blue). For abbreviations, see Additional file 1: Table S3

Mitochondrial metabolic pathways

A total of 111 mitochondrial proteins were predicted based on homology to known mitochondrial proteins (Fig. 4). Among them, the number of proteins with TargetP prediction scores and Mitoprot prediction scores > 0.5 were 42 and 50, respectively. The scores represent the probability that the protein is located in mitochondria based on mitochondrial target signals (MTSs, see Additional file 1: Table S3).

Glycolysis-derived pyruvate is imported into mitochondria via the mitochondrial pyruvate carrier (MPC). Then, pyruvate is oxidized into acetyl coenzyme A (Ace-CoA) via the pyruvate dehydrogenase complex (PDH). Ace-CoA subsequently participates in the typical tricarboxylic acid (TCA) cycle. The TCA cycle includes many different enzymes, such as citrate synthase (CS), aconitate hydratase (ACO), isocitrate dehydrogenase (IDH), 2-oxoglutarate dehydrogenase (OGDH), succinyl-CoA synthetase (SCS), succinate dehydrogenase (SDH), fumarate hydratase (FH), and malate dehydrogenase (MDH). Ace-CoA can also be formed by fatty acid oxidation. First, fatty acids are activated by the acyl-CoA synthetase long-chain family (ACSL) and then transported by the solute carrier family (SLC25) and oxidized by carnitine palmitoyl transferase (CPT). Besides, glutamate is generated from different amino acids via transaminases, such as alanine aminotransferase (ALT), aspartate aminotransferase (AST), and branched-chain amino acid aminotransferase (BCAAT). Then, glutamate is transformed to α-ketoglutarate via glutamate dehydrogenase (GDH), which can join the TCA cycle when necessary.

The mitochondrial oxidative phosphorylation (OXPHOS) is completed with the help of mitochondrial ETC. A pair of electrons from NADH are donated to the co-enzyme-ubiquinone (CoQ) through the complex I. Electrons from FADH2 are also donated to the CoQ by the complex II. Then, these electrons are donated to cytochrome c by complex III, Q-cytochrome c oxidoreductase (also called the cytochrome reductase). Finally, cytochrome c transports electrons to O2 via the complex IV cytochrome oxidase and molecular oxygen acts as the terminal electron acceptor. In this study, 11 and 5 subunits of complex I and complex IV were identified, respectively. The subunits of complex II including the succinate dehydrogenase flavoprotein subunit (SDHA) and succinate dehydrogenase iron-sulfur subunit (SDHB) were identified. The subunits of complex III include cytochrome b 5 (CYB5) and cytochrome c (CYC). ATP1, ATP2, ATP3, ATP5, ATP9, ATP16, and ATP20 of complex V, the ATP synthase, also have been identified.

Ferredoxin (YAH1), monothiol glutaredoxin (GRX5), IBA57, and NFU1 involved in the mitochondrial iron-sulfur cluster (ISC) assembly system were also identified. The translocator of the outer membrane (TOM) and translocator of the inner membrane (TIM), which are responsible for the entry and exit of proteins, were detected as well.

COX1 structure prediction

The structural prediction of COX1 in C. uncinata indicated that it has 11 transmembrane segments (Fig. 5A). Among them, segments c and g have many hydrophilic polar uncharged amino acids, while the other nine segments had more nonpolar amino acids, which indicates that they are very hydrophobic. Overall, hydrophobic amino acids accounted for a large proportion of total amino acids. The predicted 3D structure of COX1 is shown in Fig. 5B. The best template identified was the COX1 of Tetrahymena thermophila (7w5z.7.A) with a coverage of 97% and identity of 58.82%. The Global Model Quality Estimate (GMQE) score was 0.77, and the QMEAN score was 0.68. The reliability ranges of GMQE and QMEAN are both from 0 to 1, where the higher values indicate better model quality [27].

Fig. 5
figure 5

The secondary and 3D structures of COX1 of C. uncinata A There are 11 transmembrane segments in the secondary structure. Purple and blue letters represent nonpolar amino acids, while blue letters denote more hydrophobic amino acids than those in purple letters. Black letters indicate polar charged amino acids, and red letters indicate polar uncharged amino acids. B The 3D structure of COX1 predicted by Swiss-Model

Phylogenetic analysis

The COX1 gene sequence obtained in this study has been deposited in the NCBI GenBank database under the accession number OP672162. The complete COX1 gene sequence comprises 1992 nucleotides and has a GC content of 24.60%. Since there is currently (September 27, 2022) no COX1 gene sequence available for the family Chilodonellidae, the most similar sequence to C. uncinata was the orthologue of Tetrahymena nipissingi under the accession number EF070295, with 52% coverage and 76.67% identity.

The topologies generated using BI and ML algorithms were concordant, so a single consensus tree is presented herein (Fig. 6). There were no COX1 gene sequences available online of the order Chlamydodontida, so no species of this class has clustered with C. uncinata. It only clustered with the order Colpodida (represented by Maryna umbrellata and Colpoda sp.).

Fig. 6
figure 6

Phylogenetic trees (BI/ML) inferred from COX1 gene sequences. The newly obtained sequence is marked in bold. Numbers at nodes represent the Bayesian posterior probabilities (left) and bootstrap values (right)

Discussion

The mitochondria of parasitic protists evolved different shapes and sizes. The mitochondrion-related organelles (MROs) of Blastocystis hominis are predominantly elongate and range approximately from 1.2 to 3.0 μm in length and from 0.7 to 1.0 μm in width [35]. The hydrogenosomes of Tritrichomonas foetus are scattered around, and most are spherical [4]. TEM images of the marine fish parasite Cryptocaryon irritans showed that it had more typical mitochondria, which possess cristae structures [36]. In this study, we found that the mitochondria of C. uncinata are comparatively similar to those of C. irritans as both exhibit the cristae structures. Oxygen availability is considered the major driving force for the diversity of mitochondria [37]. Both C. uncinata and C. irritans parasitize the gills and fins of fish, where they have a sufficient supply of oxygen, while B. hominis and T. foetus parasitize the anaerobic intestinal and urogenital tracts of humans and cattle, respectively.

In this study, DAPI and Mito-tracker Red stainings were used to further confirm the existence and structure of C. uncinata’s mitochondria and the associated genome. The information about MROs of some common parasitic protozoa is shown in Table 1. It is worth noting that not all MROs have their own genome. Hydrogenosomes normally lack a genome, and they are entirely dependent on the nucleus for their genetic livelihood [38]. However, it has been found that the MROs of N. ovalis have a hydrogen-producing mitochondrial genome [5]. The flagellate Giardia lamblia (syn. G. intestinalis, G. duodenalis) living in the human small intestine has highly reduced mitochondria, called mitosomes [39]. Although we provide here the evidence of the presence of the mitochondrial genome in C. uncinata in this study, the complete mitochondrial genome of C. uncinata needs to be sequenced in the future.

Table 1 Mitochondria-related organelles (MROs) of parasitic protists

With the specialization of parasite mitochondria, their metabolic pathways also changed. The mitochondrial metabolic pathways of C. uncinata were inferred based on single-cell RNA-seq data. The enzymes ACSL and SLC25, which are related to the fatty acid metabolism pathway, were detected on the outer and inner membrane of C. uncinata’s mitochondria, respectively. ACSL is crucial for fatty acid uptake and activation [40]. The member A20 of the SLC25 protein family (SLC25A20) helps with the carnitine acyl-carnitine carriage, which plays an important part in the mitochondrial β-oxidation pathway [41]. The Ace-CoA produced by β-oxidation is then added to the TCA cycle to generate energy in C. uncinata. Mitochondria associated with lipid droplets have increased fat oxidation, pyruvate oxidation, and electron transport [42]. In this study, many lipid droplets surrounding the mitochondria were observed by the TEM. Moreover, a large number of annotated genes were revealed to be associated with energy production and conversion. Based on the aforementioned findings, we propose that the stored lipid droplets are the main energy source for C. uncinata when there is no food or nutrient available in the living environment. However, different parasites store different energy substances for their use. For example, the ciliate Balantidium ctenopharyngodoni, inhabiting the hindgut of grass carp, stores many starch granules in the cytoplasm [43]. The intestinal trichomonad Histomonas meleagridis uses glycogen granules as its energy storage [44]. Various energy storage forms thus also could represent differences in the main metabolic pathways of these parasites.

In mitochondria, the oxidation of pyruvate by PDH generates Ace-CoA, which can then combine with oxaloacetate (OAA) to form citrate, the first substrate of the TCA cycle [45]. The pyruvate formed in glycolysis relies on the MPC to enter the mitochondria. The MPC is encoded by three homologous genes, MPC1, MPC2, and MPC3 in Saccharomyces cerevisiae, by two genes, MPC1 and MPC2, in flies, and by three genes, MPC1, MPC1-like, and MPC2, in mammals [46, 47]. In this study, only MPC1 was detected in C. uncinata. More validation experiments are needed in the future to prove whether MPC2 is expressed or not.

All TCA cycle-related enzymes were identified in C. uncinata. The TCA cycle in the mitochondrial matrix supplies NADH and FADH2 to the ETC, each of which donates a pair of electrons to the ETC in complexes I and II, respectively [48]. All ETC complexes are located on the inner membrane of the mitochondria [49]. In aerobic mitochondria, complexes I and II have an important role in generating ubiquinol with ubiquinone. Ubiquinone can be regenerated from ubiquinol by complex II, which is important for maintaining the activity of complex I [47]. Complex I of the free-living ciliate Metopus contortus appears to lack the NAD4L and NAD6 subunits [50]; complex I of N. ovalis and Blastocystis consist of 14 and 16 subunits, respectively [5]. In this study, 11 subunits of complex I are found in C. uncinata. SDHA and SDHB of complex II were detected as well, which is consistent with reports for N. ovalis, M. contortus, and M. laminarius [5, 9, 50]. However, only SDHA of ETC was detected in the Plagiopyla cf. narasimhamurtii, which may indicate that the ETC has lost its oxidative phosphorylation function [9]. The complex III subunits, CYB5 and CYC, were detected in C. uncinata in the present study, while only the Rieske protein was detected in C. porcatum and M. contortus [51]. No genes encoding complex III, IV, and V in Metopida and Clevelandellida (Armophorea) were identified, which imply that ATP synthesis via oxidative phosphorylation was lost and substituted by substrate-level phosphorylation [8].

Finally, phylogenetic analysis based on COX1 gene sequence and 3D structure prediction of its protein product confirmed the identity of this subunit. Based on these data, it can be inferred that C. uncinata has a relatively complete ETC system. Though ETC plays an important part in contributing to energy generation in aerobic eukaryotes, multiple studies have shown that the most highly conserved function among mitochondria and MROs is the Fe-S cluster biosynthesis [52, 53]. Fe-S clusters are essential for all organisms. The most widely distributed Fe-S cluster assembly system found in bacteria and mitochondria/MROs of eukaryotes is the ISC system [37]. In this study, the iron donor YAH1, the scaffold proteins ISCU, NFU1, and IBA57, and molecular chaperones GRX5 were detected, but the necessary catalytic cysteine desulfurase was not identified.

Conclusions

Collectively, this is the first report of mitochondria of C. uncinata, studied using multiple methods. We found that C. uncinata has typical mitochondria with double membrane and cristae structures. Plenty of lipid droplets were found around the mitochondria, which may be the energy storage of C. uncinata that helps its transmission from a free-living to a parasitic lifestyle. Consistent with the morphological observations, metabolic pathways associated with fatty acid metabolism were identified based on the transcriptomic data. Except for this pathway, glycolysis, TCA cycle based on the ETC, and amino acid metabolism were also identified. The present findings identify the metabolic energy sources and cast light on the putative mechanism behind the facultative parasitism of C. uncinata and contribute to the prevention and treatment of chilodonelliasis in the future.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files. The COX1 nucleotide sequence generated in this study was deposited in the GenBank database under the accession number OP672162.

Abbreviations

TEM:

Transmission electron microscopy

COG:

Clusters of Orthologous Genes

COX1:

Cytochrome c oxidase subunit 1

TCA:

Tricarboxylic acid

ETC:

Electron transport chain

ISCs:

Iron-sulfur clusters

ORFs:

Open-reading frames

BI:

Bayesian inference

ML:

Maximum likelihood

References

  1. Stairs CW, Leger MM, Roger AJ. Diversity and origins of anaerobic metabolism in mitochondria and related organelles. Philos Trans R Soc B: Biol Sci. 2015;370:20140326.

    Article  Google Scholar 

  2. Eme L, Spang A, Lombard J, Stairs CW, Ettema TJG. Archaea and the origin of eukaryotes. Nat Rev Microbiol. 2017;15:711–23.

    Article  CAS  PubMed  Google Scholar 

  3. Müller M, Mentel M, van Hellemond JJ, Henze K, Woehle C, Gould SB, et al. Biochemistry and evolution of anaerobic energy metabolism in eukaryotes. Microbiol Mol Biol Rev. 2012;76:444–95.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Benchimol M, Johnson PJ, de Souza W. Morphogenesis of the hydrogenosome: an ultrastructural study. Biol Cell. 1996;87:197–205.

    Article  CAS  PubMed  Google Scholar 

  5. de Graaf RM, Ricard G, van Alen TA, Duarte I, Dutilh BE, Burgtorf C, et al. The organellar genome and metabolic potential of the hydrogen-producing mitochondrion of Nyctotherus ovalis. Mol Biol Evol. 2011;28:2379–91.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Mi-ichi F, Yousuf MA, Nakada-Tsukui K, Nozaki T. Mitosomes in Entamoeba histolytica contain a sulfate activation pathway. Proc Natl Acad Sci. 2009;106:21731–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Burger G, Zhu Y, Littlejohn TG, Greenwood SJ, Schnare MN, Lang BF, et al. Complete sequence of the mitochondrial genome of Tetrahymena pyriformis and comparison with Paramecium aurelia mitochondrial DNA. J Mol Biol. 2000;297:365–80.

    Article  CAS  PubMed  Google Scholar 

  8. Rotterová J, Edgcomb VP, Čepička I, Beinart R. Anaerobic ciliates as a model group for studying symbioses in oxygen‐depleted environments. J Eukaryot Microbiol. 2022;69:12912.

    Article  Google Scholar 

  9. Chen ZC, Li J, Salas-Leiva DE, Chen MY, Chen SL, Li SR, et al. Group-specific functional patterns of mitochondrion-related organelles shed light on their multiple transitions from mitochondria in ciliated protists. Mar Life Sci Tech. 2022;4:609–23.

    Article  Google Scholar 

  10. Rotterová J, Salomaki E, Pánek T, Bourland W, Žihala D, Táborský P, et al. Genomics of new ciliate lineages provides insight into the evolution of obligate anaerobiosis. Curr Biol. 2020;30:2037–50.

    Article  PubMed  Google Scholar 

  11. Bu XL, Wang RQ, Bastos Gomes GB, Ban SN, Li WX, Wu SG, et al. First record of facultative parasitism of Chilodonella uncinata based on goldfish (Carassius auratus) infection model. Aquaculture. 2021;538:1–8.

    Article  Google Scholar 

  12. Bastos Gomes GB, Jerry DR, Miller TL, Hutson KS. Current status of parasitic ciliates Chilodonella spp. (Phyllopharyngea: Chilodonellidae) in freshwater fish aquaculture. J Fish Dis. 2017;40:703–15.

    Article  CAS  PubMed  Google Scholar 

  13. Li M, Wang RQ, Bastos Gomes GB, Zou H, Li WX, Wu SG, et al. Epidemiology and identification of two species of Chilodonella affecting farmed fishes in China. Vet Parasitol. 2018;264:8–17.

    Article  PubMed  Google Scholar 

  14. Mitra AK, Haldar DP. First record of Chilodonella hexasticha (Kiernik, 1909) Kahl, 1931 (Ciliophora : Chilodonellidae) infesting a freshwater fish Nandus nandus (Hamilton) from gangetic West Bengal. India Anim Biol J. 2004;54:111–8.

    Article  Google Scholar 

  15. Bu XL, Zhao WS, Li M, Li WX, Wu SG, Zou H, et al. Transcriptomic differences between free-living and parasitic Chilodonella uncinata (Alveolata, Ciliophora). Microorganisms. 2022;10:1646.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinform. 2009;10:1–9.

    Article  Google Scholar 

  17. Taylor SW, Fahy E, Zhang B, Glenn GM, Warnock DE, Wiley S, et al. Characterization of the human heart mitochondrial proteome. Nat Biotechnol. 2003;21:281–6.

    Article  CAS  PubMed  Google Scholar 

  18. Prokisch H, Scharfe C, Camp DG, Xiao WZ, David L, Andreoli C, et al. Integrative analysis of the mitochondrial proteome in yeast. PLoS Biol. 2004;2:e160.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Smith DG, Gawryluk RM, Spencer DF, Pearlman RE, Siu KW, Gray MW. Exploring the mitochondrial proteome of the ciliate protozoon Tetrahymena thermophila: direct analysis by tandem mass spectrometry. J Mol Biol. 2007;374:837–63.

    Article  CAS  PubMed  Google Scholar 

  20. Emanuelsson O, Brunak S, Von Heijne G, Nielsen H. Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007;2:953–71.

    Article  CAS  PubMed  Google Scholar 

  21. Fukasawa Y, Tsuji J, Fu SC, Tomii K, Horton P, Imai K. MitoFates: improved prediction of mitochondrial targeting sequences and their cleavage sites*[S]. Mol Cell Proteomics. 2015;14:1113–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428:726–31.

    Article  CAS  PubMed  Google Scholar 

  23. Core Team R, R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. 2019; 14: 12–21.

  24. Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch AM. ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003;31:3784–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Buchan D, Minneci F, Nugent T, Bryson K, Jones DT. Scalable web services for the PSIPRED Protein Analysis Workbench. Nucleic Acids Res. 2013;41:349–57.

    Article  Google Scholar 

  26. Rice P, Longden I, Bleasby A. EMBOSS: the European molecular biology open software suite. Trends Genet. 2000;16:276–7.

    Article  CAS  PubMed  Google Scholar 

  27. Andrew W, Martino B, Stefan B, Gabriel S, Gerardo T, Rafal G, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:296–303.

    Article  Google Scholar 

  28. Guex N, Peitsch MC, Schwede T. Automated comparative protein structure modeling with SWISS-MODEL and Swiss-PdbViewer: a historical perspective. Electrophoresis. 2010;30:162–73.

    Article  Google Scholar 

  29. Zhang D, Gao FL, Jakovli I, Zou H, Zhang J, Li WX, et al. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol Ecol Resour. 2020;20:348–55.

    Article  PubMed  Google Scholar 

  30. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kalyaanamoorthy S, Minh BQ, Wong TKF, Von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.

    Article  CAS  PubMed  Google Scholar 

  33. Nguyen LT, Schmift HA, Von Haeseler A, Minh Bui Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.

    Article  CAS  PubMed  Google Scholar 

  34. Letunic I, Bork P. Interactive tree of life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007;23:127–8.

    Article  CAS  PubMed  Google Scholar 

  35. Nasirudeen AMA, Tan KSW. Isolation and characterization of the mitochondrion-like organelle from Blastocystis hominis. J Microbiol Methods. 2004;58:101–9.

    Article  CAS  PubMed  Google Scholar 

  36. Ma R, Fan XP, Yin F, Ni B, Gu FK. Ultrastructural features of the tomont of Cryptocaryon irritans (Ciliophora: Prostomatea), a parasitic ciliate of marine fishes. Parasitology. 2017;144:720–9.

    Article  PubMed  Google Scholar 

  37. Santos HJ, Makiuchi T, Nozaki T. Reinventing an organelle: the reduced mitochondrion in parasitic protists. Trends Parasitol. 2018;34:1038–55.

    Article  CAS  PubMed  Google Scholar 

  38. Palmer JD. Organelle genomes: going, going, gone. Science. 1997;275:790–790.

    Article  CAS  PubMed  Google Scholar 

  39. Leger MM, Kolisko M, Kamikawa R, Stairs CW, Kume K, Cepicka I, et al. Organelles that illuminate the origins of Trichomonas hydrogenosomes and Giardia mitosomes. Nat Ecol Evol. 2017;1:0092.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Mashek DG, Li LO, Coleman RA. Long-chain acyl-CoA synthetases and fatty acid channeling. Future Lipidol. 2007;2:465–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Tonazzi A, Giangregorio N, Console L, Palmieri F, Indiveri C. The mitochondrial carnitine acyl-carnitine carrier (SLC25A20): molecular mechanisms of transport, role in redox sensing and interaction with drugs. Biomolecules. 2021;11:521.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Benador IY, Veliova M, Mahdaviani K, Petcherski A, Wikstrom JD, Assali EA, et al. Mitochondria bound to lipid droplets have unique bioenergetics, composition, and dynamics that support lipid droplet expansion. Cell Metab. 2018;27:869–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Zhao WS, Li M, Xiong F, Zhang D, Wu SG, Zou H, et al. Identification of intracellular bacteria in the ciliate Balantidium ctenopharyngodoni (Ciliophora, Litostomatea). J Eukaryot Microbiol. 2020;67:417–26.

    Article  CAS  PubMed  Google Scholar 

  44. Munsch M, Lotfi A, Hafez HM, Al-Quraishy S, Mehlhorn H. Light and transmission electron microscopic studies on trophozoites and cyst-like stages of Histomonas meleagridis from cultures. Parasitol Res. 2009;104:683–9.

    Article  PubMed  Google Scholar 

  45. Zangari J, Petrelli F, Maillot B, Martinou JC. The multifaceted pyruvate metabolism: role of the mitochondrial pyruvate carrier. Biomolecules. 2020;10:1068.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Bricker DK, Taylor EB, Schell JC, Orsak T, Boutron A, Chen YC, et al. A mitochondrial pyruvate carrier required for pyruvate uptake in yeast, Drosophila, and humans. Science. 2012;337:96–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Herzig S, Raemy E, Montessuit S, Veuthey JL, Zamboni N, Westermann B, et al. Identification and functional expression of the mitochondrial pyruvate carrier. Science. 2012;337:93–6.

    Article  CAS  PubMed  Google Scholar 

  48. Nolfi-Donegan D, Braganza A, Shiva S. Mitochondrial electron transport chain: oxidative phosphorylation, oxidant production, and methods of measurement. Redox Biol. 2020;37:101674.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Guo RY, Gu J, Zong SA, Wu M, Yang MJ. Structure and mechanism of mitochondrial electron transport chain. Biomoed J. 2018;41:9–20.

    Article  Google Scholar 

  50. Lewis WH, Lind AE, Sendra KM, Onsbring H, Williams TA, Esteban GF, et al. Convergent evolution of hydrogenosomes from mitochondria by gene transfer and loss. Mol Biol Evol. 2019;37:524–39.

    Article  PubMed Central  Google Scholar 

  51. Perez-Brocal V, Clark CG. Analysis of two genomes from the mitochondrion-like organelle of the intestinal parasite Blastocystis: complete sequences, gene content, and genome organization. Mol Biol Evol. 2008;25:2475–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Lill R, Kispal G. Maturation of cellular Fe–S proteins: an essential function of mitochondria. Trends Biochem Sci. 2000;25:352–6.

    Article  CAS  PubMed  Google Scholar 

  53. Hjort K, Goldberg AV, Tsaousis AD, Hirt RP, Embley TM. Diversity and reductive evolution of mitochondria among microbial eukaryotes. Philos Trans R Soc B: Biol Sci. 2010;365:713–27.

    Article  CAS  Google Scholar 

  54. Coyne RS, Hannick L, Shanmugam D, Hostetler JB, Brami D, Joardar VS, et al. Comparative genomics of the pathogenic ciliate Ichthyophthirius multifiliis, its free-living relatives and a host species provide insights into adoption of a parasitic lifestyle and prospects for disease control. Genome Biol. 2011;12:1–16.

    Article  Google Scholar 

  55. Verner Z, Basu S, Benz C, Dixit S, Dobáková E, Faktorová D, et al. Malleable mitochondrion of Trypanosoma brucei. Int Rev Cell Mol Biol Evol. 2015;315:73–151.

    Article  CAS  Google Scholar 

  56. Tachezy J, Makki A, Hrdý I. The hydrogenosomes of Trichomonas vaginalis. J Eukaryot Microbiol. 2022;69:12922.

    Article  Google Scholar 

  57. Camacho E, Rastrojo A, Sanchiz Á, González-de la Fuente S, Aguado B, Requena JM. Leishmania mitochondrial genomes: maxicircle structure and heterogeneity of minicircles. Genes. 2019;10:758.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Tsaousis AD, Keithly JS. The mitochondrion-related organelles of Cryptosporidium species, in hydrogenosomes and mitosomes: mitochondria of anaerobic eukaryotes. Berlin: Springer; 2019. p. 231–53.

    Google Scholar 

  59. Wawrzyniak I, Roussel M, Diogon M, Couloux A, Texier C, Tan KS, et al. Complete circular DNA in the mitochondria-like organelles of Blastocystis hominis. Int J Parasitol. 2008;38:1377–82.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Yuan Xiao and Zhenfei Xing for providing the TEM service (Analysis and Testing Center, Institute of Hydrobiology, CAS). We also thank Ivan Jakovlić for improving the language.

Funding

This work was supported by the National Natural Science Foundation of China (no. 32170437), the Second Tibetan Plateau Scientific Expedition and Research Program (STEP) (no. 2019QZKK0304), the earmarked fund for CARS (No. CARS-45), the Protist 10,000 Genomics Project (P10K) Consortium, and the National Aquatic Biological Resource Center (NABRC).

Author information

Authors and Affiliations

Authors

Contributions

GW and ML designed the experiments. XB made the experiments, analyzed the experimental data, and wrote the paper. WZ helped with the data analysis. WL, HZ, and SW reviewed and revised the paper. All authors have read and agreed to the published version of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ming Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Primers used in this study. Table S2. COG annotation results. Table S3. Annotation results of mitochondria metabolism.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bu, Xl., Zhao, Ws., Li, Wx. et al. Mitochondrial metabolism of the facultative parasite Chilodonella uncinata (Alveolata, Ciliophora). Parasites Vectors 16, 92 (2023). https://doi.org/10.1186/s13071-023-05695-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-023-05695-3

Keywords