Comparing Host and Parasite Phylogenies: Gyrodactylus Flatworms Jumping from Goby to Goby
Total Page:16
File Type:pdf, Size:1020Kb
Syst. Biol. 54(5):710–718, 2005 Copyright c Society of Systematic Biologists ISSN: 1063-5157 print / 1076-836X online DOI: 10.1080/10635150500221036 Comparing Host and Parasite Phylogenies: Gyrodactylus Flatworms Jumping from Goby to Goby TINE HUYSE AND FILIP A. M. VOLCKAERT Katholieke Universiteit Leuven, Laboratory of Aquatic Ecology, Ch. De Beriotstraat´ 32, B-3000 Leuven, Belgium; E-mail: [email protected] (T.H.) Abstract.—The combination of exceptionally high species diversity, high host specificity, and a complex reproduction sys- tem raises many questions about the underlying mechanisms triggering speciation in the flatworm genus Gyrodactylus. The coevolutionary history with their goby hosts was investigated using both topology- and distance-based approaches; phylogenies were constructed of the V4 region of the 18S rRNA and the complete ITS rDNA region for the parasites, and 12S and 16S mtDNA fragments for the hosts. The overall fit between both trees was significant according to the topol- ogy-based programs (TreeMap 1.0, 2.0β and TreeFitter), but not according to the timed analysis in TreeMap 2.0β and the distance-based method (ParaFit). An absolute timing of speciation events in host and parasite ruled out the possibility of synchronous speciation for the gill parasites, favouring the distance-based result. Based on this information together with the biological background of host and parasite, the following TreeMap solution was selected. The group of gill parasites evolved from a host switch from G. arcuatus, parasitizing the three-spined stickleback onto the gobies, followed by several host-switching events among the respective goby hosts. The timing of these events is estimated to date back to the Late Pleistocene, suggesting a role for refugia-mediated mixing of parasite species. In contrast, it is suggested that cospeciation in the fin-parasites resulted in several host-associated species complexes. This illustrates that phylogenetically conserved host-switching mimics the phylogenetic signature of cospeciation, confounding topology-based programs. [Coevolution; Gobies; Gyrodactylus; host-parasite interactions; host-switching; speciation.] One of the most significant radiations of platyhelminth scribed as well (Bakke et al., 2002). As such, all speciation fish parasites is documented in the monogenean “super- modes appear plausible in this group. genus” Gyrodactylus (Cribb et al., 2002). The 409 Gyro- In spite of these interesting characteristics, together dactylus species named so far (Harris et al., 2004) are with the fact that Gyrodactylus may pose serious threats to expected to represent only 2% to 10% of the estimated aquaculture (e.g., the highly pathogenic salmon parasite diversity (Bakke et al., 2002). Gyrodactylus displays the G. salaris), quantitative host-parasite cospeciation stud- widest host range of any monogenean family (19 orders ies are lacking (Bakke et al., 2002). According to Page of bony fish), encompassing both highly specific and gen- et al. (1996), the prerequisites for coevolutionary stud- eralist species. These remarkable facts can be linked with ies are (1) the availability of a sound alpha taxonomy of their even more remarkable reproduction mode. Gyro- both host and parasite, (2) robust phylogenies of hosts dactylus spp. contain a fully grown daughter in utero, and parasites, (3) although not essential, molecular phy- which in turn encloses a developing embryo, boxed in- logenies are preferably based on homologous characters, side one another like Russian dolls; asexual, partheno- (4) wide taxon sampling, and (5) quantitative compari- genetic and sexual reproduction alternate (Cable and son of host and parasite trees by means of explicit sta- Harris, 2002). In optimal conditions, a first-born daugh- tistical tests. The first four prerequisites are fulfilled in ter can be produced within 24 hours of her parent’s birth our model consisting of the Gyrodactylus parasites living (Cable and Harris, 2002). This can result in explosive pop- on gobies of the genus Pomatoschistus. (1) The various ulation growths, and the direct life cycle (only one host Gyrodactylus taxa have been described by morphologi- involved) allows them to invade new habitats with their cal, morphometric, and molecular methods (Huyse and hosts. All these factors may accelerate population differ- Malmberg, 2004; Huyse et al., 2004a), whereas the host entiation and thus speciation. species have been extensively studied by Miller (1986). On the one hand, the direct life cycle and the high host- They occur in the Atlantic and Mediterranean coasts of specificity enforce a tight relationship of a Gyrodacty- Europe; some goby species are sympatric and hybridiza- lus species and its host, promoting coevolution (Poulin, tion has been reported (Webb,1980). (2) The phylogenetic 1992; Kearn, 1994). On the other hand, the ability to relationships of Gyrodactylus have been studied using the produce an entire population from a single individual complete ITS rDNA region and the V4 region of the 18S increases the chance for sympatric speciation and speci- rRNA (Huyse et al., 2003); the phylogeny of the host has ation by host-switching (Brooks and McLennan, 1993). been reconstructed from 12S and 16S mtDNA fragments, They may also survive for some time independently and in order to compare homologous characters. (3) The of their host and a “swimming behavior” has been ob- ITS 1 locus has been sequenced as well (Huyse et al., served (Cable et al., 2002). These features add to the 2004b). (4) In order to minimize sampling bias, sampling colonizing capability of Gyrodactylus, illustrated by sev- has been extended over time and space throughout the eral switches between fish families (Harris, 1993; Zietara natural geographical distribution of the fish hosts (see and Lumme, 2002) and even fish orders (Huyse et al., below). 2003). However, phylogenetic radiations (diversification All respective Gyrodactylus species display phylo- constrained/triggered by host phylogeny) have been de- genetic host specificity towards gobies of the genus 710 2005 HUYSE AND VOLCKAERT—COEVOLUTION: DISTANCE- VERSUS TOPOLOGY-BASED PROGRAMS 711 Pomatoschistus (Huyse et al., 2003). They can be split into position for all sequences was compared usinga5%χ 2 two groups by their genetic and ecological parameters. test on the average composition (TREE-PUZZLE v. 5.0, Group A represents a monophyletic group of closely re- Schmidt et al., 2002). Plotting transitions and transver- lated and host-specific species, mainly infecting gills, sions against divergence of the complete data set using all belonging to the subgenus G.(Mesonephrotus). Host DAMBE v. 4.0.75 (Xia and Xie, 2001) did not show sat- specificity is lower in group B, which comprises species uration. The molecular-clock hypothesis was tested as- that are dominantly found on fin and skin and belong to suming the HKY+ model with the likelihood-ratio test G. (Paranephrotus). In this article, we fulfill the last prereq- (LRT) for the clock hypothesis implemented in TREE- uisite for coevolutionary studies by comparing host and PUZZLE. ModelTest v. 3.06 (Posada and Crandall, 1998) parasite phylogenies using a combination of statistical was used to estimate the optimal model of molecular approaches. In addition, we test whether the difference evolution in a likelihood-testing framework. These pa- in niche is reflected in the evolution and distribution on rameters were used in the maximum likelihood (ML) the host species: Are exposed fin parasites more subject method using PAUP∗;trees were statistically tested by to host-switching than the more “sheltered” gill para- calculating P values for the individual branches. This sites? Did the monophyletic group of gill parasites evolve likelihood-ratio test assesses whether branch-lengths are through cospeciation with the host, as generally expected significantly different from zero by optimizing all branch for host-specific parasites (Poulin, 1992; Kearn, 1994)? lengths under the constraint that one of the branches is Alternatively, did the group of the more generalist fin zero, for each branch tested. For maximum parsimony parasites evolve through repeated host-switching events (MP) the exhaustive search method was performed us- instead of cospeciation? ing the branch and bound algorithm. In these analyses gaps were treated both as missing data and fifth character because indels might be phylogenetically informative; all MATERIAL AND METHODS sites were equally weighted. The neighbor-joining (NJ) Host and Parasite Data search was conducted (1000 replicates of tree-bissection During 1999 to 2002, 15 Gyrodactylus taxa were reconnection branch swapping) from a matrix of ML ge- collected on Gobiusculus flavescens and several Po- netic distances calculated under the optimized model. matoschistus species (Gobiidae, Teleostei) along the ML branch lengths have been generated under the spec- northeastern Atlantic continental shelf, the Baltic, Adri- ified model and topology, with and without a molec- atic, and Mediterranean Seas (see Huyse et al., 2003). ular clock assumption. Bayesian inference (BI) using Two Gyrodactylus parasites of the three-spined stickle- MrBayes v. 3.0b4 (Huelsenbeck and Ronquist, 2001) was back, Gasterosteus aculeatus, were collected as outgroup performed, using the model obtained by ModelTest. Pos- species (G. arcuatus and G. branchicus). Phylogenetic re- terior probabilities were estimated over 2.0 × 106 gen- lationships of the Gyrodactylus parasites were inferred erations,