Genomic insights into the evolutionary origin of Myxozoa within Cnidaria

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.

Fluorescence In Situ Hybridization as a Tool for Studying the Specification and Differentiation of Cell Types in Nematostella vectensis

The sea anemone Nematostella vectensis is a genetically tractable cnidarian species that has become a model organism for studying the evolution of developmental processes and genome regulation, resilience to fluctuations in environmental conditions, and the response to pollutants. Gene expression analyses are central to many of these studies, and in situ hybridization has been an important method for obtaining spatial information, in particular during embryonic development. Like other cnidarians, Nematostella embryos are of comparably low morphological complexity, but they possess many cell types that are dispersed throughout the tissue and originate from broad and overlapping areas. These features have made two-color fluorescence in situ hybridization an important method to determine potential co-expression of genes and to generate hypotheses for their functions in cell fate specification. We here share protocols for single and double fluorescence in situ hybridization in Nematostella and for the combination of fluorescence in situ hybridization and immunofluorescence.

Establishing Bilateral Symmetry in Hydrozoan Planula Larvae, a Review of Siphonophore Early Development

Siphonophores are colonial hydrozoans, characterized by complex colony organization and unparalleled zooid functional specialization. Recent genomic studies have offered an evolutionary perspective on how this morphological complexity arose, but a molecular characterization of symmetry breaking in siphonophore embryonic development is still largely missing. Here, bringing together historical data on early development with new immunohistochemical data, we review the diversity of developmental trajectories that lead to the formation of bilaterally symmetric planula larvae in siphonophores. Embryonic development, up to the planula stage, is remarkably similar across siphonophore phylogeny. Then, with the appearance of the lateral endodermal thickening (= ventral endoderm), larval development diverges between taxa, differing in the location and patterning of the primary buds, chronology of budding, establishment of growth zones, and retention of larval zooids. Our work also uncovers a number of open questions in siphonophore development, including homology of different zooids, mechanisms underlying formation and maintenance of spatially restricted growth zone(s), and molecular factors establishing a secondary dorsal-ventral axis in planulae. By discussing siphonophore development and body axes within the broader cnidarian context, we then set the framework for future work on siphonophores, which is finally achievable with the advent of culturing methods.

Myxobolus dumerilii sp. n. (Myxozoa: Myxobolidae) infecting the brain of Chinese longsnout catfish Tachysurus dumerili (Bleeker) in China

During an investigation of Myxobolus diversity in the Chinese longsnout catfish Tachysurus dumerili (Bleeker), a new species infecting the intracranial epidermis of the host was discovered. Upon opening the cranial cavity, several round whitish plasmodia measuring 0.55-0.80 mm in diameter were observed. Fresh spores (n= 50) were pyriform in the frontal view and fusiform in the sutural view, with a length of 15.4±0.6 (13.9-16.5) μm, width of 9.1±0.4 (8.3-9.8) μm, and thickness of 7.0±0.4 (6.3-7.9) μm. The spores had smooth shell surfaces and transparent membrane sheaths in the posterior. No folds, intercapsular appendix, and caudal appendages were observed. Two equal polar capsules were pyriform and measured 7.5±0.5 (6.7-8.7) μm in length and 3.2±0.3 (2.5-3.6) μm in width. The polar filaments were coiled with five to six turns and perpendicular to the polar capsule length. A BLAST search indicated M. dumerilii sp. n. was closely related to five Myxobolus species (with sequences similarities ranging from 90.54% to 96.52%) found in different organs of yellow catfish Tachysurus fulvidraco (Richardson), rather than the T. dumerili-infecting species M. branchiola Dong and Zhao, 2014 (with 90.5% sequence similarity). Phylogenetic analysis showed that M. dumerilii sp. n. didn't form sister clade with brain-infecting Myxobolus spp, but clustered with M. jianlinensis Gao et Zhao, 2020 and M. voremkhai Akhmerov, 1960 within the Siluriformes-clade with highly supported values, indicating that the host specificity may play a stronger signal than site infections during the evolution of Myxobolus species. Based on the morphological, ecological, and molecular differences observed between the newly discovered species and other available Myxobolus species, M. dumerilii sp. n., is proposed and described in this study.

Animal evolution: No cilia, no problem!

Cilia are fundamental organelles of eukaryotes with diverse functions spanning from cell motility to sensory perception. A new study presenting genomes from parasitic horsehair worms reveals highly reduced genomes that have lost the molecular machinery needed to make cilia.

Rampant loss of universal metazoan genes revealed by a chromosome-level genome assembly of the parasitic Nematomorpha

Parasites may manipulate host behavior to increase the odds of transmission or to reach the proper environment to complete their life cycle.1,2 Members of the phylum Nematomorpha (known as horsehair worms, hairworms, or Gordian worms) are large endoparasites that affect the behavior of their arthropod hosts. In terrestrial hosts, they cause erratic movements toward bodies of water,3,4,5,6 where the adult worm emerges from the host to find mates for reproduction. We present a chromosome-level genome assembly for the freshwater Acutogordius australiensis and a draft assembly for one of the few known marine species, Nectonema munidae. The assemblies span 201 Mbp and 213 Mbp in length (N50: 38 Mbp and 716 Kbp), respectively, and reveal four chromosomes in Acutogordius, which are largely rearranged compared to the inferred ancestral condition in animals. Both nematomorph genomes have a relatively low number of genes (11,114 and 8,717, respectively) and lack a high proportion (∼30%) of universal single-copy metazoan orthologs (BUSCO genes7). We demonstrate that missing genes are not an artifact of the assembly process, with the majority of missing orthologs being shared by the two independent assemblies. Missing BUSCOs are enriched for Gene Ontology (GO) terms associated with the organization of cilia and cell projections in other animals. We show that most cilium-related genes conserved across eukaryotes have been lost in Nematomorpha, providing a molecular basis for the suspected absence of ciliary structures in these animals.

The Molecular Mechanisms Employed by the Parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) from Invasion through Sporulation for Successful Proliferation in Its Fish Host

Myxozoa is a unique group of obligate endoparasites in the phylum Cnidaria that can cause emerging diseases in wild and cultured fish populations. Recently, we identified a new myxozoan species, Myxobolus bejeranoi, which infects the gills of cultured tilapia while suppressing host immunity. To uncover the molecular mechanisms underlying this successful parasitic strategy, we conducted transcriptomics analysis of M. bejeranoi throughout the infection. Our results show that histones, which are essential for accelerated cell division, are highly expressed even one day after invasion. As the infection progressed, conserved parasitic genes that are known to modulate the host immune reaction in different parasitic taxa were upregulated. These genes included energy-related glycolytic enzymes, as well as calreticulin, proteases, and miRNA biogenesis proteins. Interestingly, myxozoan calreticulin formed a distinct phylogenetic clade apart from other cnidarians, suggesting a possible function in parasite pathogenesis. Sporogenesis was in its final stages 20 days post-exposure, as spore-specific markers were highly expressed. Lastly, we provide the first catalog of transcription factors in a Myxozoa species, which is minimized compared to free-living cnidarians and is dominated by homeodomain types. Overall, these molecular insights into myxozoan infection support the concept that parasitic strategies are a result of convergent evolution.

Data-specific substitution models improve protein-based phylogenetics

Calculating amino-acid substitution models that are specific for individual protein data sets is often difficult due to the computational burden of estimating large numbers of rate parameters. In this study, we tested the computational efficiency and accuracy of five methods used to estimate substitution models, namely Codeml, FastMG, IQ-TREE, P4 (maximum likelihood), and P4 (Bayesian inference). Data-specific substitution models were estimated from simulated alignments (with different lengths) that were generated from a known simulation model and simulation tree. Each of the resulting data-specific substitution models was used to calculate the maximum likelihood score of the simulation tree and simulated data that was used to calculate the model, and compared with the maximum likelihood scores of the known simulation model and simulation tree on the same simulated data. Additionally, the commonly-used empirical models, cpREV and WAG, were assessed similarly. Data-specific models performed better than the empirical models, which under-fitted the simulated alignments, had the highest difference to the simulation model maximum-likelihood score, clustered further from the simulation model in principal component analysis ordination, and inferred less accurate trees. Data-specific models and the simulation model shared statistically indistinguishable maximum-likelihood scores, indicating that the five methods were reasonably accurate at estimating substitution models by this measure. Nevertheless, tree statistics showed differences between optimal maximum likelihood trees. Unlike other model estimating methods, trees inferred using data-specific models generated with IQ-TREE and P4 (maximum likelihood) were not significantly different from the trees derived from the simulation model in each analysis, indicating that these two methods alone were the most accurate at estimating data-specific models. To show the benefits of using data-specific protein models several published data sets were reanalysed using IQ-TREE-estimated models. These newly estimated models were a better fit to the data than the empirical models that were used by the original authors, often inferred longer trees, and resulted in different tree topologies in more than half of the re-analysed data sets. The results of this study show that software availability and high computation burden are not limitations to generating better-fitting data-specific amino-acid substitution models for phylogenetic analyses.

Expansion and collapse of VEGF diversity in major clades of the animal kingdom

Together with the platelet-derived growth factors (PDGFs), the vascular endothelial growth factors (VEGFs) form the PDGF/VEGF subgroup among cystine knot growth factors. The evolutionary relationships within this subgroup have not been examined thoroughly to date. Here, we comprehensively analyze the PDGF/VEGF growth factors throughout all animal phyla and propose a phylogenetic tree. Vertebrate whole-genome duplications play a role in expanding PDGF/VEGF diversity, but several limited duplications are necessary to account for the temporal pattern of emergence. The phylogenetically oldest PDGF/VEGF-like growth factor likely featured a C-terminus with a BR3P signature, a hallmark of the modern-day lymphangiogenic growth factors VEGF-C and VEGF-D. Some younger VEGF genes, such as VEGFB and PGF, appeared completely absent in important vertebrate clades such as birds and amphibia, respectively. In contrast, individual PDGF/VEGF gene duplications frequently occurred in fish on top of the known fish-specific whole-genome duplications. The lack of precise counterparts for human genes poses limitations but also offers opportunities for research using organisms that diverge considerably from humans. Sources for the graphical abstract: 326 MYA and older [1]; 72-240 MYA [2]; 235-65 MYA [3].

Coevolution of the Tlx homeobox gene with medusa development (Cnidaria: Medusozoa)

Cnidarians display a wide diversity of life cycles. Among the main cnidarian clades, only Medusozoa possesses a swimming life cycle stage called the medusa, alternating with a benthic polyp stage. The medusa stage was repeatedly lost during medusozoan evolution, notably in the most diverse medusozoan class, Hydrozoa. Here, we show that the presence of the homeobox gene Tlx in Cnidaria is correlated with the presence of the medusa stage, the gene having been lost in clades that ancestrally lack a medusa (anthozoans, endocnidozoans) and in medusozoans that secondarily lost the medusa stage. Our characterization of Tlx expression indicate an upregulation of Tlx during medusa development in three distantly related medusozoans, and spatially restricted expression patterns in developing medusae in two distantly related species, the hydrozoan Podocoryna carnea and the scyphozoan Pelagia noctiluca. These results suggest that Tlx plays a key role in medusa development and that the loss of this gene is likely linked to the repeated loss of the medusa life cycle stage in the evolution of Hydrozoa.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes

BACKGROUND

Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera.

RESULTS

We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species.

CONCLUSIONS

The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

Ortholinea nupchi n. sp. (Myxosporea: Ortholineidae) from the urinary bladder of the cultured olive flounder Paralichthys olivaceus, South Korea

A new myxosporean parasite, Ortholinea nupchi n. sp. (Myxozoa; Bivalvulida), was isolated from the urinary bladder of the olive flounder Paralichthys olivaceus cultured on Jeju Island, Korea. Mature spores were subspherical in the valvular and apical views and ellipsoidal in the sutural view. The spores measured 7.6 ± 0.5 μm in length, 6.7 ± 0.3 μm in thickness, and 7.3 ± 0.5 μm in width. Two pyriform polar capsules measured 3.2 ± 0.1 μm in length and 2.7 ± 0.1 μm in width and were located at the same level at the anterior half of the myxospores. The suture line was straight in the middle of the spores, and the surface ridges ranged between five and seven, forming an intricate pattern. The result of the 18S rDNA comparison showed ≤ 93.0% similarity with other Ortholinea species. The phylogenetic tree demonstrated that O. nupchi n. sp. was closest to O. auratae and clustered with oligochaete-infecting myxosporeans (OIM) having urinary system infection tropism. Based on the comparison of environmental and host factors in the phylogenetic groups of the OIM clade, we propose that the infection of O. nupchi n. sp. originated from marine oligochaetes.

GenomeFLTR: filtering reads made easy

In the last decade, advances in sequencing technology have led to an exponential increase in genomic data. These new data have dramatically changed our understanding of the evolution and function of genes and genomes. Despite improvements in sequencing technologies, identifying contaminated reads remains a complex task for many research groups. Here, we introduce GenomeFLTR, a new web server to filter contaminated reads. Reads are compared against existing sequence databases from various representative organisms to detect potential contaminants. The main features implemented in GenomeFLTR are: (i) automated updating of the relevant databases; (ii) fast comparison of each read against the database; (iii) the ability to create user-specified databases; (iv) a user-friendly interactive dashboard to investigate the origin and frequency of the contaminations; (v) the generation of a contamination-free file. Availability: https://genomefltr.tau.ac.il/.

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

Myxozoans (Cnidaria) do not Retain Key Oxygen-Sensing and Homeostasis Toolkit Genes

For aerobic organisms, both the hypoxia-inducible factor pathway and the mitochondrial genomes are key players in regulating oxygen homeostasis. Recent work has suggested that these mechanisms are not as highly conserved as previously thought, prompting more surveys across animal taxonomic levels, which would permit testing of hypotheses about the ecological conditions facilitating evolutionary loss of such genes. The Phylum Cnidaria is known to harbor wide variation in mitochondrial chromosome morphology, including an extreme example, in the Myxozoa, of mitochondrial genome loss. Because myxozoans are obligate endoparasites, frequently encountering hypoxic environments, we hypothesize that variation in environmental oxygen availability could be a key determinant in the evolution of metabolic gene networks associated with oxygen-sensing, hypoxia-response, and energy production. Here, we surveyed genomes and transcriptomes across 46 cnidarian species for the presence of HIF pathway members, as well as for an assortment of hypoxia, mitochondrial, and stress-response toolkit genes. We find that presence of the HIF pathway, as well as number of genes associated with mitochondria, hypoxia, and stress response, do not vary in parallel to mitochondrial genome morphology. More interestingly, we uncover evidence that myxozoans have lost the canonical HIF pathway repression machinery, potentially altering HIF pathway functionality to work under the specific conditions of their parasitic lifestyles. In addition, relative to other cnidarians, myxozoans show loss of large proportions of genes associated with the mitochondrion and involved in response to hypoxia and general stress. Our results provide additional evidence that the HIF regulatory machinery is evolutionarily labile and that variations in the canonical system have evolved in many animal groups.

Exploring genome gene content and morphological analysis to test recalcitrant nodes in the animal phylogeny

An accurate phylogeny of animals is needed to clarify their evolution, ecology, and impact on shaping the biosphere. Although datasets of several hundred thousand amino acids are nowadays routinely used to test phylogenetic hypotheses, key deep nodes in the metazoan tree remain unresolved: the root of animals, the root of Bilateria, and the monophyly of Deuterostomia. Instead of using the standard approach of amino acid datasets, we performed analyses of newly assembled genome gene content and morphological datasets to investigate these recalcitrant nodes in the phylogeny of animals. We explored extensively the choices for assembling the genome gene content dataset and model choices of morphological analyses. Our results are robust to these choices and provide additional insights into the early evolution of animals, they are consistent with sponges as the sister group of all the other animals, the worm-like bilaterian lineage Xenacoelomorpha as the sister group of the other Bilateria, and tentatively support monophyletic Deuterostomia.

Gene Loss may have Shaped the Cnidarian and Bilaterian Hox and ParaHox Complement

Hox and ParaHox transcription factors are important for specifying cell fates along the primary body axes during the development of most animals. Within Cnidaria, much of the research on Hox/ParaHox genes has focused on Anthozoa (anemones and corals) and Hydrozoa (hydroids) and has concentrated on the evolution and function of cnidarian Hox genes in relation to their bilaterian counterparts. Here we analyze together the full complement of Hox and ParaHox genes from species representing all four medusozoan classes (Staurozoa, Cubozoa, Hydrozoa, and Scyphozoa) and both anthozoan classes (Octocorallia and Hexacorallia). Our results show that Hox genes involved in patterning the directive axes of anthozoan polyps are absent in the stem leading to Medusozoa. For the first time, we show spatial and temporal expression patterns of Hox and ParaHox genes in the upside-down jellyfish Cassiopea xamachana (Scyphozoa), which are consistent with diversification of medusozoan Hox genes both from anthozoans and within medusozoa. Despite unprecedented taxon sampling, our phylogenetic analyses, like previous studies, are characterized by a lack of clear homology between most cnidarian and bilaterian Hox and Hox-related genes. Unlike previous studies, we propose the hypothesis that the cnidarian-bilaterian ancestor possessed a remarkably large Hox complement and that extensive loss of Hox genes was experienced by both cnidarian and bilaterian lineages.

Phylogenomics and the first higher taxonomy of Placozoa, an ancient and enigmatic animal phylum

Placozoa is an ancient phylum of extraordinarily unusual animals: miniscule, ameboid creatures that lack most fundamental animal features. Despite high genetic diversity, only recently have the second and third species been named. While prior genomic studies suffer from incomplete placozoan taxon sampling, we more than double the count with protein sequences from seven key genomes and produce the first nuclear phylogenomic reconstruction of all major placozoan lineages. This leads us to the first complete Linnaean taxonomic classification of Placozoa, over a century after its discovery: This may be the only time in the 21st century when an entire higher taxonomy for a whole animal phylum is formalized. Our classification establishes 2 new classes, 4 new orders, 3 new families, 1 new genus, and 1 new species, namely classes Polyplacotomia and Uniplacotomia; orders Polyplacotomea, Trichoplacea, Cladhexea, and Hoilungea; families Polyplacotomidae, Cladtertiidae, and Hoilungidae; and genus Cladtertia with species Cladtertia collaboinventa, nov. Our likelihood and gene content tree topologies refine the relationships determined in previous studies. Adding morphological data into our phylogenomic matrices suggests sponges (Porifera) as the sister to other animals, indicating that modest data addition shifts this node away from comb jellies (Ctenophora). Furthermore, by adding the first genomic protein data of the exceptionally distinct and branching Polyplacotoma mediterranea, we solidify its position as sister to all other placozoans; a divergence we estimate to be over 400 million years old. Yet even this deep split sits on a long branch to other animals, suggesting a bottleneck event followed by diversification. Ancestral state reconstructions indicate large shifts in gene content within Placozoa, with Hoilungia hongkongensis and its closest relatives having the most unique genetics.

Infection by the Parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) Suppresses the Immune System of Hybrid Tilapia

Myxozoa (Cnidaria) is a large group of microscopic obligate endoparasites that can cause emerging diseases, affecting wild fish populations and fisheries. Recently, the myxozoan Myxobolus bejeranoi was found to infect the gills of hybrid tilapia (Nile tilapia (Oreochromis niloticus) × Jordan/blue tilapia (O. aureus)), causing high morbidity and mortality. Here, we used comparative transcriptomics to elucidate the molecular processes occurring in the fish host following infection by M. bejeranoi. Fish were exposed to pond water containing actinospores for 24 h and the effects of minor, intermediate, and severe infections on the sporulation site, the gills, and on the hematopoietic organs, head kidney and spleen, were compared. Enrichment analysis for GO and KEGG pathways indicated immune system activation in gills at severe infection, whereas in the head kidney a broad immune suppression included deactivation of cytokines and GATA3 transcription factor responsible for T helper cell differentiation. In the spleen, the cytotoxic effector proteins perforin and granzyme B were downregulated and insulin, which may function as an immunomodulatory hormone inducing systemic immune suppression, was upregulated. These findings suggest that M. bejeranoi is a highly efficient parasite that disables the defense mechanisms of its fish host hybrid tilapia.

The Phylogenetic Position of the Enigmatic, Polypodium hydriforme (Cnidaria, Polypodiozoa): Insights from Mitochondrial Genomes

Polypodium hydriforme is an enigmatic parasite that belongs to the phylum Cnidaria. Its taxonomic position has been debated: whereas it was previously suggested to be part of Medusozoa, recent phylogenomic analyses based on nuclear genes support the view that P. hydriforme and Myxozoa form a clade called Endocnidozoa. Medusozoans have linear mitochondrial (mt) chromosomes, whereas myxozoans, as most metazoan species, have circular chromosomes. In this work, we determined the structure of the mt genome of P. hydriforme, using Illumina and Oxford Nanopore Technologies reads, and showed that it is circular. This suggests that P. hydriforme is not nested within Medusozoa, as this would entail linearization followed by recirculation. Instead, our results support the view that P. hydriforme is a sister clade to Myxozoa, and mt linearization in the lineage leading to medusozoans occurred after the divergence of Myxozoa + P. hydriforme. Detailed analyses of the assembled P. hydriforme mt genome show that: (1) it is encoded on a single circular chromosome with an estimated size of ∼93,000 base pairs, making it one of the largest metazoan mt genomes; (2) around 78% of the genome encompasses a noncoding region composed of several repeat types; (3) similar to Myxozoa, no mt tRNAs were identified; (4) the codon TGA is a stop codon and does not encode for tryptophan as in other cnidarians; (5) similar to myxozoan mt genomes, it is extremely fast evolving.

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.

Gene Manipulation in Hydractinia

The ability to regenerate lost body parts is irregularly distributed among animals, with substantial differences in regenerative potential between and within metazoan phyla. It is widely believed that regenerative animal clades inherited some aspects of their capacity to regenerate from their common ancestors but have also evolved new mechanisms that are not shared with other regenerative animals. Therefore, to gain a broad understanding of animal regenerative mechanisms and evolution, a broad sampling approach is necessary. Unfortunately, only few regenerative animals have been established as laboratory models with protocols for functional gene studies. Here, we describe the methods to establish transgenic individuals of the marine cnidarian Hydractinia. We also provide methods for transient gene expression manipulation without modifying the genome of the animals.

Molecular (co)evolution of hymenopteran social parasites and their hosts

Social parasitism describes a fascinating way of life in which species exploit the altruistic behaviour of closely related, social species. Social parasites have repeatedly evolved in the social Hymenoptera, including ants, bees, and wasps. The common ancestry and shared (social) environment with their hosts facilitates the study of molecular adaptations to the parasitic lifestyle. Moreover, when social parasites are widespread and virulent, they exert strong selection pressure on their hosts, leading to the evolution of defense mechanisms and triggering a coevolutionary arms race. Recent advances in sequencing technology now make it possible to study the molecular basis of this coevolutionary process. In addition to describing the latest developments, we highlight open research questions that could be tackled with genomic, transcriptomic, or epigenetic data.

How to build single-celled cnidarians with worm-like motility: Lessons from Myxozoa

Metazoans with worm-like morphologies across diverse and disparate groups typically demonstrate motility generated by hydrostatic skeletons involving tissue layers (muscles and epithelia). Here we present representative morphological, behavioural and molecular data for parasitic cnidarians (myxozoans) that demonstrate unprecedented variation in form and function, developing as cellular hydrostats. Motile elongate plasmodia characterise a remarkable radiation of species in the genus Ceratomyxa. The vermiform plasmodia inhabit gall bladders of a range of South American freshwater fish and exhibit undulatory motility reminiscent of nematodes but achieved at the cellular level. Collective insights from ultrastructure, confocal and light microscopy along with videos depicting movements highlight key features that we propose explain the unique motility of the plasmodia. These features include cytoskeletal elements (net forming microfilaments and microtubules), a large internal vacuole, a relatively rigid outer glycocalyx and peripherally arranged mitochondria. These constituents provide collective evidence for repurposing of the cnidarian epitheliomuscular cell to support worm-like motility at the cellular level. The apparent restriction of vermiform ceratomyxids to South American freshwaters suggests an origination via Cretaceous or Miocene marine transgressions and subsequent radiation.

Studying of Molecular Regulation of Developmental Processes of Lower Metazoans Exemplified by Cnidaria Using High-Throughput Sequencing

A unique set of features and characteristics of species of the Cnidaria phylum is the one reason that makes them a model for a various studies. The plasticity of a life cycle and the processes of cell differentiation and development of an integral multicellular organism associated with it are of a specific scientific interest. A new stage of development of molecular genetic methods, including methods for high-throughput genome, transcriptome, and epigenome sequencing, both at the level of the whole organism and at the level of individual cells, makes it possible to obtain a detailed picture of the development of these animals. This review examines some modern approaches and advances in the reconstruction of the processes of ontogenesis of cnidarians by studying the regulatory signal transduction pathways and their interactions.

A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

Background

Parasite evolution has been conceptualized as a process of genetic loss and simplification. Contrary to this model, there is evidence of expansion and conservation of gene families related to essential functions of parasitism in some parasite genomes, reminiscent of widespread mosaic evolution—where subregions of a genome have different rates of evolutionary change. We found evidence of mosaic genome evolution in the cnidarian Myxobolus honghuensis, a myxozoan parasite of fish, with extremely simple morphology.

Results

We compared M. honghuensis with other myxozoans and free-living cnidarians, and determined that it has a relatively larger myxozoan genome (206 Mb), which is less reduced and less compact due to gene retention, large introns, transposon insertion, but not polyploidy. Relative to other metazoans, the M. honghuensis genome is depleted of neural genes and has only the simplest animal immune components. Conversely, it has relatively more genes involved in stress resistance, tissue invasion, energy metabolism, and cellular processes compared to other myxozoans and free-living cnidarians. We postulate that the expansion of these gene families is the result of evolutionary adaptations to endoparasitism. M. honghuensis retains genes found in free-living Cnidaria, including a reduced nervous system, myogenic components, ANTP class Homeobox genes, and components of the Wnt and Hedgehog pathways.

Conclusions

Our analyses suggest that the M. honghuensis genome evolved as a mosaic of conservative, divergent, depleted, and enhanced genes and pathways. These findings illustrate that myxozoans are not as genetically simple as previously regarded, and the evolution of some myxozoans is driven by both genomic streamlining and expansion.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01249-8.

Proteases as Therapeutic Targets Against the Parasitic Cnidarian Ceratonova shasta: Characterization of Molecules Key to Parasite Virulence In Salmonid Hosts

Proteases and their inhibitors play critical roles in host-parasite interactions and in the outcomes of infections. Ceratonova shasta is a myxozoan pathogen that causes enteronecrosis in economically important salmonids from the Pacific Northwest of North America. This cnidarian parasite has host-specific genotypes with varying virulence, making it a powerful system to decipher virulence mechanisms in myxozoans. Using C. shasta genome and transcriptome, we identified four proteases of different catalytic types: cathepsin D (aspartic), cathepsin L and Z-like (cysteine) and aminopeptidase-N (metallo); and a stefin (cysteine protease inhibitor), which implied involvement in virulence and hence represent target molecules for the development of therapeutic strategies. We characterized, annotated and modelled their 3D protein structure using bioinformatics and computational tools. We quantified their expression in C. shasta genotype 0 (low virulence, no mortality) and IIR (high virulence and mortality) in rainbow trout Oncorhynchus mykiss, to demonstrate that there are major differences between the genotypes during infection and parasite development. High proliferation of genotype IIR was associated with high expression of the cathepsin D and the stefin, likely correlated with high nutrient demands and to regulate cell metabolism, with upregulation preceding massive proliferation and systemic dispersion. In contrast, upregulation of the cathepsin L and Z-like cysteine proteases may have roles in host immune evasion in genotype 0 infections, which are associated with low proliferation, low inflammation and non-destructive development. In contrast to the other proteases, C. shasta aminopeptidase-N appears to have a prominent role in nematocyst formation in both genotypes, but only during sporogenesis. Homology searches of C. shasta proteases against other myxozoan transcriptomes revealed a high abundance of cathepsin L and aminopeptidase homologs suggesting common gene requirements across species. Our study identified molecules of potential therapeutic significance for aquaculture and serves as a baseline for future research aimed at functional characterisation of these targets.

The cnidarian parasite Ceratonova shasta utilizes inherited and recruited venom-like compounds during infection

Background

Cnidarians are the most ancient venomous organisms. They store a cocktail of venom proteins inside unique stinging organelles called nematocysts. When a cnidarian encounters chemical and physical cues from a potential threat or prey animal, the nematocyst is triggered and fires a harpoon-like tubule to penetrate and inject venom into the prey. Nematocysts are present in all Cnidaria, including the morphologically simple Myxozoa, which are a speciose group of microscopic, spore-forming, obligate parasites of fish and invertebrates. Rather than predation or defense, myxozoans use nematocysts for adhesion to hosts, but the involvement of venom in this process is poorly understood. Recent work shows some myxozoans have a reduced repertoire of venom-like compounds (VLCs) relative to free-living cnidarians, however the function of these proteins is not known.

Methods

We searched for VLCs in the nematocyst proteome and a time-series infection transcriptome of Ceratonova shasta, a myxozoan parasite of salmonid fish. We used four parallel approaches to detect VLCs: BLAST and HMMER searches to preexisting cnidarian venom datasets, the machine learning tool ToxClassifier, and structural modeling of nematocyst proteomes. Sequences that scored positive by at least three methods were considered VLCs. We then mapped their time-series expressions in the fish host and analyzed their phylogenetic relatedness to sequences from other venomous animals.

Results

We identified eight VLCs, all of which have closely related sequences in other myxozoan datasets, suggesting a conserved venom profile across Myxozoa, and an overall reduction in venom diversity relative to free-living cnidarians. Expression of the VLCs over the 3-week fish infection varied considerably: three sequences were most expressed at one day post-exposure in the fish’s gills; whereas expression of the other five VLCs peaked at 21 days post-exposure in the intestines, coinciding with the formation of mature parasite spores with nematocysts. Expression of VLC genes early in infection, prior to the development of nematocysts, suggests venoms in C. shasta have been repurposed to facilitate parasite invasion and proliferation within the host. Molecular phylogenetics suggested some VLCs were inherited from a cnidarian ancestor, whereas others were more closely related to sequences from venomous non-Cnidarian organisms and thus may have gained qualities of venom components via convergent evolution. The presence of VLCs and their differential expression during parasite infection enrich the concept of what functions a “venom” can have and represent targets for designing therapeutics against myxozoan infections.

Quantitative Insights into the Contribution of Nematocysts to the Adaptive Success of Cnidarians Based on Proteomic Analysis

Nematocysts are secretory organelles in cnidarians that play important roles in predation, defense, locomotion, and host invasion. However, the extent to which nematocysts contribute to adaptation and the mechanisms underlying nematocyst evolution are unclear. Here, we investigated the role of the nematocyst in cnidarian evolution based on eight nematocyst proteomes and 110 cnidarian transcriptomes/genomes. We detected extensive species-specific adaptive mutations in nematocyst proteins (NEMs) and evidence for decentralized evolution, in which most evolutionary events involved non-core NEMs, reflecting the rapid diversification of NEMs in cnidarians. Moreover, there was a 33-55 million year macroevolutionary lag between nematocyst evolution and the main phases of cnidarian diversification, suggesting that the nematocyst can act as a driving force in evolution. Quantitative analysis revealed an excess of adaptive changes in NEMs and enrichment for positively selected conserved NEMs. Together, these findings suggest that nematocysts may be key to the adaptive success of cnidarians and provide a reference for quantitative analyses of the roles of phenotypic novelties in adaptation.

OUP accepted manuscript

DNA methylation, an important component of eukaryotic epigenetics, varies in pattern and function across Metazoa. Notably, bilaterian vertebrates and invertebrates differ dramatically in gene body methylation (GbM). Using the frequency of cytosine-phospho-guanines (CpGs), which are lost through mutation when methylated, we report the first broad survey of DNA methylation in Cnidaria, the ancient sister group to Bilateria. We find that: 1) GbM differentially relates to expression categories as it does in most bilaterian invertebrates, but distributions of GbM are less discretely bimodal. 2) Cnidarians generally have lower CpG frequencies on gene bodies than bilaterian invertebrates potentially suggesting a compensatory mechanism to replace CpG lost to mutation in Bilateria that is lacking in Cnidaria. 3) GbM patterns show some consistency within taxonomic groups such as the Scleractinian corals; however, GbM patterns variation across a range of taxonomic ranks in Cnidaria suggests active evolutionary change in GbM within Cnidaria. 4) Some but not all GbM variation is associated with life history change and genome expansion, whereas GbM loss is evident in endoparasitic cnidarians. 5) Cnidarian repetitive elements are less methylated than gene bodies, and methylation of both correlate with genome repeat content. 6) These observations reinforce claims that GbM evolved in stem Metazoa. Thus, this work supports overlap between DNA methylation processes in Cnidaria and Bilateria, provides a framework to compare methylation within and between Cnidaria and Bilateria, and demonstrates the previously unknown rapid evolution of cnidarian methylation.

An evolutionary genomics view on neuropeptide genes in Hydrozoa and Endocnidozoa (Myxozoa)

BACKGROUND

The animal phylum Cnidaria consists of six classes or subphyla: Hydrozoa, Scyphozoa, Cubozoa, Staurozoa, Anthozoa, and Endocnidozoa. Cnidarians have an early evolutionary origin, diverging before the emergence of the Bilateria. Extant members from this phylum, therefore, are important resources for understanding the evolution of the nervous system. Cnidarian nervous systems are strongly peptidergic. Using genomics, we have recently shown that three neuropeptide families (the X1PRX2amides, GRFamides, and GLWamides) are wide-spread in four (Scyphozoa, Cubozoa, Staurozoa, Anthozoa) out of six cnidarian classes or subphyla, suggesting that these three neuropeptide families emerged in the common cnidarian ancestor. In the current paper, we analyze the remaining cnidarian class, Hydrozoa, and the subphylum Endocnidozoa, to make firm conclusions about the evolution of neuropeptide genes in Cnidaria.

RESULTS

We analyzed sixteen hydrozoan species with a sequenced genome or transcriptome, using a recently developed software program for discovering neuropeptide genes. These species belonged to various hydrozoan subclasses and orders, among them the laboratory models Hydra, Hydractinia, and Clytia. We found that each species contained three to five neuropeptide families. A common feature for all hydrozoans was that they contained genes coding for (i) X1PRX2amide peptides, (ii) GRFamide peptides, and (iii) GLWamide peptides. These results support our previous conclusions that these three neuropeptide families evolved early in evolution. In addition to these three neuropeptide families, hydrozoans expressed up to two other neuropeptide gene families, which, however, were only occurring in certain animal groups. Endocnidozoa (Myxozoa) are microscopically small endoparasites, which are strongly reduced. For long, it was unknown to which phylum these parasites belonged, but recently they have been associated with cnidarians. We analyzed nine endocnidozoan species and found that two of them (Polypodium hydriforme and Buddenbrockia plumatellae) expressed neuropeptide genes. These genes coded for neuropeptides belonging to the GRFamide and GLWamide families with structures closely resembling them from hydrozoans.

CONCLUSIONS

We found X1PRX2amide, GRFamide, and GLWamide peptides in all species belonging to the Hydrozoa, confirming that these peptides originated in the common cnidarian ancestor. In addition, we discovered GRFamide and GLWamide peptide genes in some members of the Endocnidozoa, thereby linking these parasites to Hydrozoa.

Nematostella vectensis, an Emerging Model for Deciphering the Molecular and Cellular Mechanisms Underlying Whole-Body Regeneration

The capacity to regenerate lost or injured body parts is a widespread feature within metazoans and has intrigued scientists for centuries. One of the most extreme types of regeneration is the so-called whole body regenerative capacity, which enables regeneration of fully functional organisms from isolated body parts. While not exclusive to this habitat, whole body regeneration is widespread in aquatic/marine invertebrates. Over the past decade, new whole-body research models have emerged that complement the historical models Hydra and planarians. Among these, the sea anemone Nematostella vectensis has attracted increasing interest in regard to deciphering the cellular and molecular mechanisms underlying the whole-body regeneration process. This manuscript will present an overview of the biological features of this anthozoan cnidarian as well as the available tools and resources that have been developed by the scientific community studying Nematostella. I will further review our current understanding of the cellular and molecular mechanisms underlying whole-body regeneration in this marine organism, with emphasis on how comparing embryonic development and regeneration in the same organism provides insight into regeneration specific elements.

Morphology and molecular phylogeny of Hoferellus jutubensis n. sp. (Cnidaria: Myxozoa) found parasitising Ageneiosus inermis (Siluriformes: Auchenipteridae), in the Brazilian Amazon region

The present study describes Hoferellus jutubensis n. sp., a myxozoan parasite found in the urinary bladder of the driftwood catfish Ageneiosus inermis, captured on Jutuba Island in the state of Pará, northern Brazil. A total of 30 A. inermis specimens were examined, of which 26 (86.7%) had myxospores and polysporic plasmodia of varying shapes and sizes dispersed in the lumen of the urinary bladder, either floating freely or attached to the epithelium. In the apical view, the myxospores of Hoferellus jutubensis n. sp. are rounded, 6.1 ± 0.2 (5.7-6.3) μm long and 5.5 ± 0.3 (5.2-6.0) μm wide, with two sub-spherical polar capsules, equal in size and shape, 2.5 ± 0.2 (2.3-2.7) μm long and 1.7 ± 0.2 (1.4-2.2) μm wide. The phylogenetic analysis of a partial sequence of the SSU rDNA gene, indicated that the new species is the sister taxon of Hoferellus azevedoi, with these two species forming a Brazilian lineage of Hoferellus. The comparison of the morphological and molecular data with those of the existing members of the genus confirmed the species status of Hoferellus jutubensis n. sp., which adds one further Hoferellus taxon to the known myxosporean diversity of the Amazon basin.

Metabolic adaptations to marine environments: molecular diversity and evolution of ovothiol biosynthesis in Bacteria

Ovothiols are sulfur-containing amino acids synthesized by marine invertebrates, protozoans, and bacteria. They act as pleiotropic molecules in signaling and protection against oxidative stress. The discovery of ovothiol biosynthetic enzymes, sulfoxide synthase OvoA and β-lyase OvoB, paves the way for a systematic investigation of ovothiol distribution and molecular diversification in nature. In this work, we conducted genomic and metagenomics data mining to investigate the distribution and diversification of ovothiol biosynthetic enzymes in Bacteria. We identified the bacteria endowed with this secondary metabolic pathway, described their taxonomy, habitat and biotic interactions in order to provide insight into their adaptation to specific environments. We report that OvoA and OvoB are mostly encountered in marine aerobic Proteobacteria, some of them establishing symbiotic or parasitic relationships with other organisms. We identified a horizontal gene transfer event of OvoB from Bacteroidetes living in symbiosis with Hydrozoa. Our search within the Ocean Gene Atlas revealed the occurrence of ovothiol biosynthetic genes in Proteobacteria living in a wide range of pelagic and highly oxygenated environments. Finally, we tracked the evolutionary history of ovothiol biosynthesis from marine bacteria to unicellular eukaryotes and metazoans. Our analysis provides new conceptual elements to unravel the evolutionary and ecological significance of ovothiol biosynthesis.

Transcriptome-Wide Comparisons and Virulence Gene Polymorphisms of Host-Associated Genotypes of the Cnidarian Parasite Ceratonova shasta in Salmonids

Ceratonova shasta is an important myxozoan pathogen affecting the health of salmonid fishes in the Pacific Northwest of North America. Ceratonova shasta exists as a complex of host-specific genotypes, some with low to moderate virulence, and one that causes a profound, lethal infection in susceptible hosts. High throughput sequencing methods are powerful tools for discovering the genetic basis of these host/virulence differences, but deep sequencing of myxozoans has been challenging due to extremely fast molecular evolution of this group, yielding strongly divergent sequences that are difficult to identify, and unavoidable host contamination. We designed and optimized different bioinformatic pipelines to address these challenges. We obtained a unique set of comprehensive, host-free myxozoan RNA-seq data from C. shasta genotypes of varying virulence from different salmonid hosts. Analyses of transcriptome-wide genetic distances and maximum likelihood multigene phylogenies elucidated the evolutionary relationship between lineages and demonstrated the limited resolution of the established Internal Transcribed Spacer marker for C. shasta genotype identification, as this marker fails to differentiate between biologically distinct genotype II lineages from coho salmon and rainbow trout. We further analyzed the data sets based on polymorphisms in two gene groups related to virulence: cell migration and proteolytic enzymes including their inhibitors. The developed single-nucleotide polymorphism-calling pipeline identified polymorphisms between genotypes and demonstrated that variations in both motility and protease genes were associated with different levels of virulence of C. shasta in its salmonid hosts. The prospective use of proteolytic enzymes as promising candidates for targeted interventions against myxozoans in aquaculture is discussed. We developed host-free transcriptomes of a myxozoan model organism from strains that exhibited different degrees of virulence, as a unique source of data that will foster functional gene analyses and serve as a base for the development of potential therapeutics for efficient control of these parasites.

Morphological and molecular characterization of a novel Myxobolus species from the gastrointestinal tract of brown trout (Salmo trutta) in Spain

The genus Myxobolus Bütschli, 1882 is the largest group within the class Myxosporea and includes 905 nominal species, 18 of which have been found to infect fish belonging to the family Salmonidae. In the present study, microscopic analysis enabled detection of myxospores in 43 of 613 (7.0%) gastrointestinal tracts from brown trout (Salmo trutta) captured in several rivers in the northwest of the Iberian Peninsula. Measurement of the whole myxospores, polar capsules and other morphological characteristics, together with identification of the site of infection, has led us to propose a novel salmonid-myxobolid species, Myxobolus compostellanus n. sp. Molecular analysis of the small subunit ribosomal RNA (SSU-rRNA) gene yielded the same consensus sequence of 2039 bp in 14 fish specimens. A BLAST search indicated 97.6% similarity to Myxobolus neurobius. Phylogenetic analysis revealed that M. compostellanus n. sp. is clustered with other salmonid-infecting myxobolids. The present findings contribute to the existing knowledge about the genus Myxobolus, providing both morphological and molecular data on a novel species of Myxobolus found to infect the gastrointestinal tract of salmonids, M. compostellanus n. sp. in the brown trout (S. trutta).

Myxosporea (Myxozoa, Cnidaria) Lack DNA Cytosine Methylation

DNA cytosine methylation is central to many biological processes, including regulation of gene expression, cellular differentiation, and development. This DNA modification is conserved across animals, having been found in representatives of sponges, ctenophores, cnidarians, and bilaterians, and with very few known instances of secondary loss in animals. Myxozoans are a group of microscopic, obligate endoparasitic cnidarians that have lost many genes over the course of their evolution from free-living ancestors. Here, we investigated the evolution of the key enzymes involved in DNA cytosine methylation in 29 cnidarians and found that these enzymes were lost in an ancestor of Myxosporea (the most speciose class of Myxozoa). Additionally, using whole-genome bisulfite sequencing, we confirmed that the genomes of two distant species of myxosporeans, Ceratonova shasta and Henneguya salminicola, completely lack DNA cytosine methylation. Our results add a notable and novel taxonomic group, the Myxosporea, to the very short list of animal taxa lacking DNA cytosine methylation, further illuminating the complex evolutionary history of this epigenetic regulatory mechanism.

Relaxed selection underlies genome erosion in socially parasitic ant species

Inquiline ants are highly specialized and obligate social parasites that infiltrate and exploit colonies of closely related species. They have evolved many times convergently, are often evolutionarily young lineages, and are almost invariably rare. Focusing on the leaf-cutting ant genus Acromyrmex, we compared genomes of three inquiline social parasites with their free-living, closely-related hosts. The social parasite genomes show distinct signatures of erosion compared to the host lineages, as a consequence of relaxed selective constraints on traits associated with cooperative ant colony life and of inquilines having very small effective population sizes. We find parallel gene losses, particularly in olfactory receptors, consistent with inquiline species having highly reduced social behavioral repertoires. Many of the genomic changes that we uncover resemble those observed in the genomes of obligate non-social parasites and intracellular endosymbionts that branched off into highly specialized, host-dependent niches.

Recruitment of toxin-like proteins with ancestral venom function supports endoparasitic lifestyles of Myxozoa

Cnidarians are the oldest lineage of venomous animals and use nematocysts to discharge toxins. Whether venom toxins have been recruited to support parasitic lifestyles in the Endocnidozoa (Myxozoa + Polypodium) is, however, unknown. To examine this issue we variously employed transcriptomic, proteomic, associated molecular phylogenies, and localisation studies on representative primitive and derived myxozoans (Malacosporea and Myxosporea, respectively), Polypodium hydriforme, and the free-living staurozoan Calvadosia cruxmelitensis. Our transcriptomics and proteomics analyses provide evidence for expression and translation of venom toxin homologs in myxozoans. Phylogenetic placement of Kunitz type serine protease inhibitors and phospholipase A2 enzymes reveals modification of toxins inherited from ancestral free-living cnidarian toxins, and that venom diversity is reduced in myxozoans concordant with their reduced genome sizes. Various phylogenetic analyses of the Kunitz-type toxin family in Endocnidozoa suggested lineage-specific gene duplications, which offers a possible mechanism for enhancing toxin diversification. Toxin localisation in the malacosporean Buddenbrockia plumatellae substantiates toxin translation and thus illustrates a repurposing of toxin function for endoparasite development and interactions with hosts, rather than for prey capture or defence. Whether myxozoan venom candidates are expressed in transmission stages (e.g. in nematocysts or secretory vesicles) requires further investigation.

Phylogenetic and Selection Analysis of an Expanded Family of Putatively Pore-Forming Jellyfish Toxins (Cnidaria: Medusozoa)

Many jellyfish species are known to cause a painful sting, but box jellyfish (class Cubozoa) are a well-known danger to humans due to exceptionally potent venoms. Cubozoan toxicity has been attributed to the presence and abundance of cnidarian-specific pore-forming toxins called jellyfish toxins (JFTs), which are highly hemolytic and cardiotoxic. However, JFTs have also been found in other cnidarians outside of Cubozoa, and no comprehensive analysis of their phylogenetic distribution has been conducted to date. Here, we present a thorough annotation of JFTs from 147 cnidarian transcriptomes and document 111 novel putative JFTs from over 20 species within Medusozoa. Phylogenetic analyses show that JFTs form two distinct clades, which we call JFT-1 and JFT-2. JFT-1 includes all known potent cubozoan toxins, as well as hydrozoan and scyphozoan representatives, some of which were derived from medically relevant species. JFT-2 contains primarily uncharacterized JFTs. Although our analyses detected broad purifying selection across JFTs, we found that a subset of cubozoan JFT-1 sequences are influenced by gene-wide episodic positive selection compared with homologous toxins from other taxonomic groups. This suggests that duplication followed by neofunctionalization or subfunctionalization as a potential mechanism for the highly potent venom in cubozoans. Additionally, published RNA-seq data from several medusozoan species indicate that JFTs are differentially expressed, spatially and temporally, between functionally distinct tissues. Overall, our findings suggest a complex evolutionary history of JFTs involving duplication and selection that may have led to functional diversification, including variability in toxin potency and specificity.

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

Neural versus alternative integrative systems: molecular insights into origins of neurotransmitters

Transmitter signalling is the universal chemical language of any nervous system, but little is known about its early evolution. Here, we summarize data about the distribution and functions of neurotransmitter systems in basal metazoans as well as outline hypotheses of their origins. We explore the scenario that neurons arose from genetically different populations of secretory cells capable of volume chemical transmission and integration of behaviours without canonical synapses. The closest representation of this primordial organization is currently found in Placozoa, disk-like animals with the simplest known cell composition but complex behaviours. We propose that injury-related signalling was the evolutionary predecessor for integrative functions of early transmitters such as nitric oxide, ATP, protons, glutamate and small peptides. By contrast, acetylcholine, dopamine, noradrenaline, octopamine, serotonin and histamine were recruited as canonical neurotransmitters relatively later in animal evolution, only in bilaterians. Ligand-gated ion channels often preceded the establishment of novel neurotransmitter systems. Moreover, lineage-specific diversification of neurotransmitter receptors occurred in parallel within Cnidaria and several bilaterian lineages, including acoels. In summary, ancestral diversification of secretory signal molecules provides unique chemical microenvironments for behaviour-driven innovations that pave the way to complex brain functions and elementary cognition. This article is part of the theme issue 'Basal cognition: multicellularity, neurons and the cognitive lens'.

The myxozoan minicollagen gene repertoire was not simplified by the parasitic lifestyle: computational identification of a novel myxozoan minicollagen gene

BACKGROUND

Lineage-specific gene expansions represent one of the driving forces in the evolutionary dynamics of unique phylum traits. Myxozoa, a cnidarian subphylum of obligate parasites, are evolutionarily altered and highly reduced organisms with a simple body plan including cnidarian-specific organelles and polar capsules (a type of nematocyst). Minicollagens, a group of structural proteins, are prominent constituents of nematocysts linking Myxozoa and Cnidaria. Despite recent advances in the identification of minicollagens in Myxozoa, the evolutionary history and diversity of minicollagens in Myxozoa and Cnidaria remain elusive.

RESULTS

We generated new transcriptomes of two myxozoan species using a novel pipeline for filtering of closely related contaminant species in RNA-seq data. Mining of our transcriptomes and published omics data confirmed the existence of myxozoan Ncol-4, reported only once previously, and revealed a novel noncanonical minicollagen, Ncol-5, which is exclusive to Myxozoa. Phylogenetic analyses support a close relationship between myxozoan Ncol-1-3 with minicollagens of Polypodium hydriforme, but suggest independent evolution in the case of the myxozoan minicollagens Ncol-4 and Ncol-5. Additional genome- and transcriptome-wide searches of cnidarian minicollagens expanded the dataset to better clarify the evolutionary trajectories of minicollagen.

CONCLUSIONS

The development of a new approach for the handling of next-generation data contaminated by closely related species represents a useful tool for future applications beyond the field of myxozoan research. This data processing pipeline allowed us to expand the dataset and study the evolution and diversity of minicollagen genes in Myxozoa and Cnidaria. We identified a novel type of minicollagen in Myxozoa (Ncol-5). We suggest that the large number of minicollagen paralogs in some cnidarians is a result of several recent large gene multiplication events. We revealed close juxtaposition of minicollagens Ncol-1 and Ncol-4 in myxozoan genomes, suggesting their common evolutionary history. The unique gene structure of myxozoan Ncol-5 suggests a specific function in the myxozoan polar capsule or tubule. Despite the fact that myxozoans possess only one type of nematocyst, their gene repertoire is similar to those of other cnidarians.

Ellipsomyxa tucujuensis n. sp. (Myxozoa: Ceratomyxidae), a parasite of Satanoperca jurupari (Osteichthyes: Cichlidae) from the Brazilian Amazon

The present study describes a new coelozoic, eukaryotic microparasite of the genus Ellipsomyxa Køie, 2003 (Ceratomyxidae: Myxozoa) found parasitizing the gallbladder of Satanoperca jurupari Heckel, 1840 collected in the Curiaú River Environmental Protection Area in Macapá, Amapá state, Brazil. The fish were collected using mesh cast net. The gallbladders were examined, preserved in 80% alcohol for molecular analysis (SSU rDNA gene), and fixed in Davidson for histological slide preparation. The new parasite had a prevalence of 81% in the gallbladder, asymmetric plasmodia, irregular free spores in the bladder fluid, with no cyst formation. The spores are elliptical, with characteristics of the genus Ellipsomyxa, and they had a mean length of 10.11 (8.56-10.5) μm, mean width of 7.81 (5.96-9.56) μm, and thick walls. The polar capsules are sub-spherical in shape, slightly asymmetrical, with a mean length of 3.12 (2.31-3.99) μm and mean width of 2.5 (2.22-2.95) μm, containing polar filament with five or six coils perpendicular to the longitudinal axis of the capsule. The Bayesian Inference assigned the new species to a subclade formed by a lineage of Ellipsomyxa species from the Amazon region. Ellipsomyxa tucujuensis n. sp. is the sixth species of this genus described in fish from the Amazon region, and the first for the state of Amapá.

Gene Duplication and Gain in the Trematode Atriophallophorus winterbourni Contributes to Adaptation to Parasitism

Gene duplications and novel genes have been shown to play a major role in helminth adaptation to a parasitic lifestyle because they provide the novelty necessary for adaptation to a changing environment, such as living in multiple hosts. Here we present the de novo sequenced and annotated genome of the parasitic trematode Atriophallophorus winterbourni and its comparative genomic analysis to other major parasitic trematodes. First, we reconstructed the species phylogeny, and dated the split of A. winterbourni from the Opisthorchiata suborder to approximately 237.4 Ma (±120.4 Myr). We then addressed the question of which expanded gene families and gained genes are potentially involved in adaptation to parasitism. To do this, we used hierarchical orthologous groups to reconstruct three ancestral genomes on the phylogeny leading to A. winterbourni and performed a GO (Gene Ontology) enrichment analysis of the gene composition of each ancestral genome, allowing us to characterize the subsequent genomic changes. Out of the 11,499 genes in the A. winterbourni genome, as much as 24% have arisen through duplication events since the speciation of A. winterbourni from the Opisthorchiata, and as much as 31.9% appear to be novel, that is, newly acquired. We found 13 gene families in A. winterbourni to have had more than ten genes arising through these recent duplications; all of which have functions potentially relating to host behavioral manipulation, host tissue penetration, and hiding from host immunity through antigen presentation. We identified several families with genes evolving under positive selection. Our results provide a valuable resource for future studies on the genomic basis of adaptation to parasitism and point to specific candidate genes putatively involved in antagonistic host-parasite adaptation.

The origin of animals: an ancestral reconstruction of the unicellular-to-multicellular transition

How animals evolved from a single-celled ancestor, transitioning from a unicellular lifestyle to a coordinated multicellular entity, remains a fascinating question. Key events in this transition involved the emergence of processes related to cell adhesion, cell–cell communication and gene regulation. To understand how these capacities evolved, we need to reconstruct the features of both the last common multicellular ancestor of animals and the last unicellular ancestor of animals. In this review, we summarize recent advances in the characterization of these ancestors, inferred by comparative genomic analyses between the earliest branching animals and those radiating later, and between animals and their closest unicellular relatives. We also provide an updated hypothesis regarding the transition to animal multicellularity, which was likely gradual and involved the use of gene regulatory mechanisms in the emergence of early developmental and morphogenetic plans. Finally, we discuss some new avenues of research that will complement these studies in the coming years.

Development of Fish Parasite Vaccines in the OMICs Era: Progress and Opportunities

Globally, parasites are increasingly being recognized as catastrophic agents in both aquaculture sector and in the wild aquatic habitats leading to an estimated annual loss between 1.05 billion and 9.58 billion USD. The currently available therapeutic and control measures are accompanied by many limitations. Hence, vaccines are recommended as the "only green and effective solution" to address these concerns and protect fish from pathogens. However, vaccine development warrants a better understanding of host-parasite interaction and parasite biology. Currently, only one commercial parasite vaccine is available against the ectoparasite sea lice. Additionally, only a few trials have reported potential vaccine candidates against endoparasites. Transcriptome, genome, and proteomic data at present are available only for a limited number of aquatic parasites. Omics-based interventions can be significant in the identification of suitable vaccine candidates, finally leading to the development of multivalent vaccines for significant protection against parasitic infections in fish. The present review highlights the progress in the immunobiology of pathogenic parasites and the prospects of vaccine development. Finally, an approach for developing a multivalent vaccine for parasitic diseases is presented. Data sources to prepare this review included Pubmed, google scholar, official reports, and websites.

Advances and Discoveries in Myxozoan Genomics

Myxozoans are highly diverse and globally distributed cnidarian endoparasites in freshwater and marine habitats. They have adopted a heteroxenous life cycle, including invertebrate and fish hosts, and have been associated with diseases in aquaculture and wild fish stocks. Despite their importance, genomic resources of myxozoans have proven difficult to obtain due to their miniaturized and derived genome character and close associations with fish tissues. The first 'omic' datasets have now become the main resource for a better understanding of host-parasite interactions, virulence, and diversity, but also the evolutionary history of myxozoans. In this review, we discuss recent genomic advances in the field and outline outstanding questions to be answered with continuous and improved efforts of generating myxozoan genomic data.

Snails riding mantis shrimps: Ectoparasites evolved from ancestors living as commensals on the host's burrow wall

The molluscan class Gastropoda includes over 5,000 parasitic species whose evolutionary origins remain poorly understood. Marine snails of the genus Caledoniella (Caledoniellidae) are obligate parasites that live on the abdominal surface of the gonodactylid mantis shrimps. They have highly modified morphological characteristics specialized to the ectoparasitic lifestyle that make it difficult to infer their close relatives, thereby posing a question about their current systematic position in the superfamily Vanikoroidea. In the present study, we performed molecular phylogenetic analyses using three nuclear and three mitochondrial gene sequences to unveil the phylogenetic position of these enigmatic snails. The resulting trees recovered Caledoniella in the superfamily Truncatelloidea and within a subclade of commensal species that live on the burrow wall of marine benthic invertebrates. More specifically, Caledoniella formed the sister clade to a commensal snail species living in mantis-shrimp burrows and they collectively were sister to Sigaretornus planus (formerly in the family Tornidae or Vitrinellidae), a commensal living in echiuran burrows. This topology suggests that the species of Caledoniella achieved their ectoparasitic mode of life through the following evolutionary pathway: (1) invasion into the burrows of benthic invertebrates, (2) specialization to mantis shrimps, and (3) colonization of the host body surface from the host burrow wall with the evolution of the parasitic nature. The final step is likely to have been accompanied by the acquisition of a sucker on the metapodium, the loss of the radula and operculum, and the formation of monogamous pair bonds. The present molecular phylogeny also suggested parallel evolution of planispiral shells in a subclade of Truncatelloidea and enabled us to newly redefine the families Caledoniellidae, Elachisinidae, Teinostomatidae, Tornidae and Vitrinellidae.

Morphological Characters Can Strongly Influence Early Animal Relationships Inferred from Phylogenomic Data Sets

There are considerable phylogenetic incongruencies between morphological and phylogenomic data for the deep evolution of animals. This has contributed to a heated debate over the earliest-branching lineage of the animal kingdom: the sister to all other Metazoa (SOM). Here, we use published phylogenomic data sets ($\sim $45,000-400,000 characters in size with $\sim $15-100 taxa) that focus on early metazoan phylogeny to evaluate the impact of incorporating morphological data sets ($\sim $15-275 characters). We additionally use small exemplar data sets to quantify how increased taxon sampling can help stabilize phylogenetic inferences. We apply a plethora of common methods, that is, likelihood models and their "equivalent" under parsimony: character weighting schemes. Our results are at odds with the typical view of phylogenomics, that is, that genomic-scale data sets will swamp out inferences from morphological data. Instead, weighting morphological data 2-10$\times $ in both likelihood and parsimony can in some cases "flip" which phylum is inferred to be the SOM. This typically results in the molecular hypothesis of Ctenophora as the SOM flipping to Porifera (or occasionally Placozoa). However, greater taxon sampling improves phylogenetic stability, with some of the larger molecular data sets ($>$200,000 characters and up to $\sim $100 taxa) showing node stability even with $\geqq100\times $ upweighting of morphological data. Accordingly, our analyses have three strong messages. 1) The assumption that genomic data will automatically "swamp out" morphological data is not always true for the SOM question. Morphological data have a strong influence in our analyses of combined data sets, even when outnumbered thousands of times by molecular data. Morphology therefore should not be counted out a priori. 2) We here quantify for the first time how the stability of the SOM node improves for several genomic data sets when the taxon sampling is increased. 3) The patterns of "flipping points" (i.e., the weighting of morphological data it takes to change the inferred SOM) carry information about the phylogenetic stability of matrices. The weighting space is an innovative way to assess comparability of data sets that could be developed into a new sensitivity analysis tool. [Metazoa; Morphology; Phylogenomics; Weighting.].