A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

The Highly Repetitive Genome of Myxobolus rasmusseni, an Emerging Myxozoan Parasite of Fathead Minnows

Myxozoans are a monophyletic taxon of approximately 2,400 described species of parasites from the phylum Cnidaria. The recent focus on their negative impacts on fisheries, on their evolution from free-living ancestors, and on their emergence into new fish host populations has stressed the critical need for genomic resources for this parasitic group. Here, we describe the genome assembly and annotation of Myxobolus rasmusseni, an emerging parasite of fathead minnows in Alberta, Canada. The assembly is 174.6 Mb in size, 68% of which is made up of repetitive elements, making it one of the most repetitive animal genomes sequenced to date. Through comparisons to other myxozoans, we show that widespread gene loss, a known phenomenon of this group of parasites, is consistent with closely related species. Additionally, we assembled the M. rasmusseni mitochondrial genome, which is nearly twice the size of the typical animal mitochondrial genome yet contains only five of the canonical mitochondrial protein-coding genes and open reading frames not found in other myxozoans. These results add to our understanding of the gene- and genome-level diversity observed in myxozoans.

Evolution of myxozoan mitochondrial genomes: insights from myxobolids

BACKGROUND

Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms.

RESULTS

Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene.

CONCLUSIONS

Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

Apoptotic gene loss in Cnidaria is associated with transition to parasitism

The phylum Cnidaria consists of several morphologically diverse classes including Anthozoa, Cubozoa, Hydrozoa, Polypodiozoa, Scyphozoa, Staurozoa, and Myxozoa. Myxozoa comprises two subclasses of obligate parasites-Myxosporea and Malacosporea, which demonstrate various degrees of simplification. Myxosporea were previously reported to lack the majority of core protein domains of apoptotic proteins including caspases, Bcl-2, and APAF-1 homologs. Other sequenced Cnidaria, including the parasite Polypodium hydriforme from Polypodiozoa do not share this genetic feature. Whether this loss of core apoptotic proteins is unique to Myxosporea or also present in its sister subclass Malacosporea was not previously investigated. We show that the presence of core apoptotic proteins gradually diminishes from free-living Cnidaria to Polypodium to Malacosporea to Myxosporea. This observation does not favor the hypothesis of catastrophic simplification of Myxosporea at the genetic level, but rather supports a stepwise adaptation to parasitism that likely started from early parasitic ancestors that gave rise to Myxozoa.

Genome-Scale Analysis Reveals Extensive Diversification of Voltage-Gated K+ Channels in Stem Cnidarians

Ion channels are highly diverse in the cnidarian model organism Nematostella vectensis (Anthozoa), but little is known about the evolutionary origins of this channel diversity and its conservation across Cnidaria. Here, we examined the evolution of voltage-gated K+ channels in Cnidaria by comparing genomes and transcriptomes of diverse cnidarian species from Anthozoa and Medusozoa. We found an average of over 40 voltage-gated K+ channel genes per species, and a phylogenetic reconstruction of the Kv, KCNQ, and Ether-a-go-go (EAG) gene families identified 28 voltage-gated K+ channels present in the last common ancestor of Anthozoa and Medusozoa (23 Kv, 1 KCNQ, and 4 EAG). Thus, much of the diversification of these channels took place in the stem cnidarian lineage prior to the emergence of modern cnidarian classes. In contrast, the stem bilaterian lineage, from which humans evolved, contained no more than nine voltage-gated K+ channels. These results hint at a complexity to electrical signaling in all cnidarians that contrasts with the perceived anatomical simplicity of their neuromuscular systems. These data provide a foundation from which the function of these cnidarian channels can be investigated, which will undoubtedly provide important insights into cnidarian physiology.

Highly Dynamic Gene Family Evolution Suggests Changing Roles for PON Genes Within Metazoa

Change in gene family size has been shown to facilitate adaptation to different selective pressures. This includes gene duplication to increase dosage or diversification of enzymatic substrates and gene deletion due to relaxed selection. We recently found that the PON1 gene, an enzyme with arylesterase and lactonase activity, was lost repeatedly in different aquatic mammalian lineages, suggesting that the PON gene family is responsive to environmental change. We further investigated if these fluctuations in gene family size were restricted to mammals and approximately when this gene family was expanded within mammals. Using 112 metazoan protein models, we explored the evolutionary history of the PON family to characterize the dynamic evolution of this gene family. We found that there have been multiple, independent expansion events in tardigrades, cephalochordates, and echinoderms. In addition, there have been partial gene loss events in monotremes and sea cucumbers and what appears to be complete loss in arthropods, urochordates, platyhelminths, ctenophores, and placozoans. In addition, we show the mammalian expansion to three PON paralogs occurred in the ancestor of all mammals after the divergence of sauropsida but before the divergence of monotremes from therians. We also provide evidence of a novel PON expansion within the brushtail possum. In the face of repeated expansions and deletions in the context of changing environments, we suggest a range of selective pressures, including pathogen infection and mitigation of oxidative damage, are likely influencing the diversification of this dynamic gene family across metazoa.

Transposable element expansion and low-level piRNA silencing in grasshoppers may cause genome gigantism

Background

Transposable elements (TEs) have been likened to parasites in the genome that reproduce and move ceaselessly in the host, continuously enlarging the host genome. However, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful consequences of TE invasion by imposing small-RNA-mediated silencing. Here we compare the TE activity of two grasshopper species with different genome sizes in Acrididae (Locusta migratoria manilensis♀1C = 6.60 pg, Angaracris rhodopa♀1C = 16.36 pg) to ascertain the influence of piRNAs.

Results

We discovered that repetitive sequences accounted for 74.56% of the genome in A. rhodopa, more than 56.83% in L. migratoria, and the large-genome grasshopper contained a higher TEs proportions. The comparative analysis revealed that 41 TEs (copy number > 500) were shared in both species. The two species exhibited distinct “landscapes” of TE divergence. The TEs outbreaks in the small-genome grasshopper occurred at more ancient times, while the large-genome grasshopper maintains active transposition events in the recent past. Evolutionary history studies on TEs suggest that TEs may be subject to different dynamics and resistances in these two species. We found that TE transcript abundance was higher in the large-genome grasshopper and the TE-derived piRNAs abundance was lower than in the small-genome grasshopper. In addition, we found that the piRNA methylase HENMT, which is underexpressed in the large-genome grasshopper, impedes the piRNA silencing to a lower level.

Conclusions

Our study revealed that the abundance of piRNAs is lower in the gigantic genome grasshopper than in the small genome grasshopper. In addition, the key gene HENMT in the piRNA biogenesis pathway (Ping-Pong cycle) in the gigantic genome grasshopper is underexpressed. We hypothesize that low-level piRNA silencing unbalances the original positive correlation between TEs and piRNAs, and triggers TEs to proliferate out of control, which may be one of the reasons for the gigantism of grasshopper genomes.

Expression profiling and cellular localization of myxozoan minicollagens during nematocyst formation and sporogenesis

In free-living cnidarians, minicollagens are major structural components in the biogenesis of nematocysts. Recent sequence mining and proteomic analysis demonstrate that minicollagens are also expressed by myxozoans, a group of evolutionarily ancient cnidarian endoparasites. Nonetheless, the presence and abundance of nematocyst-associated genes/proteins in nematocyst morphogenesis have never been studied in Myxozoa. Here, we report the gene expression profiles of three myxozoan minicollagens, ncol-1, ncol-3, and the recently identified noncanonical ncol-5, during the intrapiscine development of Myxidium lieberkuehni, the myxozoan parasite of the northern pike, Esox lucius. Moreover, we localized the myxozoan-specific minicollagen Ncol-5 in the developing myxosporean stages by Western blotting, immunofluorescence, and immunogold electron microscopy. We found that expression of minicollagens was spatiotemporally restricted to developing nematocysts within the myxospores during sporogenesis. Intriguingly, Ncol-5 is localized in the walls of nematocysts and predominantly in nematocyst tubules. Overall, we demonstrate that despite being significantly reduced in morphology, myxozoans retain structural components associated with nematocyst development in free-living cnidarians. Furthermore, our findings have practical implications for future functional and comparative studies as minicollagens are useful markers of the developmental phase of myxozoan parasites.

Transcriptomic Insights into the Diversity and Evolution of Myxozoa (Cnidaria, Endocnidozoa) Toxin-like Proteins

Myxozoa is a speciose group of endoparasitic cnidarians that can cause severe ecological and economic effects. Their cnidarian affinity is affirmed by genetic relatedness and the presence of nematocysts, historically called "polar capsules". Previous studies have revealed the presence of toxin-like proteins in myxozoans; however, the diversity and evolution of venom in Myxozoa are not fully understood. Here, we performed a comparative analysis using the newly sequenced transcriptomes of five Myxobolidae species as well as some public datasets. Toxin mining revealed that myxozoans have lost most of their toxin families, while most species retained Kunitz, M12B, and CRISP, which may play a role in endoparasitism. The venom composition of Endocnidozoa (Myxozoa + Polypodium) differs from that of free-living cnidarians and may be influenced by ecological and environmental factors. Phylogenetic analyses showed that toxin families of myxozoans and free-living cnidarians were clustered into different clades. Selection analyses showed that purifying selection was the dominant evolutionary pressure in toxins, while they were still influenced by episodic adaptive selection. This suggests that the potency or specificity of a particular toxin or species might increase. Overall, our findings provide a more comprehensive framework for understanding the diversity and evolution of Myxozoa venoms.