Mitochondrial genome sequencing and analysis of scuticociliates (Uronema marinum) isolated from Takifugu rubripes

Characterization of the complete mitochondrial genome of Miamiensis avidus causing flatfish scuticociliatosis

Miamiensis avidus is a parasitic pathogen that causes the disease scuticociliatosis in teleost fish species. It is a ciliate and a free-living marine protozoan belonging to the order Philasterida, subclass Scuticociliatida, class Oligohymenophorea, and phylum Ciliophora. The complete mt-genome of M. avidus was linear and 38,695 bp in length with 47 genes, including 40 protein-coding genes, two ribosomal RNA (rRNA) genes, and five transfer RNA (tRNA) genes. Of these, 20 genes typically belong to the clusters of orthologous groups, playing roles in energy production and conversion, translation, ribosomal structure and biogenesis, and defense mechanisms. This is the first report of sequencing and characterization of the mt-genome of M. avidus, which was observed to be linear and possessing the typical ciliate mitochondrial genome organization and phylogenetic relationships. Remarkable differences were observed between M. avidus and other ciliates in the mitochondrially encoded rRNAs, extensive gene loss in ribosomal genes and tRNAs, terminal repeat sequences, and stop codon usage. A comparative and phylogenetic analysis of M. avidus and Uronema marinum of the order Hymenostomatida, which is most closely related to the order Philasterida, signified the promise of the mitogenome data of M. avidus as a valuable genetic marker in species detection and taxonomic research. The present study has potential applications in epidemiological studies and host-parasite interaction investigations facilitating disease control.

Comparison of mitochondrial genome and development of specific PCR primers for identifying two scuticociliates, Pseudocohnilembus persalinus and Uronema marinum

Background

Pseudocohnilembus persalinus and Uronema marinum (Ciliophora, Scuticociliatia), as parasitic scuticociliatid ciliates, were isolated from Scophthalmus maximus and Takifugu rubripes, respectively, in our previous studies. These ciliates are morphologically very similar; hence, it is difficult to identify specific scuticociliate species using traditional classification methods for performing taxonomic research and disease control studies.

Methods

We annotated the mitochondrial genomes of these two scuticociliates on the basis of previous sequencing, including analyses of nucleotide composition, codon usage, Ka/Ks, and p-distance. We also compared the nucleotide and amino acid similarity of the mitochondrial genomes of P. persalinus, U. marinum, and other 12 related ciliates, and a phylogenetic tree was constructed using 16 common genes. We chose the nad4 and nad7 genes to design specific PCR primers for identification.

Results

P. persalinus and U. marinum were found to have a close evolutionary relationship. Although codon preferences were similar, differences were observed in the usage of codons such as CGA, CGC, and GTC. Both Ka/Ks and p-distance were less than 1. Except for yejR, ymf57, ymf67, and ymf75, the amino acid sequence similarity between P. persalinus and U. marinum was greater than 50%.

Conclusions

The mitochondrial genomes of P. persalinus and U. marinum were thoroughly compared to provide a reference for disease prevention and control. The specific PCR primers enabled us to identify P. persalinus and U. marinum rapidly and accurately at the molecular level, thus providing a basis for classification and identification.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04821-3.

Characterization and Comparative Analyses of Mitochondrial Genomes in Single-Celled Eukaryotes to Shed Light on the Diversity and Evolution of Linear Molecular Architecture

Determination and comparisons of complete mitochondrial genomes (mitogenomes) are important to understand the origin and evolution of mitochondria. Mitogenomes of unicellular protists are particularly informative in this regard because they are gene-rich and display high structural diversity. Ciliates are a highly diverse assemblage of protists and their mitogenomes (linear structure with high A+T content in general) were amongst the first from protists to be characterized and have provided important insights into mitogenome evolution. Here, we report novel mitogenome sequences from three representatives (Strombidium sp., Strombidium cf. sulcatum, and Halteria grandinella) in two dominant ciliate lineages. Comparative and phylogenetic analyses of newly sequenced and previously published ciliate mitogenomes were performed and revealed a number of important insights. We found that the mitogenomes of these three species are linear molecules capped with telomeric repeats that differ greatly among known species. The genomes studied here are highly syntenic, but larger in size and more gene-rich than those of other groups. They also all share an AT-rich tandem repeat region which may serve as the replication origin and modulate initiation of bidirectional transcription. More generally we identified a split version of ccmf, a cytochrome c maturation-related gene that might be a derived character uniting taxa in the subclasses Hypotrichia and Euplotia. Finally, our mitogenome comparisons and phylogenetic analyses support to reclassify Halteria grandinella from the subclass Oligotrichia to the subclass Hypotrichia. These results add to the growing literature on the unique features of ciliate mitogenomes, shedding light on the diversity and evolution of their linear molecular architecture.

Incorporating mitogenome sequencing into integrative taxonomy: The multidisciplinary redescription of the ciliate Thuricola similis (Peritrichia, Vaginicolidae) provides new insights into the evolutionary relationships among Oligohymenophorea subclasses

The evolutionary relationships among Oligohymenophorea subclasses are under debate as the phylogenomic analysis using a large dataset of nuclear coding genes is significantly different to the 18S rDNA phylogeny, and it is unfortunately not stable within and across different published studies. In addition to nuclear genes, the faster-evolving mitochondrial genes have also shown the ability to solve phylogenetic problems in many ciliated taxa. However, due to the paucity of mitochondrial data, the corresponding work is scarce, let alone the phylogenomic analysis based on mitochondrial gene dataset. In this work, we presented the characterization on Thuricola similis Bock, 1963, a loricate peritrich (Oligohymenophorea), incorporating mitogenome sequencing into integrative taxonomy. As the first mitogenome for the subclass Peritrichia, it is linear, 38,802 bp long, and contains two rRNAs, 12 tRNAs, and 43 open reading frames (ORFs). As a peculiarity, it includes a central repeated region composed of tandemly repeated A-T rich units working as a bi-transcriptional start. Moreover, taking this opportunity, the phylogenomic analyses based on a set of mitochondrial genes were also performed, revealing that T. similis, as a representative of Peritrichia subclass, branches basally to other three Oligohymenophorea subclasses, namely Hymenostomatia, Peniculia, and Scuticociliatia. Evolutionary relationships among those Oligohymenophorea subclasses were discussed, also in the light of recent phylogenomic reconstructions based on a set of nuclear genes. Besides, as a little-known species, T. similis was also redescribed and neotypified based on data from two populations collected from wastewater treatment plants (WWTPs) in Brazil and Italy, by means of integrative methods (i.e., living observation, silver staining methods, scanning and transmission electron microscopy, and 18S rDNA phylogeny). After emended diagnosis, it is characterized by: (1) the sewage habitat; (2) the lorica with a single valve and small undulations; (3) the 7-22 µm-long inner stalk; and (4) the presence of only a single postciliary microtubule on the left side of the aciliferous row in the haplokinety. Among Vaginicolidae family, our 18S rRNA gene-based phylogenetic analysis revealed that Thuricola and Cothurnia are monophyletic genera, and Vaginicola could be a polyphyletic genus.

Comparative Studies on the Polymorphism and Copy Number Variation of mtSSU rDNA in Ciliates (Protista, Ciliophora): Implications for Phylogenetic, Environmental, and Ecological Research

While nuclear small subunit ribosomal DNA (nSSU rDNA) is the most commonly-used gene marker in studying phylogeny, ecology, abundance, and biodiversity of microbial eukaryotes, mitochondrial small subunit ribosomal DNA (mtSSU rDNA) provides an alternative. Recently, both copy number variation and sequence variation of nSSU rDNA have been demonstrated for diverse organisms, which can contribute to misinterpretation of microbiome data. Given this, we explore patterns for mtSSU rDNA among 13 selected ciliates (representing five classes), a major component of microbial eukaryotes, estimating copy number and sequence variation and comparing to that of nSSU rDNA. Our study reveals: (1) mtSSU rDNA copy number variation is substantially lower than that for nSSU rDNA; (2) mtSSU rDNA copy number ranges from 1.0 × 10⁴ to 8.1 × 10⁵; (3) a most common sequence of mtSSU rDNA is also found in each cell; (4) the sequence variation of mtSSU rDNA are mainly indels in poly A/T regions, and only half of species have sequence variation, which is fewer than that for nSSU rDNA; and (5) the polymorphisms between haplotypes of mtSSU rDNA would not influence the phylogenetic topology. Together, these data provide more insights into mtSSU rDNA as a powerful marker especially for microbial ecology studies.

Characterization of the complete mitochondrial genome of Barilius malabaricus and its phylogenetic implications

Bariliine fishes are important components of the ornamental as well as subsistence fishery sectors in India. Many of the species in the genus Barilius are threatened by habitat loss and therefore need to be met with conservation initiatives. Effective conservation measures, however, require clarification of species identities and resolution of the validity of many species currently treated as synonymous or sub-species. The complete mitochondrial genome data provide better insight into phylogenetic information than the short fragment or single gene based analysis. Thus, we have sequenced the complete mitochondrial genome of Barilius malabaricus, one of the important fish species in the fresh water ornamental sector, for better understanding its phylogenetic status. The 16,519 bp mitochondrial genome consists of 37 genes which classifies as 22 tRNA, 13 protein coding and 2 ribosomal RNA genes and a control region. Overall, the mitochondrial genome bears the typical gene order and composition as in other fishes. Further, the COI, cytochrome b and 16S rRNA gene sequences revealed that, B. malabaricus is genetically closer to B. canarensis and B. bakeri compared to other Barilius species. Also, the Barilius species of west flowing rivers in Western Ghats were consistently recovered as a clade distinct from other species. We therefore suggest to retain the genus name Barilius for the species from the Western Ghats until a comprehensive analysis based on both morphological and molecular markers reveals the relationship between species now variously placed in the genera Barilius and Opsarius in greater detail.

The complete mitochondrial genome of Gruberia lanceolata (Gruber, 1884) Kahl, 1932 (Ciliophora: Heterotrichea)

The ciliate Gruberia lanceolata (Gruber, 1884) Kahl, 1932 belonging to the class Heterotrichea was sampled from the coastal waters of South Korea. The complete mitogenome in its linear form and large size (∼40 kb) was obtained. It consisted of 27 protein-coding genes (PCGs), two ribosomal subunit RNA (rRNA) genes, four transfer RNAs (tRNAs), and ten unclassified open reading frames (ORFs). Their telomeric structures were capped, with repeat regions at both ends. We analyzed its phylogenetic tree using the data of its respiratory chain complex I genes. It can be suggested that the complete mitochondrial genome of G. lanceolata can be recorded as a new class of the mitogenome.

Mitochondrial genome sequencing of the monogenean Heterobothrium okamotoi isolated from the tiger puffer Takifugu rubripes in North China

In this study, monogenean Heterobothium okamotoi was isolated and identified from the gill of diseased Tiger puffer (T. rubripes) at an industrial farm in Liaoning, North China (121.3459 E, 38.9861 N). With the completion of H. okamotoi mitochondrial genome sequencing, the full-length mitochondrial genome of H. okamotoi was assembled and analyzed. All results indicate that the complete mitochondrial genome of H. okamotoi was 14,643 bp. There were 2 rRNAs, 20 tRNAs, and 12 protein-coding genes (PCGs) all located at the heavy (H) strand. Besides, the phylogenetic tree of 19 monogeneans was constructed. The results showed that H. okamotoi and Pseudochauhanea macrorchis were clustered in a clade. To sum up, our research results would further provide essential data for systematics and evolution study of H. okamotoi.

Mitochondrial genome sequencing of notorious scuticociliates (Pseudocohnilembus persalinus) isolated from Turbot (Scophthalmus maximus L.)

Previously, a pathogenic ciliate was isolated from the surface ulcer of a diseased Turbot (Scophthalmus maximus L.) at an aquaculture farm in North China. After morphological and molecular biological identification based on 18rRNA, the ciliated was identified as the notorious scuticociliates (Pseudocohnilembus persalinus). In this study, the whole sequence of the mitochondrial genomic gene of P. persalinus was carried out. The sequencing results showed that the complete sequence of P. persalinus mitogenome was 38,375 bp. There were 2 rRNAs, 4 tRNAs, and 34 protein-coding genes (PCGs), respectively, located on both the heavy strand and the light stand. 15 PCGs were on the heavy strand, and 19 PCGs on the light strand. Besides, phylogenetic trees among 11 ciliates were constructed based on the sequences of 17 PCGs located in mitogenome using BI methods. The results of clustering showed that P. persalinus and Uronema marinum was the first cluster belonging to the order Scuticociliatida. Our research results will further provide primary data for evolution and classified study of scuticociliates.