What can 18S rDNA do for bivalve phylogeny?

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

A family-level Tree of Life for bivalves based on a Sanger-sequencing approach

The systematics of the molluscan class Bivalvia are explored using a 5-gene Sanger-based approach including the largest taxon sampling to date, encompassing 219 ingroup species spanning 93 (or 82%) of the 113 currently accepted bivalve families. This study was designed to populate the bivalve Tree of Life at the family level and to place many genera into a clear phylogenetic context, but also pointing to several major clades where taxonomic work is sorely needed. Despite not recovering monophyly of Bivalvia or Protobranchia-as in most previous Sanger-based approaches to bivalve phylogeny-our study provides increased resolution in many higher-level clades, and supports the monophyly of Autobranchia, Pteriomorphia, Heteroconchia, Palaeoheterodonta, Heterodonta, Archiheterodonta, Euheterodonta, Anomalodesmata, Imparidentia, and Neoheterodontei, in addition to many other lower clades. However, deep nodes within some of these clades, especially Pteriomorphia and Imparidentia, could not be resolved with confidence. In addition, many families are not supported, and several are supported as non-monophyletic, including Malletiidae, Nuculanidae, Yoldiidae, Malleidae, Pteriidae, Arcidae, Propeamussiidae, Iridinidae, Carditidae, Myochamidae, Lyonsiidae, Pandoridae, Montacutidae, Galeommatidae, Tellinidae, Semelidae, Psammobiidae, Donacidae, Mactridae, and Cyrenidae; Veneridae is paraphyletic with respect to Chamidae, although this result appears to be an artifact. The denser sampling however allowed testing specific placement of species, showing, for example, that the unusual Australian Plebidonax deltoides is not a member of Donacidae and instead nests within Psammobiidae, suggesting that major revision of Tellinoidea may be required. We also showed that Cleidothaerus is sister group to the cementing member of Myochamidae, suggesting that Cleidothaeridae may not be a valid family and that cementation in Cleidothaerus and Myochama may have had a single origin. These results highlight the need for an integrative approach including as many genera as possible, and that the monophyly and relationships of many families require detailed reassessment. NGS approaches may be able to resolve the most recalcitrant nodes in the near future.

A phylogenetic backbone for Bivalvia: an RNA-seq approach

Bivalves are an ancient and ubiquitous group of aquatic invertebrates with an estimated 10 000–20 000 living species. They are economically significant as a human food source, and ecologically important given their biomass and effects on communities. Their phylogenetic relationships have been studied for decades, and their unparalleled fossil record extends from the Cambrian to the Recent. Nevertheless, a robustly supported phylogeny of the deepest nodes, needed to fully exploit the bivalves as a model for testing macroevolutionary theories, is lacking. Here, we present the first phylogenomic approach for this important group of molluscs, including novel transcriptomic data for 31 bivalves obtained through an RNA-seq approach, and analyse these data with published genomes and transcriptomes of other bivalves plus outgroups. Our results provide a well-resolved, robust phylogenetic backbone for Bivalvia with all major lineages delineated, addressing long-standing questions about the monophyly of Protobranchia and Heterodonta, and resolving the position of particular groups such as Palaeoheterodonta, Archiheterodonta and Anomalodesmata. This now fully resolved backbone demonstrates that genomic approaches using hundreds of genes are feasible for resolving phylogenetic questions in bivalves and other animals.

Investigating the Bivalve Tree of Life – an exemplar-based approach combining molecular and novel morphological characters

To re-evaluate the relationships of the major bivalve lineages, we amassed detailed morpho-anatomical, ultrastructural and molecular sequence data for a targeted selection of exemplar bivalves spanning the phylogenetic diversity of the class. We included molecular data for 103 bivalve species (up to five markers) and also analysed a subset of taxa with four additional nuclear protein-encoding genes. Novel as well as historically employed morphological characters were explored, and we systematically disassembled widely used descriptors such as gill and stomach ‘types’. Phylogenetic analyses, conducted using parsimony direct optimisation and probabilistic methods on static alignments (maximum likelihood and Bayesian inference) of the molecular data, both alone and in combination with morphological characters, offer a robust test of bivalve relationships. A calibrated phylogeny also provided insights into the tempo of bivalve evolution. Finally, an analysis of the informativeness of morphological characters showed that sperm ultrastructure characters are among the best morphological features to diagnose bivalve clades, followed by characters of the shell, including its microstructure. Our study found support for monophyly of most broadly recognised higher bivalve taxa, although support was not uniform for Protobranchia. However, monophyly of the bivalves with protobranchiate gills was the best-supported hypothesis with incremental morphological and/or molecular sequence data. Autobranchia, Pteriomorphia, Heteroconchia, Palaeoheterodonta, Archiheterodonta, Euheterodonta, Anomalodesmata and Imparidentia new clade ( = Euheterodonta excluding Anomalodesmata) were recovered across analyses, irrespective of data treatment or analytical framework. Another clade supported by our analyses but not formally recognised in the literature includes Palaeoheterodonta and Archiheterodonta, which emerged under multiple analytical conditions. The origin and diversification of each of these major clades is Cambrian or Ordovician, except for Archiheterodonta, which diverged from Palaeoheterodonta during the Cambrian, but diversified during the Mesozoic. Although the radiation of some lineages was shifted towards the Palaeozoic (Pteriomorphia, Anomalodesmata), or presented a gap between origin and diversification (Archiheterodonta, Unionida), Imparidentia showed steady diversification through the Palaeozoic and Mesozoic. Finally, a classification system with six major monophyletic lineages is proposed to comprise modern Bivalvia: Protobranchia, Pteriomorphia, Palaeoheterodonta, Archiheterodonta, Anomalodesmata and Imparidentia.

Phylogenetic analysis of Thecosomata Blainville, 1824 (holoplanktonic opisthobranchia) using morphological and molecular data

Thecosomata is a marine zooplankton group, which played an important role in the carbonate cycle in oceans due to their shell composition. So far, there is important discrepancy between the previous morphological-based taxonomies, and subsequently the evolutionary history of Thecosomata. In this study, the remarkable planktonic sampling of TARA Oceans expedition associated with a set of various other missions allowed us to assess the phylogenetic relationships of Thecosomata using morphological and molecular data (28 S and COI genes). The two gene trees showed incongruities (e.g. Hyalocylis, Cavolinia), and high congruence between morphological and 28S trees (e.g. monophyly of Euthecosomata). The monophyly of straight shell species led us to reviving the Orthoconcha, and the split of Limacinidae led us to the revival of Embolus inflata replacing Limacina inflata. The results also jeopardized the Euthecosomata families that are based on plesiomorphic character state as in the case for Creseidae which was not a monophyletic group. Divergence times were also estimated, and suggested that the evolutionary history of Thecosomata was characterized by four major diversifying events. By bringing the knowledge of palaeontology, we propose a new evolutionary scenario for which macro-evolution implying morphological innovations were rhythmed by climatic changes and associated species turn-over that spread from the Eocene to Miocene, and were shaped principally by predation and shell buoyancy.

Phylogenetic analysis of four nuclear protein-encoding genes largely corroborates the traditional classification of Bivalvia (Mollusca)

Revived interest in molluscan phylogeny has resulted in a torrent of molecular sequence data from phylogenetic, mitogenomic, and phylogenomic studies. Despite recent progress, basal relationships of the class Bivalvia remain contentious, owing to conflicting morphological and molecular hypotheses. Marked incongruity of phylogenetic signal in datasets heavily represented by nuclear ribosomal genes versus mitochondrial genes has also impeded consensus on the type of molecular data best suited for investigating bivalve relationships. To arbitrate conflicting phylogenetic hypotheses, we evaluated the utility of four nuclear protein-encoding genes-ATP synthase β, elongation factor-1α, myosin heavy chain type II, and RNA polymerase II-for resolving the basal relationships of Bivalvia. We sampled all five major lineages of bivalves (Archiheterodonta, Euheterodonta [including Anomalodesmata], Palaeoheterodonta, Protobranchia, and Pteriomorphia) and inferred relationships using maximum likelihood and Bayesian approaches. To investigate the robustness of the phylogenetic signal embedded in the data, we implemented additional datasets wherein length variability and/or third codon positions were eliminated. Results obtained include (a) the clade (Nuculanida+Opponobranchia), i.e., the traditionally defined Protobranchia; (b) the monophyly of Pteriomorphia; (c) the clade (Archiheterodonta+Palaeoheterodonta); (d) the monophyly of the traditionally defined Euheterodonta (including Anomalodesmata); and (e) the monophyly of Heteroconchia, i.e., (Palaeoheterodonta+Archiheterodonta+Euheterodonta). The stability of the basal tree topology to dataset manipulation is indicative of signal robustness in these four genes. The inferred tree topology corresponds closely to those obtained by datasets dominated by nuclear ribosomal genes (18S rRNA and 28S rRNA), controverting recent taxonomic actions based solely upon mitochondrial gene phylogenies.

A molecular phylogeny of bivalve mollusks: ancient radiations and divergences as revealed by mitochondrial genes

BACKGROUND

Bivalves are very ancient and successful conchiferan mollusks (both in terms of species number and geographical distribution). Despite their importance in marine biota, their deep phylogenetic relationships were scarcely investigated from a molecular perspective, whereas much valuable work has been done on taxonomy, as well as phylogeny, of lower taxa.

METHODOLOGY/PRINCIPAL FINDINGS

Here we present a class-level bivalve phylogeny with a broad sample of 122 ingroup taxa, using four mitochondrial markers (MT-RNR1, MT-RNR2, MT-CO1, MT-CYB). Rigorous techniques have been exploited to set up the dataset, analyze phylogenetic signal, and infer a single final tree. In this study, we show the basal position of Opponobranchia to all Autobranchia, as well as of Palaeoheterodonta to the remaining Autobranchia, which we here propose to call Amarsipobranchia. Anomalodesmata were retrieved as monophyletic and basal to (Heterodonta + Pteriomorphia).

CONCLUSIONS/SIGNIFICANCE

Bivalve morphological characters were traced onto the phylogenetic trees obtained from the molecular analysis; our analysis suggests that eulamellibranch gills and heterodont hinge are ancestral characters for all Autobranchia. This conclusion would entail a re-evaluation of bivalve symplesiomorphies.

Fast evolving 18S rRNA sequences from Solenogastres (Mollusca) resist standard PCR amplification and give new insights into mollusk substitution rate heterogeneity

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Comparative morphology among representatives of main taxa of Scaphopoda and basal protobranch Bivalvia (Mollusca)

This study deals with detailed morphology and anatomy of 4 species of Scaphopoda and 5 species of protobranch Bivalvia. Both classes are traditionally grouped in the taxon Diasoma, which has been questioned by different methodologies, such as molecular and developmental. This study is developed under a phylogenetic methodology with the main concern in performing it in an intelligible and testable methodology. The analyzed Scaphopoda species came from the Brazilian coast and belong to the family Dentaliidae [(1) Coccodentalium carduus; (2) Paradentalium disparile] and Gadiliidae; [(3) Polyschides noronhensis, n. sp. from Fernando de Noronha Archipelago; (4) Gadila braziliensis]. These species represent the main branches of the class Scaphopoda. From protobranch bivalves, representatives of the families Solemyidae [(5) Solemya occidentalis, from Florida; S. notialis, n. sp. from S.E. Brazil], Nuculanidae [(6) Propeleda carpentieri from Florida], and Nuculidae [(7) Ennucula puelcha, from south Brazil] are included. These species represent the main branches of the basal Bivalvia. The descriptions on the anatomy of S. occidentalis and of P. carpentieri are published elsewhere. The remaining are included here, for which a complete taxonomical treatment is performed. Beyond these species, representatives of other taxa are operationally included as part of the ingroup (indices are then shared with them), as a procedure to test the morphological monophyly of Diasoma. These taxa are: two lamellibranch bivalves [(8) Barbatia - Arcidae; (9) Serratina - Tellinidae; both published elsewhere;, and Propilidium (10) Patellogastropoda, and (11) Nautilus, basal Cephalopoda, based on basal taxa. The effective outgroups are (12) Neopilina (Monoplacophora) and (13) Hanleya (Polyplacophora). The phylogenetic analysis based on morphology revealed that the taxon Diasoma is supported by 14 synapomorphies, and is separated from Cyrtosoma (Gastropoda + Cephalopoda). Although they are not the main goal of this paper, the taxa Scaphopoda and Bivalvia are supported by 8 and by 7 synapomorphies respectively. The taxon Protobranchia resulted paraphyletic. Both scaphopod orders resulted monophyletic. The obtained cladogram is: ((((Coccodentalium carduus - Paradentalium disparile) (Polyschides noronhensis - Gadila brasiliensis)) ((Solemya occidentalis - S. notialis) (Propeleda carpenteri (Ennucula puelcha (Barbatia cancellaria - Serratina capsoides))))) (Propilidium curumim - Nautilus pompilius - Lolliguncula brevis)).

Comparative morphology among representatives of main taxa of Scaphopoda and basal protobranch Bivalvia (Mollusca)

This study deals with detailed morphology and anatomy of 4 species of Scaphopoda and 5 species of protobranch Bivalvia. Both classes are traditionally grouped in the taxon Diasoma, which has been questioned by different methodologies, such as molecular and developmental. This study is developed under a phylogenetic methodology with the main concern in performing it in an intelligible and testable methodology. The analyzed Scaphopoda species came from the Brazilian coast and belong to the family Dentaliidae [(1) Coccodentalium carduus; (2) Paradentalium disparile] and Gadiliidae; [(3) Polyschides noronhensis, n. sp. from Fernando de Noronha Archipelago; (4) Gadila braziliensis]. These species represent the main branches of the class Scaphopoda. From protobranch bivalves, representatives of the families Solemyidae [(5) Solemya occidentalis, from Florida; S. notialis, n. sp. from S.E. Brazil], Nuculanidae [(6) Propeleda carpentieri from Florida], and Nuculidae [(7) Ennucula puelcha, from south Brazil] are included. These species represent the main branches of the basal Bivalvia. The descriptions on the anatomy of S. occidentalis and of P. carpentieri are published elsewhere. The remaining are included here, for which a complete taxonomical treatment is performed. Beyond these species, representatives of other taxa are operationally included as part of the ingroup (indices are then shared with them), as a procedure to test the morphological monophyly of Diasoma. These taxa are: two lamellibranch bivalves [(8) Barbatia - Arcidae; (9) Serratina - Tellinidae; both published elsewhere;, and Propilidium (10) Patellogastropoda, and (11) Nautilus, basal Cephalopoda, based on basal taxa. The effective outgroups are (12) Neopilina (Monoplacophora) and (13) Hanleya (Polyplacophora). The phylogenetic analysis based on morphology revealed that the taxon Diasoma is supported by 14 synapomorphies, and is separated from Cyrtosoma (Gastropoda + Cephalopoda). Although they are not the main goal of this paper, the taxa Scaphopoda and Bivalvia are supported by 8 and by 7 synapomorphies respectively. The taxon Protobranchia resulted paraphyletic. Both scaphopod orders resulted monophyletic. The obtained cladogram is: ((((Coccodentalium carduus - Paradentalium disparile) (Polyschides noronhensis - Gadila brasiliensis)) ((Solemya occidentalis - S. notialis) (Propeleda carpenteri (Ennucula puelcha (Barbatia cancellaria - Serratina capsoides))))) (Propilidium curumim - Nautilus pompilius - Lolliguncula brevis)).

A simple and reliable method for DNA extraction from bivalve mantle

A simple and reliable method was developed for extracting genomic DNA from preserved mantle tissues of Pacific oyster Crassostrea gigas for reproducible PCR amplification. The method was based on destruction of the tissue using Proteinase K, Chelex 100 resin, detergents, and urea, followed by preferential capturing of genomic DNA with silica particles. Approximately 5 mg of mantle tissue provided a sufficient quality and quantity of DNA for several hundreds of PCR reactions amplifying the hypervariable mitochondrial DNA intergenic spacer, which is a useful genetic marker for population structure analysis of Pacific oyster. The method can be applied for DNA preparation from not only fresh and frozen but also ethanol-preserved mantle tissues, so this rapid and economical method can serve for investigating a large number of bivalve specimens collected in the field and next transported in ethanol at ambient temperature.

Origin and evolution of Mytilus mussel satellite DNAs

A phylogenetic reconstruction based on the amplification of 3 satellite DNAs (stDNAs) was carried out in 1 crustacean species and 15 bivalve species of the subclass Pteriomorphia (10, subfamily Mytilinae; 1, subfamily Litophaginae; 1, subfamily Modiolinae, all belonging to family Mytilidae; 1, family Arcidae; and 2, family Pectinidae). The sequences obtained showed motifs with high similarity to those of A and B boxes of tRNA promoter regions. Dot-blot hybridizations revealed that the 3 stDNAs are present mainly in high copy numbers for each species of the genus Mytilus, whereas for the other species they appear in low copy numbers. Maximum-parsimony trees evidenced a tendency to group Mytilus clones together, and species containing these sequences as a single copy were distributed among the different mytilids. Finally, the possible origin and evolution of these stDNAs is discussed.

Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences

The Mollusca represent one of the most morphologically diverse animal phyla, prompting a variety of hypotheses on relationships between the major lineages within the phylum based upon morphological, developmental, and paleontological data. Analyses of small-ribosomal RNA (SSU rRNA) gene sequence have provided limited resolution of higher-level relationships within the Mollusca. Recent analyses suggest large-subunit (LSU) rRNA gene sequences are useful in resolving deep-level metazoan relationships, particularly when combined with SSU sequence. To this end, LSU (approximately 3.5 kb in length) and SSU (approximately 2 kb) sequences were collected for 33 taxa representing the major lineages within the Mollusca to improve resolution of intraphyletic relationships. Although the LSU and combined LSU+SSU datasets appear to hold potential for resolving branching order within the recognized molluscan classes, low bootstrap support was found for relationships between the major lineages within the Mollusca. LSU+SSU sequences also showed significant levels of rate heterogeneity between molluscan lineages. The Polyplacophora, Gastropoda, and Cephalopoda were each recovered as monophyletic clades with the LSU+SSU dataset. While the Bivalvia were not recovered as monophyletic clade in analyses of the SSU, LSU, or LSU+SSU, the Shimodaira-Hasegawa test showed that likelihood scores for these results did not differ significantly from topologies where the Bivalvia were monophyletic. Analyses of LSU sequences strongly contradict the widely accepted Diasoma hypotheses that bivalves and scaphopods are closely related to one another. The data are consistent with recent morphological and SSU analyses suggesting scaphopods are more closely related to gastropods and cephalopods than to bivalves. The dataset also presents the first published DNA sequences from a neomeniomorph aplacophoran, a group considered critical to our understanding of the origin and early radiation of the Mollusca.

Molecular phylogeny of the family Pectinidae (Mollusca: Bivalvia) based on mitochondrial 16S and 12S rRNA genes

Pectinidae is a large bivalve family characterised by almost circular, flat shells. Species are distributed worldwide and fall into three life-styles: swimming, byssally attached to hard substrates, and cemented to rocks with one valve. Despite these very different life strategies, pectinid shells are highly conservative in shape and offer few clues for the unravelling of phylogenetic issues. Consequently, phylogenetic studies based on morphological features have not yielded conclusive results. We thus set out to analyse partial sequences of mitochondrial 12S and 16S rRNA genes from 23 species of 16 genera with molecular techniques. The results are largely in contrast, both at the genus and the subfamily level, with the systematic classifications based on adult morphological characters, whereas they agree with the morphological classifications based on the more conserved non-adaptive features.

Bilaterian phylogeny based on analyses of a region of the sodium-potassium ATPase beta-subunit gene

Molecular investigations of deep-level relationships within and among the animal phyla have been hampered by a lack of slowly evolving genes that are amenable to study by molecular systematists. To provide new data for use in deep-level metazoan phylogenetic studies, primers were developed to amplify a 1.3-kb region of the alpha subunit of the nuclear-encoded sodium-potassium ATPase gene from 31 bilaterians representing several phyla. Maximum parsimony, maximum likelihood, and Bayesian analyses of these sequences (combined with ATPase sequences for 23 taxa downloaded from GenBank) yield congruent trees that corroborate recent findings based on analyses of other data sets (e.g., the 18S ribosomal RNA gene). The ATPase-based trees support monophyly for several clades (including Lophotrochozoa, a form of Ecdysozoa, Vertebrata, Mollusca, Bivalvia, Gastropoda, Arachnida, Hexapoda, Coleoptera, and Diptera) but do not support monophyly for Deuterostomia, Arthropoda, or Nemertea. Parametric bootstrapping tests reject monophyly for Arthropoda and Nemertea but are unable to reject deuterostome monophyly. Overall, the sodium-potassium ATPase alpha-subunit gene appears to be useful for deep-level studies of metazoan phylogeny.

Complete sequences of the highly rearranged molluscan mitochondrial genomes of the Scaphopod Graptacme eborea and the bivalve Mytilus edulis

We have determined the complete sequence of the mitochondrial genome of the scaphopod mollusk Graptacme eborea (14,492 nts) and completed the sequence of the mitochondrial genome of the bivalve mollusk Mytilus edulis (16,740 nts). (The name Graptacme eborea is a revision of the species formerly known as Dentalium eboreum.) G. eborea mtDNA contains the 37 genes that are typically found and has the genes divided about evenly between the two strands, but M. edulis contains an extra trnM and is missing atp8, and it has all genes on the same strand. Each has a highly rearranged gene order relative to each other and to all other studied mtDNAs. G. eborea mtDNA has almost no strand skew, but the coding strand of M. edulis mtDNA is very rich in G and T. This is reflected in differential codon usage patterns and even in amino acid compositions. G. eborea mtDNA has fewer noncoding nucleotides than any other mtDNA studied to date, with the largest noncoding region only 24 nt long. Phylogenetic analysis using 2,420 aligned amino acid positions of concatenated proteins weakly supports an association of the scaphopod with gastropods to the exclusion of Bivalvia, Cephalopoda, and Polyplacophora, but it is generally unable to convincingly resolve the relationships among major groups of the Lophotrochozoa, in contrast to the good resolution seen for several other major metazoan groups.

Searching factors causing implausible non-monophyly: ssu rDNA phylogeny of Isopoda Asellota (Crustacea: Peracarida) and faster evolution in marine than in freshwater habitats

This contribution addresses two questions: which alignment patterns are causing non-monophyly of the Asellota and what is the phylogenetic history of this group? The Asellota are small benthic crustaceans occurring in most aquatic habitats. In view of the complex morphological apomorphies known for this group, monophyly of the Asellota has never been questioned. Using ssu rDNA sequences of outgroups and of 16 asellote species from fresh water, littoral marine habitats and from deep-sea localities, the early divergence between the lineages in fresh water and in the ocean, and the monophyly of the deep-sea taxon Munnopsidae are confirmed. Relative substitution rates of freshwater species are much lower than in other isopod species, rates being highest in some littoral marine genera (Carpias and Jaera). Furthermore, more sequence sites are variable in marine than in freshwater species, the latter conserve outgroup character states. Monophyly is recovered with parsimony methods, but not with distance and maximum likelihood analyses, which tear apart the marine from the freshwater species. The information content of alignments was studied with spectra of supporting positions. The scarcity of signal (=apomorphic nucleotides) supporting monophyly of the Asellota is attributed to a short stem-line of this group or to erosion of signal in fast evolving marine species. Parametric boostrapping in combination with spectra indicates that a tree model cannot explain the data and that monophyly of the Asellota should not be rejected even though many topologies do not recover this taxon.

A molecular phylogeny of Heterodonta (Bivalvia) based on small ribosomal subunit RNA sequences

Within Heterodonta, phylogenesis has so far been studied almost exclusively on the basis of morphological data. Results have often been discordant, and an exhaustive molecular approach has not yet been attempted. The present study was undertaken to clarify the phylogenetic relationships obtaining among Heterodonta families through the analysis of 18S rRNA gene. To do this, the whole sequence of this gene was analyzed in 29 species of eight superfamilies of the order of Veneroida (Arcticoidea, Cardioidea, Galeommatoidea, Mactroidea, Solenoidea, Tellinoidea, Tridacnoidea, and Veneroidea) and in two superfamilies of Myoida (Pholaloidea and Myoidea). The study was extended by constructing phylogenetic trees using partial sequences. This strategy made it possible to include 11 additional species by introducing three further superfamilies: Chamoidea, Corbiculoidea, and Hiatellinoidea. At variance with the conclusions reached on the basis of morphological features, the molecular data clearly show that the Myoida species included in this study belong to Veneroida, thus undermining the legitimacy of the division of Heterodonta into two orders, and that considerable differences in the phylogenetic relationships obtain among superfamilies.

Gastropod arginine kinases from Cellana grata and Aplysia kurodai. Isolation and cDNA-derived amino acid sequences

Arginine kinase (AK) was isolated from the radular muscle of the gastropod molluscs Cellana grata (subclass Prosobranchia) and Aplysia kurodai (subclass Opisthobranchia), respectively, by ammonium sulfate fractionation, Sephadex G-75 gel filtration and DEAE-ion exchange chromatography. The denatured relative molecular mass values were estimated to be 40 kDa by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The isolated enzyme from Aplysia gave a Km value of 0.6 mM for arginine and a Vmax value of 13 micromole Pi min(-1) mg protein(-1) for the forward reaction. These values are comparable to other molluscan AKs. The cDNAs encoding Cellana and Aplysia AKs were amplified by polymerase chain reaction, and the nucleotide sequences of 1,608 and 1,239 bp, respectively, were determined. The open reading frame for Cellana AK is 1044 nucleotides in length and encodes a protein with 347 amino acid residues, and that for A. kurodai is 1077 nucleotides and 354 residues. The cDNA-derived amino acid sequences were validated by chemical sequencing of internal lysyl endopeptidase peptides. The amino acid sequences of Cellana and Aplysia AKs showed the highest percent identity (66-73%) with those of the abalone Nordotis and turbanshell Battilus belonging to the same class Gastropoda. These AK sequences still have a strong homology (63-71%) with that of the chiton Liolophura (class Polyplacophora), which is believed to be one of the most primitive molluscs. On the other hand, these AK sequences are less homologous (55-57%) with that of the clam Pseudocardium (class Bivalvia), suggesting that the biological position of the class Polyplacophora should be reconsidered.

Ribosomal DNA and resolution of branching order among the ascomycota: how many nucleotides are enough?

Molecular phylogenies for the fungi in the Ascomycota rely heavily on 18S rRNA gene sequences but this gene alone does not answer all questions about relationships. Particularly problematical are the relationships among the first ascomycetes to diverge, the Archiascomycetes, and the branching order among the basal filamentous ascomycetes, the Euascomycetes. Would more data resolve branching order? We used the jackknife and bootstrapping resampling approach that constitutes the "pattern of resolved nodes" method to address the relationship between number of variable sites in a DNA sequence alignment and support for taxonomic clusters. We graphed the effect of increasing sizes of subsamples of the 18S rRNA gene sequences on bootstrap support for nodes in the Ascomycota tree. Nodes responded differently to increasing data. Some nodes, those uniting the filamentous ascomycetes for example, would still have been well supported with only two thirds of the 18S rRNA gene. Other nodes, like the one uniting the Archiascomycetes as a monophyletic group, would require about double the number of variable sites available in the 18S gene for 95% neighbor-joining bootstrap support. Of the several groups emerging at the base of the filamentous ascomycetes, the Pezizales receive the most support as the first to diverge. Our analysis suggests that we would also need almost three times as much sequence data as that provided by the 18S gene to confirm the basal position for the Pezizales and more than seven times as much data to resolve the next group to diverge. If more data from other genes show the same pattern, the lack of resolution for the filamentous ascomycetes may indicate rapid radiation within this clade.

Arginine kinase from Nautilus pompilius, a living fossil. Site-directed mutagenesis studies on the role of amino acid residues in the Guanidino specificity region

Arginine kinases were isolated from the cephalopods Nautilus pompilius, Octopus vulgaris, and Sepioteuthis lessoniana, and the cDNA-derived amino acid sequences have been determined. Although the origin and evolution of cephalopods have long been obscure, this work provides the first molecular evidence for the phylogenetic position of Cephalopoda in molluscan evolution. A crystal structure for Limulus arginine kinase showed that four amino acid residues (Ser(63), Gly(64), Val(65), and Tyr(68)) are hydrogen-bonded with the substrate arginine. We introduced three independent mutations, Ser(63) --> Gly, Ser(63) --> Thr, and Tyr(68) --> Ser, in Nautilus arginine kinase. One of the mutants had a considerably reduced substrate affinity, accompanied by a decreased V(max). In other mutants, the activity was lost almost completely. It is known that substantial conformational changes take place upon substrate binding in arginine kinase. We hypothesize that the hydrogen bond between Asp(62) and Arg(193) stabilizes the closed, substrate-bound state. Site-directed mutagenesis studies strongly support this hypothesis. The mutant (Asp(62) --> Gly or Arg(193) --> Gly), which destabilizes the maintenance of the closed state and/or perhaps disrupts the unique topology of the catalytic pocket, showed only a very weak activity (0.6-1.5% to the wild-type).

Phylogenetic relationships among mytilidae (Bivalvia): 18S rRNA data suggest convergence in mytilid body plans

Nearly complete sequences were determined for small-subunit (18S) rRNA genes from seven species representative of four subfamilies of Mytilidae: Modiolus modiolus and M. auriculatus (Modiolinae); Lithophaga lithophaga and L. nigra (Lithophaginae); Musculus senhousie and M. discors (Crenellinae); and Hormomya domingensis (Mytilinae). Small-subunit rRNA gene sequences were also determined for Solemya reidi (Subclass Protobranchia), Mya arenaria (Subclass Heterodonta), and Elliptio complanata (Subclass Paleoheterodonta) as outgroup taxa. Phylogenetic analyses including these and other nearly complete bivalve small-subunit rRNA sequences demonstrate support for the monophyly of the family Mytilidae and the subfamilies Crenellinae and Lithophaginae. However, the subfamilies Mytilinae and Modiolinae appear polyphyletic. Likelihood, parsimony, and distance analyses support the placement of H. domingensis (Mytilinae) in a clade with G. demissa (Modiolinae). This clade is distinct from those containing other species traditionally assigned to these two subfamilies. Kishino-Hasegawa tests support these nontraditional relationships, suggesting that the mytiliform and/or modioliform body plans have evolved independently in at least two mytilid lineages.

Phylogenetic analysis of molluscan mitochondrial LSU rDNA sequences and secondary structures

Mollusks are an extraordinarily diverse group of animals with an estimated 200,000 species, second only to the phylum Arthropoda. We conducted a comparative analysis of complete mitochondrial ribosomal large subunit sequences (LSU) of a chiton, two bivalves, six gastropods, and a cephalopod. In addition, we determined secondary structure models for each of them. Comparative analyses of nucleotide variation revealed substantial length variation among the taxa, with stylommatophoran gastropods possessing the shortest lengths. Phylogenetic analyses of the nucleotide sequence data supported the monophyly of Albinaria, Euhadra herklotsi + Cepaea nemoralis, Stylommatophora, Cerithioidea, and when only transversions are included, the Bivalvia. The phylogenetic limits of the mitochondrial LSU rRNA gene within mollusks appear to be up to 400 million years, although this estimate will have to be tested further with additional taxa. Our most novel finding was the discovery of phylogenetic signal in the secondary structure of rRNA of mollusks. The absence of entire stem/loop structures in Domains II, III, and V can be viewed as three shared derived characters uniting the stylommatophoran gastropods. The absence of the aforementioned stem/loop structure explains much of the observed length variation of the mitochondrial LSU rRNA found within mollusks. The distribution of these unique secondary structure characters within mollusks should be examined.

Molecular evidence on the evolution of the Bivalvia

Despite widespread agreement on the monophyly of several major taxa of bivalves, others remain uncertain and the relationships among them are debated. The present study compares new and published morphological phylogenies with new analyses based on 18S gene sequences.

All but one family and all superfamilies in the Bivalvia were monophyletic in all the analyses. Several higher taxa, including most subclasses and orders, were also resolved as monophyletic. Only Myoida shows strong evidence for polyphyly, with at least two origins from Veneroida. Autobranchia was supported as monophyletic in the parsimony analyses. Within Pteriomorphia, Ostreoida is the sister taxon of Pterioida, if not derived from within it, rather than closest to Pectinoida. The numerous points of agreement with morphology based analyses suggests that both types of evidence are converging on a common phylogeny; however, differences remain to be resolved by further study.

Molecular phylogeny of the Bivalvia inferred from 18S rDNA sequences with particular reference to the Pteriomorphia

The classification of the Pteriomorphia, a major division of the Bivalvia with Ordovician origins, is controversial both among palaeontologists and neontologists. To elucidate phylogenetic relationships new near-complete 18S rDNA sequences of 26 Pteriomorphia, three Protobranchia, three Heterodonta, one Anomalodesmata and three Scaphopoda were obtained, aligned with 71 other published molluscan sequences, and analysed with parsimony, maximum likelihood and spectral analysis. Although Bivalvia appear diphyletic due to heterogeneity of substitution rates among lineages, monophyly of Protobranchia, Heteroconchia and Pteriomorphia is supported. The heteroconch Lucinida, Myoida and Venerida are not monophyletic, and Anomalodesmata arise from within Heteroconchia. The basal nodes of Pteriomorphia have little support but two major clades, [Pinnoidea (Ostreoidea + Pterioidea)] and [(Anomioidea + Plicatuloidea) + (Limoidea + Pectinoidea)], are resolved with more confidence. The strongly supported clade of Anomioidea + Plicatuloidea, the separation of Pinnoidea from Pterioidea and most of the intrafamiliar relationships are not in accordance with morphological classifications. Combining these results with selected morphological characters, a phylogenetic hypothesis is proposed showing Mytiloidea and Arcoidea as the basal pteriomorph groups, the latter giving rise to the clade uniting the pinnoid-ostreoid-pterioid and the anomioid-limoid-pectinoid lines.

Comparative sperm ultrastructure in pteriomorphian bivalves with special reference to phylogenetic and taxonomic implications

Pteriomorphian spermatozoa, like those of most other bivalves, are of the classic aquasperm type (conical acrosomal vesicle, short to rod-shaped nucleus, short midpiece composed of two centrioles and a ring of spherical mitochondria, a simple flagellum). Whereas most other bivalve subclasses show at least some defining acrosomal feature(s), this does not appear to be the case within the Pteriomorphia. While this could indicate non-monophyletic status, it also correlates with the fact that the Pteriomorphia are a very old and very successful group of bivalves. Acrosomal similarities suggest a close link between the Ostreoidea and Limoidea (acrosomal vesicle with wedge-shaped apical zone; radiating plates present but not well developed); and between the Pterioidea, Pinnoidea and Pectinoidea (dense anterior layer; very well developed radiating plates). For supposedly closely related taxa, the Arcoidea and Limopsoidea (both Arcoida) differ markedly from each other in acrosomal shape and substructure. The affinities of the Anomioidea and even more so the Mytiloida remain uncertain, the latter possibly connected with the Pterioidea or, more likely, removed from the rest of the Pteriomorphia (mytiloid acrosomes show concentric lamellae). A very close relationship between the Pectinidae and Spondylidae of the Pectinoidea is demonstrated (dense anterior layer of acrosome recurved). Within the Mytilidae (Mytiloidea) there is substantial variation in sperm morphology between supraspecific taxa especially at the subfamial level.

Initial results on the molecular phylogeny of the Nudibranchia (Gastropoda, Opisthobranchia) based on 18S rDNA data

This study investigated nudibranch phylogeny on the basis of 18S rDNA sequence data. 18S rDNA sequence data of 19 taxa representing the major living orders and families of the Nudibranchia were analyzed. Representatives of the Cephalaspidea, Anaspidea, Gymnomorpha, Prosobranchia, and Pulmonata were also sequenced and used as outgroups. An additional 28 gastropod sequences taken from GenBank were also included in our analyses. Phylogenetic analyses of these more than 50 gastropod taxa provide strong evidence for support of the monophyly of the Nudibranchia. The monophyly of the Doridoidea, Cladobranchia, and Aeolidoidea within the Nudibranchia are also strongly supported. Phylogenetic utility and information content of the 18S rDNA sequences for Nudibranchia, and Opisthobranchia in general, are examined using the program SplitsTree as well as phylogenetic reconstructions using distance and parsimony approaches. 0Results based on these molecular data are compared with hypotheses about nudibranch phylogeny inferred from morphological data.

Molecular phylogeny of brachiopods and phoronids based on nuclear-encoded small subunit ribosomal RNA gene sequences

Brachiopod and phoronid phylogeny is inferred from SSU rDNA sequences of 28 articulate and nine inarticulate brachiopods, three phoronids, two ectoprocts and various outgroups, using gene trees reconstructed by weighted parsimony, distance and maximum likelihood methods. Of these sequences, 33 from brachiopods, two from phoronids and one each from an ectoproct and a priapulan are newly determined. The brachiopod sequences belong to 31 different genera and thus survey about 10% of extant genus-level diversity. Sequences determined in different laboratories and those from closely related taxa agree well, but evidence is presented suggesting that one published phoronid sequence (GenBank accession UO12648) is a brachiopod-phoronid chimaera, and this sequence is excluded from the analyses. The chiton, Acanthopleura , is identified as the phenetically proximal outgroup; other selected outgroups were chosen to allow comparison with recent, non-molecular analyses of brachiopod phylogeny. The different outgroups and methods of phylogenetic reconstruction lead to similar results, with differences mainly in the resolution of weakly supported ancient and recent nodes, including the divergence of inarticulate brachiopod sub-phyla, the position of the rhynchonellids in relation to long- and short-looped articulate brachiopod clades and the relationships of some articulate brachiopod genera and species. Attention is drawn to the problem presented by nodes that are strongly supported by non-molecular evidence but receive only low bootstrap resampling support. Overall, the gene trees agree with morphology-based brachiopod taxonomy, but novel relationships are tentatively suggested for thecideidine and megathyrid brachiopods. Articulate brachiopods are found to be monophyletic in all reconstructions, but monophyly of inarticulate brachiopods and the possible inclusion of phoronids in the inarticulate brachiopod clade are less strongly established. Phoronids are clearly excluded from a sister-group relationship with articulate brachiopods, this proposed relationship being due to the rejected, chimaeric sequence (GenBank UO12648). Lineage relative rate tests show no heterogeneity of evolutionary rate among articulate brachiopod sequences, but indicate that inarticulate brachiopod plus phoronid sequences evolve somewhat more slowly. Both brachiopods and phoronids evolve slowly by comparison with other invertebrates. A number of palaeontologically dated times of earliest appearance are used to make upper and lower estimates of the global rate of brachiopod SSU rDNA evolution, and these estimates are used to infer the likely divergence times of other nodes in the gene tree. There is reasonable agreement between most inferred molecular and palaeontological ages. The estimated rates of SSU rDNA sequence evolution suggest that the last common ancestor of brachiopods, chitons and other protostome invertebrates (Lophotrochozoa and Ecdysozoa) lived deep in Precambrian time. Results of this first DNA-based, taxonomically representative analysis of brachiopod phylogeny are in broad agreement with current morphology-based classification and systematics and are largely consistent with the hypothesis that brachiopod shell ontogeny and morphology are a good guide to phylogeny.

A priori estimation of phylogenetic information conserved in aligned sequences

A new phenomenological approach to explorative data analysis, the estimation of spectra of supporting positions, allows the search for conserved tracks left by phylogeny in DNA sequences. Spectra of supporting positions can be generated without reference to a tree topology or a model of sequence evolution and are therefore an ideal tool for a priori estimation of information content of data sets. Analysis of published 18S rDNA alignments shows that signal to noise relationship varies greatly in a way not detected by conventional tree-construction methods.