The secondary structure of the unusually long 18S ribosomal RNA of the myxozoan Sphaerospora truttae and structural evolutionary trends in the Myxozoa

https://doi.org/10.1016/j.ijpara.2007.03.014Get rights and content

Abstract

The nearly complete 18S rRNA sequence of the myxozoan parasite Sphaerospora truttae shows an extraordinary length (2,552 bp) in comparison with other myxozoans and with metazoans in general (average 1,800–1,900 bp). The sequence shows nucleotide insertions in most variable regions of the 18S rRNA (V2, V4, V5 and V7), with especially large expansion segments in V4 and V7. In the myxozoans, nucleotide insertions and specific secondary structures in these regions of the gene were found to be strongly related to large scale phylogenetic clustering and thus with the invertebrate host type. Whereas expansion segments were generally found to be absent in the malacasporeans and the clade of primary marine myxozoan species, they occur in all taxa of the clade containing freshwater species, where they showed a consistent secondary structure throughout. The longest expansion segments occur in S. truttae, Sphaerospora elegans and Leptotheca ranae, which represent a clade that has emerged after the malacosporeans and before the radiation of all other myxozoan genera. These three species demonstrate structural links to the malacosporeans as well as other unique features. A smaller number of nucleotide insertions in different subhelices and specific secondary structures appear to have evolved independently in two marine genera, i.e. Ceratomyxa and Parvicapsula. The secondary structural elements of V4 and V7 of the myxozoan 18S rRNAs were found to be highly informative and revealed evolutionary trends of various regions of the gene hitherto unknown, since previous analyses have been based on primary sequence data excluding these regions. Furthermore, the unique features of the V4 region in S. truttae allowed for the design of a highly specific PCR assay for this species.

Introduction

The myxozoan Sphaerospora truttae Fischer-Scherl, El-Matbouli & Hoffmann 1986 is a significant pathogen in some Scottish freshwater salmon hatcheries (McGeorge et al., 1997, Özer and Wootten, 2000, Holzer et al., 2003). Study of the 18S rRNA gene region of S. truttae revealed its unusual length (Holzer et al., 2003) and resulted in phylogenetic positioning of the species in a clade of marine myxozoans (Holzer et al., 2004). Recently, Fiala (2006) conducted a comprehensive phylogenetic analysis of a considerably larger dataset of myxozoan 18S rRNA sequences which allocated S. truttae, with high bootstrap support, to a separate clade, basal to the marine myxozoans. The outcome of phylogenetic analyses is highly dependent on the alignment of primary sequence data, and rRNAs are often hard to align due to frequent base changes and extensive insertions in the so-called variable regions (Neefs and De Wachter, 1990). This characteristic of the rRNA suggests that, for phylogenetic reconstruction, these molecules should benefit from a higher order structure of rRNA sequences, i.e. their secondary and, if possible, tertiary structure. The general secondary structure model of the 18S rRNA in eukaryotes (Fig. 1) contains six variable regions whose secondary structure can aid considerably in alignment and improve the reliability of subsequent phylogenetic analyses (e.g. Kjer, 1995, Kjer, 2004, Titus and Frost, 1996). Robust phylogenies are especially difficult in the phylum Myxozoa, which shows extremely high sequence variability between most species, probably due to the “fast clock” characteristic (Kent et al., 1996) of the myxozoan 18S rRNA chronometer. Thus, information on 18S rRNA secondary structure might be of special importance in determining myxozoan phylogeny. However, even if primary sequence data based on secondary structure cannot be compared due to large differences between some species, common structural evolutionary changes in the gene might reveal phylogenetic relationships between these taxa. Taxon-specific rRNA morphologies have been reported from various groups of the animal kingdom and can be phylogenetically indicative (e.g. Crease and Taylor, 1998, Lydeard et al., 2000, Gillespie et al., 2005, Misof et al., 2006). The current study aimed to identify regions in the 18S rRNA of S. truttae and other myxozoans that show frequent base changes and nucleotide insertions, to determine their secondary structure and to compare these structural elements in order to reconstruct the evolution of the 18S rRNA structure in the phylum. This could then be used to confirm the phylogenetic positioning of taxa based on primary sequence data only and to optimise future alignment procedures. This is of particular importance for the positioning of unusual sequences such as that of S. truttae.

Furthermore, the identification of unique inserts in the 18S rRNA sequence of S. truttae has allowed the development of a specific PCR assay, which enables early detection of low levels of the parasite in the fish host and offers the possibility of identifying extrapiscine life cycle stages and thus the alternate host of S. truttae, which is so far unknown.

Section snippets

DNA extraction, cloning and sequencing

In order to avoid inaccuracy in the primary sequence of S. truttae, the 18S rDNA gene of S. truttae was re-sequenced from a cloned insert. A previously obtained sample of DNA from isolated S. truttae spores (Holzer et al., 2003) was used. These spores were derived from infected 1-year-old Atlantic salmon reared in a freshwater hatchery on the north coast of Scotland, which shows annually recurring S. truttae infections with 100% prevalence in the farmed fish and which was used in previous

18S rRNA primary and secondary structure analysis and phylogeny

The nearly complete sequence (10 bp missing) of the cloned 18S rDNA of S. truttae is 2,552 bp long and was submitted to GenBank under the Accession No. AM410773. It was slightly longer than the sequence previously obtained for S. truttae without cloning (AJ581915) but was otherwise identical.

The alignment on which the phylogenetic analysis was based consisted of 1,584 nucleotide positions. A total of 968 ambiguously aligned positions located in V4 and V7 were excluded from the analysis. The MP

Discussion

Molecular studies of the 18S rDNA of the myxozoan S. truttae previously showed that the mature gene is of extraordinary length (Holzer et al., 2003, Holzer et al., 2004). In the current study, the 18S rDNA was sequenced from a cloned PCR product with an almost complete length of 2,552 bp. It is thus considerably larger than the 18S rRNA of other myxozoan species sequenced to date, with the exception of S. elegans and L. ranae whose sequences have similarly long inserts and secondary structure

Acknowledgements

This study was funded by a doctoral fellowship to A.S. Holzer by the Austrian Academy of Sciences and by the External Facilities of the University of Stirling. We thank M. Jirků for providing the L. ranae sequence before its official publication date.

References (52)

  • C.L. Anderson et al.

    A triploblast origin for Myxozoa?

    Nature

    (1998)
  • K.B. Andree et al.

    Relationships among members of the genus Myxobolus, Myxozoa: Bilvalvidae based on small subunit ribosomal DNA sequences

    J. Parasitol.

    (1999)
  • J.R. Barta et al.

    Phylogenetic relationships among eight Eimeria species infecting domestic fowl inferred using complete small subunit ribosomal DNA sequences

    J. Parasitol.

    (1997)
  • J.L. Bartholomew et al.

    The life cycle of Ceratomyxa shasta, a myxosporean parasite of salmonids, requires a freshwater polychaete as an alternate host

    J. Parasitol.

    (1997)
  • J.L. Bartholomew et al.

    Involvement of Manayunkia speciosa (Annelida: Polychaeta: Sabellidae) in the life cycle of Parvicapsula minibicornis, a myxozoan parasite of Pacific salmon

    J. Parasitol.

    (2006)
  • B.C. Campbell et al.

    Paraphyly of Homoptera and Auchenorrhyncha inferred from 18S rDNA nucleotide sequences

    Syst. Entomol.

    (1995)
  • E.U. Canning et al.

    Ecology, development and pathogenicity of Buddenbrockia plumatellae Schröder, 1910 (Myxozoa, Malacosporea) (syn. Tetracapsula bryozoides) and establishment of Tetracapsuloides n. gen. for Tetracapsula bryosalmonae

    J. Eukaryot. Microbiol.

    (2002)
  • N. Chalwatzis et al.

    Strongly expanded 18S rRNA genes correlated with a peculiar morphology in the insect order of Strepsiptera

    Zoology

    (1995)
  • C.P. Choe et al.

    Analysis of the primary sequence and secondary structure of the unusually long SSU rRNA of the soil bug, Armadillidium vulgare

    J. Mol. Evol.

    (1999)
  • T.J. Crease et al.

    The unusually long small-subunit ribosomal RNA of the crustacean Daphnia pulex: sequence and predicted secondary structure

    J. Mol. Evol.

    (1998)
  • T.J. Crease et al.

    The origin and evolution of variable-region helices in V4 and V7 of the small-subunit ribosomal RNA of branchipod crustaceans

    Mol. Biol. Evol.

    (1998)
  • P. De Rijk et al.

    DCSE, an interactive tool for sequence alignment and secondary structure research

    Comput. Appl. Biosci.

    (1993)
  • P. De Rijk et al.

    RnaViz 2: an improved representation of RNA secondary structure

    Bioinformatics

    (2003)
  • E. Eszterbauer

    Genetic relationship among gill-infecting Myxobolus species (Myxosporea) of cyprinids: molecular evidence of importance of tissue-specificity

    Dis. Aquat. Org.

    (2004)
  • S.W. Feist et al.

    Proliferative kidney disease and renal myxosporidiosis in juvenile salmonids from rivers in England and Wales

    J. Fish Dis.

    (2002)
  • I. Fiala et al.

    The phylogeny of marine and freshwater species of the genus Chloromyxum Mingazzini, 1910 (Myxosporea: Bivalvulida) based on small subunit ribosomal RNA gene sequences

    Folia Parasitol.

    (2004)
  • Cited by (68)

    • The SSU rRNA secondary structures of the Plagiorchiida species (Digenea), its applications in systematics and evolutionary inferences

      2020, Infection, Genetics and Evolution
      Citation Excerpt :

      Indeed, the interaction sites of two eukaryotic ribosomal subunits on h44 (the bridge between rpS17e (40S) and rpS24e (60S)) are classified as areas with a low degree of variability. It is also known that ES12 adjacent to h44 in different animal species are relatively variable in nucleotide sequences, but with the common pattern of the secondary structures (Johnston, 2006; Armache et al., 2010; Mallatt et al., 2010). In this regard, it is extremely interesting to find out why there are some variations in the ES12-h44 secondary structure.

    • Novel Myxobolus and Ellipsomyxa species (Cnidaria: Myxozoa) parasiting Brachyplatystoma rousseauxii (Siluriformes: Pimelodidae) in the Amazon basin, Brazil

      2018, Parasitology International
      Citation Excerpt :

      However, as more myxozoans have been described and sequenced, the division of species by host habitat has been shown to be inaccurate, with freshwater species falling in the “marine” clade and vice versa. Growing evidence supports the alternate hypothesis, that these two lineages may be better defined by their definitive, invertebrate hosts, with “marine” myxozoans having polychaete hosts and “freshwater” having oligochaete hosts [55, 59]. Our present work illustrates the problem of the view of “marine” and “freshwater” lineages: although Ellipsomyxa species and closely-related myxozoans all fall within the “marine” lineage, the actual taxa in these clades infect marine, estuarine or freshwater fish hosts (see Fig. 5).

    View all citing articles on Scopus

    Nucleotide sequence data reported in this paper is available in the GenBank[TM], EMBL and DDBJ databases under the Accession No. AM410773.

    View full text