Fast evolving 18S rRNA sequences from Solenogastres (Mollusca) resist standard PCR amplification and give new insights into mollusk substitution rate heterogeneity
© Meyer et al; licensee BioMed Central Ltd. 2010
Received: 4 June 2009
Accepted: 9 March 2010
Published: 9 March 2010
The 18S rRNA gene is one of the most important molecular markers, used in diverse applications such as molecular phylogenetic analyses and biodiversity screening. The Mollusca is the second largest phylum within the animal kingdom and mollusks show an outstanding high diversity in body plans and ecological adaptations. Although an enormous amount of 18S data is available for higher mollusks, data on some early branching lineages are still limited. Despite of some partial success in obtaining these data from Solenogastres, by some regarded to be the most "basal" mollusks, this taxon still remained problematic due to contamination with food organisms and general amplification difficulties.
We report here the first authentic 18S genes of three Solenogastres species (Mollusca), each possessing a unique sequence composition with regions conspicuously rich in guanine and cytosine. For these GC-rich regions we calculated strong secondary structures. The observed high intra-molecular forces hamper standard amplification and appear to increase formation of chimerical sequences caused by contaminating foreign DNAs from potential prey organisms. In our analyses, contamination was avoided by using RNA as a template. Indication for contamination of previously published Solenogastres sequences is presented. Detailed phylogenetic analyses were conducted using RNA specific models that account for compensatory substitutions in stem regions.
The extreme morphological diversity of mollusks is mirrored in the molecular 18S data and shows elevated substitution rates mainly in three higher taxa: true limpets (Patellogastropoda), Cephalopoda and Solenogastres. Our phylogenetic tree based on 123 species, including representatives of all mollusk classes, shows limited resolution at the class level but illustrates the pitfalls of artificial groupings formed due to shared biased sequence composition.
The small subunit (SSU) 18S rRNA gene is one of the most frequently used genes in phylogenetic studies (see below) and an important marker for random target PCR in environmental biodiversity screening . In general, rRNA gene sequences are easy to access due to highly conserved flanking regions allowing for the use of universal primers . Their repetitive arrangement within the genome provides excessive amounts of template DNA for PCR, even in smallest organisms. The 18S gene is part of the ribosomal functional core and is exposed to similar selective forces in all living beings . Thus, when the first large-scale phylogenetic studies based on 18S sequences were published - first and foremost Field et al.'s  phylogeny of the animal kingdom - the gene was celebrated as the prime candidate for reconstructing the metazoan tree of life. And in fact, 18S sequences later provided evidence for the splitting of Ecdysozoa and Lophotrochozoa [, see also ], thus contributing to the most recent revolutionary change in our understanding of metazoan relationships. Methodological innovation within the last years came from the incorporation of secondary structure into phylogenetic analyses. In particular RNA specific substitution models considering paired sites in rRNA genes have been shown to outperform standard DNA models [7–12].
During recent years and with increased numbers of taxa included into molecular phylogenies, however, two problems became apparent. First, there are prevailing sequencing impediments in representatives of certain taxa, such as the mollusk classes Solenogastres and Tryblidia [13, 14], selected bivalve taxa (pers. comment H. Dreyer), and the enigmatic crustacean class Remipedia (pers. comment H. Glenner). Failure to obtain 18S sequences of single taxa is considered a common phenomenon but is rarely ever reported. Secondly, in contrast to initially high hopes, 18S cannot resolve nodes at all taxonomic levels and its efficacy varies considerably among clades. This has been discussed as an effect of rapid ancient radiation within short periods . Multigene analyses are currently thought to give more reliable results for tracing deep branching events in Metazoa but 18S still is extensively used in phylogenetic analyses.
Considering the wide range of studies based on the 18S gene as a molecular marker, both sequencing problems and the applicability of 18S for phylogenetic inferences need to be scrutinized. To address these questions, we focus on the Mollusca, the second largest animal phylum. There are eight higher taxa (classes) defined within Mollusca: the aplacophoran Solenogastres (= Neomeniomorpha) and Caudofoveata (= Chaetodermomorpha), two small clades with about 250 and 150 currently described species; Polyplacophora (ca. 920 species); and the conchiferan clades Tryblidia (= Monoplacophora; 29 species described), Scaphopoda (ca. 520 species), Cephalopoda (ca. 1,000 species), Bivalvia (ca. 30,000 species), and Gastropoda (40,000-150,000 species). The monophyly of the phylum is well established based on morphological characters, but 18S phylogenies often show the Mollusca as polyphyletic or paraphyletic with low resolution of the deeper nodes [e.g: [16–19]]. One of the main deficiencies of all published studies on mollusk phylogeny is the underrepresentation of the minor taxa Solenogastres, Caudofoveata, and Monoplacophora. The aplacophoran Solenogastres and Caudofoveata together with Polyplacophora (chitons) are traditionally thought to represent early branching taxa within Mollusca, but their relative position is still controversially discussed [for review see [20, 21]]. Despite of their key position in mollusk phylogeny and evolution, the number of molecular phylogenetic studies including any aplacophoran mollusk species is low: Winnipeninckx et al. ; Okusu et al. ; Lindgren et al. ; Passamaneck et al ; Giribet et al. ; Dunn et al. ; Wägele et al. ; Wilson et al. . To date there are no more than four caudofoveate and five solenogaster (partial) 18S sequences published in Genbank, the taxon sampling spanning no more than two genera of Caudofoveata and Solenogastres, respectively.
In solenogasters, Okusu and Giribet  described severe contamination issues caused by cnidarian prey. Here, we specifically address contamination issues and point out technical problems hampering amplification during PCR. This is important not only for prompting representative taxonomic sampling for phylogenetic analyses, but also for avoiding under- or overestimation of biodiversity in environmental screening programs. Moreover, we evaluate the usefulness of the 18S gene for phylogenetic inferences by combining our new authentic solenogaster 18S data with published molluscan sequences and - based on an extended taxon sampling - analyze sequence divergence and compositional heterogeneity within the phylum.
Initial experiments using standard PCR protocols and genomic DNA from starved specimens resulted in sequences from prey organisms and epibionts, or in chimerical PCR products. At best, cDNA templates led to shortened 18S sequences. Finally, authentic solenogaster 18S sequences for three species were obtained via isolation of total RNA from starved specimens, followed by reverse transcription and utilizing additives for GC-rich templates. 10% DMSO was applied in sequencing reactions.
Sequence composition and alignment
Sequence composition within selected taxa (no = number of specimens) characterized using the disparity index (I D ), compositional distance (Dc), percent GC content, gene length in base pairs (bp) and average branchlengths (brls) with standard deviation (σ).
Heterodonta (ex. Anom.)
Results from the disparity analyses on the 123 taxon alignment with the disparity index between taxa depicted in the lower triangle and the mcmc-test results for the homogeneity of substitution pattern in the upper triangle (% of taxa).
Heterodonta s.A. (14)
GC rich 18S and the chimera problem
Amplification problems are a frequent phenomenon, even if using standard markers such as the 18S rRNA gene. Secondary structures and GC-rich sections as pronounced as for the solenogaster 18S RNAs shown here, have not been reported earlier. May GC-rich sequences cause hampered PCR reactions in other taxa, too? Based on our results, we assume this to be likely, but equally startling as the failure to obtain gene sequences is the danger to produce chimerical products.
GC-rich sequences demand higher melting temperatures to be converted into and kept as single stranded DNA molecules. Under standard amplification conditions, sequence sections exclusively composed of GC residues will terminate elongation steps in PCR by forming highly stable stems or GC clamps causing amplification to be refractory. Incompletely extended primers can anneal to heterologous 18S sequences and, despite of some degree of nucleotide mismatching, they will often be completed in the subsequent polymerization step resulting in chimeras . Alternatively, compatible priming sequences from genes with lower GC contents are favored and successfully amplified . Single stranded sequences with pronounced GC-stem regions also may lead the Taq polymerase to skip the 'locked' sections, which leads to shortened PCR products lacking the stem-loop regions [e.g.: ]. This potentially explains the lack of data for the GC-rich regions in the previously published solenogaster sequences, in contrast to the Wirenia argentea, Simrothiella margaritacea, and Micromenia fodiens sequences described herein. This interpretation is corroborated by the severe amplification problems we experienced when using recombinant plasmids as templates under standard PCR conditions (denaturing step at 94°C). The cloned Wirenia argentea 18S inserts could not be amplified using conventional PCR protocols but resulted in blank agarose gels (not documented), although the number of templates accessible from clone amplification usually exceeds the gene copy number of genomic DNA extractions by far. Failure to amplify solenogaster 18S under such conditions demonstrates that clean starting material is just the first step. In addition to ours, a number of protocols for difficult or GC rich templates are published [36, 37], supplemented by several commercial kits.
The solenogaster midgut is extremely voluminous and can hold undigested food that may provide considerable amounts of DNA templates leading to non-target or chimerical amplification products (see above). To avoid contamination with prey organisms, Okusu and Giribet  suggested the use of gonad tissue, larval material or of species that are not predators of metazoan organisms. The first two options are feasible where material is available, but the third option can be misleading when the diet is unknown. A prey-predator relationship between Helicoradomenia and non-cnidarian metazoan organisms, probably polychaetes, has been proposed earlier based on transmission electron microscopy data .
Starving animals in the laboratory for several days can reduce the amount of contaminating prey tissue considerably , but exogenous DNA was amplified even after starvation times of up to three months when using standard PCR protocols. This may be due to residues of prey cnidocysts held back in the midgut epithelial cells. Thus, we found isolation of total RNA from starved specimens, followed by reverse transcription, to be most effective to avoid contamination. Isolation of total RNA followed by DNase digestion makes up to 95% of rRNA available for RT-PCR and destroys contaminating non-target DNA templates. Due to the ubiquitous presence of RNases causing instability of RNA outside of living cells, the presence of considerable amounts of exogenous RNA within an isolate is highly unusual, but one always has to bear in mind the possibility of contamination with parasites, epibionts, and undigested prey tissue.
Considering the pitfalls of phylogenetic analyses based on a single gene and taking into account well-established knowledge on mollusk relationships, a number of groupings in our tree are ambiguous. This concerns for example the clade composed of Cephalopoda and Solenogastres, where molecular inferences may reflect a shared bias in base composition rather than a sistergroup relationship. This assumption is corroborated by the two experiments testing the influence of taxon sampling by manipulating the gastropod dataset (figure 4). In Patellogastropoda and Heterobranchia there are species with substitution patterns similar to Solenogastres (table 2). For both taxa, the exclusion of species breaking up their long branches resulted in a placement next to species with a similar sequence composition. The branching pattern including all major gastropod lineages is robust in all analyses both within the Gastropoda and the clade (((Solenogastres+Cephalopoda)+Caudofoveata)+Scaphopoda). When we add the heterobranch taxa from the original analysis to the three extremely fast evolving Aeolidina, the latter taxon is correctly placed inside the heterobranch radiation (data not shown). Within Solenogastres our species selection aimed to reflect the taxon diversity by selecting representatives from two different orders (three families), but the trimming of the alignment adapted to the broad taxon sampling diminished the phylogenetic signal at the species level.
Our new and more representative dataset of 123 taxa covered for the first time all mollusk classes, but our analyses across the Mollusca resulted in short internal branches combined with single long branching clades. This scenario is known to be critical in phylogenetic tree reconstructions . The structure aware and non-stationary models greatly improved the results by detecting "problematic" taxa and by resulting in increased node support for well established subgroups. Thus we were able to recover nearly all mollusk classes, even though problematic fast evolving taxa, such as Patellogastropoda, were included. In rejecting the major mollusk clade Conchifera (mollusks with a single or bivalved shell), however, our phylogenetic tree does not reflect traditional views. We detected a sister group of Polyplacophora and Monoplacophora as proposed by Giribet et al.  and Wilson et al. , but not with full support. Both taxa show only minor differences in sequence composition (table 2), which might mislead analyses of class level relationships.
The failure to detect monophyletic Bivalvia probably best demonstrates the limits of 18S analyses in Mollusca using contemporary methods [see also [19, 22, 41]]. Bivalvia are known to show a number of exceptional morphological characters as well as several unusual structural features in the mitochondrial genome . Thus, despite not obviously fast evolving, here the alignment methods may have failed to identify homologous positions correctly. Alternatively, the phylogenetic signal in bivalves may be largely eroded. The close relationship between Scaphopoda and Protobranchia is worth mentioning because it (partly) supports the so-called Diasoma hypothesis, recently rejected by a number of molecular and morphological analyses, but here again proposing a sistergroup relationship between Bivalvia and Scaphopoda ,
Sequence divergence in Mollusca
Elevated substitution rates are known to be gene specific but also characteristic for certain lineages across different genes [48, 49]. High substitution rates in the generally well-conserved 18S gene thus may point to fast evolving taxa. To allow for general conclusions across the Mollusca, we inferred substitution rates for all published molluscan 18S sequences longer than 1600 bp. The solenogaster 18S sequences obtained in this study are among the fastest evolving 18S sequences within the Mollusca. Solenogastres bear a number of assumed ancestral mollusk features, but substitution rates in the 18S gene are nearly doubled in our selected species compared to other presumably "basal" mollusks, such as Caudofoveata and Polyplacophora. If the assumed plesiomorphic morphological characters in Solenogastres are in fact conserved ancestral features, then molecular and morphological rates of evolution are unlinked, at least for the 18S gene. Similar cases of possibly ancestral morphology combined with exceptionally high substitution rates as shown in table 1 are Nautilus, the "living fossil" cephalopods, and Patellogastropoda, the true limpets [50, 51]. A number of factors, for example functionality of proteins and RNAs, generation time, metabolic rate, population size, and life histories, are thought to influence substitution rates [52, 53].
We show that the solenogaster 18S gene has an exceptional base composition, resulting in a number of technical difficulties. Amplification problems are a common phenomenon but can be overcome by combining known methods in a new framework and by employing alternative strategies. We suggest to include assays with modified PCR methods and to creatively vary PCR conditions if amplification fails or if it leads to 'peculiar' results [see also ].
We show that the practical amplification issues of 18S are conquerable, whereas the future of mollusk class level phylogeny appears not to lie in this gene. The incorporation of additional data, such as structural information to identify homologous positions between sequences, is a promising approach to improve alignments and phylogenetic analyses. Within the Mollusca, where the sequence composition is highly variable between clades, both alignment methods and models of evolution used in phylogenetic analyses at date still remain the main bottlenecks in tracing deep phylogeny. Thus, multigene analyses are more effective to resolve such ancient splits, as recently demonstrated .
Animals and PCR conditions
All Solenogastres specimens were collected off Bergen, Norway and starved up to three months at 4°C in natural seawater. RNA was extracted from six to twenty animals of each species using Trizol (Invitrogen, Germany) and followed by DNase digestion on NucleoSpin II columns (Macherey-Nagel, Germany). RT was performed for 45 min at 55°C using 100 ng of RNA as template. Either random or the gene specific R1843 primer  were applied using Superscript III (Invitrogen, Germany) or BcaBEST™ RNA PCR Kit Ver. 1.1 (TaKaRa Bio Inc., Japan). The GC-rich PCR system (Roche) was used to amplify cDNA adding the supplied GC-rich resolution solution to a final concentration of 0.5 M. Three non overlapping fragments were amplified in all species. The thermal profile for the primers 500F (5'GCGGCGCGACGATCGAAATGAGTCGG3') and 2000R (5'GCCTTATCCCGAGCACGCGCGGGGTTCG3'), annealing at 58°C was setup as time incremental PCR with 1'35" + 5"/cycle at 72°C for 25 times, and a final elongation at 68°C for 7'. The two smaller, neighboring 18S regions were annealed at 53°C [primer F19  and 500R (5'CCGACTCATTTCGATCGTCGCGCCGC3')] or 56°C [primer 2000F (5'CGAACCCCGCGCGTGCTCGG3')/R1843 ]. All PCR products were excised from 0.8% TAE agarose gels and isolated using the Agarose-Out kit (EURx, Poland). All large central 18S PCR fragments (500F/2000R) were cloned using the TOPO-TA vector (Invitrogen, Germany) and electrocompetent cells DH10B pulsed at 1800 V. Recombinant plasmids from over night grown liquid cultures (1.5 ml) were isolated using the GeneMatrix Miniprep purification kit (EURx, Poland). Sequencing was performed using a M13 primer, gene specific primer or Wiar900F (CCGCGGCCGCCTCG) and R427 . Cycle sequencing was done applying the Taq Dye Terminator system (Big Dye 3.1, Applied Biosystems) including 10% DMSO in each reaction. All sequences were submitted to Genbank: Wirenia argentea FJ649599, Simrothiella margaritacea FJ649600 and Micromenia fodiens FJ649601.
Melting curves of double stranded DNA products were analyzed with Stitch Profiles . We determined the critical temperature to melt half of the DNA (helicity = 0.5) in Wirenia argentea and than applied 90.7° using the Blossey and Carlon  parameter set. The profiles were aligned to the public available sequences for Epimenia babai (AY212107, AY212106), Epimenia sp. (AY377657) and Helicoramenia sp. (AY145377, AY212108). Secondary structures were calculated at 72°C and 65°C with MFold  at the Institut Pasteur webserver http://bioweb2.pasteur.fr.
Sequence composition (table 1)
18S data (>1.6 kb) from 834 mollusks and 40 non-mollusk outgroup species including all sequence data available for Sipunculida (11 species) and Kamptozoa (5 species) were collected from GenBank and aligned using Mafft version 6.7 [61, 62] in the E-INS-i mode. GC content, gene length, base composition distance and the disparity index were calculated from this alignment using Mega version 4.0 . The secondary structure from Mimachlamys varia L49051 was extracted from the European Ribosomal Database and used to fit a consensus secondary structure on the alignment using RNA-Salsa  in combination with the Alicut perl script http://zfmk.de/web/Forschung/Abteilungen/AG_Wgele/Software/Utilities/index.en.html. Ambiguously aligned hypervariable loop regions within the 870 taxa alignment were identified by eye and coded as binary absence (gap) and present (any base) data  before fed into RAxML version 7.2.2 . The S16 model was used for paired sites, the GTR + Γ for non paired sites and the two states BIN model for ambiguously aligned regions. The resulting ML tree was used to infer the branchlengths from each species tip to the most basal node (last common ancestor of Priapulus caudatus and all remaining lophotrochozoans) using the tool Branchlength at the HIV database http://www.hiv.lanl.gov/content/sequence/HIV/HIVTools.html.
Alignment and phylogenetic inferences
The species chosen for phylogenetic inferences aimed to represent the molecular diversity observed in the 870 taxa alignment. If possible, short branching species were preferred. An alignment of 123 species was completed using Mafft (see above) and subsequently trimmed using Aliscore . A guide tree was constructed from this alignment in RAxML 7.2.2 using the GTR + Γ model and a constrain tree for the molluscan classes. The resulting ML tree was used as guide tree in Prank+F  with two subsequent alignment steps using the flags -once, -realbranches and maxbranches = 0.15 respectively. A consensus secondary structure was provided as described in the previous section. The same alignment procedure was repeated for a second alignment including three Aeolidina (Heterobranchia). The final alignment was edited by eye and trimmed using a python script (available from the authors) in combination with Alicut to discard columns with a threshold of 90% missing data while maintaining a valid secondary structure file. Prank is very conservative in detecting homologous positions and thus this relatively low cut off already provides dense high quality alignments. The trimmed 123 taxa alignment was used to calculate the disparity index and to test the substitution pattern with 1,000 Monte Carlo replicates as implemented in Mega 4.0.
ML trees for all 14 doublet models implemented in RAxML were calculated to select the best fitting substitution model from the 6-, 7-, and 16-state category of models, respectively, applying the AIC as described in Tsagkogeorga et al. . For each of these models we conducted 200 inferences on the original alignment to search for the best ML tree and 1,000 bootstrap replicates enabling the bootstopping option for the frequencies based criterion . The inferences were repeated with RY (G, A = R; C, T = Y) coded loop regions. The wRF distance between the three ML trees was calculated using RAxML 7.2.2.
Phase version 2.0 http://www.bioinf.manchester.ac.uk/resources/phase/ was used for the bayesian inference. The burnin for each setting was determined conducting fully sampled 3,000,000 generation pre-runs (burnin = 0). We adapted the settings from the example control file tol40.mcmc according our dataset in changing the models to GTR + Γ6 and S7B + Γ6 with the proposal priorities for model 1 to = 7 and model 2 to = 24. All chains sampled 10,000,000 generations, discarding 2,000,000 generations as burnin. File conversion was performed using phase2tracer.pl http://hymenoptera.tamu.edu/rna/download.php. The best model found was optimized using the phase tool Optimizer to generate a model block with fixed substitution rates for a time heterogeneous run. The time-heterogeneous run estimated the base frequencies separately for three independent (GTR + Γ6, S7B + Γ6) models. Four chains run analogous as described for the mixed model above and additionally a fifth chain sampled 50,000,000 generations, discarding 10,000,000 generations, to determine the consensus tree and Bayesian posterior probabilities. Convergence of likelihood values and parameters was monitored using Tracer .
We gratefully acknowledge Prof A. Krause for the access to the computer cluster at the FH-Bingen (Bioinformatik) and Prof J. Markl for general support. Thanks to Christoffer Schander for inviting AM to University of Bergen for sampling and laboratory work (funding provided by Bergens Forskningsstiftelse, Norway, project nr. 481103). We are grateful to Mario Dejung for providing two Python scripts and to Karen Meusemann for advice using the Phase package. We thank Nerida Wilson for providing the 18S sequence from Laevipilina hyalina prior to release in Genbank and Eivind Tøstesen for providing high quality graphics from the Stitch Profile analyses. Thanks to two anonymous reviewers for their valuable comments and especially to one of the reviewers for her/his detailed suggestions to improve the manuscript. Funding has been provided by the University of Bergen free researcher-initiated program (project nr. 226270 Friforsk_08_Todt) and the Deutsche Forschungsgemeinschaft within the priority program Deep Metazoan Phylogeny (Li 998/3-1).
- Chenuil A: Choosing the right molecular genetic markers for studying biodiversity: from molecular evolution to practical aspects. Genetica. 2006, 127 (1): 101-120. 10.1007/s10709-005-2485-1.PubMedGoogle Scholar
- Hillis DM, Dixon MT: Ribosomal DNA: molecular evolution and phylogenetic inference. The Quarterly Review of Biology. 1991, 66 (4): 411-453. 10.1086/417338.View ArticlePubMedGoogle Scholar
- Moore PB, Steitz TA: The involvement of RNA in ribosome function. Nature. 2002, 418 (6894): 229-235. 10.1038/418229a.View ArticlePubMedGoogle Scholar
- Field KG, Olsen GJ, Lane DJ, Giovannoni SJ, Ghiselin MT, Raff EC, Pace NR, Raff RA: Molecular phylogeny of the animal kingdom. Science. 1988, 239 (4841): 748-753. 10.1126/science.3277277.View ArticlePubMedGoogle Scholar
- Aguinaldo AMA, Turbeville JM, Linford LS, Rivera MC, Garey JR, Raff RA, Lake JA: Evidence for a clade of nematodes, arthropods and other moulting animals. Nature. 1997, 387 (6632): 489-493. 10.1038/387489a0.View ArticlePubMedGoogle Scholar
- Halanych KM: The new view of animal phylogeny. Annual Review of Ecology Evolution and Systematics. 2004, 35: 229-256. 10.1146/annurev.ecolsys.35.112202.130124.View ArticleGoogle Scholar
- Telford M, Wise M, Gowri-Shankar V: Consideration of RNA secondary structure significantly improves likelihood-based estimates of phylogeny: Examples from the Bilateria. Molecular Biology and Evolution. 2005, 22 (4): 1129-1136. 10.1093/molbev/msi099.View ArticlePubMedGoogle Scholar
- von Reumont B, Meusemann K, Szucsich N, Dell'Ampio E, Gowri-Shankar V, Bartel D, Simon S, Letsch H, Stocsits R, Luan Y-x, et al: Can comprehensive background knowledge be incorporated into substitution models to improve phylogenetic analyses? A case study on major arthropod relationships. BMC Evolutionary Biology. 2009, 9 (1): 119-10.1186/1471-2148-9-119.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsagkogeorga G, Turon X, Hopcroft RR, Tilak MK, Feldstein T, Shenkar N, Loya Y, Huchon D, Douzery EJ, Delsuc F: An updated 18S rRNA phylogeny of tunicates based on mixture and secondary structure models. BMC Evolutionary Biology. 2009, 9: 187-10.1186/1471-2148-9-187.PubMed CentralView ArticlePubMedGoogle Scholar
- Jow H, Hudelot C, Rattray M, Higgs P: Bayesian phylogenetics using an RNA substitution model applied to early mammalian evolution. Molecular Biology and Evolution. 2002, 19 (9): 1591-1601.View ArticlePubMedGoogle Scholar
- Dohrmann M, Janussen D, Reitner J, Collins AG, Worheide G: Phylogeny and Evolution of Glass Sponges (Porifera, Hexactinellida). Systematic Biology. 2008, 57 (3): 388-405. 10.1080/10635150802161088.View ArticlePubMedGoogle Scholar
- Mallatt J, Craig CW, Yoder MJ: Nearly complete rRNA genes assembled from across the metazoan animals: Effects of more taxa, a structure-based alignment, and paired-sites evolutionary models on phylogeny reconstruction. Molecular Phylogenetics and Evolution. 2009.Google Scholar
- Giribet G, Okusu A, Lindgren AR, Huff SW, Schroedl M, Nishiguchi MK: Evidence for a clade composed of molluscs with serially repeated structures: Monoplacophorans are related to chitons. Proceedings of the National Academy of Sciences. 2006, 103 (20): 7723-7728. 10.1073/pnas.0602578103.View ArticleGoogle Scholar
- Okusu A, Giribet G: New 18S rRNA sequences from neomenioid aplacophorans and the possible origin of persistent exogenous contamination. Journal of Molluscan Studies. 2003, 69 (4): 385-387. 10.1093/mollus/69.4.385.View ArticleGoogle Scholar
- Abouheif E, Zardoya R, Meyer A: Limitations of Metazoan 18S rRNA Sequence Data: Implications for Reconstructing a Phylogeny of the Animal Kingdom and Inferring the Reality of the Cambrian Explosion. Journal of Molecular Evolution. 1998, 47 (4): 394-405. 10.1007/PL00006397.View ArticlePubMedGoogle Scholar
- Adoutte A, Balavoine G, Lartillot N, Lespinet O, Prud'homme B, de Rosa R: The new animal phylogeny: Reliability and implications. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (9): 4453-4456. 10.1073/pnas.97.9.4453.PubMed CentralView ArticlePubMedGoogle Scholar
- Peterson KJ, Eernisse DJ: Animal phylogeny and the ancestry of bilaterians: inferences from morphology and 18S rDNA gene sequences. Evolution & Development. 2001, 3 (3): 170-205. 10.1046/j.1525-142x.2001.003003170.x.View ArticleGoogle Scholar
- Passamaneck Y, Halanych KM: Lophotrochozoan phylogeny assessed with LSU and SSU data: Evidence of lophophorate polyphyly. Molecular Phylogenetics and Evolution. 2006, 40 (1): 20-28. 10.1016/j.ympev.2006.02.001.View ArticlePubMedGoogle Scholar
- Passamaneck YJ, Schander C, Halanych KM: Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences. Molecular Phylogenetics and Evolution. 2004, 32 (1): 25-38. 10.1016/j.ympev.2003.12.016.View ArticlePubMedGoogle Scholar
- Haszprunar G, Schander C, Halanych KM: Relationships of higher molluscan taxa. 2008, Berkley, Los Angeles, London: University of California PressView ArticleGoogle Scholar
- Todt C, Okusu A, Schander C, Schwabe E: Phylogeny and Evolution of the Mollusca. 2008, Berkley, Los Angeles, London: University of California PressGoogle Scholar
- Winnepenninckx B, Backeljau T, De Wachter R: Investigation of molluscan phylogeny on the basis of 18S rRNA sequences. Molecular Biology and Evolution. 1996, 13 (10): 1306-1317.View ArticlePubMedGoogle Scholar
- Lindgren AR, Giribet G, Nishiguchi MK: A combined approach to the phylogeny of Cephalopoda (Mollusca). Cladistics. 2004, 20 (5): 454-486. 10.1111/j.1096-0031.2004.00032.x.View ArticleGoogle Scholar
- Dunn C, Hejnol A, Matus D, Pang K, Browne W, Smith S, Seaver E, Rouse G, Obst M, Edgecombe G: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452 (7188): 745-749. 10.1038/nature06614.View ArticlePubMedGoogle Scholar
- Wägele JW, Letsch H, Klussmann-Kolb A, Mayer C, Misof B, Wägele H: Phylogenetic support values are not necessarily informative: the case of the Serialia hypothesis (a mollusk phylogeny). Frontiers in Zoology. 2009, 6 (1): 12-10.1186/1742-9994-6-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson NG, Rouse GW, Giribet G: Assessing the molluscan hypothesis Serialia (Monoplacophora + Polyplacophora) using novel molecular data. Molecular Phylogenetics and Evolution. 2009, 54 (1): 187-193. 10.1016/j.ympev.2009.07.028.View ArticlePubMedGoogle Scholar
- Tostesen E: A stitch in time: Efficient computation of genomic DNA melting bubbles. Algorithms for Molecular Biology. 2008, 3 (10):Google Scholar
- Kumar S, Gadagkar SR: Disparity Index: A Simple Statistic to Measure and Test the Homogeneity of Substitution Patterns Between Molecular Sequences. Genetics. 2001, 158 (3): 1321-1327.PubMed CentralPubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA 4: Molecular evolutionary genetics analysis (MEGA) software version 4. Molecular Biology and Evolution. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Gowri-Shankar V, Rattray M: A reversible jump method for bayesian phylogenetic inference with a nonhomogeneous substitution model. Molecular Biology and Evolution. 2007, 24 (6): 1286-1299. 10.1093/molbev/msm046.View ArticlePubMedGoogle Scholar
- Robinson D, Foulds L: Comparison of weighted labelled trees. Combinatorial Mathematics VI. 1979, 119-126. full_text.View ArticleGoogle Scholar
- Posada D, Buckley TR: Model Selection and Model Averaging in Phylogenetics: Advantages of Akaike Information Criterion and Bayesian Approaches Over Likelihood Ratio Tests. Systematic Biology. 2004, 53 (5): 793-808. 10.1080/10635150490522304.View ArticlePubMedGoogle Scholar
- Wang G, Wang Y: Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes. Applied and Environmental Microbiology. 1997, 63 (12): 4645-4650.PubMed CentralPubMedGoogle Scholar
- Meyerhans A, Vartanian J-P, Wain-Hobson S: DNA recombination during PCR. Nucleic Acids Research. 1990, 18 (7): 1687-1691. 10.1093/nar/18.7.1687.PubMed CentralView ArticlePubMedGoogle Scholar
- Musso M, Bocciardi R, Parodi S, Ravazzolo R, Ceccherini I: Betaine, Dimethyl Sulfoxide, and 7-Deaza-dGTP, a Powerful Mixture for Amplification of GC-Rich DNA Sequences. Journal of Molecular Diagnostics. 2006, 8 (5): 544-550. 10.2353/jmoldx.2006.060058.PubMed CentralView ArticlePubMedGoogle Scholar
- Hube F, Reverdiau P, Iochmann S, Gruel Y: Improved PCR method for amplification of GC-Rich DNA sequences. Molecular Biotechnology. 2005, 31 (1): 81-84. 10.1385/MB:31:1:081.View ArticlePubMedGoogle Scholar
- Frey UH, Bachmann HS, Peters J, Siffert W: PCR-amplification of GC-rich regions: 'slowdown PCR'. Nature Protocols. 2008, 3 (8): 1312-1317. 10.1038/nprot.2008.112.View ArticlePubMedGoogle Scholar
- Todt C, von Salvini-Plawen L: The digestive tract of Helicoradomenia (Solenogastres, Mollusca), aplacophoran molluscs from the hydrothermal vents of the East Pacific Rise. Invertebrate Biology. 2005, 124 (3): 230-253. 10.1111/j.1744-7410.2005.00023.x.View ArticleGoogle Scholar
- Todt C, von Salvini-Plawen L: Ultrastructure of the midgut epithelium of Wirenia argentea (Mollusca: Solenogastres). Journal of Molluscan Studies. 2004, 70 (3): 213-224. 10.1093/mollus/70.3.213.View ArticleGoogle Scholar
- Felsenstein J: Cases in which parsimony or compatibility methods will be positively misleading. Systematic Zoology. 1978, 27: 401-410. 10.2307/2412923.View ArticleGoogle Scholar
- Steiner G, Müller M: What can 18S rDNA do for bivalve phylogeny?. Journal of Molecular Evolution. 1978, 43 (1): 58-70. 10.1007/BF02352300.View ArticleGoogle Scholar
- Gissi C, Iannelli F, Pesole G: Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species. Heredity. 2008, 101 (4): 301-320. 10.1038/hdy.2008.62.View ArticlePubMedGoogle Scholar
- Runnegar B, Pojeta J: Molluscan Phylogeny: The Paleontological Viewpoint. Science. 1974, 186 (4161): 311-317. 10.1126/science.186.4161.311.View ArticlePubMedGoogle Scholar
- Sørensen MV, Giribet G: A modern approach to rotiferan phylogeny: combining morphological and molecular data. Molecular Phylogenetics and Evolution. 2006, 40 (2): 585-608. 10.1016/j.ympev.2006.04.001.View ArticlePubMedGoogle Scholar
- Giribet G, Distel DL, Polz M, Sterrer W, Wheeler WC: Triploblastic Relationships with Emphasis on the Acoelomates and the Position of Gnathostomulida, Cycliophora, Plathelminthes, and Chaetognatha: A Combined Approach of 18S rDNA Sequences and Morphology. Systematic Biology. 2000, 49 (3): 539-562. 10.1080/10635159950127385.View ArticlePubMedGoogle Scholar
- Baguñà J, Martinez P, Paps J, Riutort M: Back in time: a new systematic proposal for the Bilateria. Philosophical Transactions of the Royal Society B : Biological Sciences. 2008, 363 (1496): 1481-1491. 10.1098/rstb.2007.2238.View ArticleGoogle Scholar
- Paps J, Baguñà J, Riutort M: Lophotrochozoa internal phylogeny: new insights from an up-to-date analysis of nuclear ribosomal genes. Proceedings of the Royal Society B : Biological Sciences. 2009, 276 (1660): 1245-1254.PubMed CentralView ArticlePubMedGoogle Scholar
- Hedges BS, Kumar S: Genomic clocks and evolutionary timescales. Trends in Genetics. 2003, 19 (4): 200-206. 10.1016/S0168-9525(03)00053-2.View ArticleGoogle Scholar
- Welch JJ, Bromham L: Molecular dating when rates vary. Trends in Ecology & Evolution. 2005, 20 (6): 320-327. 10.1016/j.tree.2005.02.007.View ArticleGoogle Scholar
- Ponder WF, Lindberg DR: Towards a phylogeny of gastropod molluscs: an analysis using morphological characters. Zoological Journal of the Linnean Society. 1997, 119 (2): 83-265. 10.1111/j.1096-3642.1997.tb00137.x.View ArticleGoogle Scholar
- Harasewych MG, McArthur AG: A molecular phylogeny of the Patellogastropoda (Mollusca : Gastropoda). Marine Biology. 2000, 137 (2): 183-194. 10.1007/s002270000332.View ArticleGoogle Scholar
- Thomas JA, Welch JJ, Woolfit M, Bromham L: There is no universal molecular clock for invertebrates, but rate variation does not scale with body size. Proceedings of the National Academy of Sciences of the United States of America. 2006, 103 (19): 7366-7371. 10.1073/pnas.0510251103.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith SA, Donoghue MJ: Rates of Molecular Evolution Are Linked to Life History in Flowering Plants. Science. 2008, 322 (5898): 86-89. 10.1126/science.1163197.View ArticlePubMedGoogle Scholar
- Wintzingerode Fv, Goebel UB, Stackebrandt E: Determination of microbial diversity in environmental samples: pitfalls of PCR-based rRNA analysis. FEMS Microbiology Reviews. 1997, 21 (3): 213-229. 10.1111/j.1574-6976.1997.tb00351.x.View ArticleGoogle Scholar
- Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, et al: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452 (7188): 745-749. 10.1038/nature06614.View ArticlePubMedGoogle Scholar
- Elwood H, Olsen G, Sogin M: The small-subunit ribosomal RNA gene sequences from the hypotrichous ciliates Oxytricha nova and Stylonychia pustulata . Molecular Biology and Evolution. 1985, 2 (5): 399-410.PubMedGoogle Scholar
- Turbeville J, Schulz J, Raff R: Deuterostome phylogeny and the sister group of the chordates: evidence from molecules and morphology. Molecular Biology and Evolution. 1994, 11 (4): 648-655.PubMedGoogle Scholar
- Bleidorn C: Phylogenetic relationships and evolution of Orbiniidae (Annelida, Polychaeta) based on molecular data. Zoological Journal of the Linnean Society. 2005, 144 (1): 59-73. 10.1111/j.1096-3642.2005.00160.x.View ArticleGoogle Scholar
- Blossey R, Carlon E: Reparametrizing the loop entropy weights: effect on DNA melting curves. Physical review E, Statistical, nonlinear, and soft matter physics. 2003, 68: 0619111-0619118.View ArticleGoogle Scholar
- Zuker M, Mathews DH, Turner DH, (eds): Turner Algorithms and Thermodynamics for RNA Secondary Structure Prediction: A Practical Guide in RNA Biochemistry and Biotechnology. 1999, Kluwer Academic PublishersGoogle Scholar
- Katoh K, Misawa K, Kuma K-i, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Research. 2002, 30 (14): 3059-3066. 10.1093/nar/gkf436.PubMed CentralView ArticlePubMedGoogle Scholar
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Briefings in Bioinformatics. 2008, 9 (4): 286-298. 10.1093/bib/bbn013.View ArticlePubMedGoogle Scholar
- Stocsits RR, Letsch H, Hertel J, Misof B, Stadler PF: Accurate and efficient reconstruction of deep phylogenies from structured RNAs. Nucleic Acids Research. 2009, 37 (18): 6184-6193. 10.1093/nar/gkp600.PubMed CentralView ArticlePubMedGoogle Scholar
- González D, Cubeta MA, Vilgalys R: Phylogenetic utility of indels within ribosomal DNA and [beta]-tubulin sequences from fungi in the Rhizoctonia solani species complex. Molecular Phylogenetics and Evolution. 2006, 40 (2): 459-470. 10.1016/j.ympev.2006.03.022.View ArticlePubMedGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22 (21): 2688-2690. 10.1093/bioinformatics/btl446.View ArticlePubMedGoogle Scholar
- Misof B, Misof K: A Monte Carlo Approach Successfully Identifies Randomness in Multiple Sequence Alignments : A More Objective Means of Data Exclusion. Systematic Biology. 2009, 58 (1): 21-34. 10.1093/sysbio/syp006.View ArticlePubMedGoogle Scholar
- Loytynoja A, Goldman N: Phylogeny-aware gap placement prevents errors in sequence alignment and evolutionary analysis. Science. 2008, 320 (5883): 1632-1635. 10.1126/science.1158395.View ArticlePubMedGoogle Scholar
- Pattengale ND, Alipour M, Bininda-Emonds ORP, Moret BME, Stamatakis A: How Many Bootstrap Replicates are Necessary?. Proceedings of RECOMB. 2009Google Scholar
- Drummond A, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007, 7 (1): 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.