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.

Centro de Investigación Aplicada y Transferencia Tecnológica en Recursos Marinos Almirante Storni CCT CONICET CENPAT San Antonio Oeste Argentina

Department of Integrative Biology Oregon State University

Department of Microbiology Oregon State University

Department of Molecular Biology and Genetics Faculty of Science University of South Bohemia České Budějovice Czech Republic

Institute of Parasitology Biology Centre of the Czech Academy of Sciences České Budějovice Czech Republic

Zobrazit více v PubMed

Alama-Bermejo G, Holzer AS, Bartholomew JL.. 2019. Myxozoan adhesion and virulence: Ceratonova shasta on the move. Microorganisms 7(10):397. PubMed PMC

Ali IKM, et al.2007. Evidence for a link between parasite genotype and outcome of infection with Entamoeba histolytica. J Clin Microbiol. 45(2):285–289. PubMed PMC

Atkinson SD, Bartholomew JL.. 2010. a. Disparate infection patterns of Ceratomyxa shasta (Myxozoa) in rainbow trout Oncorhynchus mykiss and Chinook salmon Oncorhynchus tshawytscha correlate with Internal Transcribed Spacer-1 sequence variation in the parasite. Int J Parasitol. 40(5):599–604. PubMed

Atkinson SD, Bartholomew JL.. 2010. b. Spatial, temporal and host factors structure the Ceratomyxa shasta (Myxozoa) population in the Klamath River basin. Infect Genet Evol. 10(7):1019–1026. PubMed

Atkinson SD, Bartholomew JL, Lotan T.. 2018. Myxozoans: ancient metazoan parasites find a home in phylum Cnidaria. Zoology 129:66–68. PubMed

Atkinson SD, Hallett SL, Bartholomew JL.. 2018. Genotyping of individual Ceratonova shasta (Cnidaria: Myxosporea) myxospores reveals intra-spore ITS-1 variation and invalidates the distinction of genotypes II and III. Parasitology 145(12):1588–1586. PubMed

Atkinson SD, Jones SR, Adlard RD, Bartholomew JL.. 2011. Geographical and host distribution patterns of Parvicapsula minibicornis (Myxozoa) small subunit ribosomal RNA genetic types. Parasitology 138(8):969–977. PubMed

Bane KS, et al.2016. The actin filament-binding protein coronin regulates motility in Plasmodium sporozoites. PLoS Pathog. 12(7):e1005710. PubMed PMC

Barragan A, Sibley LD.. 2002. Transepithelial migration of Toxoplasma gondii is linked to parasite motility and virulence. J Exp Med. 195(12):1625–1633. PubMed PMC

Bartholomew JL, Whipple MJ, Stevens DG, Fryer JL.. 1997. The life cycle of Ceratomyxa shasta, a myxosporean parasite of salmonids, requires a freshwater polychaete as an alternate host. J Parasitol. 83(5):859–868. PubMed

Berthelot C, et al.2014. The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates. Nat Commun. 5:3657. PubMed PMC

Bjork SJ, Bartholomew JL.. 2009. a. Effects of Ceratomyxa shasta dose on a susceptible strain of rainbow trout and comparatively resistant Chinook and coho salmon. Dis Aquat Org. 86:29–37. PubMed

Bjork SJ, Bartholomew JL.. 2009. b. The effects of water velocity on the Ceratomyxa shasta infectious cycle. J Fish Dis. 32(2):131–142. PubMed PMC

Bjork SJ, Bartholomew JL.. 2010. Invasion of Ceratomyxa shasta (Myxozoa) and comparison of migration to the intestine between susceptible and resistant fish hosts. Int J Parasitol. 40(9):1087–1095. PubMed

Bookwalter CS, et al.2017. Reconstitution of the core of the malaria parasite glideosome with recombinant Plasmodium class XIV myosin A and Plasmodium actin. J Biol Chem. 292(47):19290–19303. PubMed PMC

Bouzid M, Hunter PR, Chalmers RM, Tyler KM.. 2013. Cryptosporidium pathogenicity and virulence. Clin Microbiol Rev. 26(1):115–134. PubMed PMC

Bowcock AM, et al.1994. High resolution of human evolutionary trees with polymorphic microsatellites. Nature 368(6470):455–457. PubMed

Cardoso J, et al.2011. Analysis of proteasomal proteolysis during the in vitro metacyclogenesis of Trypanosoma cruzi. PLoS One 6(6):e21027. PubMed PMC

Casadevall A, Pirofski L.. 1999. Host–pathogen interactions: redefining the basic concepts of virulence and pathogenicity. Infect Immun. 67(8):3703–3713. PubMed PMC

Chang E. 2013. Transcriptomic evidence that enigmatic parasites Polypodium hydriforme and Myxozoa are cnidarians [dissertation/PhD’s thesis]. [Kansas (KS)]: University of Kansas.

Chang ES, et al.2015. Genomic insights into the evolutionary origin of Myxozoa within Cnidaria. Proc Natl Acad Sci U S A. 112(48):14912–14917. PubMed PMC

Chen X, et al.2006. Inhibitors of Plasmodium falciparum methionine aminopeptidase 1b possess antimalarial activity. Proc Natl Acad Sci U S A. 103(39):14548–14553. PubMed PMC

Choe JE, Welch MD.. 2016. Actin-based motility of bacterial pathogens: mechanistic diversity and its impact on virulence. Pathog Dis. 74(8):ftw099. PubMed PMC

Clark KF, Greenwood SJ.. 2016. Next-generation sequencing and the crustacean immune system: the need for alternatives in immune gene annotation. Integr Comp Biol. 56(6):1113–1130. PubMed

Cogliati M, Barchiesi F, Spreghini E, Tortorano AM.. 2012. Heterozygosis and pathogenicity of Cryptococcus neoformans AD-hybrid isolates. Mycopathologia 173(5–6):347–357. PubMed

Coletta-Filho HD, Bittleston LS, Lopes JRS, Daugherty MP, Almeida RPP.. 2015. Genetic distance may underlie virulence differences among isolates of a bacterial plant pathogen. J Plant Pathol. 97:465–469.

Crête-Lafrenière A, Weir LK, Bernatchez L.. 2012. Framing the Salmonidae family phylogenetic portrait: a more complete picture from increased taxon sampling. PLoS One 7(10):e46662. PubMed PMC

Daly GM, et al.2015. Host subtraction, filtering and assembly validations for novel viral discovery using next generation sequencing data. PLoS One 10(6):e0129059. PubMed PMC

Dardé ML. 2008. Toxoplasma gondii, “new” genotypes and virulence. Parasite 15(3):366–371. PubMed

David M, Dzamba M, Lister D, Ilie L, Brudno M.. 2011. SHRiMP2: sensitive yet practical short read mapping. Bioinformatics 27(7):1011–1012. PubMed

Dittman AH, et al.2010. Homing and spawning site selection by supplemented hatchery- and natural-origin Yakima River spring Chinook salmon. Trans Am Fish Soc. 139(4):1014–1028.

Domanico MJ, Phillips RB, Oakley TH.. 1997. Phylogenetic analysis of Pacific salmon (genus Oncorhynchus) using nuclear and mitochondrial DNA sequences. Can J Fish Aquat Sci. 54(8):1865–1872.

Eichenberger RM, Ramakrishnan C, Russo G, Deplazes P, Hehl AB.. 2017. Genome-wide analysis of gene expression and protein secretion of Babesia canis during virulent infection identifies potential pathogenicity factors. Sci Rep. 7(1):3357. PubMed PMC

Foox J, Ringuette M, Desser SS, Siddall ME.. 2015. In silico hybridization enables transcriptomic illumination of the nature and evolution of Myxozoa. BMC Genomics. 16(1):840. PubMed PMC

González-Rodríguez VE, Garrido C, Cantoral JM, Schumacher J.. 2016. The F-actin capping protein is required for hyphal growth and full virulence but is dispensable for septum formation in Botrytis cinerea. Fungal Biol. 120(10):1225–1235. PubMed

Gupta I, Aggarwal S, Singh K, Yadav A, Khan S.. 2018. Ubiquitin Proteasome pathway proteins as potential drug targets in parasite. Sci Rep. 8(1):8399. PubMed PMC

Haas BJ, et al.2013. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 8(8):1494–1512. PubMed PMC

Hallett SL, Bartholomew JL.. 2012. Myxobolus cerebralis and Ceratomyxa shasta In: Woo PTK, Buchmann K, editors. Fish parasites: pathobiology and protection. Oxfordshire, UK: CABI. p. 141–172.

Hallett SL, et al.2012. Density of the waterborne parasite Ceratomyxa shasta and its biological effects on salmon. Appl Environ Microbiol. 78(10):3724–3731. PubMed PMC

Hartigan A, Kosakyan A, Pecková H, Eszterbauer E, Holzer A. 2020. Transcriptome of Sphaerospora molnari (Cnidaria, Myxosporea) blood stages provides proteolytic arsenal as potential therapeutic targets against sphaerosporosis in common carp. BMC Genomics 21(1):404. PubMed PMC

Hartigan A, et al.2016. New cell motility model observed in parasitic cnidarian Sphaerospora molnari (Myxozoa: Myxosporea) blood stages in fish. Sci Rep. 6:39093. PubMed PMC

Hartikainen H, et al.2016. Assessing myxozoan presence and diversity using environmental DNA. Int J Parasitol. 46(12):781–792. PubMed

Holland JW, Okamura B, Hartikainen H, Secombes CJ.. 2011. A novel minicollagen gene links cnidarians and myxozoans. Proc R Soc B Biol Sci. 278(1705):546–553. PubMed PMC

Holzer AS, et al.2018. The joint evolution of the Myxozoa and their alternate hosts: a cnidarian recipe for success and vast biodiversity. Mol Ecol. 27(7):1651–1666. PubMed

Hurst CN, Bartholomew JL.. 2012. Ceratomyxa shasta genotypes cause differential mortality in their salmonid hosts. J Fish Dis. 35(10):725–732. PubMed

Hurst CN, Holt RA, Bartholomew JL.. 2012. Dam removal and implications for fish health: Ceratomyxa shasta in the Williamson River, Oregon, USA. N Am J Fish Manage. 32(1):14–23.

Jiménez-Guri E, Philippe H, Okamura B, Holland PWH.. 2007. Buddenbrockia is a cnidarian worm. Science 317(5834):116–118. PubMed

Jones CE, Brown AL, Baumann U.. 2007. Estimating the annotation error rate of curated go database sequence annotations. BMC Bioinf. 8:170. PubMed PMC

Katoh K, Standley DM.. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 30(4):772–780. PubMed PMC

Kelley GO, Adkison MA, Leutenegger CM, Hedrick RP.. 2003. Myxobolus cerebralis: identification of a cathepsin Z-like protease gene (MyxCP-1) expressed during parasite development in rainbow trout, Oncorhynchus mykiss. Exp Parasitol. 105(3–4):201–210. PubMed

Kenkel CD, Bay LK.. 2017. Novel transcriptome resources for three scleractinian coral species from the Indo-Pacific. GigaScience 6(9):1–4. PubMed PMC

Konstantinova IM, Tsimokha AS, Mittenberg AG.. 2008. Role of proteasomes in cellular regulation. Int Rev Cell Mol Biol. 267:59–124. PubMed

Kurath G, et al.2003. Phylogeography of infectious haematopoietic necrosis virus in North America. J Gen Virol. 84(4):803–814. PubMed

Lentini G, et al.2015. Identification and characterization of Toxoplasma SIP, a conserved apicomplexan cytoskeleton protein involved in maintaining the shape, motility and virulence of the parasite. Cell Microbiol. 17(1):62–78. PubMed PMC

Li B, Dewey CN.. 2011. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf. 12:323. PubMed PMC

Li H. 2011. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27(21):2987–2993. PubMed PMC

Li H, et al.2009. The Sequence Alignment/Map (SAM) format and SAMtools. Bioinformatics 25(16):2078–2079. PubMed PMC

Lien S, et al.2016. The Atlantic salmon genome provides insights into rediploidization. Nature 533(7602):200–205. PubMed PMC

Lin MF, et al.2017. Analyses of corallimorpharian transcriptomes provide new perspectives on the evolution of calcification in the Scleractinia (Corals). Genome Biol Evol. 9(1):150–160. PubMed PMC

Manasanch EE, Orlowski RZ.. 2017. Proteasome inhibitors in cancer therapy. Nat Rev Clin Oncol. 14(7):417–433. PubMed PMC

McKerrow J, Caffrey C, Kelly B, Loke P, Sajid M.. 2006. Proteases in parasitic diseases. Annu Rev Pathol Mech Dis. 1(1):497–536. PubMed

Mennerat A, Nilsen F, Ebert D, Skorping A.. 2010. Intensive farming: evolutionary implications for parasites and pathogens. Evol Biol. 37(2–3):59–67. PubMed PMC

Meyer E, et al.2009. Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics. 10(1):219. PubMed PMC

Minakawa N, Kraft G.. 2005. Homing behaviour of juvenile coho salmon (Oncorhynchus kisutch) within an off-channel habitat. Ecol Freshw Fish. 14(2):197–201.

Morais ER, et al.2017. Effects of proteasome inhibitor MG-132 on the parasite Schistosoma mansoni. PLoS One 12(9):e0184192. PubMed PMC

Muñoz C, Francisco JS, Gutiérrez B, González J.. 2015. Role of the ubiquitin-proteasome systems in the biology and virulence of protozoan parasites. BioMed Res Int. 2015:1–13. PubMed PMC

Nesnidal MP, Helmkampf M, Bruchhaus I, El-Matbouli M, Hausdorf B.. 2013. Agent of whirling disease meets orphan worm: phylogenomic analyses firmly place Myxozoa in Cnidaria. PLoS One 8(1):e54576. PubMed PMC

Parra G, Bradnam K, Korf I.. 2007. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23(9):1061–1067. PubMed

Pearson WR. 2013. An introduction to sequence similarity (‘homology’) searching. Curr Protoc Bioniformatics Chapter 3:Unit 3. PubMed PMC

Rawlings ND, Barrett AJ, Finn RD.. 2016. Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 44(D1):D343–D350. PubMed PMC

Ray RA, Holt RA, Bartholomew JL.. 2012. Relationship between temperature and Ceratomyxa shasta-induced mortality in Klamath River salmonids. J Parasitol. 98(3):520–526. PubMed

Rodrigues M, et al.2016. Profiling of adhesive-related genes in the freshwater cnidarian Hydra magnipapillata by transcriptomics and proteomics. Biofouling 32(9):1115–1129. PubMed PMC

Schnoes AM, Brown SD, Dodevski I, Babbitt PC.. 2009. Annotation error in public databases: misannotation of molecular function in enzyme superfamilies. PLoS Comput Biol. 5(12):e1000605. PubMed PMC

Seib KL, Dougan G, Rappuoli R.. 2009. The key role of genomics in modern vaccine and drug design for emerging infectious diseases. PLoS Genet. 5(10):e1000612. PubMed PMC

Shpirer E, Diamant A, Cartwright P, Huchon D.. 2018. A genome wide survey reveals multiple nematocyst-specific genes in Myxozoa. BMC Evol Biol. 18(1):138. PubMed PMC

Stewart ZK, Pavasovic A, Hock DH, Prentis PJ.. 2017. Transcriptomic investigation of wound healing and regeneration in the cnidarian Calliactis polypus. Sci Rep. 7:41458. PubMed PMC

Stinson MET, Atkinson SD, Bartholomew JL.. 2018. Widespread distribution of Ceratonova shasta (Cnidaria: Myxosporea) genotypes indicates both evolutionary adaptation to its salmonid fish hosts. J Parasitol. 104(6):645–650. PubMed

Stinson ME, Bartholomew JL.. 2012. Predicted redistribution of Ceratomyxa shasta genotypes with salmonid passage in the Deschutes River, Oregon. J Aquat Anim Health 24(4):274–280. PubMed

Tanaka K. 2009. The proteasome: overview of structure and functions. Proc Jpn Acad Ser B 85(1):12–36. PubMed PMC

Tine M, et al.2014. European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat Commun. 5:5770. PubMed PMC

Tirloni L, et al.2014. Proteomic analysis of cattle tick Rhipicephalus (Boophilus) microplus saliva: a comparison between partially and fully engorged females. PLoS One 9(4):e94831. PubMed PMC

Tyndall JDA, Nall T, Fairlie DP.. 2005. Proteases universally recognize beta strands in their active sites. Chem Rev. 105(3):973–1000. PubMed

Varga G, et al.2017. Proteasome subunit beta type 1 P11A polymorphism is a new prognostic marker in multiple myeloma. Clin Lymphoma Myeloma Leuk. 17(11):734–742. PubMed

Williams TH, et al.2016. Viability assessment for Pacific salmon and steelhead listed under the Endangered Species Act: Southwest. California: U.S. Department of Commerce, NOAA Technical Memorandum. NMFS-SWFSC-564.

Xu P, et al.2014. Genome sequence and genetic diversity of the common carp, Cyprinus carpio. Nat Genet. 46(11):1212–1219. PubMed

Yahalomi D, et al. . 2020. A cnidarian parasite of salmon (Myxozoa: Henneguya ) lacks a mitochondrial genome. Proc Natl Acad Sci USA. 117(10):5358–5363. PubMed PMC

Yang Y, et al.2014. The genome of the myxosporean Thelohanellus kitauei shows adaptations to nutrient acquisition within its fish host. Genome Biol Evol. 6(12):3182–3198. PubMed PMC

Zardoya R, Garrido-Pertierra A, Bautista JM.. 1995. The complete nucleotide sequence of the mitochondrial DNA genome of the rainbow trout, Oncorhynchus mykiss. J Mol Evol. 41(6):942–951. PubMed

Zheng J, Jia H, Zheng Y.. 2015. Knockout of leucine aminopeptidase in Toxoplasma gondii using CRISPR/Cas9. Int J Parasitol. 45(2–3):141–148. PubMed

Comparative transcriptomics reveal stage-dependent parasitic adaptations in the myxozoan Sphaerospora molnari

Ceratonova shasta: a cnidarian parasite of annelids and salmonids

Method for Isolation of Myxozoan Proliferative Stages from Fish at High Yield and Purity: An Essential Prerequisite for In Vitro, In Vivo and Genomics-Based Research Developments

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

To React or Not to React: The Dilemma of Fish Immune Systems Facing Myxozoan Infections

Evolutionary Analysis of Cystatins of Early-Emerging Metazoans Reveals a Novel Subtype in Parasitic Cnidarians

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

Differences in inflammatory responses of rainbow trout infected by two genotypes of the myxozoan parasite Ceratonova shasta

Genetic Diversity of Serine Protease Inhibitors in Myxozoan (Cnidaria, Myxozoa) Fish Parasites

Zobrazit více v PubMed

Dryad
10.5061/dryad.tx95x69tt