Evaluation of genome skimming to detect and characterise human and livestock helminths

The identification of gastrointestinal helminth infections of humans and livestock almost exclusively relies on the detection of eggs or larvae in faeces, followed by manual counting and morphological characterisation to differentiate species using microscopy-based techniques. However, molecular approaches based on the detection and quantification of parasite DNA are becoming more prevalent, increasing the sensitivity, specificity and throughput of diagnostic assays. High-throughput sequencing, from single PCR targets through to the analysis of whole genomes, offers significant promise towards providing information-rich data that may add value beyond traditional and conventional molecular approaches; however, thus far, its utility has not been fully explored to detect helminths in faecal samples. In this study, low-depth whole genome sequencing, i.e. genome skimming, has been applied to detect and characterise helminth diversity in a set of helminth-infected human and livestock faecal material. The strengths and limitations of this approach are evaluated using three methods to characterise and differentiate metagenomic sequencing data based on (i) mapping to whole mitochondrial genomes, (ii) whole genome assemblies, and (iii) a comprehensive internal transcribed spacer 2 (ITS2) database, together with validation using quantitative PCR (qPCR). Our analyses suggest that genome skimming can successfully identify most single and multi-species infections reported by qPCR and can provide sufficient coverage within some samples to resolve consensus mitochondrial genomes, thus facilitating phylogenetic analyses of selected genera, e.g. Ascaris spp. Key to this approach is both the availability and integrity of helminth reference genomes, some of which are currently contaminated with bacterial and host sequences. The success of genome skimming of faecal DNA is dependent on the availability of vouchered sequences of helminths spanning both taxonomic and geographic diversity, together with methods to detect or amplify minute quantities of parasite nucleic acids in mixed samples.

The increased sensitivity of qPCR in comparison to Kato-Katz is required for the accurate assessment of the prevalence of soil-transmitted helminth infection in settings that have received multiple rounds of mass drug administration

Background

The most commonly used diagnostic tool for soil-transmitted helminths (STH) is the Kato-Katz (KK) thick smear technique. However, numerous studies have suggested that the sensitivity of KK can be problematic, especially in low prevalence and low intensity settings. An emerging alternative is quantitative polymerase chain reaction (qPCR).

Methods

In this study, both KK and qPCR were conducted on stool samples from 648 participants in an STH epidemiology study conducted in the delta region of Myanmar in June 2016.

Results

Prevalence of any STH was 20.68% by KK and 45.06% by qPCR. Prevalence of each individual STH was also higher by qPCR than KK, the biggest difference was for hookworm with an approximately 4-fold increase between the two diagnostic techniques. Prevalence of Ancylostoma ceylanicum, a parasite predominately found in dogs, was 4.63%, indicating that there is the possibility of zoonotic transmission in the study setting. In individuals with moderate to high intensity infections there is evidence for a linear relationship between eggs per gram (EPG) of faeces, derived from KK, and DNA copy number, derived from qPCR which is particularly strong for Ascaris lumbricoides.

Conclusions

The use of qPCR in low prevalence settings is important to accurately assess the epidemiological situation and plan control strategies for the ‘end game’. However, more work is required to accurately assess STH intensity from qPCR results and to reduce the cost of qPCR so that is widely accessible in STH endemic countries.

The first mitochondrial genomes of endosymbiotic rhabdocoels illustrate evolutionary relaxation of atp8 and genome plasticity in flatworms

The first three mitochondrial (mt) genomes of endosymbiotic turbellarian flatworms are characterised for the rhabdocoels Graffilla buccinicola, Syndesmis echinorum and S. kurakaikina. Interspecific comparison of the three newly obtained sequences and the only previously characterised rhabdocoel, the free-living species Bothromesostoma personatum, reveals high mt genomic variability, including numerous rearrangements. The first intrageneric comparison within rhabdocoels shows that gene order is not fully conserved even between congeneric species. Atp8, until recently assumed absent in flatworms, was putatively annotated in two sequences. Selection pressure was tested in a phylogenetic framework and is shown to be significantly relaxed in this and another protein-coding gene: cox1. If present, atp8 appears highly derived in platyhelminths and its functionality needs to be addressed in future research. Our findings for the first time allude to a large degree of undiscovered (mt) genomic plasticity in rhabdocoels. It merits further attention whether this variation is correlated with a symbiotic lifestyle. Our results illustrate that this phenomenon is widespread in flatworms as a whole and not exclusive to the better-studied neodermatans.

Quantitative PCR-Based Diagnosis of Soil-Transmitted Helminth Infections: Faecal or Fickle?

Treatment and control programmes tackling soil-transmitted helminth (STH) infections require sensitive, reliable, and accurate diagnostic tools. There is a growing need for measures of infection intensity as programmes approach STH control. Quantitative real-time PCR (qPCR) is well suited to the detection of DNA targets present in stool, even in low-prevalence settings. Detecting low levels of infection becomes increasingly important when the breakpoint of transmission is approached, and is vital when monitoring for recrudescence once control, or possibly 'elimination', is achieved. We address key challenges and questions that remain as barriers to incorporating qPCR as a cornerstone diagnostic tool for STH infections.

Phylogenetic relationships within Dicrocoeliidae (Platyhelminthes: Digenea) from birds from the Czech Republic using partial 28S rDNA sequences

Partial (D1-D3) 28S rRNA gene sequences from 16 isolates of digenean parasites of the family Dicrocoeliidae recovered from 16 bird species from the Czech Republic were used for phylogenetic reconstruction. Comparison with sequences available from GenBank suggests that the genus Brachylecithum is paraphyletic, requiring further validation and possible systematic revision. Although partial 28S rDNA is relatively conserved, analyses suggest that the following taxa are synonymous: Lutztrema attenuatum = L. monenteron = L. microstomum, Brachylecithum lobatum = B. glareoli. Zonorchis petiolatus is reassigned back to the genus Lyperosomum with L. collurionis as a junior synonym. The study revealed how complicated the systematics of the family Dicrocoeliidae is currently. The morphology of the group is variable, and the current distinguishing characters at species and even generic level are not sufficiently distinctive; it is difficult to identify the specimens correctly and identification of GenBank isolates is not reliable. Extensive sampling of isolates for both molecular and morphological studies is necessary to resolve the relationships within the family.

The first next-generation sequencing approach to the mitochondrial phylogeny of African monogenean parasites (Platyhelminthes: Gyrodactylidae and Dactylogyridae)

BACKGROUND

Monogenean flatworms are the main ectoparasites of fishes. Representatives of the species-rich families Gyrodactylidae and Dactylogyridae, especially those infecting cichlid fishes and clariid catfishes, are important parasites in African aquaculture, even more so due to the massive anthropogenic translocation of their hosts worldwide. Several questions on their evolution, such as the phylogenetic position of Macrogyrodactylus and the highly speciose Gyrodactylus, remain unresolved with available molecular markers. Also, diagnostics and population-level research would benefit from the development of higher-resolution genetic markers. We aim to offer genetic resources for work on African monogeneans by providing mitogenomic data of four species (two belonging to Gyrodactylidae, two to Dactylogyridae), and analysing their gene sequences and gene order from a phylogenetic perspective.

RESULTS

Using Illumina technology, the first four mitochondrial genomes of African monogeneans were assembled and annotated for the cichlid parasites Gyrodactylus nyanzae, Cichlidogyrus halli, Cichlidogyrus mbirizei (near-complete mitogenome) and the catfish parasite Macrogyrodactylus karibae (near-complete mitogenome). Complete nuclear ribosomal operons were also retrieved, as molecular vouchers. The start codon TTG is new for Gyrodactylus and for Dactylogyridae, as is the incomplete stop codon TA for Dactylogyridae. Especially the nad2 gene is promising for primer development. Gene order was identical for protein-coding genes and differed between the African representatives of these families only in a tRNA gene transposition. A mitochondrial phylogeny based on an alignment of nearly 12,500 bp including 12 protein-coding and two ribosomal RNA genes confirms that the Neotropical oviparous Aglaiogyrodactylus forficulatus takes a sister group position with respect to the other gyrodactylids, instead of the supposedly 'primitive' African Macrogyrodactylus. Inclusion of the African Gyrodactylus nyanzae confirms the paraphyly of Gyrodactylus. The position of the African dactylogyrid Cichlidogyrus is unresolved, although gene order suggests it is closely related to marine ancyrocephalines.

CONCLUSIONS

The amount of mitogenomic data available for gyrodactylids and dactylogyrids is increased by roughly one-third. Our study underscores the potential of mitochondrial genes and gene order in flatworm phylogenetics, and of next-generation sequencing for marker development for these non-model helminths for which few primers are available.

Curious bivalves: Systematic utility and unusual properties of anomalodesmatan mitochondrial genomes

Mitogenomic trees for Bivalvia have proved problematic in the past, but several highly divergent lineages were missing from these analyses and increased representation of these groups may yet improve resolution. Here, we add seven new sequences from the Anomalodesmata and one unidentified semelid species (Bryopa lata, Euciroa cf. queenslandica, Laternula elliptica, Laternula truncata, Lyonsia norwegica, Myadora brevis, Tropidomya abbreviata, "Abra" sp.). We show that relationships in a mitogenomic tree for the Class are improved by the addition of seven anomalodesmatans from this highly divergent clade, but are still not completely consistent with relationships recovered in studies of nuclear genes. We suggest that some anomalous relationships (for instance the non-monophyly of Bivalvia) may be partially explained by compositional heterogeneity in the mitogenome and suggest that the addition of more taxa may help resolve both this effect and possible instances of long branch attraction. We also identify several curious features about anomalodesmatan mitogenomes. For example, many protein-coding gene boundaries are poorly defined in marine bivalves, but particularly so in anomalodesmatans, primarily due to non-conserved boundary sequences. The use of transcriptomic and genomic data together enabled better definition of gene boundaries, the identification of possible pseudogenes and suggests that most genes are translated monocistronically, which contrasts with many other studies. We also identified a possible case of gene duplication of ND5 in Myadora brevis (Myochamidae). Mitogenome size in the Anomalodesmata ranges from very small compact molecules, with the smallest for Laternula elliptica (Laternulidae) only 14,622bp, to Bryopa lata (Clavagellidae) which is at least 31,969bp long and may be >40,000bp. Finally, sampled species show a high degree of sequence divergence and variable gene order, although intraspecific variation in Laternula elliptica is very low.

Next-Generation Mitogenomics: A Comparison of Approaches Applied to Caecilian Amphibian Phylogeny

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a ‘traditional’ Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina’s HiSeq and MiSeq, Roche’s 454 GS FLX, and Life Technologies’ Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case.

The mitochondrial genome and ribosomal operon of Brachycladium goliath (Digenea: Brachycladiidae) recovered from a stranded minke whale

Members of the Brachycladiidae are known to cause pathologies implicated in cetacean strandings and it is important to develop accurate diagnostic markers to differentiate these and other helminths found in cetaceans. Brachycladium goliath (van Beneden, 1858) is a large trematode found, as adults, usually in the hepatic (bile) and pancreatic ducts of various cetaceans. Complete sequences were determined for the entire mitochondrial genome, and phylogenetically informative nuclear genes contained within the ribosomal operon, from a small piece of an individual worm taken from a common minke whale Balaenoptera acutorostrata Lacépède, 1804. Genomic DNA was sequenced using an Illumina MiSeq platform. The mtDNA is 15,229 bp in length consisting of 12 protein-coding genes, 22 tRNA genes, 2 rRNA genes and 2 non-coding regions of which the larger is comprised of 4 tandemly repeated units (260 bp each). The ribosomal RNA operon is 9297 bp long. These data provide a rich resource of molecular markers for diagnostics, phylogenetics and population genetics in order to better understand the role, and associated pathology of helminth infections in cetaceans.

Whipworms in humans and pigs: origins and demography

Background

Trichuris suis and T. trichiura are two different whipworm species that infect pigs and humans, respectively. T. suis is found in pigs worldwide while T. trichiura is responsible for nearly 460 million infections in people, mainly in areas of poor sanitation in tropical and subtropical areas. The evolutionary relationship and the historical factors responsible for this worldwide distribution are poorly understood. In this study, we aimed to reconstruct the demographic history of Trichuris in humans and pigs, the evolutionary origin of Trichuris in these hosts and factors responsible for parasite dispersal globally.

Methods

Parts of the mitochondrial nad1 and rrnL genes were sequenced followed by population genetic and phylogenetic analyses. Populations of Trichuris examined were recovered from humans (n = 31), pigs (n = 58) and non-human primates (n = 49) in different countries on different continents, namely Denmark, USA, Uganda, Ecuador, China and St. Kitts (Caribbean). Additional sequences available from GenBank were incorporated into the analyses.

Results

We found no differentiation between human-derived Trichuris in Uganda and the majority of the Trichuris samples from non-human primates suggesting a common African origin of the parasite, which then was transmitted to Asia and further to South America. On the other hand, there was no differentiation between pig-derived Trichuris from Europe and the New World suggesting dispersal relates to human activities by transporting pigs and their parasites through colonisation and trade. Evidence for recent pig transport from China to Ecuador and from Europe to Uganda was also observed from their parasites. In contrast, there was high genetic differentiation between the pig Trichuris in Denmark and China in concordance with the host genetics.

Conclusions

We found evidence for an African origin of T. trichiura which were then transmitted with human ancestors to Asia and further to South America. A host shift to pigs may have occurred in Asia from where T. suis seems to have been transmitted globally by a combination of natural host dispersal and anthropogenic factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1325-8) contains supplementary material, which is available to authorized users.

New digeneans (Opecoelidae) from hydrothermal vent fishes in the south eastern Pacific Ocean, including one new genus and five new species

A new genus and five new species of digeneans are reported from fishes at hydrothermal vent sites in the South East Pacific Rise region. Biospeedotrema n. gen. (Opecoelidae: Stenakrinae) is distinguished from other stenakrines by the more or less symmetrical testicular configuration, with the uterus passing between the testes, sometimes distinctly into the post-testicular region. Biospeedotrema jolliveti n. gen., n. sp. from Ventichthys biospeedoi (Ophidiidae) is distinguished by the vitelline fields which extend only slightly into the post-testicular region, the intestinal bifurcation is dorsal to the ventral sucker, the genital pore is slightly dextrally submedian or median, the cirrus sac is short and the caeca are broad and overlap the testes, usually reaching into the post-testicular region. Biospeedotrema parajolliveti n. sp. from Thermichthys hollisi differs from Biospeedotrema jolliveti in being squat, always just wider than long, the tegument is wrinkled, the testes are lobate, and the caeca only just reach to the testes. Biospeedotrema biospeedoi n. sp. from T. hollisi differs from its congeners in its body-shape, uterine extent posterior to the testes and the small vitellarium. Caudotestis ventichthysi n. sp. (Opecoelidae: Stenakrinae) from V. biospeedoi is distinguished from its five congeners in various combinations of caecal length, cirrus sac length, internal seminal vesicle shape, vitelline extent and distribution, forebody length and egg-size. Buticulotrema thermichthysi n. sp. (Opecoelidae: Opecoelininae) from T. hollisi (Bythitidae) is distinguished from its only congener by its very long, very strongly muscular oesophagus, bifurcating dorsally to the posterior part of the ventral sucker, the long, narrow pars prostatica and distal male duct and the sinistral genital pore at the level of the pharynx. The phylogenetic position for three of these species, Buticulotrema thermichthysi, Biospeedotrema jolliveti and Biospeedotrema biospeedoi, is assessed based on ssrDNA and lsrDNA sequences, which verify the position of these species in the Opecoelidae. 

Problematic barcoding in flatworms: A case-study on monogeneans and rhabdocoels (Platyhelminthes)

Some taxonomic groups are less amenable to mitochondrial DNA barcoding than others. Due to the paucity of molecular information of understudied groups and the huge molecular diversity within flatworms, primer design has been hampered. Indeed, all attempts to develop universal flatworm-specific COI markers have failed so far. We demonstrate how high molecular variability and contamination problems limit the possibilities for barcoding using standard COI-based protocols in flatworms. As a consequence, molecular identification methods often rely on other widely applicable markers. In the case of Monogenea, a very diverse group of platyhelminth parasites, and Rhabdocoela, representing one-fourth of all free-living flatworm taxa, this has led to a relatively high availability of nuclear ITS and 18S/28S rDNA sequences on GenBank. In a comparison of the effectiveness in species assignment we conclude that mitochondrial and nuclear ribosomal markers perform equally well. In case intraspecific information is needed, rDNA sequences can guide the selection of the appropriate (i.e. taxon-specific) COI primers if available.

Adding resolution to ordinal level relationships of tapeworms (Platyhelminthes: Cestoda) with large fragments of mtDNA

The construction of a stable phylogeny for the Cestoda, indicating the interrelationships of recognised orders and other major lineages, has proceeded iteratively since the group first received attention from phylogenetic systematists. Molecular analyses using nuclear ribosomal RNA gene fragments from the small (ssrDNA) and large (lsrDNA) subunits have been used to test competing evolutionary scenarios based on morphological data but could not arbitrate between some key conflicting hypotheses. To the ribosomal data, we have added a contiguous fragment of mitochondrial (mt) genome data (mtDNA) of partial nad1-trnN-trnP-trnI-trnK-nad3-trnS-trnW-cox1-trnT-rrnL-trnC-partial rrnS, spanning 4034-4447 bp, where new data for this region were generated for 18 species. Bayesian analysis of mtDNA and rDNA as nucleotides, and where appropriate as amino acids, demonstrated that these two classes of genes provide complementary signal across the phylogeny. In all analyses, except when using mt amino acids only, the Gyrocotylidea is sister group to all other Cestoda (Nephroposticophora), and Amphilinidea forms the sister group to the Eucestoda. However, an earliest-diverging position of Amphilinidea is strongly supported in the mt amino acid analysis. Amphilinidea exhibit a unique tRNA arrangement (nad1-trnI-trnL2-trnP-trnK-trnV-trnA-trnN-nad3), whereas Gyrocotylidea shares that of the derived lineages, providing additional evidence of the uniqueness of amphilinid genes and genomes. The addition of mtDNA to the rDNA genes supported the Caryophyllidea as the sister group to (Spathebothriidea+remaining Eucestoda), a hypothesis consistently supported by morphology. This relationship suggests a history of step-wise evolutionary transitions from simple monozoic, unsegmented tapeworms to the more familiar polyzoic, externally segmented (strobilate) forms. All our data partitions recovered Haplobothriidea as the sister group to Diphyllobothriidae. The sister-group relationship between Diphyllidea and Trypanorhyncha, as previously established using rDNA, is not supported by the mt data, although it is supported by the combined mt and rDNA analysis. With regards to the more derived taxa, in all except the mt amino acid analysis, the following topology is supported: (Bothriocephalidea (Litobothriidea (Lecanicephalidea (Rhinebothriidea (Tetraphyllidea, (Acanthobothrium, Proteocephalidea), (Nippotaeniidea, Mesocestoididae, Tetrabothriidea, Cyclophyllidea)))))), where the Tetraphyllidea are paraphyletic. Evidence from the mt data provides strong (nucleotides) to moderate (amino acids) support for Tetraphyllidea forming a group to the inclusion of Proteocephalidea, with the latter consistently forming the sister group to Acanthobothrium. The interrelationships among Nippotaeniidea, Mesocestoididae, Tetrabothriidea and Cyclophyllidea remain ambiguous and require further systematic attention. Mitochondrial and nuclear rDNA data provide conflicting signal for certain parts of the cestode tree. In some cases mt data offer results in line with morphological evidence, such as the interrelationships of the early divergent lineages. Also, Tetraphyllidea, although remaining paraphyletic with the inclusion of the Proteocephalidea, does not include the most derived cestodes; a result which has consistently been obtained with rDNA.

A molecular phylogeny of bryozoans

We present the most comprehensive molecular phylogeny of bryozoans to date. Our concatenated alignment of two nuclear ribosomal and five mitochondrial genes includes 95 taxa and 13,292 nucleotide sites, of which 8297 were included. The number of new sequences generated during this project are for each gene:ssrDNA (32), lsrDNA (22), rrnL (38), rrnS (35), cox1 (37), cox3 (34), and cytb (44). Our multi-gene analysis provides a largely stable topology across the phylum. The major groups were unambiguously resolved as (Phylactolaemata (Cyclostomata (Ctenostomata, Cheilostomata))), with Ctenostomata paraphyletic. Within Phylactolaemata, (Stephanellidae, Lophopodidae) form the earliest divergent clade. Fredericellidae is not resolved as a monophyletic family and forms a clade together with Plumatellidae, Cristatellidae and Pectinatellidae, with the latter two as sister taxa. Hyalinella and Gelatinella nest within the genus Plumatella. Cyclostome taxa fall into three major clades: i. (Favosipora (Plagioecia, Rectangulata)); ii. (Entalophoroecia ((Diplosolen, Cardioecia) (Frondipora, Cancellata))); and iii. (Articulata ((Annectocyma, Heteroporidae) (Tubulipora (Tennysonia, Idmidronea)))), with suborders Tubuliporina and Cerioporina, and family Plagioeciidae each being polyphyletic. Ctenostomata is composed of three paraphyletic clades to the inclusion of Cheilostomata: ((Alcyonidium, Flustrellidra) (Paludicella (Anguinella, Triticella)) (Hislopia (Bowerbankia, Amathia)) Cheilostomata); Flustrellidra nests within the genus Alcyonidium, and Amathia nests within the genus Bowerbankia. Suborders Carnosa and Stolonifera are not monophyletic. Within the cheilostomes, Malacostega is paraphyletic to the inclusion of all other cheilostomes. Conopeum is the most early divergent cheilostome, forming the sister group to ((Malacostega, Scrupariina, Inovicellina) ((Hippothoomorpha, Flustrina) (Lepraliomorpha, Umbonulomorpha))); Flustrina is paraphyletic to the inclusion of the hippothoomorphs; neither Lepraliomorpha nor Umbonulomorpha is monophyletic. Ascophorans are polyphyletic, with hippothoomorphs grouping separately from lepraliomorphs and umbonulomorphs; no cribrimorphs were included in the analysis. Results are discussed in the light of molecular and morphological evidence. Ancestral state reconstruction of larval strategy in Gymnolaemata revealed planktotrophy and lecithotrophy as equally parsimonious solutions for the ancestral condition. More comprehensive taxon sampling is expected to clarify this result. We discuss the extent of non-bryozoan contaminant sequences deposited in GenBank and their impact on the reconstruction of metazoan phylogenies and those of bryozoan interrelationships.

The mitochondrial genome of Gyrodactylus derjavinoides (Platyhelminthes: Monogenea)--a mitogenomic approach for Gyrodactylus species and strain identification

Systematists and evolutionary biologists are constantly on the lookout for new sources of characters to discriminate amongst taxa and estimate interrelationships within and between taxa. Entire mitochondrial genomes provide a wealth of data, both at the nucleotide and amino acid level. Molecular markers are of particular utility when applied to small, morphologically conserved taxa, as is the case for many monogenean ectoparasites of fish. Gyrodactylus species display a considerable degree of anatomical conservatism, complicating diagnostics based solely on morphology, and some are significant pests of wild and cultured fish. Here we sequenced the complete mitochondrial genome of Gyrodactylus derjavinoides Malmberg, Collins, Cunningham & Behiar 2007, one of the most frequently found gyrodactylid species on salmonids in Scandinavia, and compared it with the recently published genomes of Gyrodactylus salaris Malmberg, 1957 and Gyrodactylus thymalli Zitnan 1960. Through comparative sliding window analysis we identified regions of high sequence variability and designed new primer sequences. In total, 6 new primer pairs have been developed, amplifying fragments of cox1, cox3, nad1, nad2, nad4, nad5 and atp6. Together, they amplify regions capturing almost half the nucleotide variability present in the complete mitochondrial genome. These degenerate primers should also work for other Gyrodactylus species parasitizing salmonids. In addition, we developed a multiplex assay that simultaneously amplifies four fragments in a single PCR reaction. Besides the diagnostic value, these fragments can be used for studying the transmission dynamics of Gyrodactylus, providing crucial information for an improved understanding of the spread and epidemiology of these important fish pathogens.

The complete mitochondrial DNA sequence of the monogenean Gyrodactylus thymalli (Platyhelminthes: Monogenea), a parasite of grayling (Thymallus thymallus)

We present the complete mitochondrial (mt) genome of Gyrodactylus thymalli, a monogenean ectoparasite on grayling (Thymallus thymallus). The circular genome is 14788 bp in size and includes all 35 genes recognized from other flatworm mt genomes. The overall A+T content of the mt genome is 62.8%. Twenty regions of non-coding DNA ranging from 1 to 111 bp in length were identified in addition to 2 highly conserved large non-coding regions 799 bp and 767 bp in size. Compared to the recently described mt DNA of the closely related G. salaris from Atlantic salmon from Signaldalselva, Norway, the mitochondrial genome of G. thymalli from Hnilec, Slovakia, differs on average by 2.2%.

The complete mitochondrial genome of the bird schistosome Trichobilharzia regenti (Platyhelminthes: Digenea), causative agent of cercarial dermatitis

The complete mitochondrial (mt) genome of the neuropathogenic bird schistosome Trichobilharzia regenti was fully sequenced in order to develop molecular markers for future diagnostic, molecular ecological, population, and phylogenetic studies. The genome was 14,838 bp in length, with a 68.4% AT bias in protein coding regions. A repeat element (3 x 184 bp) between trnV and trnW distinguished a single short noncoding region. As 9 of 14 genera of schistosomes parasitize birds, future characterization of their mt genomes is desirable for species-specific and strain- or population-specific diagnostic markers; this concerns not only the nasal representatives, e.g., T. regenti characterized in this study, but also numerous species within the predominant group of visceral (blood dwelling) bird schistosomes.

Bird schistosomes of wildfowl in the Czech Republic and Poland

In 2005, we dissected 102 wildfowl from the Czech Republic and 73 wildfowl from Poland including representatives of Anseriformes, Gruiformes and Gaviiformes. Schistosome infection was found in a total of 21 (29%) and 23 (23%) birds from Poland and the Czech Republic, respectively. All infected birds belonged to the order Anseriformes. The prevalences of nasal and visceral species were, respectively, 22% and 16% in Poland and 6% and 19% in the Czech Republic. Four species of schistosomes were found: Bilharziella polonica Kowalewski, 1895, Trichobilharzia regenti Horák, Koláfová et Dvorák, 1998, T. szidati Neuhaus, 1952, and an undetermined schistosome from the intestinal wall of Anas penelope L. The finding of T. szidati represents the first record of the parasite from natural final host since the species description.

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.

Added resolution among ordinal level relationships of tapeworms (Platyhelminthes: Cestoda) with complete small and large subunit nuclear ribosomal RNA genes

The addition of large subunit ribosomal DNA (lsrDNA) to small subunit ribosomal DNA (ssrDNA) has been shown to add resolution to phylogenies at various taxonomic levels for a diversity of phyla. We added nearly complete lsrDNA (4057-4593bp) sequences to ssrDNA (1940-2228bp) for 26 ingroup and 3 outgroup taxa in an attempt to provide an improved ordinal phylogeny for the Cestoda. Ten lsrDNA and seven ssrDNA sequences were generated from new taxa and 13 existing partial lsrDNA sequences were sequenced to completion. The majority of phylogenetic signal in the combined analysis came from lsrDNA (69.6% of parsimonious informative sites, as opposed to 30.4% obtained from ssrDNA), resulting in almost identical topologies for lsrDNA and lsr+ssrDNA (pairwise symmetric distance=6) in model-based analyses. Topology testing found trees based on partial lsrDNA (domains D1-D3)+ssrDNA and complete lsr+ssrDNA to differ significantly; the addition of lsrDNA domains D4-D12 had a significant effect on topology. Overall nodal support was greatest in the combined analysis and weakest for ssrDNA only. Our molecular phylogenies differed significantly from those based on morphology alone. Acetabulate lineages form a monophyletic group, with the Tetraphyllidea being paraphyletic. Support for the combined data was high for the following topology: (Litobothriidea (Lecanicephalidea (Rhinebothrium/Rhodobothrium (Clistobothrium (Pachybothrium(Acanthobothrium Proteocephalidea) (Mesocestoididae, Nippotaeniidea, Cyclophyllidea, Tetrabothriidea)))))); all genus names refer to tetraphyllidean lineages. Although the interrelationships among the four most derived taxa remain uncertain, overall ambiguity of the acetabulate interrelationships was reduced. The Pseudophyllidea were recovered as polyphyletic, with support for a sister-group relationship between Diphyllobothriidae and Haplobothriidea. The monophyly of the Trypanorhyncha was recovered for the first time based on molecular data. The positions of the Trypanorhyncha, Diphyllidea and "Bothriocephaliidea" in relation to other orders remains ambiguous. Higher congruence was found between trees based on model-based phylogenetic methods than with those constructed under the parsimony criterion. Although some uncertainties remain, the addition of lsrDNA D4-D12 has provided an overall more resolved and better supported cestode phylogeny, which further promotes the utility of complete lsrDNA as phylogenetic marker where ssrDNA alone proves inadequate.

A common origin of complex life cycles in parasitic flatworms: evidence from the complete mitochondrial genome of Microcotyle sebastis (Monogenea: Platyhelminthes)

BACKGROUND

The parasitic Platyhelminthes (Neodermata) contains three parasitic groups of flatworms, each having a unique morphology, and life style: Monogenea (primarily ectoparasitic), Trematoda (endoparasitic flukes), and Cestoda (endoparasitic tapeworms). The evolutionary origin of complex life cyles (multiple obligate hosts, as found in Trematoda and Cestoda) and of endo-/ecto-parasitism in these groups is still under debate and these questions can be resolved, only if the phylogenetic position of the Monogenea within the Neodermata clade is correctly estimated.

RESULTS

To test the interrelationships of the major parasitic flatworm groups, we estimated the phylogeny of the Neodermata using complete available mitochondrial genome sequences and a newly characterized sequence of a polyopisthocotylean monogenean Microcotyle sebastis. Comparisons of inferred amino acid sequences and gene arrangement patterns with other published flatworm mtDNAs indicate Monogenea are sister group to a clade of Trematoda+Cestoda.

CONCLUSION

Results confirm that vertebrates were the first host for stem group neodermatans and that the addition of a second, invertebrate, host was a single event occurring in the Trematoda+Cestoda lineage. In other words, the move from direct life cycles with one host to complex life cycles with multiple hosts was a single evolutionary event. In association with the evolution of life cycle patterns, our result supports the hypothesis that the most recent common ancestor of the Neodermata giving rise to the Monogenea adopted vertebrate ectoparasitism as its initial life cycle pattern and that the intermediate hosts of the Trematoda (molluscs) and Cestoda (crustaceans) were subsequently added into the endoparasitic life cycles of the Trematoda+Cestoda clade after the common ancestor of these branched off from the monogenean lineage. Complex life cycles, involving one or more intermediate hosts, arose through the addition of intermediate hosts and not the addition of a vertebrate definitive host. Additional evidence is required from monopisthocotylean monogeneans in order to confirm the monophyly of the group.

Insight into the role of cetaceans in the life cycle of the tetraphyllideans (Platyhelminthes: Cestoda)

Four types of tetraphyllidean larvae infect cetaceans worldwide: two plerocercoids differing in size, 'small' (SP) and 'large' (LP), and two merocercoids referred to as Phyllobothrium delphini and Monorygma grimaldii. The latter merocercoid larvae parasitize marine mammals exclusively and exhibit a specialised cystic structure. Adult stages are unknown for any of the larvae and thus the role of cetaceans in the life cycle of these species has been a long-standing problem. The SP and LP forms are thought to be earlier stages of P. delphini and M. grimaldii that are presumed to infect large pelagic sharks that feed on cetaceans. A molecular analysis of the D2 variable region of the large subunit ribosomal DNA gene based on several individuals of each larval type collected from three Mediterranean species of cetaceans showed consistent and unique molecular signatures for each type regardless of host species or site of infection. The degree of divergence suggested that LP, P. delphini and M. grimaldii larvae may represent separate species, whereas SP may be conspecific with M. grimaldii. In all host species, individuals of SP accumulated in the gut areas in which the lymphoid tissue was especially developed. We suggest therefore that these larvae use the lymphatic system to migrate to the abdominal peritoneum and mesenteries where they develop into forms recognizable as M. grimaldii. The plerocercoid stage of P. delphini remains unknown. In a partial phylogenetic tree of the Tetraphyllidea, all larvae formed a clade that included a representative of the genus Clistobothrium, some species of which parasitize sharks such as the great white which is known to feed on cetaceans. A bibliographic examination of tetraphyllidean infections in marine mammals indicated that these larvae are acquired mostly offshore. In summary, the evidence suggests that cetaceans play a significant role in the life cycle of these larvae. In addition, it seems clear that cetaceans act as natural intermediate hosts for P. delphini and M. grimaldii, as within these hosts they undergo development from the plerocercoid stage to the merocercoid stage. Because tetraphyllidean species use fish, cephalopods and other marine invertebrates as intermediate hosts, the inclusion of cetaceans in the life cycle would have facilitated their transmission to apex predators such as the large, lamnid sharks. The biological significance of infections of LP in cetaceans is unclear, but infections do not seem to be accidental as such larvae show high prevalence and abundance as well as a high degree of site specificity, particularly in the anal crypts and bile ducts.

Genetic diversity in the trypanorhynch cestode Tentacularia coryphaenae Bosc, 1797: evidence for a cosmopolitan distribution and low host specificity in the teleost intermediate host

Partial large subunit (28S) rRNA gene (LSU) sequences were studied from Tentacularia coryphaenae (Cestoda: Trypanorhyncha) plerocercoids from the southern Java coast, Indonesia, collected from two different localities and five different host species. The teleost hosts belonged to four fish families with an overlapping depth range of 0-885 m. The LSU sequences were identical, demonstrating that all specimens belonged to the same species. They also corresponded to a sequence of T. coryphaenae from the Blue shark Prionace glauca in the North Atlantic, giving genetic evidence for the cosmopolitan distribution of the species. A 1,851 bp region of mitochondrial (mt) DNA (coding for partial cytochrome c oxidase subunit I (cox1), complete trnT and partial 16S ribosomal RNA) showed a very low level of intra-specific variation of 1%. Pairwise comparisons of published sequences for partial LSU rDNA and the same region of mtDNA demonstrated that the same regions varied by 8% in the mtDNA for two genotypes (G1 and G4) of Echinococcus granulosus (order Cyclophyllidea), at 16% in newly sequenced Kotorella pronosoma from the same trypanorhynch family and at 23% in Grillotia pristiophori from a different superfamily. The high genetic homogeneity in T. coryphaenae is explained by a constant gene flow between different regions and hosts along the Indonesian coast caused by extensive migrations of the second intermediate/paratenic and also the final hosts. Implications for the zoogeographical distribution, host specificity of the species and future research are discussed.

Characterization of the mitochondrial genome ofDiphyllobothrium latum(Cestoda: Pseudophyllidea) – implications for the phylogeny of eucestodes

The complete nucleotide sequence of the mitochondrial genome was determined for the fish tapewormDiphyllobothrium latum. This genome is 13 608 bp in length and encodes 12 protein-coding genes (but lacks theatp8), 22 transfer RNA (tRNA) and 2 ribosomal RNA (rRNA) genes, corresponding to the gene complement found thus far in other flatworm mitochondrial (mt) DNAs. The gene arrangement of this pseudophyllidean cestode is the same as the 6 cyclophyllidean cestodes characterized to date, with only minor variation in structure among these other genomes; the relative position oftrnS2andtrnL1is switched inHymenolepis diminuta. Phylogenetic analyses of the concatenated amino acid sequences for 12 protein-coding genes of all complete cestode mtDNAs confirmed taxonomic and previous phylogenetic assessments, withD. latumbeing a sister taxon to the cyclophyllideans. High nodal support and phylogenetic congruence between different methods suggest that mt genomes may be of utility in resolving ordinal relationships within the cestodes. All species ofDiphyllobothriuminfect fish-eating vertebrates, andD. latumcommonly infects humans through the ingestion of raw, poorly cooked or pickled fish. The complete mitochondrial genome provides a wealth of genetic markers which could be useful for identifying different life-cycle stages and for investigating their population genetics, ecology and epidemiology.

The mitochondrial genome of Gyrodactylus salaris (Platyhelminthes: Monogenea), a pathogen of Atlantic salmon (Salmo salar)

In the present study, we describe the complete mitochondrial (mt) genome of the Atlantic salmon parasite Gyrodactylus salaris, the first for any monogenean species. The circular genome is 14 790 bp in size. All of the 35 genes recognized from other flatworm mitochondrial genomes were identified, and they are transcribed from the same strand. The protein-coding and ribosomal RNA (rRNA) genes share the same gene arrangement as those published previously for neodermatan mt genomes (representing cestodes and digeneans only), and the genome has an overall A+T content of 65%. Three transfer RNA (tRNA) genes overlap with other genes, whereas the secondary structure of 3 tRNA genes lack the DHU arm and 1 tRNA gene lacks the TΨC arm. Eighteen regions of non-coding DNA ranging from 4 to 112 bp in length, totalling 278 bp, were identified as well as 2 large non-coding regions (799 bp and 768 bp) that were almost identical to each other. The completion of the mt genome offers the opportunity of defining new molecular markers for studying evolutionary relationships within and among gyrodactylid species.

Mitogenomics and phylogenomics reveal priapulid worms as extant models of the ancestral Ecdysozoan

Research into arthropod evolution is hampered by the derived nature and rapid evolution of the best-studied out-group: the nematodes. We consider priapulids as an alternative out-group. Priapulids are a small phylum of bottom-dwelling marine worms; their tubular body with spiny proboscis or introvert has changed little over 520 million years and recognizable priapulids are common among exceptionally preserved Cambrian fossils. Using the complete mitochondrial genome and 42 nuclear genes from Priapulus caudatus, we show that priapulids are slowly evolving ecdysozoans; almost all these priapulid genes have evolved more slowly than nematode orthologs and the priapulid mitochondrial gene order may be unchanged since the Cambrian. Considering their primitive bodyplan and embryology and the great conservation of both nuclear and mitochondrial genomes, priapulids may deserve the popular epithet of "living fossil." Their study is likely to yield significant new insights into the early evolution of the Ecdysozoa and the origins of the arthropods and their kin as well as aiding inference of the morphology of ancestral Ecdysozoa and Bilateria and their genomes.

Robinia aurata n. g., n. sp. (Digenea: Hemiuridae) from the mugilid Liza aurata with a molecular confirmation of its position within the Hemiuroidea

Robinia aurata n. g., n. sp. is described from Liza aurata (Mugilidae), the golden grey mullet, from the Ebro Delta, Spanish Mediterranean. The new genus differs from all other hemiurid genera in the combined possession of muscular flanges and a vestigial ecsoma. Within the Bunocotylinae, which currently accommodates 2 genera, Bunocotyle and Saturnius, the new genus exhibits a unique combination of blind caeca, Juel's organ, post-ovarian bulk of the uterus in the hind-body, and tegumental papillae surrounding the oral and ventral sucker apertures. Furthermore, Robinia n. g. differs from both Bunocotyle and Saturnius in the nature of the muscular extensions around the oral sucker, with the shape of a muscular belt in the latter and numerous muscular papillae in the former. The phylogenetic hypothesis for the Bunocotylinae developed from sequence data analyses based on partial lsrDNA and complete ssrDNA combined (22 species) and V4 domain of the ssrRNA gene (37 species) supports the erection of the new genus and confirms its position within the Hemiuroidea. Both molecular analyses confirmed the monophyly of the Hemiuroidea, its division into 2 major clades and the polyphyly of the Derogenidae, as in previous studies, and suggest that the Gonocercinae (with 2 genera, Gonocerca and Hemipera), may require a distinct familial status. Finally, there was poor support for the distinct status of the Lecithasteridae and Hemiuridae, following previous suggestions based on different sequence data sets. A key to genera of the Bunocotylinae is presented.

Validity reassessment of Trichobilharzia species using Lymnaea stagnalis as the intermediate host

The systematics within the genus Trichobilharzia is complicated. After the description of the type species Trichobilharzia ocellata, the name was routinely used for nearly all European findings of ocellate furcocercariae. T. ocellata was also described from North America and Japan. However, the identity of T. ocellata remains questionable. Comparison of data from the literature showed differences among various T. ocellata isolates and led us to the conclusion that the North American and the Japanese findings are not identical with European T. ocellata. In addition, the description of T. szidati corresponds with the recently reported European T. ocellata isolates. Sequence analysis of the ITS region confirmed that they are identical.

Morphological and molecular characterization of tetraphyllidean merocercoids (Platyhelminthes: Cestoda) of striped dolphins (Stenella coeruleoalba) from the Western Mediterranean

Two types of tetraphyllidean merocercoids, Phyllobothrium delphini and Monorygma grimaldii, are well known from most cetaceans world-wide. The role of cetaceans in the life-cycle of these merocercoids is unclear because their specific identity is as yet unknown. The problem is compounded by poor descriptions of both merocercoids. We used light and scanning electron microscopy, and histological techniques to provide a thorough description of merocercoids collected from 11 striped dolphins, Stenella coeruleoalba, from the Spanish Mediterranean. We also described, for the first time, specimens of P. delphini with immature proglottides. Our merocercoids were morphologically similar to those described previously, except in the structure of the apical organ. Intra- and inter-sample variability in the morphology of the apical organ suggested that it degenerates during larval development. A subsample of 16 specimens of P. delphini and M. grimaldii was characterized for the D2 variable region of the large subunit ribosomal RNA gene (LSU) and compared with published tetraphyllidean cestode LSU sequences. P. delphini showed 2 unique signatures that differed from one another by a single base, whereas all sequences of M. grimaldii were identical. This suggests that each type may represent a single species, contrary to previous speculations based on morphological data. All merocercoid specimens formed a clade together with Clistobothrium montaukensis. Based on the low degree of divergence, all specimens of this clade are predicted to be congeneric.

Molecular phylogenetic analysis of the genus Gyrodactylus (Platyhelminthes: Monogenea) inferred from rDNA ITS region: subgenera versus species groups

Analyses of small subunit ribosomal RNA gene sequences of representatives of major taxa of Monopisthocotylea were performed to identify the sister group of Gyrodactylus. Nuclear ribosomal DNA sequences from the complete internal transcribed spacer (ITS) region were used to infer phylogeny of 37 Gyrodactylus species and Gyrodactyloides bychowskii, Macrogyrodactylus polypteri and Gyrdicotylus gallieni, using maximum likelihood, parsimony and Bayesian inference. The genus Gyrodactylus appeared to be a monophyletic group in all analyses, based on the present data set. Within the genus, there were 3 major groups recognized by high bootstrap values and posterior probabilities. None of the 6 subgenera appeared to be monophyletic, and the most basal subgenus G. (Gyrodactylus) was paraphyletic. Characteristics of the excretory system of Gyrodactylus do not seem to be conservative enough to reveal subgenera within Gyrodactylus and we suggest abandoning existing subgenera as indicators of phylogeny. The grouping of species based on the morphology of the ventral bar and marginal hooks seems to have sufficient power to infer relationships between the Gyrodactylus species.

A phylogeny based on three mitochondrial genes supports the division of Schistosoma intercalatum into two separate species

Two recognized strains of Schistosoma intercalatum, one from the Democratic Republic of Congo (DRC), formerly Zaire, and the other from Cameroon, have been investigated using DNA sequences from 3 mitochondrial genes, cytochrome oxidase subunit 1 (cox1), NADH dehydrogenase subunit 6 (nad6) and the small ribosomal RNA gene (rrnS). In addition, partial DNA sequences from the nuclear large subunit ribosomal RNA gene (lsrDNA) were included within the study. Although partial lsrDNA alone reveals little taxonomic information, phylogenetic analysis of the mitochondrial data demonstrates a clear dichotomy between the 2 purported strains and it is proposed that they should be treated as distinct taxa. The 'original' S. intercalatum now falls relatively basal in the S. haematobium group, while the proposed new species is more derived and sister taxon to S. bovis and S. curassoni.

Phylogeny and evolution of developmental mode in temnopleurid echinoids

The phylogenetic relationships of 24 nominal species of temnopleurid echinoid were established using molecular and morphological data sets. The analysis combined sequence data from mitochondrial 16S rRNA and cytochrome c oxidase subunit I genes and the nuclear 18S-like small subunit rRNA gene with morphological data concerning coronal, lantern, spine, and pedicellarial traits. All four data sets contain similar phylogenetic information, although each provides support at a different taxonomic level. Two data congruence tests (Templeton's test and the incongruence length difference test) suggested no significant heterogeneity between the data sets, and all data were combined in a total evidence analysis. The resulting well-resolved phylogeny suggests that Microcyphus, Amblypneustes, and Holopneustes are not monophyletic genera, and that Temnopleurus (Temnopleurus) and Temnopleurus (Toreumatica) are not closely related and should not be regarded as subgenera. In contrast to previous morphological analyses, Mespilia is found to be more closely related to Temnotrema and Toreumatica than it is to Microcyphus. The phylogeny was used to test a series of hypotheses about the evolution of developmental patterns. All species of Amblypneustes, Holopneustes, and Microcyphus are lecithotrophic, and many of these taxa are restricted to southern Australia. Planktotrophy is the ancestral condition for the temnopleurids, and the 11 instances of lecithotrophic nonplanktotrophy in this clade can be accounted for by a single developmental transition that occurred an estimated 4.4-7.4 million years ago, apparently before the migration of Microcyphus to southern Australia. The switch to a nonplanktotrophic mode of development is unidirectional with no evidence of reversals.

A molecular phylogeny of the Littorininae (Gastropoda: Littorinidae): unequal evolutionary rates, morphological parallelism, and biogeography of the Southern Ocean

A molecular phylogeny is presented for the subfamily Littorininae (including representatives of all subgeneric taxa and all members of a group of southern-temperate species formerly classified as 'Nodilittorina'), based on sequence data from two nuclear (18S rRNA, 28S rRNA) and two mitochondrial (12S rRNA, CO1) genes. The phylogeny shows considerable disagreement with earlier hypotheses derived from morphological data. In particular, 'Nodilittorina' is polyphyletic and is here divided into four genera (Echinolittorina, Austrolittorina, Afrolittorina new genus, and the monotypic Nodilittorina s.s.). The phylogenetic relationships of 'Littorina' striata have been controversial and it is here transferred to the genus Tectarius, a surprising relationship for which there is little morphological support. The relationships of the enigmatic Mainwaringia remain poorly resolved, but it is not a basal member of the subfamily. The two living species of Mainwaringia are remarkable for a greatly elevated rate of evolution in all four genes examined; it is suggested that this may be connected with their protandrous hermaphroditism, which is unique in the family. The molecular phylogeny provides a new framework for the adaptive radiation of the Littorininae, showing more frequent shifts between habitats and climatic regimes than previously suspected, and striking parallelism of morphological characters. The fossil record of littorinids is poor, but ages of clades are estimated using a calibration based on a Lower Eocene age of the genus Littoraria. Using these estimates, the antitropical distribution of Littorina and Afrolittorina is an ancient pattern of possibly Cretaceous age. The five members of Austrolittorina show a Gondwanan distribution in Australia, New Zealand, and South America. Based on the morphological uniformity within this clade, relatively recent (Plio-Pleistocene) trans-Pacific dispersal events seemed a likely explanation, as proposed for numerous other congeneric marine taxa. However, molecular estimation of ages of divergence suggest an initial vicariance between Australian and South American lineages at 40-73Ma, contemporary with the later stages of fragmentation of the Gondwanan supercontinent, followed by more recent (but still mid-Cenozoic) dispersal events across the Tasman Sea and the Pacific Ocean. Afrolittorina is another Cretaceous clade, now restricted to southern Africa and southern Australia, but divergence between these lineages (29-55Ma) post-dates Gondwanan fragmentation. Within both Austrolittorina and Afrolittorina all sister-species divergences are estimated to fall in the range 10-47Ma, so that there is no evidence for speciation events in the Plio-Pleistocene.

Phylogeny and classification of the Digenea (Platyhelminthes: Trematoda)

Complete small subunit ribosomal RNA gene (ssrDNA) and partial (D1-D3) large subunit ribosomal RNA gene (lsrDNA) sequences were used to estimate the phylogeny of the Digenea via maximum parsimony and Bayesian inference. Here we contribute 80 new ssrDNA and 124 new lsrDNA sequences. Fully complementary data sets of the two genes were assembled from newly generated and previously published sequences and comprised 163 digenean taxa representing 77 nominal families and seven aspidogastrean outgroup taxa representing three families. Analyses were conducted on the genes independently as well as combined and separate analyses including only the higher plagiorchiidan taxa were performed using a reduced-taxon alignment including additional characters that could not be otherwise unambiguously aligned. The combined data analyses yielded the most strongly supported results and differences between the two methods of analysis were primarily in their degree of resolution. The Bayesian analysis including all taxa and characters, and incorporating a model of nucleotide substitution (general-time-reversible with among-site rate heterogeneity), was considered the best estimate of the phylogeny and was used to evaluate their classification and evolution. In broad terms, the Digenea forms a dichotomy that is split between a lineage leading to the Brachylaimoidea, Diplostomoidea and Schistosomatoidea (collectively the Diplostomida nomen novum (nom. nov.)) and the remainder of the Digenea (the Plagiorchiida), in which the Bivesiculata nom. nov. and Transversotremata nom. nov. form the two most basal lineages, followed by the Hemiurata. The remainder of the Plagiorchiida forms a large number of independent lineages leading to the crown clade Xiphidiata nom. nov. that comprises the Allocreadioidea, Gorgoderoidea, Microphalloidea and Plagiorchioidea, which are united by the presence of a penetrating stylet in their cercariae. Although a majority of families and to a lesser degree, superfamilies are supported as currently defined, the traditional divisions of the Echinostomida, Plagiorchiida and Strigeida were found to comprise non-natural assemblages. Therefore, the membership of established higher taxa are emended, new taxa erected and a revised, phylogenetically based classification proposed and discussed in light of ontogeny, morphology and taxonomic history.

The phylogeny of the Schistosomatidae based on three genes with emphasis on the interrelationships of Schistosoma Weinland, 1858

Schistosomes are digenean flukes, parasitic of birds, mammals and crocodiles. The family Schistosomatidae contains species of considerable medical and veterinary importance, which cause the disease schistosomiasis. Previous studies, both morphological and molecular, which have provided a good deal of information on the phylogenetics of this group, have been limited in the number of species investigated or the type or extent of molecular data used. This paper presents the most comprehensive phylogeny to date, based on the sequences of 3 genes, complete ribosomal small subunit rRNA and large ribosomal subunit rRNA, and mitochondrial cytochrome oxidase 1, sequenced from 30 taxa including at least 1 representative from 10 of the 13 known genera of the Schistosomatidae and 17 of the 20 recognized Schistosoma species. The phylogeny is examined using morphological characters, intermediate and definitive host associations and biogeography. Theories as to the origins and spread of Schistosoma are also explored. The principal findings are that Ornithobilharzia and Austrobilharzia form a sister group to the Schistosoma; mammalian schistosomes appear paraphyletic and 2 Trichobilharzia species, T. ocellata and T. szidati, seem to be synonymous. The position of Orientobilharzia within the Schistosoma is confirmed, as is an Asian origin for the Schistosoma, followed by subsequent dispersal through India and Africa.

Phylogenetic relationships of Echinostoma Rudolphi, 1809 (Digenea: Echinostomatidae) and related genera re-assessed via DNA and morphological analyses

In order to investigate the relationships within the Echinostomatidae two data sets of gene sequences were analysed. The first consisted of all previously published ND1 sequences (20) together with 17 new sequences. The latter represented six species from the cosmopolitan genera Echinostoma, Echinoparyphium, Hypoderaeum and Isthmiophora. The second data-set of ITS sequences again included all previously published sequences (12) and three new sequences from species of Echinostoma, Echinoparyphium and Isthmiophora. All new isolates, as well as voucher material from five previously sequenced isolates, were identified on the basis of morphological characters. The phylogenetic trees inferred from the ND1 data set helped to clarify the generic affiliation of all isolates and confirmed the morphological identifications. The only exception was Echinoparyphium aconiatum, whose current position in the genus Echinoparyphium was not supported by the sequence data. Although the ITS data provided insufficient resolution for an unequivocal solution to the relationships within the genus Echinostoma, it supported the identification of Echinoparyphium ellisi and the distinct species status of three isolates of Echinostoma revolutum as predicted from the ND1 data.

Evidence for the co-existence of separate strains or species of Ligula in Lough Neagh, Northern Ireland

LIGULA (Cestoda: Pseudophyllidea) infections in gudgeon (Gobio gobio) and roach (Rutilus rutilus) differ markedly in the pathology that is observed in the host, particularly with respect to a tissue response and the extent of inhibition of gonadal development. The entire internal transcribed spacer (ITS) region (ITS-1, 5.8S and ITS-2) and the large subunit domains D1-D3 were sequenced and compared in parasites from these fish from Lough Neagh, Northern Ireland, together with a single specimen from minnow (Phoxinus phoxinus) from Wales. Sufficient differences were observed between parasites from R. rutilus and G. gobio to support the suggestion that they may represent different strains/species. In contrast, Ligula from P. phoxinus closely resembled those from R. rutilus. Ligula infections in G. gobio were recorded prior to the introduction of R. rutilus. The co-existence of separate strains or species of Ligula in Lough Neagh probably resulted from the introduction of R. rutilus to these waters, correlated with an increase in the number of great crested grebes (Podiceps cristatus).

Phylogenetics of the Monogenea--evidence from a medley of molecules

Nuclear ribosomal DNA sequences of Monogenea from both complete small and partial large (D1-D2) subunits were determined and added to previously published sequences in order to best estimate the molecular phylogeny of the group. A total of 35 ssrDNA, 100 D1 lsrDNA and 51 D2 lsrDNA monogenean sequences were used, representing a total of 27 families. From these sequences different data sets were assembled and analysed to make the best use of all available molecular phylogenetic information from the taxa. Maximum parsimony and minimum evolution trees for each data partition were rooted against published sequences from the Cestoda, forcing the Monogenea to appear monophyletic. There was broad agreement between tree topologies estimated by both methods and between genes. Well-supported nodes were restricted to deeply diverging major groupings and more derived taxa with the lsrDNA data but were at most nodes with ssrDNA. The Polyonchoinea showed the greatest resolution with a general pattern of ((Monocotylidae(Capsalidae(Udonellidae+Gyrodactylidea)))((Anoplodiscidae+Sundanonchidae)(Pseudomurraytrematidae+Dactylogyridae))). The Heteronchoinea readily split into the Polystomatoinea+Oligonchoinea, and Chimaericolidae and Hexabothriidae were successively the most basal of oligonchoinean taxa. Relationships within the Mazocraeidea, comprising 27 families of which 15 were sampled here, were largely unresolved and appear to reflect a rapid radiation of this group that is reflected in very short internal branches for ssrDNA and D1 lsrDNA, and highly divergent D2 lsrDNA. A reduced morphological matrix, employing only those families represented by molecules, contrasted sharply with respect to polyonchoinean interrelationships. Deep branches of the Heteronchoinea were similar for both classes of data but also showed that the interrelationships of the mazocraeidean families are labile and susceptible to sampling.

Gastropod Egg Capsules and Their Contents From Deep-Sea Hydrothermal Vent Environments

Egg capsules from three different prosobranch gastropods were retrieved from the Galapagos Rift and Juan de Fuca Ridge deep-sea hydrothermal vent fields. The morphology of these capsules and their excapsulated embryos and larvae are described and illustrated. Based on their capsule type and the protoconch morphology of their contained larvae, 29 lenticular capsules from the Galapagos Rift could be attributed to a provisionally described neogastropod turrid, Phymorhynchus sp. But 3 inflated, triangular capsules from the Galapagos Rift, and 56 different egg capsules from the Juan de Fuca Ridge, each shaped like an inflated pouch, could not be unambiguously assigned to a member of the known vent gastropod fauna. The mode of development and potential for dispersal is inferred from egg capsule type, the number of embryos per capsule, and protoconch characters comparable to those of confamilial shallow-water gastropods for which the type of development is known. These criteria and a comparison to the known juvenile shell morphology of Phymorhynchus sp., suggest that, after encapsulation, this species develops planktotrophically and is capable of long-range dispersal. Similar evidence suggests that the larvae contained in the inflated triangular capsules from the Galapagos Rift may also develop planktotrophically after hatching; but the larvae in the pouch-like egg capsules from the Juan de Fuca Ridge probably develop non-planktotrophically without a dispersal stage. These developmental patterns are characteristic of shallow-water members of the systematic groups to which these species belong, indicating, as previous studies have shown, that vent gastropods can persist in these patchy, ephemeral environments in the absence of unique adaptations allowing dispersal between active hydrothermal sites.