Evolutionary Relationships of Spirurina (Nematoda: Chromadorea: Rhabditida) with Special Emphasis on Dracunculoid Nematodes Inferred from SSU Rrna Gene Sequences Q
Total Page:16
File Type:pdf, Size:1020Kb
International Journal for Parasitology 36 (2006) 1067–1075 www.elsevier.com/locate/ijpara Evolutionary relationships of Spirurina (Nematoda: Chromadorea: Rhabditida) with special emphasis on dracunculoid nematodes inferred from SSU rRNA gene sequences q Martina Wijova´ a,b, Frantisˇek Moravec a, Alesˇ Hora´k a,b, Julius Lukesˇ a,b,* a Institute of Parasitology, Biology Centre of the Czech Academy of Sciences, Cˇ eske´ Budeˇjovice, Czech Republic b Faculty of Biological Sciences, University of South Bohemia, Cˇ eske´ Budeˇjovice, Czech Republic Received 12 January 2006; received in revised form 23 March 2006; accepted 19 April 2006 Abstract The analysis of 26 new small subunit rRNA sequences obtained from helminths that primarily parasitize fishes sampled from five continents provided well-supported trees, allowing us to study the phylogenetic relationships among spirurid nematodes. The analyses have shown that Dracunculoidea is a paraphyletic taxon and Anguillicolidae and Gnathostomatidae constitute the basal branch of the suborder Spirurina. The genera Philometra and Philometroides appear to be paraphyletic, while on the higher taxonomic level, good correlation between the morphology-based system and molecular data was observed. Neither co-evolution of the studied helminths with their hosts, nor phylogeographic pattern, are apparent in our dataset. Ó 2006 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved. Keywords: Nematoda; Spirurina; Philometridae; SSU rRNA; Phylogeny; Taxonomy 1. Introduction mented by the sequencing of other conserved genes, such as the large subunit rRNA, ITS region or elongation factor 1 Nematodes, which are typically considered a phylum, rep- (Gasser and Newton, 2000; Chilton et al., 2001). Although resent an extremely species-rich group. However, from the the recent dataset of nematode SSU rRNA genes comprises morphological perspective, they constitute a trait-poor more than 300 taxa (Blaxter, 2003), it remains strongly group due to a life-style that seems to preclude appendage biased towards some groups, while sampling remains poor evolution and cephalisation in both free-living and parasitic for several other important groups. One of the underrepre- species (Blaxter, 2003). The current nematode taxonomy is sented groups is the spirurine nematodes, an order/suborder thus primarily based on the morphology of the oesophagus, of obligatory parasites that includes helminths of cold- and male and female reproductive organs and life-cycle patterns warm-blooded terrestrial and aquatic vertebrates, including (Chabaud, 1975a,b; Anderson and Bain, 1976). An initial humans (Chabaud, 1975a,b; Anderson and Bain, 1976). framework for the molecular phylogeny of nematodes, based Opinions about the phylogenetic relationships among on the small subunit (SSU) rRNA dataset (Blaxter et al., the main groups of Spirurida in Anderson’s conception 1998; Dorris et al., 1999), was for a handful of taxa comple- (2000) have been developing gradually, which is reflected in taxonomic systems proposed by different authors (for q Note: The nucleotide sequences reported in this paper have been a recent review see De Ley and Blaxter, 2002). Spirurida deposited in the GenBank under the accession numbers AY852267– was assigned the rank of an order by Chitwood (1933) AY852269, DQ442659–DQ442679, DQ490223, DQ494195. * (subclass Secernentea, class Nematoda in his system), being Corresponding author. Address: Institute of Parasitology, Czech later subdivided into the suborders Camallanina and Spir- Academy of Sciences, Branisˇovska 31, 37005 Cˇ eske´ Budeˇjovice, Czech Republic. Tel.: +420 38 7775416; fax: +420 38 5310388. urina (Chabaud, 1974). The former suborder is comprised E-mail address: [email protected] (J. Lukesˇ). of the superfamilies Camallanoidea and Dracunculoidea, 0020-7519/$30.00 Ó 2006 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved. doi:10.1016/j.ijpara.2006.04.005 1068 M. Wijova´ et al. / International Journal for Parasitology 36 (2006) 1067–1075 whereas the suborder Spirurina represented a group of 10 for 1 min, 44 °C for 1 min and 72 °C for 2 min followed superfamilies containing hundreds of named species (Cha- by 24 cycles with the annealing temperature increased to baud, 1974). Several classifications, all based on morphol- 48 °C. The amplicons of expected size were gel-purified ogy and life-cycles, have been created by Chitwood using the Jetquick gel extraction kit (Genomed) and cloned (1933, 1950), Yamaguti (1961), Ivashkin et al. (1971) and into the Topo TA Cloning Vector (Invitrogen). Both Chabaud (1974, 1975a). In the most recent classification strands were sequenced using a Beckman Coulter automat- systems, based on morphology and biology, Spirurida ic sequencer. Internal oligonucleotides designed to match appeared as an order consisting of 28 families (Moravec the conserved regions were used to complete the sequences, et al., 1998; Anderson, 2000). However, these nematodes which have been assembled with Seqman (DNAStar). are classified within a newly defined suborder Spirurina that is placed in Rhabditida, Chromadorea, in the recent 2.2. Phylogenetic analysis classification system based on molecular data (De Ley and Blaxter, 2002). In contrast to previous systems, Spir- Dataset A was created from 29 SSU rRNA sequences of urina in this conception (subdivided into infraorders Asc- species listed in Tables 1 and 2. Plectus aquatilis (Plectidae) aridomorpha, Gnathostomatomorpha, Oxyuridomorpha, and Teratocephalus lirellus (Teratocephalidae) were chosen Rhigonematomorpha and Spiruromorpha, with Dracuncu- as outgroups according to De Ley and Blaxter (2002). For loidea as incertae sedis) is much broader, also including datasets B and C, enriched for sequences of Camallanidae former Ascaridida, Oxyurida and Rhigonematida. and Philometridae, respectively, different outgroups select- In the SSU rRNA-based trees published so far (Blaxter ed based on the analyses of dataset A have been used (see et al., 1998; Blaxter, 2003), the spirurine clade is robustly Section 3 for details). In order to include the partial philo- supported. This finding is complicated by the lack of any metrid SSU rRNA sequences available in GenBank for a sequences from several key representatives, which obscures number of philometrids, an additional 900 bp-long data- the relationships within the clade and with other clades. set D was constructed containing 18 specimens of Philo- Therefore, high expectations are placed on molecular phy- metridae and was rooted with Dracunculus spp. Multiple logeny. Despite the paramount importance of some spiru- alignments of different transition/transversion (Ts/Tv) rines for human health, this approach has not been used weights (1:1/2/3/5) were created with Clustal X 1.83 because of the apparent lack of interest among molecular (Thompson et al., 1997) and further edited using BioEdit biologists in these organisms (Lukesˇ et al., 2005). We have 7.0.4.1 (Hall, 1999) for all datasets. Ambiguously aligned undertaken sequence analysis of the SSU rRNA gene from regions and gaps were removed either with Bioedit or members of major groups within this zooparasitic taxon. Gblocks software (Castresana, 2000) prior to the analyses. The obtained dataset enabled us not only to examine phy- Maximum likelihood trees were calculated from all the logenetic relationships among spirurines and other nema- datasets and Ts/Tv setting under the GTR+C4+I model todes but also to address interesting problems, such as: of evolution using PHYML 2.4.4 (Guindon and Gascuel, (i) what is the ancestral invasion strategy of the spirurines? 2003); gamma shape parameter and proportion of invari- (ii) Did taxa with two-host life-cycles evolve from the base ants (PINVAR) were estimated from the dataset (see figure of the group? (iii) What is the pattern of definitive and/or captions for actual values). The model of evolution was intermediate host usage? (iv) Is there a phylogeographic chosen according to the Akaike criterion as implemented correlation or co-evolution with the host? Furthermore, in Modeltest 3.7 (Posada and Crandall, 1998). Maximum we propose the incorporation of the SSU rRNA-based likelihood (ML) and maximum parsimony (MP) bootstrap phylogenies into a comprehensive classification. support values were calculated by PHYML 2.4.4 and PAUP*4.0b10 (Swofford, 2002) from 1000 replications. 2. Materials and methods Bayesian posterior probabilities were assessed under the above described model with MrBayes 3.0b4 (Huelsenbeck 2.1. Organisms, DNA extraction and PCR and Ronquist, 2001) software (Bayesian inference—BI), where the Markov chain was set to 2·106 generations, Taxa sampled for phylogenetic analysis are listed in every 100th tree was sampled and the first 105 generations Table 1 and taxa for which sequences have been retrieved were omitted from phylogeny reconstruction. Additional from GenBank are listed in Table 2. Most specimens were LogDet analysis of all alignments of dataset A, run using stored in 70% ethanol for between several days up to a dec- PAUP*4.0b10, was performed to unmask possible cases ade prior to DNA isolation (Table 1). Genomic DNA was of the long-branch attraction phenomenon in topologies isolated using the Jetquick Tissue DNA kit (Genomed). obtained with the above-mentioned methods. About 10–50 ng of genomic DNA was used for PCR amplification of the SSU rRNA gene using oligonucleotide 3. Results primers D-1F (GCCTATAATGGTGAAACCGCGAAC) and D-1R (CCGGTTCA