- Research article
- Open Access
Phylogenetic placement of the enigmatic parasite, Polypodium hydriforme, within the Phylum Cnidaria
BMC Evolutionary Biology volume 8, Article number: 139 (2008)
The Erratum to this article has been published in BMC Evolutionary Biology 2009 9:165
Polypodium hydriforme is a parasite with an unusual life cycle and peculiar morphology, both of which have made its systematic position uncertain. Polypodium has traditionally been considered a cnidarian because it possesses nematocysts, the stinging structures characteristic of this phylum. However, recent molecular phylogenetic studies using 18S rDNA sequence data have challenged this interpretation, and have shown that Polypodium is a close relative to myxozoans and together they share a closer affinity to bilaterians than cnidarians. Due to the variable rates of 18S rDNA sequences, these results have been suggested to be an artifact of long-branch attraction (LBA). A recent study, using multiple protein coding markers, shows that the myxozoan Buddenbrockia, is nested within cnidarians. Polypodium was not included in this study. To further investigate the phylogenetic placement of Polypodium, we have performed phylogenetic analyses of metazoans with 18S and partial 28S rDNA sequences in a large dataset that includes Polypodium and a comprehensive sampling of cnidarian taxa.
Analyses of a combined dataset of 18S and partial 28S sequences, and partial 28S alone, support the placement of Polypodium within Cnidaria. Removal of the long-branched myxozoans from the 18S dataset also results in Polypodium being nested within Cnidaria. These results suggest that previous reports showing that Polypodium and Myxozoa form a sister group to Bilateria were an artifact of long-branch attraction.
By including 28S rDNA sequences and a comprehensive sampling of cnidarian taxa, we demonstrate that previously conflicting hypotheses concerning the phylogenetic placement of Polypodium can be reconciled. Specifically, the data presented provide evidence that Polypodium is indeed a cnidarian and is either the sister taxon to Hydrozoa, or part of the hydrozoan clade, Leptothecata. The former hypothesis is consistent with the traditional view that Polypodium should be placed in its own cnidarian class, Polypodiozoa.
Polypodium hydriforme is an endocellular parasite whose unusual life cycle, peculiar morphology, and high rates of DNA evolution, have led to much controversy regarding its phylogenetic position within metazoans [1–5]. Polypodium spends most of its life inside the oocytes of acipenseriform fishes (sturgeons and paddlefish). During this time, Polypodium develops from a binucleate cell into an inside-out planuliform larva and then into an elongate inside-out stolon; the epidermal cell layer is located internal to the body and the gastrodermis is located externally [6–8]. The embryo, larva and stolon are surrounded by a protective polyploid cell, which also functions in digestion . Just prior to host spawning, Polypodium everts to the normal position of cell layers, revealing tentacles scattered along the stolon. During eversion, the yolk of the host oocyte fills the gastral cavities of the parasite, supplying the future free-living stage with nutrients [6, 7]. Finally, upon emerging from the host egg in fresh water, the free-living stolon (Figure 1A) fragments into individual medusoid-like forms (Figure 1B) that go on to multiply by means of longitudinal fission, form sexual organs, and ultimately infect host fish with their gametophores [6–9].
Two conflicting hypotheses have been proposed regarding the phylogenetic placement of Polypodium. The first, more traditional hypothesis is that Polypodium is a cnidarian. Some have suggested it is nested within a derived group of hydrozoans, the Narcomedusae [10–13] or the cnidarian class Scyphozoa ; while others have suggested it belongs to a separate cnidarian class, Polypodiozoa [1, 15, 16]. The assignment of Polypodium to Cnidaria is based primarily on morphological evidence, most notably the fact that Polypodium possesses nematocysts [17, 18], the stinging structures characteristic of all cnidarians. In addition, the presence of tentacles and overall body-plan organization of Polypodium are reminiscent of cnidarians, although it is unclear if the adult free-living stage is homologous to a polyp or medusa stage. This hypothesis is supported by a cladistic analysis of small subunit nuclear ribosomal DNA (18S rDNA) sequences in conjunction with morphological characters (including nematocysts) . In this study, Polypodium falls within the medusozoan clade of cnidarians, although the non-cnidarian placozoan, Trichoplax [19, 20], also fell within this clade, rendering Cnidaria paraphyletic.
The second hypothesis is that Polypodium is the sister taxon to Myxozoa, a diverse group of parasites in aquatic animals, and that Polypodium + Myxozoa is the sister group to Bilateria [2–4]. This hypothesis is derived from cladistic analyses utilizing 18S rDNA sequences [2–4]. However, because Polypodium and myxozoans have unusually high divergence rates in their 18S rDNA sequences, these cladistic analyses have been criticized by a number of authors who suggest that the data might be unduly affected by long-branch attraction (LBA) [5, 21, 22]. Despite some attempts to overcome the effects of LBA through the use of a maximum likelihood (ML) approach [21–23] and pruning long branches [5, 22], these results have been largely silent on the placement of Polypodium. For instance, Kim et al.  applied a maximum likelihood approach to 18S rDNA sequence data and found that myxozoans and Polypodium did not group together. Instead, Polypodium was part of an unresolved polytomy that included several cnidarian lineages and Trichoplax, as well as myxozoans + Bilateria. Most recently, Jimenez-Guri et al.  utilized multiple protein-coding gene sequences in a ML analysis and found the myxozoan, Buddenbrockia plumatellae nested within cnidarians. Unfortunately, this study had relatively limited sampling of cnidarians and did not include Polypodium.
In an attempt to resolve this controversy, we sequenced an additional marker in Polypodium, a partial gene sequence of the large nuclear ribosomal unit (28S rDNA), and greatly expanded the taxonomic sampling of cnidarian sequences. Using this approach, we provide evidence that Polypodium is nested within Cnidaria and does not group with myxozoans.
All taxa used in this study are arranged taxonomically in Table 1. 155 sequences were obtained from GenBank. 45 new cnidarian sequences for 18S and 59 for 28S (including 2 new 18S and 2 new partial 28S from Polypodium taxa) were generated for this study and deposited in GenBank (see Table 1 for accession numbers). Polypodium hydriforme sequences were obtained from both North American and Eurasian hosts. Eurasian samples were collected from two individuals of Acipenser ruthenus. North American samples were collected from Polyodon spathula and Scaphirhynchus platorynchus. This is the first reported presence of Polypodium infection in Scaphirhynchinae. While Polypodium was recovered from the oocytes of S. platorynchus, the sample from which we extracted sequence data was found externally attached to its presumed host. More specific collection data for Polypodium specimens are associated with each sequence submitted to GenBank (see Table 1 for accession numbers).
All Polypodium sequences were newly generated for this study. We did not include the previously published 18S Polypodium sequence (GenBank accession number U37526) because of concern over the quality of the sequence which included a number of ambiguities. Furthermore, while the two new Polypodium 18S sequences (from hosts Acipenser ruthensus and Polyodon spathula) differed from each other by a total of 8 sites they differed from #U37526 by 77 and 83 sites respectively. These differences included a large number of insertions and deletions. The two new 28S sequences (from hosts Acipenser ruthensus and Scaphirhynchus platorynchus) only differed from each other by 2 sites.
Position of Polypodium
The complete combined dataset of 18S rDNA and partial 28S rDNA contains 4842 characters, 2901 of which are variable and 2124 parsimony informative. Both the ML and parsimony topologies reconstructed from the combined dataset suggest that Polypodium is nested within a monophyletic Cnidaria, and myxozoans are the sister taxon to bilaterians (Figure 2). The ML bootstrap values supporting a monophyletic Cnidaria (including Polypodium), a monophyletic Medusozoa (including Polypodium) and the Polypodium + hydrozoan clade are 73, 67 and 73 respectively (Figure 2A, and Additional file 1). Parsimony analysis of the combined dataset differs from that of ML in that Polypodium is nested within a group of hydrozoans, the leptothecates (Figure 2B). The parsimony bootstrap values supporting a monophyletic Cnidaria and Hydrozoa, with Polypodium nested within these clades are 50 and 51 respectively (Figure 2B). The clade nested within hydrozoans, that includes Polypodium + leptothecates is weakly supported in the sub-sampling tests with a bootstrap value of less than 50.
The analyses using partial 28S rDNA sequences alone (129 sampled taxa) contains 1756 characters, 1196 of which are parsimony informative. The ML topology using this dataset reveals Polypodium nested within Cnidaria, specifically within leptothecate hydrozoans, (Additional file 2). This analysis however fails to recover a monophyletic Cnidaria, as the anthozoans are placed outside the Cnidaria + Bilateria clade. Analysis of the 18S rDNA dataset alone (132 taxa, 3038 characters, 1469 parsimony informative) under both optimality criteria conflicts with the combined and partial 28S topologies. The 18S rDNA topology for both criteria place Polypodium at the base of Bilateria (Figure 3A, Additional files 3, 4 and 5). However, the ML topology also reflects a sister relationship between Polypodium and myxozoans (Figure 3A and Additional file 3A) while the parsimony topology does not (Additional files 4 and 5). Moreover, under parsimony criteria the position of myxozoans is dependent upon how gaps are coded: if gaps are coded as a fifth character state, myxozoans are placed as a highly derived clade of bilaterians (Additional file 4); if gaps are coded as missing, myxozoans are placed as sister to all metazoans (Additional file 5). The 18S analysis showing placement of Polypodium with Bilateria, and more specifically as sister to myxozoans, is consistent with previously reported studies using the same marker [2–4], but raises similar concerns of long-branch attraction .
Test of long-branch attraction
Myxozoans and Polypodium have unusually high rates of evolution in their 18S and 28S rDNA sequences relative to the other sampled taxa. To investigate the influence of myxozoans on the placement of Polypodium, we removed the myxozoans from our three datasets and re-ran each analysis. Under the ML analysis of 18S rDNA, the removal of myxozoans results in the placement of Polypodium nested within Cnidaria (Figure 3B and Additional file 3B). This result suggests that the placement of Polypodium at the base of bilaterians in the 18S analysis (Figure 3A) was indeed an artifact of LBA. The placement of Polypodium within Cnidaria was not effected by the removal of myxozoans in the 28S (Additional file 6) and combined datasets (Additional file 7).
To investigate the possible role of LBA on myxozoan placement, we removed Polypodium from the combined ML analyses and found that it did not affect the position of Myxozoa at the base of the Bilateria (not shown). Given that bilaterians also form long branches, we tried removing all bilaterian sequences in the combined ML analysis. This resulted in a Myxozoa + Polypodium clade nested within Cnidaria (not shown). However, when Polypodium and bilaterians were removed, myxozoans fell outside the cnidarians (not shown). Similar effects of myxozoan placement to long-branches were also found in parsimony analyses of the combined dataset (not shown).
Polypodiumis a cnidarian
Our metazoan dataset of 18S and partial 28S rDNA sequences, with a large taxonomic sample of cnidarians, places Polypodium within a monophyletic Cnidaria. This accords with the fact that Polypodium possesses nematocysts [17, 18] and a cnidarian-like body plan [7–9, 12]. The precise placement of Polypodium within Cnidaria is less certain. The ML combined analysis places Polypodium as sister to Hydrozoa (Figure 2A), a hypothesis consistent with the suggestion that Polypodium be considered a separate class of cnidarians, Polypodiozoa . By contrast, the combined parsimony analysis (Figure 2B) and the ML analyses of 28S alone (Additional file 2 and 6) place Polypodium within the hydrozoan clade Leptothecata. Given that leptothecates have relatively high rates of evolution within hydrozoans, one possible explanation for the conflicting hypotheses is that the placement of Polypodium within leptothecates is an artifact of LBA and that the combined data, in conjunction with the ML approach (Figure 2A), overcame this localized LBA artifact.
Evolution of Polypodiumlife-history characters
Although the fresh water habitat of Polypodium is unusual for cnidarians, it is not unheard of, especially within hydrozoans. For instance, the model organism Hydra and the jellyfish Craspedacusta are both exclusively fresh-water hydrozoans. Hydra and Craspedacusta are distantly related  and our analyses do not indicate a close phylogenetic affinity of Polypodium to either of the clades containing these taxa. Thus, it appears that in the evolution of cnidarians, invasion to fresh-water habitats has happened at least three separate times.
Although Polypodium is the only known intracellular cnidarian parasite, other cnidarians have adopted parasitic life-styles [11, 26–29]. For example, parasites belonging to the Narcomedusae (Hydrozoa) have been reported to live in the stomach cavities of other narcomedusae [11, 27] and anthomedusae . In addition, the sea anemone Edwardsiella lineata parasitizes the stomach cavity of the ctenophore Mnemiopsis leidyi  and the anemone Peachia quinquecapitata is reported to parasitize the stomachs of hydromedusa .
Effects of long-branch attraction
The well-documented effects of long-branch attraction artifacts (reviewed in Bergsten ) are particularly concerning when investigating relationships amongst early-diverging metazoans, where rates between lineages vary greatly . Suggestions for avoiding LBA artifacts include choice of appropriate markers [31, 32], increased taxonomic sampling to effectively break up long branches [33, 34] and utilization of best-fit models that incorporate rate variation [21–23]. Previous conflicting reports that show Polypodium and myxozoans form a sister taxon to Bilateria [2–4] can be explained by limited taxon sampling and an inadequate number of informative characters in their analyses, both of which confound long-branch problems. In this study, the increased taxonomic sampling of cnidarians and the addition of 28S rDNA sequence data proved critical to placing the highly divergent Polypodium taxon within Cnidaria. The choice of optimality criteria (ML vs. parsimony) both supported Polypodium as a cnidarian but did affect the placement within Cnidaria.
Our analyses are inconclusive in the placement of Myxozoa within metazoans. We found that myxozoans consistently grouped with long-branched taxa and that removal of long-branches resulted in myxozoans being placed to the next longest branch. For example myxozoans group with Polypodium in the absence of Bilateria and group with Bilateria in the absence of Polypodium (not shown).
Jimenez-Guri et al.  sampled the myxozoan, Buddenbrockia, and found it to fall within Cnidaria, as the sister group to two hydrozoan representatives and a single scyphozoan. Previous studies have suggested a sister group relationship between cnidarians and myxozoans [2–4], and some morphological evidence has been used to support this view . Although our present study does not support this relationship, further investigation is merited. Myxozoans are a highly diverse group (reviewed in Kent et al. ) that comprise two clades, the Myxosporea and the Malacosporea . We were only able to include 28S rDNA sequences from myxosporeans, although the malacosporean Buddenbrockia was included in our 18S analysis and found to group with other myxozoans and outside of Cnidaria. Future studies with a comprehensive sampling of myxozoans together with Polypodium, in a dataset that includes a large taxonomic sampling of cnidarians, should shed further light on the relationships between myxozoans and Polypodium.
Although previous molecular phylogenetic hypotheses conflicted with the traditional interpretation of cnidarian affinity for Polypodium, the molecular evidence we present, using an augmented dataset, ultimately confirms and reconciles this traditional hypothesis and suggests that Polypodium is indeed a cnidarian. This study also reaffirms the importance to large taxonomic sampling and inclusion of additional informative characters for avoiding long-branch attraction artifacts.
DNA isolation, amplification and sequencing
Genomic DNA was extracted using Qiagen DNeasy kits according to manufacturer's protocol (QIAGEN Inc., Mississauga, ON) or a standard phenol/chloroform protocol. The latter method involved tissue digestion with proteinase K (20 mg/ml) in a lysis buffer (20 mM Tris-CL pH 8.0, 5 mM EDTA pH 8.0, 400 mM NaCl, 2%SDS), extraction with phenol/chloroform (1:1), precipitation with 2.5 vol. 95% EtOH, and elution in TE or H2O.
An approximately 1.8 kb portion of the gene coding for 18S was amplified and sequenced with universal eukaryotic primers as described by Medlin et al. , with the annealing temperature modified to 57°C. With the exception of Polypodium samples, a nearly complete, roughly 3 kb portion of the gene coding for 28S was amplified and sequenced with an approach modified from that reported in Collins et al. . 28S was directly amplified in two fragments with combinations of primers F63mod+R2077sq and F1379+R3264 from Medina et al.  or newly developed medusozoan specific primers F97+R2084 and F1383+R3238 (F97: CCYYAGTAACGGCGAGT, R2084: AGAGCCAATCCTTTTCC, F1383: GGACGGTGGCCATGGAAGT, and R3238: SWACAGATGGTAGCTTCG). Amplifications of 28S were conducted with the following thermal profile: 4 minutes at 94°C; 30 cycles of 30 seconds at 94°C, 1 minute at 45°C, and 3 minutes at 72°C; and 10 minutes at 72°C. For Polypodium, a portion of the 5' end of 28S (approx. 0.8–1.0 kbps) was amplified using two universal metazoan primers (fw1and rev2) as reported by Sonnenberg et al. . Sequencing was carried out using amplification primers and F635sq and R635sq from Medina et al. .
All gene fragments were purified and sequenced by Cogenics, Inc. (Houston, TX) and assembled and edited using Sequencher v4.5 (Gene Code Co., 2005). Sequences for each marker were aligned using the program MUSCLE . The 28S sequence alignment was then trimmed to reflect only that region which included sequence data for Polypodium. This trimmed 28S dataset was analyzed separately and used in conjunction with the complete 18S sequences to create the combined dataset.
Phylogenetic analyses were performed using both maximum likelihood (ML) and parsimony criteria. ML searches were performed using GARLI v0.951.OsX-GUI  under an assumed GTR model with rates estimated from the data. The assumed model of nucleotide substitution was selected by using the Akaike Information Criterion (AIC) as implemented in ModelTest . Each run was repeated 10 times from random starting trees using default termination conditions. Each run gave identical topologies and similar likelihood scores. 100 bootstrap replications were run in GARLI v0.951.0sX-GUI  under the same parameters.
To assess the effect that omitting length-variable regions has on topology, we removed these regions from the combined dataset, using the less stringent settings of Gblocks . This dataset contained 126 metazoan taxa, 2415 characters, 1391 of which are parsimony informative. We found that removal of length-variable regions had no effect on the placement of Polypodium and minimal effect on overall topology in our combined ML analyses (Additional file 8). Therefore we performed all other analyses with the complete datasets, including the more variable regions.
Parsimony analyses were performed using TNTv.1.1 . Separate tree searches were performed with gaps coded as missing and gaps coded as a fifth state. However, with one exception (see results for myxozoan placement with 18S data) there was no significant difference in topology. Numerous search methods available in TNT were utilized to search the tree space but the following approach was found to consistently recover trees with minimum lengths from our datasets. The implemented search was a driven new technology search with a random seed of 0 (where 0 = time). Default settings for sectorial searches (RSS and CSS) and tree fusing were used , with 5 replicates per repetition, and a requirement that the global optimum be found 20 times. TBR branch swapping was performed on the resulting trees and a strict consensus was calculated. TNT was used to calculate standard bootstrap values (1000 replicates). Alignments and trees for 18S, 28S and combined datasets have been submitted to TreeBASE http://www.treebase.org/treebase/index.html.
Raikova EV: On the systematic position of Polypodium hydriforme Ussov, (Coelenterata). Porifera and Cnidaria Contemporary state and perspectives of investigations. Edited by: Koltum VM, Stepanjants SD. 1988, Leningrad: Zoological Institute of Academy of Sciences of USSR, 116-122.
Siddall ME, Martin DS, Bridge D, Cone DM, Desser SS: The demise of a phylum of protests: Phylogeny of Myxozoa and other parasitic Cnidaria. J Parasitol. 1995, 81: 961-967. 10.2307/3284049.
Siddall ME, Whiting MF: Long-branch abstractions. Cladistics. 1999, 15: 9-24. 10.1111/j.1096-0031.1999.tb00391.x.
Zrzavý J, Hypša V: Myxozoa, Polypodium, and the origin of the Bilateria: The phylogenetic position of "Endocnidozoa" in the light of the rediscovery of Buddenbrockia. Cladistics. 2003, 19 (3): 164-169. 10.1016/S0748-3007(03)00007-0.
Hanelt B, Van Schyndel D, Adema CM, Lewis L, Loker ES: The phylogenetic position of Rhopalura ophiocomae (Orthonectida) based on 18S ribosomal DNA sequence analysis. Mol Biol Evol. 1996, 13: 1187-1191.
Raikova EV: Life cycle and systematic position of Polypodium hydriforme Ussov (Coelenterata), a cnidarian parasite of the eggs of Acipenseridae. Publ Seto Mar Biol Lab. 1973, 20: 165-173.
Raikova EV: Morphology, ultrastructure and development of the parasitic larva and its surrounding trophamnion of Polypodium hydriforme Ussov (Coelenterata). Cell Tissue Res. 1980, 206 (3): 487-500. 10.1007/BF00237977.
Raikova EV: Life cycle, cytology, and morphology of Polypodium hydriforme, a coelenterate parasite of the eggs of acipenseriform fishes. J Parasitol. 1994, 80 (1): 1-22. 10.2307/3283338.
Raikova EV, Suppes VC, Hoffmann GL: The parasitic coelenterate, Polypodium hydriforme Ussov, from the eggs of the American acipenseriform Polyodon spathula. J Parasitol. 1979, 65 (5): 804-810. 10.2307/3280366.
Berrill NJ: Development and medusa-bud formation in the Hydromedusae. Quart Rev Biol. 1950, 25: 292-316. 10.1086/397707.
Bouillon J: Considérations sur les dévéloppement des narcoméduses et sur leur position phylogenétique. Indo-Malayan Zool. 1987, 4: 189-278.
Hyman L: The invertebrates. I. Protozoa through Ctenophora. New York and London: McGraw-Hill, 1940
Schimkevitch VM: Essay on contemporary state of the problem of the development of Hydrozoa. Vestnik Estestvoznanija. 1890, 4: 171-176.
Lipin A: Geschlechtliche Form, Phylogenie und systematische Stellung von Polypodium hydriforme Ussov. Zool Jahrb Anat. 1925, 47: 541-635.
Bouillon J, Medel MD, Pagès F, Gili JM, Boero F, Gravili C: Fauna of the Mediterranean Hydrozoa. Scientia Marina. 2004, 68 (Suppl 2): 5-438.
Bouillon J, Gravili c, Pagès F, Gili JM, Boero F: An introduction to Hydrozoa. 2006, Paris: Publications Scientifiques du Muséum, Paris
Raikova EV: Fine structure of the nematocytes of Polypodium hydriforme Ussov (Cnidaria). Zoologica Scripta. 1990, 19 (1): 1-11. 10.1111/j.1463-6409.1990.tb00236.x.
Ibragimov A, Raikova E: Nematocysts of Polypodium hydriforme, a cnidarian parasite of acipenseriform fishes. Hydrobiologia. 2004, 531: 165-171. 10.1007/s10750-004-2651-y.
Dellaporta SL, Xu A, Sagasser S, Jakob W, Moreno MA, Buss LW, Schierwater B: Mitochondrial genome of Trichoplax adhaerens supports Placozoa as the basal lower metazoan phylum. Proc Natl Acad Sci USA. 2006, 103 (23): 8751-8756. 10.1073/pnas.0602076103.
Collins AG: Phylogeny of Medusozoa and the evolution of cnidarian life cycles. J Evol Biol. 2002, 15 (3): 418-432. 10.1046/j.1420-9101.2002.00403.x.
Huelsenbeck JP: Is the Felsenstein zone a fly trap?. Syst Biol. 1997, 46: 247-264. 10.2307/2413636.
Kim JH, Kim W, Cunningham CW: A new perspective on lower metazoan relationships from 18S rDNA sequences. Mol Biol Evol. 1999, 16: 423-427.
Cunningham CW, Zhu H, Hillis DM: Best-fit maximum likelihood models for phylogenetic inference: Empirical tests with known phylogenies. Evolution. 1998, 52: 978-987. 10.2307/2411230.
Jimenez-Guri E, Philippe H, Okamura B, Holland PWH: Buddenbrockia is a cnidarian worm. Science. 2007, 317 (5834): 116-118. 10.1126/science.1142024.
Collins AG, Schuchert P, Marques AC, Jankowski T, Medina M, Schierwater B: Medusozoan phylogeny and character evolution clarified by new large and small subunit rDNA data and an assessment of the utility of phylogenetic mixture models. Syst Biol. 2006, 55 (1): 97-115. 10.1080/10635150500433615.
Osborn DA: Cnidarian "parasites" on Solmissus incisa, a Narcomedusa. Sci Mar. 2000, 64: 157-163.
Pagès F, Corbera J, Lindsay D: Piggybacking pycongonids and parasitic narcomedusae on Pandea rubra (Anthomedusae, Pandeidae). Plankton Benthos Res. 2007, 2 (2): 83-90. 10.3800/pbr.2.83.
Bumann D, Puls G: Infestation with larvae of the sea anemone Edwardsia lineata affects nutrition and growth of the ctenophore Mnemiopsis leidyi. Parasitology. 1996, 113: 123-128.
Spaulding JG: The life cycle of Peachia quinquecapitata, an anemone parasitic on medusae during its larval development. Biol Bull. 1972, 143 (2): 440-453. 10.2307/1540065.
Bergsten J: A review of long-branch attraction. Cladistics. 2005, 21: 163-193. 10.1111/j.1096-0031.2005.00059.x.
Rokas A, King N, Finnerty J, Carroll SB: Conflicting phylogenetic signals at the base of the metazoan tree. Evol Dev. 2003, 5 (4): 346-359. 10.1046/j.1525-142X.2003.03042.x.
Xiang QY, Moody ML, Soltis DE, Fan C, Soltis PS: Relationships within Cornales and circumscription of Cornaceae – matK and rbcL sequence data and effects of outgroups and long branches. Mol Phyl Evol. 2002, 24: 35-57. 10.1016/S1055-7903(02)00267-1.
Hillis DM: Taxonomic sampling, phylogenetic accuracy, and investigator bias. Syst Biol. 1998, 47: 3-8. 10.1080/106351598260987.
Zwickl DJ, Hillis DM: Increased taxon sampling greatly reduces phylogenetic error. Syst Biol. 2002, 51: 588-598. 10.1080/10635150290102339.
Raikova EV: Cytomorphological characters of Polypodium hydriforme and problems of myxozoan and cnidarian phylogeny. Tsitologiia. 2005, 47 (10): 933-939.
Kent ML, Andree KB, Bartholomew JL, El-Matbouli M, Desser SS, Devlin RH, Feist SW, Hedrick RP, Hoffmann RW, Khattra J, Hallett SL, Lester RJG, Longshaw M, Oswaldo P, Siddall ME, Xiao C: Recent advances in our knowledge of the Myxozoa. J Eukaryot Microbiol. 2001, 48 (4): 395-413. 10.1111/j.1550-7408.2001.tb00173.x.
Canning EU, Curry A, Feist SW, Longshaw M, Okamura B: A new class and order of myxozoans to accommodate parasites of bryozoans with ultrastructural observations on Tetracapsula bryosalmonae (PKX Organism). J Eukaryot Microbiol. 2000, 47 (5): 456-468. 10.1111/j.1550-7408.2000.tb00075.x.
Medlin LH, Elwood HJ, Stickel S, Sogin ML: The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene. 1988, 71: 491-499. 10.1016/0378-1119(88)90066-2.
Medina M, Collins AG, Silberman JD, Sogin ML: Evaluating hypotheses of basal animal phylogeny using complete sequences of large and small subunit rRNA. Proc Natl Acad Sci USA. 2001, 98 (18): 9707-9712. 10.1073/pnas.171316998.
Sonnenberg R, Nolte A, Tautz D: An evaluation of LSU rDNA D1-D2 sequences for their use in species identification. Front Zool. 2007, 4 (1): 6-10.1186/1742-9994-4-6.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. 2006, Austin: The University of Texas
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 2000, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-552.
Goloboff P, Farris J, Nixon K: TNT: Tree Analysis Using New Technology. 2003, Tucuman, Argentina: Published by the authors, 1.1
Goloboff P: Analyzing large data sets in reasonable times: Solutions for composite optima. Cladistics. 1999, 14: 415-428. 10.1111/j.1096-0031.1999.tb00278.x.
We thank S. Ash Bullard, Jan Dean, Bobby Reed and Ron Bruch for contributions of Polypodium specimens. We also thank Peter Schuchert, Casey Dunn and Michael Dawson for contributions of other cnidarian specimens. NME acknowledges the instructors at the 2007 Ohio State Cladistics workshop, especially P. Goloboff, for help with TNT. We thank M. Holder for discussions, B. Bentlage, B. Lieberman, A. Nawrocki, and five anonymous reviewers for helpful comments on earlier versions of this manuscript. This work was supported by grants from NSF AToL EF-0531779 (to PC and AGC) and support for AL from NSF (PEET DEB-9978086) and FAPESP (06/02960-8/05821-9/60327-0).
This study was inspired by the work of EVR and originally conceived by EVR, AL and AGC. NME performed most of the data collection and submission of the new sequences. NME and PC performed the analyses. PC took the lead in organizing the study and drafting the manuscript with contributions from NME. All other authors provided comments and suggestions and approved the final manuscript.
An erratum to this article is available at http://dx.doi.org/10.1186/1471-2148-9-165.