Vintage helminths yield valuable molecules

Museum and archival collections of parasites are available throughout the world but, although they represent a huge diversity of species and forms, they tend to be used solely for reference to morphology, if at all. As biochemical techniques begin to overcome the problems associated with ancient, degraded and formalin-fixed tissues, the value of such collections increases. Molecular data are now available for rare, elusive and extinct species, as well as those densely sampled for epidemiological, biogeographical or clinical collections. Here, Elisabeth Herniou, Auriol Pearce and Tim Littlewood describe some of the advances and pitfalls associated with retrieving DNA from formalin-fixed helminth material and suggest just some of the new ways that parasitologists can tap into these resources.

Mitochondrial gene content, arrangement and composition compared in African and Asian schistosomes

Complete sequences were obtained for the coding portions of the mitochondrial (mt) genomes of Schistosoma mansoni (NMRI strain, Puerto Rico; 14 415 bp), S. japonicum (Anhui strain, China; 14 085 bp) and S. mekongi (Khong Island, Laos; 14 072 bp). Each comprises 36 genes: 12 protein-encoding genes (cox1-3, nad1-6, nad4L, atp6 and cob); two ribosomal RNAs, rrnL (large subunit rRNA or 16S) and rrnS (small subunit rRNA or 12S); as well as 22 transfer RNA (tRNA) genes. The atp8 gene is absent. A large segment (9.6 kb) of the coding region (comprising 14 tRNAs, eight complete and two incomplete protein-encoding genes) for S. malayensis (Baling, Malaysian Peninsula) was also obtained. Each genome also possesses a long non-coding region that is divided into two parts (a small and a large non-coding region, the latter not fully sequenced in any species) by one or more tRNAs. The protein-encoding genes are similar in size, composition and codon usage in all species except for cox1 in S. mansoni (609 aa) and cox2 in S. mekongi (219 aa), both of which are longer than homologues in other species. An unexpected finding in all the Schistosoma species was the presence of a leucine zipper motif in the nad4L gene. The gene order in S. mansoni is strikingly different from that seen in the S. japonicum group and other flatworms. There is a high level of identity (87-94% at both the nucleotide and amino acid levels) for all protein-encoding genes of S. mekongi and S. malayensis. The identity between genes of these two species and those of S. japonicum is less (56-83% for amino acids and 73-79% for nucleotides). The identity between the genes of S. mansoni and the Asian schistosomes is far less (33-66% for amino acids and 54-68% for nucleotides), an observation consistent with the known phylogenetic distance between S. mansoni and the other species.

The nature and evolution of the association among digeneans, molluscs and fishes

Patterns of association of digenean families and their mollusc and vertebrate hosts are assessed by way of a new database containing information on over 1000 species of digeneans for life-cycles and over 5000 species from fishes. Analysis of the distribution of digenean families in molluscs suggests that the group was associated primitively with gastropods and that infection of polychaetes, bivalves and scaphopods are all the results of host-switching. For the vertebrates, infections of agnathans and chondrichthyans are apparently the result of host-switching from teleosts. For digenean families the ratio of orders of fishes infected to superfamilies of molluscs infected ranges from 0.5 (Mesometridae) to 16 (Bivesiculidae) and has a mean of 5.6. Individual patterns of host association of 13 digenean families and superfamilies are reviewed. Two, Bucephalidae and Sanguinicolidae, are exceptional in infecting a range of first intermediate hosts qualitatively as broad as their range of definitive hosts. No well-studied taxon shows narrower association with vertebrate than with mollusc clades. The range of definitive hosts of digeneans is characteristically defined by eco-physiological similarity rather than phylogenetic relationship. The range of associations of digenean families with mollusc taxa is generally much narrower. These data are considered in the light of ideas about the significance of different forms of host association. If Manter's Second Rule (the longer the association with a host group, the more pronounced the specificity exhibited by the parasite group) is invoked, then the data may suggest that the Digenea first parasitised molluscs before adopting vertebrate hosts. This interpretation is consistent with most previous ideas about the evolution of the Digenea but contrary to current interpretations based on the monophyly of the Neodermata. The basis of Manter's Second Rule is, however, considered too flimsy for this interpretation to be robust. Problems of the inference of the evolution of patterns of parasitism in the Neodermata are discussed and considered so intractable that the truth may be presently unknowable.

Digenean parasites of deep-sea teleosts: a review and case studies of intrageneric phylogenies

Studies on the digenean parasites of deep-sea (> 200 m depth) teleosts are reviewed and two case study generic phylogenies are presented based on LSU rDNA and ND1 mtDNA sequences. The phylogeny of the lepocreadiid genus Lepidapedon, the most common deep-sea digenean genus, is not clearly resolved as the two gene trees are not compatible. It can be inferred, however, that the genus has radiated in the deeper waters off the continental shelf, mainly in fishes of the gadiform family Macrouridae. Steringophorus, a fellodistomid genus, is better resolved. In this case a deep-sea radiation is also indicated, but the pattern of host-specificity is not clear, with evidence of much host-switching. Results of studies of the parasites of the macrourid fish Coryphaenoides (Nematonurus) armatus from various depths have reinforced recent views on the lack of zoned depth-related communities in the deep-sea. The diversity of deep-sea digeneans is relatively low with only 18 families (of about 60) reported. Little, or nothing, is known from most deep-sea areas and nothing from trenches and mid-ocean ridge systems.

Interrelationships and evolution of the tapeworms (Platyhelminthes: Cestoda)

Interrelationships of the tapeworms (Platyhelminthes: Cestoda) were examined by use of small (SSU) and large (LSU) subunit ribosomal DNA sequences and morphological characters. Fifty new complete SSU sequences were added to 21 sequences previously determined, and 71 new LSU (D1-D3) sequences were determined for the complementary set of taxa representing each of the major lineages of cestodes as currently understood. New sequences were determined for three amphilinidean taxa, but were removed from both alignments due to their excessively high degree of divergence from other cestode sequences. A morphological character matrix coded for supraspecific taxa was constructed by the modification of matrices from recently published studies. Maximum-parsimony (MP) analyses were performed on the LSU, SSU, LSU+SSU, and morphological data partitions, and minimum-evolution (ME) analyses utilizing a general time reversible model of nucleotide substitution including estimates of among-site rate heterogeneity were performed on the molecular data partitions. Resulting topologies were rooted at the node separating the Gyrocotylidea from the Eucestoda. The LSU data were found to be more informative than the SSU data and were more consistent with inferences from morphology, although nodal support was generally weak for most basal nodes. One class of transitions was found to be saturated for comparisons between the most distantly related taxa (gyrocotylideans vs cyclophyllideans and tetrabothriideans). Differences in the topologies resulting from MP and ME analyses were not statistically significant. Nonstrobilate orders formed the basal lineages of trees resulting from analysis of LSU data and morphology. Difossate orders were basal to tetrafossate orders, the latter of which formed a strongly supported clade. A clade including the orders Cyclophyllidea, Nippotaeniidea, and Tetrabothriidea was supported by all data partitions and methods of analysis. Paraphyly of the orders Pseudophyllidea, Tetraphyllidea, and Trypanorhyncha was consistent among the molecular data partitions. Inferences are made regarding a monozoic (nonsegmented) origin of the Eucestoda as represented by the Caryophyllidea and for the evolution of the strobilate and acetabulate/tetrafossate conditions having evolved in a stepwise pattern.

Changes in mitochondrial genetic codes as phylogenetic characters: two examples from the flatworms

Shared molecular genetic characteristics other than DNA and protein sequences can provide excellent sources of phylogenetic information, particularly if they are complex and rare and are consequently unlikely to have arisen by chance convergence. We have used two such characters, arising from changes in mitochondrial genetic code, to define a clade within the Platyhelminthes (flatworms), the Rhabditophora. We have sampled 10 distinct classes within the Rhabditophora and find that all have the codon AAA coding for the amino acid Asn rather than the usual Lys and AUA for Ile rather than the usual Met. We find no evidence to support claims that the codon UAA codes for Tyr in the Platyhelminthes rather than the standard stop codon. The Rhabditophora are a very diverse group comprising the majority of the free-living turbellarian taxa and the parasitic Neodermata. In contrast, three other classes of turbellarian flatworm, the Acoela, Nemertodermatida, and Catenulida, have the standard invertebrate assignments for these codons and so are convincingly excluded from the rhabditophoran clade. We have developed a rapid computerized method for analyzing genetic codes and demonstrate the wide phylogenetic distribution of the standard invertebrate code as well as confirming already known metazoan deviations from it (ascidian, vertebrate, echinoderm/hemichordate).

The interrelationships of proseriata (Platyhelminthes: seriata) tested with molecules and morphology

Proseriate flatworms are common members of the interstitial benthic fauna worldwide, predominantly occupying marine environments. As minute animals, having relatively few characters useful for cladistic analysis, they have been difficult to present in a phylogenetic framework using morphology alone. Here we present a new morphological matrix consisting of 16 putatively homologous characters and two molecular data sets to investigate further this major group of free-living members of the Platyhelminthes. Complete 18S rDNA (representing 277 parsimony-informative characters) from 17 ingroup taxa and partial 28S rDNA spanning variable expansion regions D1 to D3 and D1 to D6 (representing 219 and 361 parsimony-informative characters, respectively) from 27 and 14 ingroup taxa, respectively, were determined and aligned as complementary data sets. Morphological and molecular data sets were analyzed separately and together to determine underlying phylogenetic patterns and to resolve conflict between published scenarios based on morphology alone. The monophyly of the Proseriata cannot be confirmed categorically with any of these data sets. However, the constituent taxa are confirmed as basal members of the Neoophora, and a sister group relationship with Tricladida is rejected. Similarly, the monophyly of one of the two subtaxa of the Proseriata, the Lithophora, could not be confirmed with molecules. Concerning intragroup relationships, we could reject one of the two phylogenetic trees formerly proposed, as well as the clade Otoplanidae + Coelogynoporidae. However, a clade Otoplanidae + Archimonocelididae + Monocelididae (to which the Monotoplanidae belong) was supported, and the position of the genus Calviria shifted from the Archimonocelididae to the Coelogynoporidae.

Acoel flatworms: earliest extant bilaterian Metazoans, not members of Platyhelminthes

Because of their simple organization the Acoela have been considered to be either primitive bilaterians or descendants of coelomates through secondary loss of derived features. Sequence data of 18S ribosomal DNA genes from non-fast evolving species of acoels and other metazoans reveal that this group does not belong to the Platyhelminthes but represents the extant members of the earliest divergent Bilateria, an interpretation that is supported by recent studies on the embryonic cleavage pattern and nervous system of acoels. This study has implications for understanding the evolution of major body plans, and for perceptions of the Cambrian evolutionary explosion.

Phylogenetic relationships of the family Campulidae (Trematoda) based on 18S rRNA sequences

Traditionally, the family Campulidae has been associated either with the family Fasciolidae, parasites of ruminants, or the Acanthocolpidae, parasites of fishes, based on morphological similarities. Since morphology does not seem to resolve clearly the problem of the relationships of campulids, we have used the sequences of the 18S rRNA gene of the campulids Zalophotrema hepaticum, Campula oblonga and Nasitrema globicephalae, the fasciolid Fasciola hepatica, the acanthocolpid Stephanostomum baccatum and the outgroup Schistosoma mansoni to infer a phylogeny. Maximum parsimony and neighbour-joining methods were applied. Both methods indicated that campulids are closer to acanthocolpids than fasciolids. In order to confirm this relationship, we generated a second phylogeny using all the partial sequences of the 18S published for trematodes: Lobatostoma manteri, Echinostoma caproni, Calicophoron calicophorum, Tetracerasta blepta, Gyliauchen sp. and Opistorchis viverrini, plus those mentioned above, and Dicrocoelium dendriticum. The aspidogastrean L. manteri was used as the outgroup. Results were identical to the first analysis. According to this and the most recent Digenean phylogeny, which considers campulids and acanthocolpids as sister groups, we suggest that a common origin for these 2 groups would imply a host-switching process. The life-cycle of acanthocolpids includes marine gastropods as first intermediate hosts, and fishes as second intermediate and definitive hosts. In this context, the hypothesis would be that trematodes whose cycle ended in fishes were able to switch to mammalian hosts.

The phylogenetic position of Udonella (Platyhelminthes)

Phylogenetic analysis of molecular data from complete 18S rRNA and partial 28S rRNA genes, of a variety of platyhelminths, places the enigmatic Udonella caligorum firmly as a monopisthocotylean monogenean. Both maximum parsimony and a modified distance measure, operating under a maximum likelihood model, gave identical solutions for each data set. These data further support morphological evidence from ultrastructural studies indicating the neodermatan affinities of Udonella, namely shared features in sensory receptors, surface tegument, sperm structure and spermiogenesis. The molecular data reject the class Udonellidea and the placement of udonellids as sister-group to the Neodermata. As shown previously with molecular data, the monogeneans appear as a paraphyletic assemblage comprising strongly monophyletic Monopisthocotylea and Polyopisthocotylea. Their relationships with the trematodes and cestodes are not resolved with 28S rDNA or 18S rDNA alone.

A robust molecular phylogeny of the Tricladida (Platyhelminthes: Seriata) with a discussion on morphological synapomorphies

The suborder Tricladida (Platyhelminthes: Turbellaria, Seriata) comprises most well-known species of free-living flatworms. Four infraorders are recognized: (i) the Maricola (marine planarians); (ii) the Cavernicola (a group of primarily cavernicolan planarians); (iii) the Paludicola (freshwater planarians); and (iv) the Terricola (land planarians). The phylogenetic relationships among these infraorders have been analysed using morphological characters, but they remain uncertain. Here we analyse the phylogeny and classification of the Tricladida, with additional, independent, molecular data from complete sequences of 18S rDNA and 18S rRNA. We use maximum parsimony and neighbour-joining methods and the characterization of a unique gene duplication event involving the Terricola and the dugesiids to reconstruct the phylogeny. The results show that the Maricola is monophyletic and is the primitive sister group to the rest of the Tricladida (the Paludicola plus the Terricola). The Paludicola are paraphyletic since the Terricola and one paludicolan family, the Dugesiidae, share a more recent common ancestor than the dugesiids with other paludicolans (dendrocoelids and planariids). A reassessment of morphological evidence may confirm the apparent redundancy of the existing infraorders Paludicola and Terricola. In the meantime, we suggest replacing the Paludicola and Terricola with a new clade, the Continenticola, which comprises the families Dugesiidae, Planariidae, Dendrocoelidae and the Terricola.

Gnathostomulida--an enigmatic metazoan phylum from both morphological and molecular perspectives

On the basis of few and contentious morphological characters Gnathostomulids have been thought to be the sister-group of either the Platyhelminthes or the Syndermata (Rotifera + Acanthocephala). We provide a full 18S rDNA sequence for a species of Gnathostomula and attempt to resolve its position among the Metazoa, on the basis of molecular evidence. Sixty sequences, representing 30 nominal phyla and including new entoproct and gastrotrich sequences, were used to reconstruct phylogenies using maximum-parsimony, neighbor-joining, and minimum evolution models. We were unable to support either of the morphological hypotheses outright and, moreover, our data supported more strongly a third possible relationship with the gnathostomulids as a member of the Nematoda + Chaetognatha clade. Superficially, as active benthic, vermiform creatures with sclerotized cuticular jaws, they fit a predicted ancestral form of the Nematoda + Chaetognatha clade and, as such, would arguably be members of the Ecdysozoa. The molecular data at least call for a reevaluation of the morphological data and a denser sampling of the lesser phyla. Data from morphology and molecules act synergistically in estimating phylogeny; morphology alone provided limited phylogenetic signal and alternative phylogenetic hypotheses, whereas the molecular solution suggested an alternative topology which, when interpreted in the light of comparative anatomy, may suggest previously unconsidered possibilities.

Parasite speciation within or between host species?--phylogenetic evidence from site-specific polystome monogeneans

Partial nuclear 28S ribosomal RNA and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences (953 and 385 nucleotides, respectively) of one fish monogenean (outgroup) and six polystome monogeneans (four Polystomoides spp. from the oral cavities and urinary bladders of freshwater turtles in Australia and Malaya, two Neopolystoma spp. from the urinary bladder and conjunctival sac of a freshwater turtle in Australia) were used to examine the question of whether congeneric species infecting different sites in the same host species have speciated in that host by adapting to different sites, or whether species infecting a particular site in one host have given rise to species infecting the same site in different hosts. Results show unequivocally that congeneric species infecting the same site, even of host species belonging to different suborders and occurring on different continents, are more closely related than congeneric species infecting different sites of the same host species. This is interpreted as meaning that speciation has not occurred in one host. Morphological evolution of polystomes has been very slow: few differences between species and even genera have evolved over a period of at least 150 Myr, and this is matched by low substitution rates of nucleotides, and the ambiguous position of species of different genera, depending on whether COI or 28S rDNA sequences are used.

Molecular phylogenetics of the four Schistosoma species groups determined with partial 28S ribosomal RNA gene sequences

Partial 28S ribosomal RNA (rRNA) gene sequences, including the variable domains D1, D2 and D3, were determined for representative species from the 4 Schistosoma species groups. On an alignment of 1345 bp from S. mansoni, S. haematobium, S. spindale and S. japonicum (with Heterobilharzia americana chosen as an outgroup), both maximum likelihood and maximum parsimony analyses provide a robust molecular phylogeny for the genus; ((((S. haematobium, S. spindale), S. mansoni), S. japonicum), H. americana). When analysed separately, both domain D1 and domain D2 yielded similarly informative data whereas D3 failed to resolve the phylogeny. These results confirm a phylogeny previously suggested by 18S rRNA gene sequences, corroborating the status of S. spindale as a sister taxon to S. haematobium, and demonstrate the utility of 28S rRNA gene sequence data for resolving phylogenies within the Schistosomatidae.

A combined morphological and molecular phylogeny for sea urchins (Echinoidea: Echinodermata)

Phylogenetic relationships of higher taxa of echinoids have been investigated using a 163 character morphological data base and molecular sequences from large and small subunit (LSU and SSU) ribosomal RNA (rRNA) genes. The complete SSU rRNA gene has been sequenced for 21 taxa, with representatives from nine of the 14 extant orders of Echinoidea. Partial LSU sequences, representing the first 400 base pairs (b.p.) from the 5' end were also sequenced for three taxa to complement an existing data base of ten taxa. The two molecular sequences provided a total of 371 variable sites, of which 143 were phylogenetically informative (compared to 145 phylogenetically informative sites from morphological data). Morphological, LSU and SSU data have been analysed separately and together. Morphological and SSU sequence data generate topologies that are not significantly in conflict (under Templeton's test), but the strong signal pairing arbaciids with clypeasteroids in the LSU derived tree marks the LSU sequence data as anomalous for this taxon. A 'total evidence' approach derived a tree very similar in topology to that derived from morphological data. Rooted on the stem group echinoid Archaeocidaris, our total evidence tree suggested relationships of higher taxa as follows: Cidaroida (Phormosomatidae (Echinothuriidae (Diadematidae ((Spatangoida (Clypeasteroida, Cassiduloida)) ((Calycina, Arbacioida) (Stomopneustidae (Glyphocidaridae (Temnopleuridae (Echinometridae (Echinidae, Stronglyocentridae)))))))))). Phylogenetic analyses run both with and without key fossil taxa yielded slightly different topologies. It is important to include fossil taxa in a phylogenetic analysis where there are long stem-group branches or where the crown group is highly derived.

Molecular phylogenetics of cupped oysters based on partial 28S rRNA gene sequences

Partial sequences of 28S-like rDNA were amplified using PCR and sequenced for eight species of oyster and one species of mussel. Phylogenetic relationships among seven species of Crassostreinid oyster were inferred from aligned sequences by parsimony and maximum-likelihood methods. Of the 315 sites that varied, 90 were phylogenetically informative in parsimony analysis. Inference by maximum parsimony (MP) is consistent with maximum-likelihood (ML) analysis for the major lineages, yielding a tree with the topology (Mytilus edulis (Ostrea edulis ((Crassostrea rivularis (C. belcheri, C. gigas))(C. virginica, C. rhizophorae, Saccostrea cuccullata, S. commercialis)))). MP and ML analyses resolved the systematic relationships of the Saccostrea and Atlantic Crassostrea differently such that a polytomy linking these four taxa is preferred with the data available. Molecular data support a later divergence of the tropical Pacific Saccostrea from a common ancestor of the Atlantic Crassostrea species. Molecular data from domains D1, D2, and partial D3 of the 28S rDNA supply sufficient phylogenetic information to determine systematic relationships among the extant oyster taxa, from the major species groups to the family level, thus providing valuable characters that are able to supplement the paucity of morphological characters so far recognized.

Small subunit ribosomal RNA gene sequence of the oyster parasite Perkinsus marinus

The small subunit rRNA gene of the oyster pathogen Perkinsus marinus was characterized from cells of infected oyster hemolymph by polymerase chain reaction and molecular cloning. The gene, 1,793 nucleotides in size, has 77.2% sequence similarity to that of its host, the eastern oyster Crassostrea virginica. The sequence was confirmed using recently available in vitro cultures of P. marinus. DNA from pure P. marinus culture was amplified with specific primers synthesized according to the sequence from infected oyster hemolymph, and predicted size fragments were obtained. Furthermore, restriction digests yielded fragments of expected size in amplified rDNA from in vitro cultures. The P. marinus sequence has 97.5% similarity to the Perkinsus sp. sequence from the Australian mollusc Anadara trapezia.

Cellular responses of oysters infected with Haplosporidium nelsoni: changes in circulating and tissue-infiltrating hemocytes

The protozoan parasite Haplosporidium nelsoni (MSX) elicits an inflammatory-type response in oysters, Crassostrea virginica. We assayed circulating hemocytes of oysters exposed to H. nelsoni to quantify the effects of parasitism, selection for resistance, and season on total and differential counts. All factors had a significant (P < 0.02) effect on total counts, but explained relatively little of the overall variance [season (12%) > selection (4%) > infection (2%)]. Circulating hemocyte densities increased with infection intensity in resistant animals, but were depressed in susceptible oysters with light infections. Counts in both groups were lowest in August and highest in May, at which times densities in susceptible oysters were reduced (P < 0.02) compared to those in resistant animals. No differences existed in November. Susceptible oysters may be more debilitated by infections than resistant animals, resulting in impaired circulation and depressed cell counts. The number of tissue-infiltrating hemocytes increased with infection intensity (P < 0.0001), but showed no statistically significant (P > 0.05) association with numbers of circulating hemocytes in individual oysters. Size-frequency analysis with a Coulter counter indicated that the proportion of large cells (presumed to be granular hemocytes) was lower (P < 0.0001) in susceptible oysters, which were also heavily infected, compared to resistant oysters, which had very few infections. The loss of granular hemocytes may stem from degranulation associated with tissue damage and inflammation. Present evidence suggests that the principal role of the hemocyte response in MSX disease is to plug lesions, remove debris, and repair tissues and that these functions may help oysters survive infection.(ABSTRACT TRUNCATED AT 250 WORDS)