Skip Navigation

Systematic Biology 2006 55(1):97-115; doi:10.1080/10635150500433615
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Collins, A. G.
Right arrow Articles by Schierwater, B.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Collins, A. G.
Right arrow Articles by Schierwater, B.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2006 Society of Systematic Biologists

Medusozoan Phylogeny and Character Evolution Clarified by New Large and Small Subunit rDNA Data and an Assessment of the Utility of Phylogenetic Mixture Models

Edited by Tim Collins: Associate Editor

Allen G. Collins1, Peter Schuchert2, Antonio C. Marques3, Thomas Jankowski4, Mónica Medina5 and Bernd Schierwater6

1 NMFS, National Systematics Laboratory, National Museum of Natural History, MRC-153, Smithsonian Institution P.O. Box 37012, Washington, DC, 20013–7012, USA; E-mail: CollinsA{at}SI.edu (A.G.C.)
2 Muséum d'Histoire Naturelle 1 route Malagnou, CH-1211, Genève, Switzerland
3 Department of Zoology, Institute of Biosciences, University of São Paulo São, Paulo, Brazil
4 W & T, Swiss Federal Institute of Aquatic Science and Technology (EAWAG-ETH) Ueberlandstrasse 133, CH-8600, Duebendorf, Switzerland
5 DOE Joint Genome Institute 2800 Mitchell Drive, Walnut Creek, California, 94598, USA
6 ITZ, Ecology & Evolution, TiHo Hannover Bünteweg 17d, D-30559, Hannover, Germany


    Abstract
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 
A newly compiled data set of nearly complete sequences of the large subunit of the nuclear ribosome (LSU or 28S) sampled from 31 diverse medusozoans greatly clarifies the phylogenetic history of Cnidaria. These data have substantial power to discern among many of the competing hypotheses of relationship derived from prior work. Moreover, LSU data provide strong support at key nodes that were equivocal based on other molecular markers. Combining LSU sequences with those of the small subunit of the nuclear ribosome (SSU or 18S), we present a detailed working hypothesis of medusozoan relationships and discuss character evolution within this diverse clade. Stauromedusae, comprising the benthic, so-called stalked jellyfish, appears to be the sister group of all other medusozoans, implying that the free-swimming medusa stage, the motor nerve net, and statocysts of ecto-endodermal origin are features derived within Medusozoa. Cubozoans, which have had uncertain phylogenetic affinities since the elucidation of their life cycles, form a clade—named Acraspeda—with the scyphozoan groups Coronatae, Rhizostomeae, and Semaeostomeae. The polyps of both cubozoans and hydrozoans appear to be secondarily simplified. Hydrozoa is comprised by two well-supported clades, Trachylina and Hydroidolina. The position of Limnomedusae within Trachylina indicates that the ancestral hydrozoan had a biphasic life cycle and that the medusa was formed via an entocodon. Recently hypothesized homologies between the entocodon and bilaterian mesoderm are therefore suspect. Laingiomedusae, which has often been viewed as a close ally of the trachyline group Narcomedusae, is instead shown to be unambiguously a member of Hydroidolina. The important model organisms of the Hydra species complex are part of a clade, Aplanulata, with other hydrozoans possessing direct development not involving a ciliated planula stage. Finally, applying phylogenetic mixture models to our data proved to be of little additional value over a more traditional phylogenetic approach involving explicit hypothesis testing and bootstrap analyses under multiple optimality criteria.

Keywords: 18S; 28S; Cubozoa; Hydrozoa; medusa; molecular systematics; polyp; Scyphozoa; Staurozoa

Received February 2, 2005; Revised April 25, 2005; Accepted September 12, 2005


Cnidaria has been the subject of several recent molecular studies addressing the evolutionary relationships among its major component groups (Bridge et al., 1992, 1995; Odorico and Miller, 1997; Kim et al., 1999; Collins, 2000, 2002; Medina et al., 2001). This research has demonstrated nearly beyond doubt that Cnidaria consists of two large clades, Anthozoa and Medusozoa, the latter name being a reference to the typical, though far from universal, adult pelagic medusa stage of the group. More importantly, these studies have invigorated interest in medusozoan character evolution and are at least partly responsible for spurring a number of contemporary investigations of medusozoan relationships based on morphological and life history features (Schuchert, 1993; Bouillon and Boero, 2000; Marques and Collins, 2004). These phylogenetic endeavors, both molecular and morphological, continue the efforts going back well over a century (Haeckel, 1879; Brooks, 1886) to provide a historical context for understanding how the tremendous diversity of life-history strategies, morphology at the individual and colony levels, development, and species richness of Medusozoa originated.

Somewhat ironically, the bounty of recent and older analyses leaves us faced with an unprecedented number of alternative and somewhat contradictory ideas about the phylogeny of Medusozoa. A cursory glance at Figure 1 might suggest that a great deal of uncertainty about medusozoan phylogeny exists. However, many specific hypotheses of cnidarian relationships are present in all or nearly all of the proposed relationships shown in Figure 1. Taken together, these studies provide strong evidence that their common elements are true because they are based upon different types of data assessed by diverse researchers. Moreover, these studies encompass methodologies ranging from scenario approaches, i.e., synthetic, narrative discussions of plausible histories, to more explicit cladistic and likelihood analyses. Each study also contains various caveats and focal points—e.g., Hydrozoa, Scyphozoa, etc.—and should be consulted directly for analytical details.


Figure 1
View larger version (286K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1 Comparison of hypotheses of phylogenetic relationships among medusozoan cnidarian groups published since 1940. Note that some precladistic, scenario-based hypotheses involve paraphyletic taxa, as indicated by parentheses. Taxa used in cladistic studies are marked (tt), for terminal taxon, in order to indicate that the monophyly of these groups was not explicitly tested in the analysis (though note that in some cases, as for example by Schuchert [1993], evidence supporting monophyly of the groups was presented.) For molecular studies, numbers in parentheses following taxon names indicate the number of species sampled. The cross next to the taxon Conulatae indicates that the group is extinct.

 
One of several examples of consensus (Fig. 1) involves the scyphozoan groups Semaeostomeae and Rhizostomeae, which are universally viewed as close allies forming the clade Discomedusae. In fact, precladistic scenarios, being readily compatible with the concept of paraphyletic taxa, often viewed rhizostomes as direct descendants from semaeostomes (Mayer, 1910; Hyman, 1940; Thiel, 1966), presaging results later obtained from small subunit (SSU) rDNA data (Collins, 2002). Within Hydrozoa, the holopelagic groups Trachymedusae and Narcomedusae form either a clade or a paraphyletic assemblage according to the vast majority of prior studies (Fig. 1). Similarly, nearly all analyses have suggested that the hydrozoan groups typified by highly polymorphic polyp colonies—Anthoathecata, Leptothecata, and Siphonophorae—are closely allied in a clade known as Hydroidolina (a name introduced by Collins, 2000). Finally, as noted above, all other considerations of cnidarian relationships since 1940, despite a couple of notable exceptions (i.e., Hyman, 1940; Brusca and Brusca, 1990), have concluded that Anthozoa is the sister group of Medusozoa.

If the aforementioned relationships are accepted as a reasonably certain backbone representing our present understanding of big-picture medusozoan phylogeny, then great progress can be made by clarifying the phylogenetic positions of a relatively small number of groups. For example, within Hydrozoa, Limnomedusae and two species-poor groups, Actinulida and Laingiomedusae, have uncertain positions, potentially being allied with either of the two hydrozoan clades about which there is little disagreement (Narcomedusae + Trachymedusae or Anthoathecata + Leptothecata + Siphonophorae; see Fig. 1). Similarly, among nonhydrozoan medusozoans, the class Cubozoa, and the groups Stauromedusae and Coronatae (traditionally classified as orders in the class Scyphozoa), have occupied a number of different hypothesized positions relative to each other and Discomedusae (Fig. 1). In order to reduce the uncertainty of the evolutionary relationships of these groups, additional evidence must be gathered and analyzed until a preponderance of results indicates relationships that can be accepted with some measure of confidence.

To further evaluate competing phylogenetic hypotheses involving medusozoan groups, we analyze sequences of the nearly complete (~ 3300 bp) large subunit of the nuclear ribosome (LSU or 28S) from 33 cnidarians, including 28 newly generated medusozoan sequences. In addition, we examine the most complete set (121 sequences, of which 25 medusozoan representatives are new) of cnidarian SSU (or 18S) data yet compiled. We complement more standard analyses of these data using phylogenetic mixture models in a Bayesian context (Pagel and Meade, 2004) in order to determine if explicitly accounting for heterogeneous patterns in our rDNA data helps us obtain the best possible inference of cnidarian phylogeny. These data, particularly the new LSU sequences, add further weight to some previously proposed hypotheses by strongly contradicting others and suggest a novel hypothesis for the phylogeny of Medusozoa as a whole. Our results establish a more robust working hypothesis of cnidarian phylogeny that can be used to examine character evolution within Cnidaria and should prove useful as a framework for comparative studies employing a growing number of different cnidarian species as model organisms.


    Materials
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 
Genomic DNA was extracted from specimens using the DNAzol (Chomczynski et al., 1997) and Invisorb Spin Tissue (Invitek, Berlin) kits. Nearly complete genes coding for the large subunit of the nuclear ribosome (LSU or 28S) were amplified from 28 medusozoans using the primers F63mod and R3264 from Medina et al. (2001). In order to yield sufficient template quantities for sequencing, two internal primers were matched with terminal primers (F63sq+R2077sq and F1379+R3264 from Medina et al., 2001) in a second PCR to amplify approximately the first and last two-thirds of the entire fragment using diluted template obtained from the initial PCR. Cycle sequence reactions for LSU PCR products were carried out using primers from Medina et al. (2001), as well as the primers F2800 and R2800 (from Voigt et al., 2004). Sequences of the gene coding for SSU were obtained using standard PCR and sequencing primers (Medlin et al., 1988). We used DYEnamicTM E.T.-Terminator cycle sequencing kit (Amersham Biosciences) for all cycle sequencing reactions, which were subsequently visualized with a Megabace 500 Sequencer (Amersham Biosciences). Sequences and collection data for the source specimens are deposited in GenBank (see Appendix 1 for accession numbers).

An initial LSU alignment was obtained by employing Clustal W on a set of roughly 10 LSU sequences, including five obtained from GenBank (Appendix 1). Using Seaview (Galtier et al., 1996), this coarse alignment was adjusted with the aim of placing putatively homologous sites into columns. As additional LSU sequences were incorporated into the alignment, refinements guided by visual inspection were made by adjusting both single sequences and blocks of sequences. By a similar method, new SSU sequences were aligned to a preexisting data set (Collins, 2002). Positions from both the LSU and SSU alignments that were subjectively determined to be difficult to align were identified. In addition, the program Gblocks (Castresana, 2000)—using the default settings except that positions with gaps were allowed to be chosen—was used to determine nonconserved regions of the alignment. Those regions judged as nonconserved by Gblocks were excluded from phylogenetic analysis. Aligned data sets, including all exclusion sets, are available from TreeBase (S1383).


    Phylogenetic Methods
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 
Our general approach to phylogenetic inquiry is to first examine whether our new data provide substantial evidence against already existing hypotheses. Then we move on to explicitly assess those hypotheses favored by our data. In all phylogenetic analyses, anthozoans were used as outgroup taxa. Under the criterion of maximum likelihood (ML), searches (10 replicates with taxa added randomly to the starting tree) were conducted using PAUP*4.0 (Swofford, 2000) for optimal trees conforming to constraint topologies representing 13 different hypotheses of cnidarian relationships (Table 1, Fig. 1), given the LSU and combined data. The SSU data were not examined in this way because prior work already indicated the hypotheses likely to be favored by these data (Collins, 2002). In addition to the constrained analyses, searches in the absence of topological constraints were carried out. For these analyses, an assumed model of nucleotide evolution was obtained by using the Akaike Information Criterion (AIC) as implemented in ModelTest (Posada and Crandall, 1998).


View this table:
[in this window]
[in a new window]

 
Table 1 Results of approximately unbiased (AU) and Kishino-Hasegawa (KH) tests of different hypotheses identified a priori, given the LSU and combined data.

 
In order to gauge the extent to which the LSU and combined data conflict with prior hypotheses, we performed Approximately Unbiased (AU) and Kishino-Hasegawa (KH) tests (Hasegawa and Kishino, 1989; Shimodaira, 2002) using the software Consel (Shimodaira and Hasegawa, 2001), run with 10 sets of 100,000 bootstrap replicates. The KH test is somewhat controversial, having often been used inappropriately to discern among hypotheses chosen after preferences in the data have already been identified (Goldman et al., 2000), in contrast to how it is used here. Two alternative tests, the SH (Shimodaira and Hasegawa, 1999) and SOWH (Swofford et al., 1996) tests, were not employed because the latter may have a tendency to reject hypotheses that are true, whereas the former tends not to falsify unlikely alternatives (Buckley, 2002; Shimodaira, 2002).

Because the use of different optimality criteria may lead to incorrect results under some conditions, for instance, maximum parsimony (MP) when evolutionary rates vary strongly across taxa (Felsenstein, 1978) or ML when rates of evolution vary over time (Kolaczkowski and Thornton, 2004), we analyzed our data in a variety of ways. In addition to ML, we conducted MP (200 replicate searches with taxa added randomly to the starting tree) and minimum evolution (ME) searches for optimal topologies based on the LSU and combined data sets. For ME analyses, we used the same assumed model as was used for the ML searches. Bootstrap analyses under each criterion (500, 500, and 200 replicates for MP, ME, and ML, respectively) were conducted in order to assess node support deriving from the LSU and combined data. Using the SSU data set, tree searches were conducted for the optimal MP (1000 replicates searches keeping a maximum of 500 unique trees during each replicate) and ME topologies. Node support was assessed under MP and ME by bootstrap analysis (200 random replicate data sets, with 5 random replicate heuristic searches under MP). ML analyses were not conducted with the SSU data sets due to computational limitations.

We also used phylogenetic mixture models (Pagel and Meade, 2004) in an attempt to obtain the best possible inference of cnidarian phylogeny using our rDNA data. Considerable variability in substitution rate is evident across both SSU and LSU (Van de Peer et al., 1997; Ben Ali et al., 1999), and typically a gamma parameter is used to take into account differences in the rate of evolution across sites. However, the evolution of different regions of ribosomal genes are likely to differ not just in rate but in kind, e.g., as a result of the complementary base-pairing involved in ribosomal structure. Even within categories of characters (such as stem or loop), there is likely to be significant variation in substitution patterns (Pagel and Meade, 2004). This observation led Pagel and Meade (2004) to develop the software BayesPhylogenies (Pagel and Meade, 2004), which can be used to fit two or more different models to a set of data without specifying which sites are best fit by which models. This approach, which is only implemented in BayesPhylogenies, appears to provide a powerful means for taking into account data heterogeneity.

Using BayesPhylogenies, we fit from one to six general time-reversible (GTR) models including a gamma parameter to our data sets—LSU, combined, and SSU—until we were reasonably sure of over-parameterization. We did this because there are uncertainties about how best to determine the most appropriate number of models for a given set of data (Pagel and Meade, 2004) and we wanted to observe how the number of models used impacts inferences of phylogeny. We analyzed the data by Markov chain Monte Carlo (MCMC) method, allowing the Markov chain to proceed for 500,000 iterations for the LSU and combined data sets and 1,000,000 times for the SSU data set. For each data set and model combination, we ran a minimum of four independent analyses to ensure that runs were converging to the same parameter space. Posterior probabilities (pp) were obtained by constructing a majority rule consensus of 100 trees sampled every 10,000 from 410,000 to 500,000 of the converged LSU and combined runs and 910,000 to 1,000,000 from the converged SSU runs.


    Results
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 
Of the 3619 characters in our LSU data set, 3247 were designated as conserved by Gblocks. Of these, 2130 are constant and 865 are parsimony-informative over the taxa we have examined. Our SSU data set contains 2139 positions, of which 1898 are conserved. Of these, 1142 are constant and 483 are parsimony-informative. In the combined data set, 1164 of the 5145 characters are parsimony-informative and 3620 are constant. For the LSU, SSU, and combined data sets, ModelTest indicated that the most appropriate model of nucleotide evolution is one that has six substitution rates, an assumed proportion of invariant sites, and a gamma-shape parameter (GTR+I+G), though estimates of the model parameters differ (see legends for Figs. 2, 3, 4).


Figure 2
View larger version (194K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2 Unconstrained phylogenetic hypothesis of relationships among 31 medusozoan cnidarians, rooted with two anthozoans, based on ML analysis of LSU data. Numbers in parentheses indicate the number of species sampled for that taxon. Bootstrap values under ML, ME, and MP criteria are shown at the nodes. Nodes with bootstrap indices of 100 under all three criteria are indicated with a bold 100. Bootstrap values less than 50 are represented by <. The assumed model (GTR+I+G) of nucleotide evolution for ML and ME tree searches has six substitution rates (A-C, 0.8456; A-G, 2.6613; A-T, 1.0709; C-G, 0.8204; C-T, 4.9158; and G-T, 1.0000), an assumed proportion of invariant sites (0.4776), and a gamma-shape parameter (0.5251). The length of the bar indicates 0.05 substitutions per site.

 


Figure 3
View larger version (343K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3 Phylogenetic hypothesis of relationships among 84 medusozoan cnidarians, rooted with 37 anthozoans, based on MP analysis of SSU data, which yielded this strict consensus of 68,528 MP trees of length 4062. Numbers in parentheses indicate the number of species sampled for that taxon. Bootstrap values under ME and MP criteria are shown at the nodes. Nodes with bootstrap indices of 100 under both criteria are indicated with a bold 100. Nodes with both bootstrap values above 95 are indicated with the lower of the two scores and +. Bootstrap values less than 50 are represented by <. The assumed model (GTR+I+G) of nucleotide evolution for ML and ME tree searches has six substitution rates (A-C, 1.2183; A-G, 3.1905; A-T, 1.3619; C-G, 1.1074; C-T, 5.5110; and G-T, 1.0000), an assumed proportion of invariant sites (0.4294), and a gamma-shape parameter (0.5314).

 


Figure 4
View larger version (190K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4 Phylogenetic hypothesis of relationships among 29 medusozoan cnidarians, rooted with two anthozoans, based on ML analysis of LSU and SSU data combined. Numbers in parentheses indicate the number of species sampled for that taxon. Bootstrap values under ML, ME, and MP criteria are shown at the nodes. Nodes with bootstrap indices of 100 under all three criteria are indicated with a bold 100. Bootstrap values less than 50 are represented by <. The assumed model (GTR+I+G) of nucleotide evolution for ML and ME tree searches has six substitution rates (A-C, 0.9431; A-G, 2.8946; A-T, 1.1569; C-G, 0.9034; C-T, 5.1750; and G-T, 1.0000), an assumed proportion of invariant sites (0.5333), and a gamma-shape parameter (0.5605). The length of the bar indicates 0.05 substitutions per site.

 
Results of hypothesis testing are summarized in Table 1. In all cases, significance or non-significance of the AU and KH tests matched. Using the 95 percent level as the conventional arbitrary cutoff, five hypothesized clades (3, 4, 8, 9, and 10) are rejected by both the LSU and combined data sets. By these tests, one hypothesis (5) is significantly suboptimal based on the LSU data, but is not rejected by the combined data.

Four hypothesized clades listed in Table 1 (6, 7, 11, and 13) are present in the most optimal (ML) trees (Figs. 2 and 4). Gegenbaur's (1856) "Acraspeda" most closely approximates clade 6; this name refers to the absence of a velum typical of hydrozoan medusae. We call clade 7 Scyphozoa following Marques and Collins (2004), clade 11 Trachylina, and clade 13 Hydroidolina. Finally, three hypothetical clades (1, 2, and 12) are contradicted by the most optimal trees (Figs. 2, 4), but are not found to be so suboptimal that they are rejected by the AU and KH tests employed here.

Figure 2 shows the relationships for which our compiled LSU data are most probable. Tree searches using all three optimality criteria, ME, ML, and MP, are largely congruent (ME and MP not shown), differing only in the relationships among representatives of the hydroidolinan groups Aplanulata, Other Capitata, Filifera, Laingiomedusae, Leptothecata, and Siphonophorae. With the exception of nodes involving these taxa, bootstrap support indices under the three criteria are generally high and consistent. Figure 3 shows the topology that requires the fewest assumed changes in our SSU data. In terms of the larger groups, the strict consensus of the MP trees matches the LSU hypothesis (Fig. 2) closely. However, several key clades (e.g., Stauromedusae, Cubozoa, Scyphozoa, and Hydrozoa) form a polytomy based on the SSU data. A phylogenetic hypothesis for Cnidaria based on ML analysis of combined LSU and SSU data is presented in Figure 4. The combined results mirror those of the LSU and SSU data sets examined separately (Figs. 2, 3).

In addition to those listed by Daly et al. (2002) and Collins and Daly (2005), i.e., Z92904 [GenBank] , AF099103 [GenBank] , and AF099104 [GenBank] , two cnidarian SSU sequences in GenBank are potentially misidentified. An octocoral sequence labeled Virgularia gustaviana (GenBank accession no. Z86106 [GenBank] ) consistently branches with the ceriantharian tube anemones rather than within Octocorallia, which otherwise is a highly supported clade based on SSU data (Fig. 3; and see Song and Won, 1997; Berntson et al., 1999; Won et al., 2001). Another potentially erroneous sequence is that of Cassiopea sp. (GenBank accession no. AF099675 [GenBank] ), a rhizostome jellyfish. This sequence tends to branch ambiguously near the base of Hydrozoa or Medusozoa and far away from the two other rhizostomes that had been sampled for the SSU gene. Here, however, we present a new sequence for Cassiopea xamachana, a culture of which is maintained in the Schierwater lab in Hannover, Germany. This sequence differs greatly from AF099675 [GenBank] and falls in a well supported clade with all other rhizostomes sampled (Fig. 3). Z92904 [GenBank] , AF099103 [GenBank] , AF099104 [GenBank] , Z86106 [GenBank] , and AF099675 [GenBank] are excluded from the analyses presented here.

Summarizing the results of the phylogenetic mixture models is somewhat difficult because of the large number of different analyses. We fit from one to six models with a gamma parameter to each data set resulting in 18 different estimates of cnidarian phylogeny. In each case, at least three out of the four independent runs of BayesPhylogenies converged to the same parameter space. Those runs that reached a suboptimal region of tree space, as judged by log-likelihoods, were discarded from calculations of posterior probabilities (pp). A comparison of all the average log-likelihoods of the resulting trees under the different number of models is shown Table 2. The use of an additional model always results in a better mean tree score. We always incorporated a gamma parameter because preliminary analyses indicated that use of a gamma parameter with n models (numbers of estimated parameters 6n + 1) results in an average score that exceeded that of n + 1 models without a gamma parameter (number of estimated parameters 6n + 6) for our data sets.


View this table:
[in this window]
[in a new window]

 
Table 2 Comparison of results under different numbers of phylogenetic mixture models.

 
Picking the correct number of models that best fit the data is not a simple exercise. The most obvious potential indicator would be a positive log-Bayes factor (Gelman et al., 1995; Raftery, 1996), but these are dependent on the assumption that characters are independent, an assumption that is clearly violated by rDNA data. Nevertheless, approximate log-Bayes factors, calculated as twice the increase in log-likelihood score when going from n models to n + 1 less 32 (Pagel and Meade, 2005), are shown in Table 2. This approximation can be made because BayesPhylogenies MCMC runs assign uniform priors to trees and model parameters (Pagel and Meade, 2004). Two other possible subjective indicators of overparameterization are a marked increase in the average standard deviation of the rate parameters and/or a precipitous decline in the improvements of the overall log-likelihood scores as models are added (Pagel and Meade, 2004, 2005). Both these measures are presented in Table 2.

By these three measures, there is no clear indication for what the most appropriate number of models is for any of our data sets (Table 2). The approximate Bayes factors are negative when going from Q4 to Q5, Q5 to Q6, and Q4 to Q5 for the LSU, combined, and SSU data sets, respectively. The average standard deviation of the rate parameters increases most markedly when going from Q3 to Q4, Q2 to Q3, and Q2 to Q3 for the three data sets, whereas there appear to be relatively steady declines in the improvements of the scores (Table 2).

In general, we found that the number of different models did not have a great impact on the resultant topologies for the LSU, combined, and SSU data sets. Figure 5 shows the preferred topology when three models (Q3) are applied to the LSU data. At each node, Bayesian pp's under all numbers of models fit to both the LSU or combined data illustrate that the majority of hypothesized relationships are equivalent when additional models are fit to both data sets. Most nodes have pp's equal to 100, and discrepancies arise only for those nodes that had low bootstrap support (Figs. 2, 4, and 5).


Figure 5
View larger version (50K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 5 Preferred topology based on the Bayesian analysis with three GTR models and a gamma parameter (Q3 in Table 2) fit to the LSU data. Posterior probabilities under all six models for both the LSU (above) and combined data (below) are shown at all nodes. A single value of 100 is shown at those nodes receiving high support no matter how many models are fit to either the LSU or combined data sets. Note that several nodes have been collapsed, e.g., Acraspeda contains all scyphozoan and cubozoan representatives. Every hidden node received Bayesian posterior probabilities of 100 for every combination of models and data.

 
The placement of Stauromedusae and the monophyly of Acraspeda (Cubozoa plus Scyphozoa) are of particular importance to the ensuing discussion. As additional models are fit to the LSU and combined data, support for Stauromedusae being the sister group of all other medusozoans is eroded (Fig. 5). Instead, the pp for a clade including Stauromedusae and Hydrozoa exceeds 50 when four or more models are fit to the LSU data. When the combined data are analyzed, the probability of Stauromedusae having a basal position within Medusozoa increases and the probability that Scyphozoa and Cubozoa form the clade Acraspeda decreases slightly (Fig. 5).

Similarly, fitting additional models to the SSU data has little impact on hypothesized relationships (Fig. 6). When three models and a gamma parameter are fit to the SSU data (Fig. 6), there are 68 resolved nodes within Medusozoa. Of these, 45 have pp's exceeding 90 (33 nodes with pp = 100) no matter how many models are fit to the data. Of the remaining 23 nodes, 15 are favored (pp > 50) irrespective of the number of models. Thus, just eight nodes are impacted by the number of models being fit to the data and the vast majority of these (seven) are present in most instances. In comparison to the strict consensus of most parsimonious trees (Fig. 3), 44 of the 51 nodes within Medusozoa have high support (pp > 90) no matter how many numbers of models are fit to the data. Of the remaining seven medusozoan nodes in the MP tree, five have pp's that exceed 70 under all models. Of the two remaining nodes, each one has a pp that is below 50 for just a single instance. In sum, applying phylogenetic mixture models to the SSU data results in topologies that are highly consistent with each other (Fig. 6). These topologies are also consistent with the strict consensus of MP trees (Fig. 3), but have greater resolution. The question then arises whether the increased resolution resulting from the Bayesian analyses represents true relationships. Because we cannot know the true phylogeny, these hypothetical nodes can only be judged in light of other information and future analyses.


Figure 6
View larger version (323K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 6 Preferred topology based on the Bayesian analysis with three GTR models and a gamma parameter (Q3 in Table 2) fit to the SSU data. Posterior probabilities when one to six different models are fit to the data are shown at each node, with the exception of internal nodes within the major anthozoan groups. A single value of 100 or a range of high values are shown at those nodes receiving high support (> 90) no matter how many models are fit to the data. Posterior probabilities less than 50 are represented by <. Seventeen nodes not present in the MP analysis (Fig. 3) are indicated with*.

 

    Discussion
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 
Quality of the LSU Marker
The LSU marker appears to be of great use for phylogenetic analysis of this diverse group of metazoans. Hypotheses based on LSU data (Fig. 2) share many similarities with those based on the analysis of nonmolecular characters (Fig. 1), provide support for nodes that were equivocal based on SSU data alone, and are highly congruent with those derived from SSU data (Fig. 3; Collins, 2002). More importantly, there are no strong contradictions between conclusions drawn from the LSU and SSU data. One might suspect that both markers share similar biases, given that LSU and SSU are part of the complex of ribosomal genes that evolve by concerted evolution. In fact, a simple visual comparison of the phylogram in figure 3 of Collins (2002), to Figure 2 here reveals that SSU and LSU exhibit similar lineage-specific variation in evolutionary rates, e.g., cubozoans and especially leptothecates are relatively long-branched taxa.

On the other hand, SSU data are also congruent with the generally faster evolving data of partial mitochondrial 16S sequences in their indications of relationships among less inclusive groups of cnidarians, e.g., within Hydrozoa (Collins et al., 2005), Hexacorallia (Daly et al., 2003), and Stauromedusae (Collins and Daly, 2005). Although 16S is another ribosomal gene, there is no mechanism known that explains why this mitochondrial gene would be biased in the same way as nuclear SSU and LSU. The important conclusion that can be drawn is that combining data from all readily available rDNA markers—including complete mitochondrial 16S and 12S and complete nuclear SSU, LSU, and 5.8S—might be an effective and efficient strategy for pursuing large-scale studies aimed at elucidating cnidarian relationships at all taxonomic scales. Of course, this does not eliminate the need to develop and use alternative markers to test ribosomal-based hypotheses and seek resolution where it is presently lacking.

Utility of Phylogenetic Mixture Models
When employing phylogenetic mixture models, following a simple rule to determine the number of models to fit to the data is unwise, given the subjective nature of the decision. Instead, changes (if any) in topology that occur as models are added should be examined. Figures 5 and 6 show that this choice has relatively little influence on our phylogenetic inferences based on our data sets. Furthermore, applying phylogenetic mixture models to our data resulted in topologies that were overwhelmingly consistent with those derived from ML (Figs. 2, 4) and MP (Fig. 3) analyses. Not surprisingly, taxa that group with low bootstrap support are those most likely to have alternative relationships in the Bayesian results and whose inferred relationships are most influenced by the application of different numbers of models. In general, we conclude that little insight was gained by applying phylogenetic mixture models to these data sets because there seems little justification to favor any one particular set of relationships over another.

Nevertheless, our SSU data set was too large to analyze using maximum likelihood and the Bayesian approach allowed us to examine our data using a method that takes advantage of the likelihood function. The Bayesian topology for the SSU data (Fig. 6) contains all of the medusozoan nodes revealed in the strict consensus of MP trees (Fig. 3), as well as 17 additional nodes (denoted with asterisks in Fig. 6). Whether these nodes accurately reflect phylogenetic history or not can only be judged in light of other evidence. For most of these hypothesized relationships, there is no relevant independent evidence, and so their veracity can only be judged through future phylogenetic analyses.

Evidence can, however, be brought to bear on a few of these hypothetical clades. First, Figure 6 suggests the possibility of a relationship between Stauromedusae and Cubozoa (Fig. 6). This hypothesis accords with the results of a cladistic analysis of morphological and life history characters (Marques and Collins, 2004), but a reanalysis of these same data, except that scoring of characters dealing with the fossil group Conulatae was altered, found that Stauromedusae was the sister group to all other medusozoans (Van Iten et al., in press). Moreover, although LSU data do not have sufficient strength to statistically falsify a close relationship between Stauromedusae and Cubozoa (Table 1), this hypothesis has a pp of 0 no matter how many models are fit to the LSU data. Within the rhizostomes sampled here, a clade containing Cassiopea and Cotylorhiza agrees with the views of Stiasny (1921) and Thiel (1970). In addition, many of the relationships revealed among the capitate taxa sampled here were also found through an analysis of mitochondrial 16S data (Collins et al., 2005). In contrast, relationships among the members of Aplanulata contradict those derived from 16S data (Collins et al., 2005). On the whole, it seems premature to conclude that the application of phylogenetic mixture models has allowed us to efficiently extract more phylogenetic information from this relatively taxon-rich data set.

Relationships among the Cnidarian Classes
Our working hypothesis of relationships among the cnidarian classes (Fig. 7) has Anthozoa as the earliest diverging class, Staurozoa as the sister group of all other medusozoans, and Cubozoa and Scyphozoa forming a clade that is the sister group of Hydrozoa. In order to communicate that Stauromedusae has a discrete history and that it possesses features that are correspondingly distinct from those of the other cnidarian classes, Marques and Collins (2004) recognized the new class Staurozoa. Our molecular data contradict the placement of Stauromedusae derived from their cladistic analysis of morphological and life-history characters. However, it reaffirms that Stauromedusae is a primary clade of Cnidaria and that it does not form a clade with the scyphozoan groups with which it has traditionally been classified. Instead, a growing body of both molecular (Collins and Daly, 2005) and morphological (Van Iten et al., in press) evidence suggests that Stauromedusae is the sister group to all other medusozoans (Fig. 7). For classification purposes, using Staurozoa as a separate class in Cnidaria conveys that the phylum is more diverse at a fundamental level than previously recognized.


Figure 7
View larger version (211K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 7 Working hypothesis for cnidarian relationships with selected hypothesized ancestral characters mapped at nodes. Diagrams representing taxa are after figures in Mayer (1910), with the exception of representatives of Octocorallia and Siphonophorae, which are modified from Hyman (1940), and Hexacorallia, which was drawn by Crissy Huffard.

 
Cubozoa was designated as a class (Werner, 1973) because its members have a life cycle and polyp morphology that are rather distinct from those of scyphozoans (Werner et al., 1971). Rather than producing ephyrae through strobilation, cubozoan medusae result from a total metamorphosis of the polyp (Werner et al., 1971; Studebaker, 1972; Arneson and Cutress, 1976; Yamaguchi and Hartwick, 1980). This difference is qualified by exceptions: in some species, part of the polyp sometimes remains after metamorphosis into a juvenile medusa and is capable of developing into a new polyp and producing a subsequent medusa (Stangl et al., 2002; Straehler-Prohl and Jarms, 2005). Some of the starkest differences between cubozoans and scyphozoans involve polyp morphology. Whereas scyphozoan polyps exhibit strong fourfold symmetry due to the presence of four gastric septa and four intramesogleal muscle strands, cubozoan polyp organization is essentially radial. These facts were stressed by Werner (1973), who championed the idea that Cubozoa is more closely related to Hydrozoa than it is to Scyphozoa.

Since Werner's 1973 work, the position of Cubozoa within Medusozoa has been rather contentious, with some favoring a hydrozoan affinity (Arneson and Cutress, 1976; Petersen, 1979; Leonard, 1980) and others a closer relationship to scyphozoans (Larson, 1976; Salvini-Plawen, 1978; Satterlie and Spencer, 1980; Schuchert, 1993; Ax, 1996; Marques and Collins, 2004). The two prior molecular analyses that contained samples of the relevant groups to test these competing alternatives (Bridge et al., 1995; Collins, 2002) provided no compelling evidence for a definitive placement of Cubozoa. Because LSU data provide strong evidence that Cubozoa forms a clade with Scyphozoa (Acraspeda), it seems likely that the absence of four gastric septae and interradial muscle fibers are due to evolutionary loss in the ancestry of Cubozoa (Fig. 7). This accords well with observations that cubopolyps contain intramesogleal muscles, though not concentrated in four fibers (Chapman, 1978) and that numerous similarities exist between cubozoan and scyphozoan medusae, including neuromuscular and nerve net organization (Satterlie, 1979, 2002) and the presence of rhopalia (Fig. 7).

Our working hypothesis suggests that hydrozoan polyps, which are radial and lack gastric pockets, were also derived from an ancestor that was more complex in these respects (Fig. 7). Furthermore, the ancestral hydrozoan likely had a biphasic life cycle with the medusa formed via lateral budding and development through a tissue mass known as the entocodon (Fig. 7). This complex trait involving an ontogenetic process as well as morphology is typical of Limnomedusae and the hydroidolinan groups (Anthoathecata, Leptothecata, and Siphonophorae), but it is absent in Actinulida, Narcomedusae, and Trachymedusae, and not yet determined for Laingiomedusae. Because Limnomedusae is at the base of Trachylina (Fig. 7; Collins, 2002) rather than among the hydroidolinan groups (Petersen, 1979; Bouillon and Boero, 2000), the last common ancestor of extant hydrozoans likely produced medusae by a developmental process involving an entocodon and lateral budding from polyps. The polyp stage and the attendant mode of medusa development were apparently lost in the lineage leading to Trachymedusae and Narcomedusae.

Origin of the Medusa Stage: A Link to Bilaterian Mesoderm?
Whether the adult form of the ancestral cnidarian was a polyp or a medusa has been a classic dispute. Brooks (1886) contended that a pelagic ancestor directly gave rise to the holopelagic trachyline hydrozoans, implying that even anthozoans were derived from hydrozoan ancestors. This opinion found its way into important textbooks (e.g., Hyman, 1940; Brusca and Brusca, 1990), but was probably a minority view. Instead, many argued that sessile polyps represent the ancestral cnidarian form and that the fundamental phylogenetic divergence within Cnidaria separated Anthozoa from the medusa-bearing groups (Haeckel, 1879; Hadzi, 1953; Werner, 1973; Salvini-Plawen, 1978). This latter phylogenetic hypothesis is widely accepted today because of new data, most notably that the mitochondrial genomes of medusozoan cnidarians are linear (Bridge et al., 1992). Despite the consensus phylogenetic view, however, the nature of the ancestral cnidarian remains somewhat controversial, with some authors supporting the idea that a pelagic medusa stage has probably been lost in the lineage leading to Anthozoa (Schuchert, 1993; Spring et al., 2002; Muller et al., 2003; Scholtz, 2004).

Most recent arguments championing the idea that the ancestral cnidarian possessed a medusa stage are built on the observation that several genes involved in specifying mesoderm and differentiating striated muscle in bilaterian animals are shared by the hydrozoan Podocoryna carnea during the ontogeny of its medusa as it develops from the entocodon (Spring et al., 2000, 2002; Muller et al., 2003). Because the medusae of cnidarians have striated muscles lining their subumbrellar surfaces and hydrozoan polyps do not, these authors conclude that a swimming medusa stage was present in the ancestral cnidarian. These conclusions are based on observations of a hydrozoan and on the hypothesized homologies between this hydrozoan and bilaterians, implicitly suggesting that the ancestral process of medusa production for Cnidaria involved lateral budding of a medusa that grew from the entocodon.

Our working hypothesis, however, suggests that medusa production via lateral budding and growth of the entocodon is a synapomorphy of Hydrozoa, nested deeply within Medusozoa (Fig. 7). Therefore, hypothesized homologies between the entocodon and bilaterian mesoderm (Spring et al., 2002; Mueller et al., 2003; Seipel and Schmid, 2005) are suspect. Nevertheless, the intriguing similarities between the genes involved in myogenesis in the hydrozoan Podocoryna carnea and bilaterians still need to be explained. One potential key to understanding these similarities is that, in contrast with the lack of striated muscles in Podocoryna polyps, diverse cnidarian polyps do possess striated myofibers, including those of scyphozoans (Chia et al., 1984; Matsuno and Kawaguti, 1991), cubozoans (Chapman, 1978; Golz, 1993), and anthozoans (Amerongen and Peteya, 1980). Thus, the proposed homologies between the hydrozoan Podocoryna and bilaterians in myogenesis seem quite plausible. Study of additional model cnidarians may reveal that homologous genetic pathways involved in differentiation and specification of the striated muscles are broadly distributed within Cnidaria.

By providing reasonably strong evidence that Stauromedusae—comprising the sessile, so-called stalked jellyfishes—is the earliest diverging lineage within Medusozoa, our data further weaken the argument that there is a pelagic medusa stage in the ancestry of Anthozoa. Instead, the medusa appears to be derived within Medusozoa (Fig. 7). Furthermore, the motor nerve net (also known as giant fiber nerve net) and statocysts of ecto-endodermal origin likely developed in association with the pelagic, motile medusa stage. Despite these shared features, differences in the medusae of Hydrozoa and Acraspeda (e.g., lateral budding versus metamorphosis at the oral end of a polyp, presence versus absence of velum, absence versus presence of gastric filaments and coronal muscle, etc.; see Fig. 7) have led some to argue that the two are not homologous (Thiel, 1966; Salvini-Plawen, 1978).

A plausible, though speculative, scenario would involve a "reinvented" pelagic medusa in the ancestry of Hydrozoa. It is conceivable that the lateral budding of medusae in hydrozoans arose through the transformation of a laterally budding polyp into an entire medusa. The polyps of numerous species of cubozoans and scyphozoans laterally bud secondary polyps or vermiform stages known as frustules, which creep some distance before differentiating into polyps (Hofmann and Crow, 2002; Fischer and Hofmann, 2004). An interesting observation from the hydrozoan species Podocoryna carnea demonstrates that there may be some plasticity with respect to the morphologies that can be produced at budding sites. Polyps of Podocoryna carnea, if cultivated under stressful conditions (e.g., high temperature, suboptimal sea water, etc.), will sometimes start to develop little hydranths at the positions where medusae would otherwise arise (Peter Schuchert, personal observation).

The preceding discussion highlights the difficulties in assuming putative homology between entire life phases of the different cnidarian groups. For instance, the so-called "stalked jellyfishes" have frequently been viewed as possessing a mix of polyp and medusa characters that indicate a state of "degeneracy." In other words, stauromedusans have been called stalked jellyfishes because they were thought to be descended from an ancestor with a pelagic medusa phase (Thiel, 1966; Uchida, 1972; Werner, 1973). Given the likely phylogenetic position of Stauromedusae (Fig. 7), it might actually be more accurate to think of medusae as free-swimming stauromedusans. However, attempting to homologize entire life stages across Cnidaria is almost certainly an oversimplification because some medusozoan characters may have arisen prior to the origin of a medusa and a benthic-pelagic life cycle. It follows that separating features into groups specific to polyp or medusa (as seen for instance in the cladistic scoring of Marques and Collins, 2004) is not entirely appropriate, though difficult to get around. Focusing on the potential homologies between characters rather than entire life stages may help us move beyond the old debate about which came first, polyp or medusa. As it stands, the possibility that anthozoan groups have or had characters that are presently thought of as features of the medusa is largely unexplored.

Relationships and Evolution within Scyphozoa
Ephyrae and polydisk strobilation are likely synapomorphies for Scyphozoa (Fig. 7). Within this clade, rhizostome jellyfishes have long been considered to be direct descendants from semaeostome ancestors because of similarities in their radial canal systems (Mayer, 1910; Hyman, 1940; Thiel, 1966). This hypothesis is confirmed by both SSU (Collins, 2002; Figs. 3, 6) and LSU (Figs. 2, 4, 5) data. Of the three semaeostome groups (Cyaneidae, Pelagiidae, and Ulmaridae), Pelagiidae (Chrysaora spp.) appears to be the earliest diverging. Ulmaridae (Aurelia and Phacellophora) most likely gave rise to Rhizostomeae (Figs. 3, 6). Within rhizostomes there are two main alliances of taxa, Cepheida and Rhizostomida (Stiasny, 1921; Thiel, 1970), which Thiel (1970) envisioned as independently derived from within Semaeostomeae. However, members of Cepheida (Cassiopea and Cotylorhiza) and Rhizostomida (Catostylus and Stomolophus) form a well-supported monophyletic group (Figs. 3, 6). Thus, the monodisk strobilation and fused mouth arms of rhizostomes likely have a single origin (Fig. 7).

Relationships and Evolution within Hydrozoa
The new LSU and SSU data presented here greatly clarify the scope and likely phylogenetic position of Limnomedusae. Of the various groups that have been classified within Limnomedusae (Moerisiidae, Monobrachiidae, Olindiasidae, and Proboscidactylidae) only Monobrachiidae and Olindiasidae (= Olindiidae) are supported by molecular data as a paraphyletic assemblage within the clade Trachylina (Figs. 3, 6). Our data reinforce hypotheses based on morphology that suggest that Moerisiidae and Proboscidactylidae are members of Hydroidolina (Rees, 1958; Edwards, 1973; Petersen, 1990; Schuchert, 1996). Limnomedusae appears to be most easily characterized by two symplesiomorphies, i.e., ecto-endodermal statocysts and a life cycle that includes a polyp stage, shared with other members of Trachylina (Trachymedusae and Narcomedusae) and Hydroidolina, respectively.

The freshwater jellyfishes, here represented by the genera Astrohydra (LSU only), Craspedacusta (SSU only), and Limnocnida, appear to have a single origin. Polyps of these genera are very similar (Bouillon, 1957; Hashimoto, 1981), but Craspedacusta and Limnocnida have adult medusae that differ markedly. Little is known about the medusa stage of Astrohydra, which has only once been observed and then only in a juvenile state (Hashimoto, 1985). As with other members of Limnomedusae, Craspedacusta has medusae with four-parted mouths and gonads distributed along the radial canals. In contrast, the jellyfish of Limnocnida has a large circular mouth and gonads distributed around the edge of its gut. Indeed, prior to knowledge of its polyp stage, Limnocnida proved difficult to classify (Mayer, 1910). Nevertheless, both species possess the ecto-endodermal statocysts that are typical of other trachylines. The unusual medusa form of Limnocnida appears to be relatively recently derived from a more typical limnomedusan form. The only brackish water species within Limnomedusae (Maeotias marginata) is always well supported as the sister group to the freshwater species (Figs. 3, 6).

LSU and SSU data analyzed separately and together indicate that the single representative of Laingiomedusae sampled here is an unambiguous member of the clade Hydroidolina (Figs. 2, 3, 4, 5, 6). Laingiomedusae was erected by Bouillon (1978) as a new hydrozoan subclass (equivalent in his ranking system to Limnomedusae, Trachymedusae, etc.) and contains just four known species in three genera. The complete life cycle is unknown for any representative of Laingiomedusae, as only the medusa stages have been observed. They have tentacles that emerge on the exumbrellar side of a scalloped margin, characters reminiscent of Narcomedusae, and indeed some analyses of morphology have suggested that the two groups share a close alliance (Brusca and Brusca, 1990; Marques and Collins, 2004). In contrast, Schuchert (1996) argued that Laingiomedusae, in particular the genus sampled here, Fabienna (placed in Laingiomedusae by Bouillon and Barnett, 1999), may actually be closely related to anthoathecate hydrozoans of the family Proboscidactylidae; both possess a solid radial canal and macrobasic euryteles. This hypothesis is strongly supported by SSU data (Figs. 3, 6). Furthermore, there is a strong similarity between Laingiomedusae and another filiferan, Thecocodium quadratum, which has a markedly lobed umbrella margin and tentacles that issue on the exumbrella (see Jarms, 1987). We predict that, upon discovery, the polyps of Fabienna will be filiferan hydroids. If future sampling shows that the name-giving species of Laingiomedusae, Laingia jaumotti, forms a clade with Fabienna and the filiferan group Proboscidactylidae, then the group should be demoted to family level and relegated to Anthoathecata.

Within Hydroidolina, relationships are uncertain (Figs. 2, 3, 4, 5, 6). We have no new data for Siphonophorae and there is little to add to the discussion of SSU data by Collins (2002). Similarly, the position of Leptothecata within Hydroidolina remains uncertain. We have added an SSU sequence for Melicertum octocostatum, a species that has a typical leptothecate medusa stage but has a hydroid stage lacking a theca. This species, though of uncertain position, appears to have diverged early in the history of Leptothecata (Figs. 3, 6), suggesting the possibility that the typical leptothecate medusa feature, gonads distributed along the radial canals, arose prior to the origin of thecae on polyps. Better sampling of Leptothecata and Hydroidolina is needed in order to pinpoint the origins of these characters.

Unlike Siphonophorae and Leptothecata, the phyletic status of Anthoathecata (= Capitata plus Filifera), a group for which we know of no synapomorphy (Schuchert, 1996), is unresolved by these data. Indeed, there is no molecular evidence for monophyly of either of the two anthoathecate subtaxa, Capitata or Filifera (Collins, 2000, 2002; Collins et al., 2005; and Figs. 2, 3, 4, 5, 6). It would not be surprising if both Leptothecata and Siphonophorae are found to be derived from within Anthoathecata. For Siphonophorae, this scenario has certainly been entertained (Haeckel, 1888; Garstang, 1946; Leloup, 1955; Totton, 1965). It is equally true that our data provide no strong evidence against monophyly of Anthoathecata, Capitata, or Filifera (with the exception of Laingiomedusae). These phylogenetic questions must be regarded as those that are most in doubt within Medusozoa, and their importance is underscored by the fact that Hydroidolina represents the majority of diversity within Medusozoa. It may be that an especially rapid evolutionary radiation separated the main groups within Hydroidolina. Regardless of the root cause for uncertainty about hydroidolinan phylogeny, the best strategy will be to continue sampling new taxa and markers and to compile a general phylogenetic hypothesis by identifying stable alliances among its component groups.

As an example, Collins et al. (2005) proposed, based on mitochondrial 16S data, that the important model organisms in the Hydra species group, Hydridae, form a clade with the anthoathecate families Candelabridae, Corymorphidae, and Tubulariidae. Species of Hydridae completely lack medusae, and different authors have come to conflicting conclusions concerning the group's phylogenetic position within Hydrozoa (e.g., Naumov, 1960; Petersen, 1990; Stepan'yjants et al., 2000). Species of Hydridae have direct development that does not involve the ciliated planula typical of Cnidaria; this character is shared by the anthoathecate groups Acaulidae, Candelabridae, Corymorphidae, Margelopsidae, Paracorynidae, Tricyclusidae, and Tubulariidae. In reference to this putative synapomorphy, the name Aplanulata has been given to the putative clade uniting Hydridae with Candelabridae, Corymorphidae, and Tubulariidae (Collins et al., 2005). The SSU data presented here (Figs. 3, 6), and to a lesser extent the LSU data (Fig. 2), affirm the validity of Aplanulata. Acaulidae, Margelopsidae, Paracorynidae, and Tricyclusidae have yet to be sampled for molecular data, but it seems likely that these groups are also part of Aplanulata.

Exploring Characters of Evolutionary Interest
Because Cnidaria is an early diverging lineage of Metazoa, it is likely that many characters of evolutionary interest have yet to be adequately examined. Such characters will need to be studied intensively in a small number of model organisms that are chosen with phylogenetic relationships in mind. Consider the case of gap junctions and the underlying proteins that form them. Prior work using antibodies to rat connexins suggested that connexins, the proteins underlying vertebrate gap junctions, are present in the hydrozoan Hydra (Fraser et al., 1987). However, sequencing of expressed sequence tags of Hydra have failed to reveal connexins, but instead have turned up innexins (Alexopoulos et al., 2004), proteins that form gap junctions in non-vertebrate bilaterians. This result prompted Alexopolous et al. (2004) to conclude that innexin-based gap junctions arose in the common ancestor of cnidarians and bilaterians. However, Hydrozoa is nested relatively deeply in Cnidaria, gap junctions are not apparent in the other cnidarian groups (Mackie et al., 1984), and connexin-like proteins have been documented in an anthozoan (Germain and Anctil, 1996). Thus, further investigations of the genomes of distantly related cnidarians are necessary in order to evaluate the generality of the character states exhibited by Hydra.

Fortunately, recent efforts have established the anthozoan Nematostella as a new model organism (e.g., Darling et al., 2005). The synergistic impact of having two distantly related cnidarians developed as model organisms is likely to lead to many important insights into the evolution of form, development, and genomes. Nevertheless, when differences arise between results obtained from Hydra and those from Nematostella, there will be little justification, in the absence of relevant information obtained from outgroup taxa, for assuming one set of results is more likely to represent ancestral states in Cnidaria. Instead, other model organisms that are most likely to reveal ancestral character states for Anthozoa and Medusozoa will have to be developed in order to properly interpret the differences.


View this table:
[in this window]
[in a new window]

 
Appendix 1 Sequences analyzed.

 

    Acknowledgment
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 
We thank the Smithsonian Tropical Research Institute and R. Collin, who provided support allowing for the collection of specimens; C. Cunningham for allowing the use of three previously unpublished SSU sequences; J. Collins for assistance collecting samples; C. Widmer for samples of Aglauropsis aeora; G. Mackie for samples of Aglantha digitale, Proboscidactyla flavicirrata, and Melicertum octocostatum; I. Kossevitch for samples of Monobrachium parasitum; A. De Wever for a sample of Limnocnida tanganyicae; G. Jarms and I. Straehler-Pohl for providing a DNA sample of H. octoradiatus; M. Dawson for samples of Depastromorpha africana; J. Collins, N. Stoletzki, B. Bentlage, and other students at ITZ, TiHo-Hannover, for assistance deriving sequences; M. Pagel for answering simplistic questions about Bayes factors; the Smithsonian's Laboratory of Analytical Biology for the use of its computer cluster to conduct revised analyses; E. Begovic, A. Lindner, S. Nichols, V. Pearse, and J. Valentine for commenting on earlier versions of this manuscript; and D. Bridge, M. Daly, T. Collins, and R. Page for many helpful suggestions during the review process. Grants to BS from the Deutsche Forschungsgemeinschaft (Schi-277/10-2, Schi-277/20-1) are gratefully acknowledged, as are grants to ACM from Fundação de Amparo à Pesquisa do Estado de São Paulo and Conselho Nacional de Desenvolvimento Científico e Tecnológico.


    References
 Top
 Abstract
 Materials
 Phylogenetic Methods
 Results
 Discussion
 Acknowledgment
 References
 

    Alexopoulos H., Böttger A., Fischer S., Levin A., Wolf A., Fujisawa T., Hayakawa S., Gojobori T., Davies J. A., David C. N., Bacon J. P. Evolution of gap junctions: the missing link? Curr. Biol. (2004) 14:R879–R880.[CrossRef][Web of Science][Medline]

    Amerongen H. M., Peteya D. J. Ultrastructural study of 2 kinds of muscle in sea-anemones Stomphia and Coccinea, the existence of fast and slow muscles. J. Morph (1980) 166:145–154.[CrossRef]

    Arneson A. C., Cutress C. E. Life history of Carybdea alata Reynaud, 1831 (Cubomedusae). In: Coelenterate ecology and behavior—Mackie G. O., ed. (1976) New York: Plenum Press. 227–236.

    Ax P. Multicellular animals: A new approach to the phylogenetic order in nature (1996) Berlin: Springer Verlag.

    Berntson E. A., France S. C., Mullineaux L. S. Phylogenetic relationships within the Class Anthozoa (Phylum Cnidaria) based on nuclear 18S rDNA sequences. Mol. Phylog. Evol. (1999) 13:417–433.[CrossRef][Web of Science][Medline]

    Bouillon J. Limnocnida congoensis nouvelle espède Limnoméduse du bassin du Congo. Rev. Zool. Bot. Afr (1957) 56:388–395.

    Bouillon J. Hydromèduses de la mer de Bismarck (Papousie, Nouvelle-Guinee). II.—Limnomedusa, Narcomedusa, Trachymedusa et Laingiomedusa. Cahiers Biol. Mar. (1978) 19:473–487.

    Bouillon J., Barnett T. J. The Marine Fauna of New Zealand: Hydromedusae (Cnidaria: Hydrozoa). NIWA Biodiversity Memoir (1999) 113:1–136.

    Bouillon J., Boero F. The Hydrozoa: a new classification in the light of old knowledge. Thalass. Salent. (2000) 24:3–45.

    Bridge D., Cunningham C. W., Desalle R., Buss L. W. Class-level relationships in the phylum Cnidaria—molecular and morphological evidence. Mol. Biol. Evol. (1995) 12:679–689.[Abstract]

    Bridge D., Cunningham C. W., Schierwater B., Desalle R., Buss L. W. Class-level relationships in the phylum Cnidaria—evidence from mitochondrial genome structure. Proc. Natl. Acad. Sci. USA (1992) 89:8750–8753.[Abstract/Free Full Text]

    Brooks W. K. Life history of the Hydromedusae. A discussion of the medusae and of the significance of metagenesis. Mem. Boston Soc. Nat. Hist. (1886) 3:359–430.

    Brusca R. C., Brusca G. J. Invertebrates (1990) Sunderland, Massachusetts: Sinauer Associates.

    Buckley T.-R. Model misspecification and probabilistic tests of topology: Evidence from empirical data sets. Syst. Biol. (2002) 51:509–523.[Abstract/Free Full Text]

    Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. (2000) 17:540–552.[Abstract/Free Full Text]

    Chapman D. M. Microanatomy of the cubopolyp, Tripedalia cystophora (Class Cubozoa). Helgol. wissen. Meeresunters (1978) 31:128–168.[CrossRef]

    Chia F.-S., Amerongen H. M., Petaya D. G. Ultrastructure of the neuromuscular system of the polyp of Aurelia aurita L. 1758 (Cnidaria, Scyphozoa). J. Morph. (1984) 180:69–79.[CrossRef]

    Chomczynski P., Mackey K., Drews R., Wilfinger W. DNAzol: A reagent for the rapid isolation of genomic DNA. Biotechniques (1997) 22:550–553.[Web of Science][Medline]

    Collins A. G. Towards understanding the phylogenetic history of Hydrozoa: Hypothesis testing with 18S gene sequence data. Scient. Mar. (2000) 64:5–22.

    Collins A. G. Phylogeny of Medusozoa and the evolution of cnidarian life cycles. J. Evol. Biol. (2002) 15:418–432.[CrossRef][Web of Science]

    Collins A. G., Daly M. A new deepwater species of Stauromedusae, Lucernaria janetae (Cnidaria, Staurozoa, Lucernariidae), and a preliminary investigation of stauromedusan phylogeny based on nuclear and mitochondrial rDNA data. Biol. Bull (2005) 208:221–230.[Abstract/Free Full Text]

    Collins A. G., Winkelman S., Hadrys H., Schierwater B. Phylogeny of Capitata and Corynidae (Cnidaria, Hydrozoa) in light of mitochondrial 16S rDNA data. Zool. Scripta (2005) 34:91–99.[CrossRef][Web of Science]

    Daly M., Fautin D. G., Cappola V. A. Systematics of the Hexacorallia (Cnidaria: Anthozoa). Zool. J. Linn. Soc. (2003) 139:419–437.[CrossRef]

    Darling J. A., Reitzel A. R., Burton P. M., Mazza M. E., Ryan J. F., Sullivan J. C., Finnerty J. R. Rising starlet: The starlet sea anemone, Nematostella vectensis. Bioessays (2005) 27:211–221.[CrossRef][Web of Science][Medline]

    Edwards C. Contributory thoughts on form, function, habitat and classification of hydroids and hydromedusae. Publ. Seto Mar. Biol. Lab. (1973) 20:11–22.

    Felsenstein J. Cases in which parsimony or compatibility methods will be positively misleading. Syst. Zool. (1978) 27:401–410.[Abstract/Free Full Text]

    Fischer A. B., Hofmann D. K. Budding, bud morphogenesis, and regeneration in Carybdea marsupialis Linnaeus, 1758 (Cnidaria: Cubozoa). Hydrobiologia (2004) 530/531:331–337.

    Fraser S. E., Green C. R., Bode H. R., Gilula N. B. Selective disruption of gap junctional communication interferes with a patterning process in hydra. Science (1987) 237:49–55.[Abstract/Free Full Text]

    Galtier N., Gouy M., Gautier C. SeaView and Phylo_win, two graphic tools for sequence alignment and molecular phylogeny. Comput. Applic. Biosci. (1996) 12:543–548.[Abstract/Free Full Text]

    Garstang W. The morphology and relations of the Siphonophora. Q. J. Microsc. Sci. (1946) 87:103–193.

    Gegenbaur C. Versuch eines Systemes der Medusen mit Beschreibung neuer oder wenig gekannter Formen; zugleich ein Beitrag zur Kenntniss der Fauna des Mittelmeeres. Z. wiss. Zool. (1856) 8:202–273.

    Gelman A., Carlin J. B., Stern H. S., Rubin D. B. Bayesian data analysis (1995) London: Chapman and Hall. 420–438.

    Germain G., Anctil M. Evidence for intercellular coupling and connexin-like protein in the luminescent endoderm of Renilla koellikeri(Cnidaria, Anthozoa). Biol. Bull. (1996) 191:353–366.[Abstract]

    Goldman N., Anderson J. P., Rodrigo A. G. Likelihood-based tests of topologies in phylogenetics. Syst. Biol. (2000) 49:652–670.[Abstract/Free Full Text]

    Golz R. Anchorage and retraction of nematocytes in the tentacles of the cubopolyp Carybdea marsupialis are mediated by a species-specific mesogleal support. Cell Tiss. Res. (1993) 274:173–180.[CrossRef][Web of Science]

    Hadzi J. An attempt to reconstruct the system of animal classification. Syst. Zool. (1953) 2:145–154.[Free Full Text]

    Haeckel E. H. Das System der Medusen: Erster Theil einer Monographie der Medusen (1879) Jena.: G. Fischer.

    Haeckel E. H. Report on the Siphonophorae collected by H.M.S. Challenger during the years 1873–76 (1888) Edinburgh: Her Majesty's Stationery Office.

    Hasegawa M., Kishino H. Heterogeneity of tempo and mode of mitochondrial DNA evolution among mammalian orders. Japn. J. Genet. (1989) 64:243–258.[CrossRef]

    Hashimoto H. The fresh water hydroid Astrohydra japonica new genus new species. Annotat. Zool. Japn (1981) 54:207–212.

    Hashimoto H. Medusa of fresh-water hydroid Astrohydra japonica. Zool. Sci (1985) 2:761–766.[Web of Science]

    Hofmann D. K., Crow G. Induction of larval metamorphosis in the tropical scyphozoan Mastigias papua: striking similarity with upside-down jellyfish Cassiopea spp. (with notes on related species). Vie Mil (2002) 52:141–147.

    Hyman L. H. The invertebrates, Protozoa through Ctenophora (1940) Volume 1. New York: McGraw-Hill.

    Jarms G. Thecocodium quadratum (Werner 1965) redescribed, T. penicillatum sp. nov. and a method for rearing hydrozoans. In: Modern trends in the systematics, ecology, and evolution of Hydroids and Hydromedusae—Bouillon J., Boero F., Cicogna F., Cornelius P. F. S., eds. (1987) Oxford, UK: Clarendon Press. 57–66.

    Kim J. H., Kim W., Cunningham C. W. A new perspective on lower metazoan relationships from 18S rDNA sequences. Mol. Biol. Evol. (1999) 16:423–427.[Web of Science][Medline]

    Kolaczkowski B., Thornton J. W. Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous. Nature (2004) 431:980–984.[CrossRef][Medline]

    Larson R. J. Cubomedusae: Feeding, functional morphology, behavior and phylogenetic position. In: Coelenterate Ecology and Behavior—Mackie G. O., ed. (1976) New York: Plenum Publishing. 237–245.

    Leloup E. Siphonophores. Rep. Scient. Results Michael Sars N. Atlant. Deep Sea Exped. (1955) 5:1–24.

    Leonard J. L. Cubomedusae belong to the class Cubozoa, not Scyphozoa. Nature (1980) 284:377.

    Mackie G. O., Anderson P. A. V., Singia C. L. Apparent absence of gap junctions in 2 classes of Cnidaria. Biol. Bull. (1984) 167:120–123.[Abstract/Free Full Text]

    Marques A. C., Collins A. G. Cladistic analysis of Medusozoa and cnidarian evolution. Invert. Biol. (2004) 123:23–42.

    Matsuno A., Kawaguti S. An ultrastructural study of the polyp and strobila of Atorella japonica (Cnidaria, Coronatae) with respect to muscles and nerves. Hydrobiologia (1991) 216/217:39–43.

    Mayer A. G. Medusae of the World (1910) Washington, D.C.: Carnegie Institution.

    Medina M., Collins A. G., Silberman J. D., Sogin M. L. Evaluating hypotheses of basal animal phylogeny using complete sequences of large and small subunit rRNA. Proc. Natl. Acad. Sci. USA (2001) 98:9707–9712.[Abstract/Free Full Text]

    Muller P., Seipel K., Yanze N., Reber-Muller S., Streitwolf-Engel R., Stierwald M., Spring J., Schmid V. Evolutionary aspects of developmentally regulated helix-loop-helix transcription factors in striated muscle of jellyfish. Dev. Biol. (2003) 255:216–229.[CrossRef][Web of Science][Medline]

    Naumov D. V. Hydroids and Hydromedusae of the USSR. Translated from Russian (1969) (1960) Jerusalem. Israel Program for Scientific Translations.

    Odorico D. M., Miller D. J. Internal and external relationships of the Cnidaria: Implications of primary and predicted secondary structure of the 5'-end of the 23S-like rDNA. Proc. Roy. Soc. B (1997) 264:77–82.[Medline]

    Pagel M., Meade A. A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data. Syst. Biol. (2004) 53:571–581.[Abstract/Free Full Text]

    Pagel M., Meade A. Mixture models in phylogenetic inference. In: Mathematics of evolution and phylogeny—Gascuel O., ed. (2005) Oxford, UK: Clarendon Press. 121–142.

    Petersen K. W. Development of coloniality in Hydrozoa. In: Biology and systematics of colonial organisms—Larwood G., Rosen B. R., eds. (1979) New York: Academic Press. 105–139.

    Petersen K. W. Evolution and taxonomy in capitate hydroids and medusae (Cnidaria: Hydrozoa). Zool. J. Linn. Soc. (1990) 100:101–231.[CrossRef]

    Posada D., Crandall K. A. Model Test: Testing the model of DNA substitution. Bioinformatics (1998) 14:817–818.[Abstract/Free Full Text]

    Raftery A. E. Hypothesis testing and model selection. In: Markov chain Monte Carlo in practice—Gilks W. R., Richardson S., Spiegelhalter D. J., eds. (1996) London: Chapman and Hall. 163–188.

    Rees W. J. The relationships of Moerisia lyonsi Boulenger and the family Moerisiidae, with capitate hydroids. Proc. Zool. Soc. Lond (1958) 130:537–545.

    Salvini-Plawen. On the origin and evolution of the lower Metazoa. Z. Zool. Syst. Evolutions. (1978) 16:40–88.

    Satterlie R. A. Neuronal control of swimming in jellyfish: A comparative story. Can. J. Zool. (2002) 80:1654–1669.[CrossRef]

    Satterlie R. A. Central control of swimming in the cubomedusan jellyfish Carybdea rastonii. J. Comp. Physiol (1979) 133:357–367.[CrossRef]

    Satterlie R. A., Spencer A. N. Reply to Leonard's paper Cubomedusae belong to the class Cubozoa, not Scyphozoa. Nature (1980) 284:377.

    Scholtz G. Coelenterata versus Acrosomata—zur Position der Rippenquallen (Ctenophora) im phylogenetischen System der Metazoa. Sber. Ges. Naturf. Freunde Berlin (2004) 43:15–33.

    Schuchert P. Phylogenetic analysis of the Cnidaria. Z. Zool. Sys. Evolutionsf. (1993) 31:161–173.

    Schuchert P. The marine fauna of New Zealand: Athecate hydroids and their Medusae (Cnidaria: Hydrozoa) (1996) Wellington: New Zealand Oceanographic Institute.

    Seipel K., Schmid V. Evolution of striated muscle: Jellyfish and the origin of triploblasty. Dev. Biol. (2005) 282:14–26.[CrossRef][Web of Science][Medline]

    Shimodaira H. An approximately unbiased test of phylogenetic tree selection. Syst. Biol. (2002) 51:492–508.[Abstract/Free Full Text]

    Shimodaira H., Hasegawa M. Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol. Biol. Evol. (1999) 16:1114–1116.[Web of Science]

    Shimodaira H., Hasegawa M. CONSEL: For assessing the confidence of phylogenetic tree selection. Bioinformatics (2001) 17:1246–1247.[Abstract/Free Full Text]

    Song J.-I., Won J. H. Systematic relationship of the anthozoan orders based on the partial nuclear 18S rDNA sequences. Kor. J. Biol. Sci. (1997) 1:43–52.

    Spring J., Yanze N., Josch C., Middel A. M., Winninger B., Schmid V. Conservation of Brachyury, Mef2, and Snail in the myogenic lineage of jellyfish: A connection to the mesoderm of Bilateria. Dev. Biol. (2002) 244:372–384.[CrossRef][Web of Science][Medline]

    Spring J., Yanze N., Middel A. M., Stierwald M., Groger H., Schmid V. The mesoderm specification factor twist in the life cycle of jellyfish. Dev. Biol. (2000) 228:363–375.[CrossRef][Web of Science][Medline]

    Stangl K., Salvini-Plawen L. V., Holstein T. W. Staging and induction of medusa metamorphosis in Carybdea marsupialis(Cnidaria, Cubozoa). Vie Mil (2002) 52:131–140.

    Stepan'yjants S. D., Anokhin B. A., Kuznetsova V. G. Hydrida composition and place in the system of Hydroidea (Cnidaria: Hydrozoa). Trudy Zool. Inst. Ross. Akad. Nauk. (2000) 286:155–162.

    Straehler-Pohl I., Jarms G. Life cycle of Carybdea marsupialis Linnaeus, 1758 (Cubozoa, Carybdeidae) reveals metamorphosis to be a modified strobilation. Mar. Biol. (2005) 147(6):1271–1277.[CrossRef]

    Studebaker J. P. Development of the cubomedusa. In: Carybdea marsupialis (1972) Mayaguez: University of Puerto Rico. M.S. thesis.

    Swofford D. L. PAUP*. 4.0 Phylogenetic analysis using parsimony (*and other methods), version Ver. 4B10 (2000) Sunderland, Massachusetts: Sinauer Associates.

    Swofford D. L., Olsen G. J., Waddell P. J., Hillis D. M. Phylogenetic inference. In: Molecular systematics—Hillis D. M., Moritz C., Mable B. K., eds. (1996) Sunderland, Massachusetts: Sinauer Associates. 407–514.

    Thiel H. The evolution of the Scyphozoa, a review. In: Cnidaria and their Evolution—Rees W. J., ed. (1966) London: Academic Press. 77–117.

    Totton A. K. A synopsis of the Siphonophora (1965) London: British Museum of Natural History.

    Uchida T. The systematic position of the Stauromedusae. Publ. Seto Mar. Biol. Lab. (1972) 20:133–139.

    Van de, Peer Y., Jansen J., De Rijk P., De Wachter R. Database on the structure of small ribosomal subunit RNA. Nucleic Acids Res. (1997) 25:111–116.[Abstract/Free Full Text]

    Van Iten H., De Moraes Leme J., Simões M. G., Marques A. C., Collins A. G. Reassessment of the phylogenetic position of conulariids (?Ediacaran-Triassic) within the subphylum Medusozoa (phylum Cnidaria). J. Syst. Palaeont.

    Voigt O., Collins A. G., Pearse V. B., Pearse J. S., Ender A., Hadrys H., Schierwater B. Placozoa - no longer a phylum of one. Curr. Biol. (2004) 14:R944–R945.[CrossRef][Web of Science][Medline]

    Werner B. New investigations on systematics and evolution of the class Scyphozoa and the phylum Cnidaria. Pub. Seto Mar. Biol. Lab. (1973) 20:35–61.

    Werner B., Cutress C. E., Studebaker J. P. Life cycle of Tripedalia cystophora Conant (Cubomedusae). Nature (1971) 232:582–583.[CrossRef][Medline]

    Won J. H., Rho B. J., Song J. I. A phylogenetic study of the Anthozoa (phylum Cnidaria) based on morphological and molecular characters. Cor. Reefs (2001) 20:39–50.[CrossRef]

    Yamaguchi M., Hartwick R. Early life history of the sea wasp, Chironex fleckeri (Class Cubozoa). In: Developmental and cellular biology of Coelenterates—Tardent P., Tardent R., eds. (1980) Amsterdam: Elsevier/North Holland Biomedical Press. 11–16.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Syst BiolHome page
L. Leclere, P. Schuchert, C. Cruaud, A. Couloux, and M. Manuel
Molecular Phylogenetics of Thecata (Hydrozoa, Cnidaria) Reveals Long-Term Maintenance of Life History Traits despite High Frequency of Recent Character Changes
Syst Biol, October 1, 2009; 58(5): 509 - 526.
[Abstract] [Full Text] [PDF]


Home page
Biol. Bull.Home page
C. Nielsen
How Did Indirect Development With Planktotrophic Larvae Evolve?
Biol. Bull., June 1, 2009; 216(3): 203 - 215.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Collins, A. G.
Right arrow Articles by Schierwater, B.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Collins, A. G.
Right arrow Articles by Schierwater, B.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?