Phylogenetic support values are not necessarily informative: the case of the Serialia hypothesis (a mollusk phylogeny)

Neuronal segmentation in cephalopod arms

The prehensile arms of the cephalopod are among these animals' most remarkable features, but little is known about the neural circuitry governing arm and sucker movements1,2. Here, we investigated the cellular and molecular organization of the arm nervous system, focusing on the massive axial nerve cords (ANCs) in the octopus arms which collectively harbor four times as many neurons as the central brain3. We found that the ANC is segmented. In transverse cross sections, the ANC cell body layer wraps around the neuropil with no apparent segregation of sensory and motor neurons. In longitudinal sections, however, ANC neurons form segments, setting up a modular organization to the adjoining ANC neuropil. The septa between each segment are, in contrast, neuron-poor but contain nerve exits, vasculature and abundant collagen. Surprisingly, nerves exiting from neighboring septa differ in their fiber trajectories indicating that multiple adjoining segments must cooperate to innervate the arm musculature fully. The nerves for each sucker also exit through septa and set up a spatial "suckerotopy" in the ANC. A strong link between ANC segmentation and flexible sucker-laden arms was confirmed by comparative study of squid arms and tentacles. The ANC segmental modules represent a new template for understanding the motor control of octopus soft tissues. They also provide the first example of nervous system segmentation in a mollusc4.

Neuronal segmentation in cephalopod arms

The prehensile arms of the cephalopod are among these animals most remarkable features, but the neural circuitry governing arm and sucker movements remains largely unknown. We studied the neuronal organization of the adult axial nerve cord (ANC) of Octopus bimaculoides with molecular and cellular methods. The ANCs, which lie in the center of every arm, are the largest neuronal structures in the octopus, containing four times as many neurons as found in the central brain. In transverse cross section, the cell body layer (CBL) of the ANC wraps around its neuropil (NP) with little apparent segregation of sensory and motor neurons or nerve exits. Strikingly, when studied in longitudinal sections, the ANC is segmented. ANC neuronal cell bodies form columns separated by septa, with 15 segments overlying each pair of suckers. The segments underlie a modular organization to the ANC neuropil: neuronal cell bodies within each segment send the bulk of their processes directly into the adjoining neuropil, with some reaching the contralateral side. In addition, some nerve processes branch upon entering the NP, forming short-range projections to neighboring segments and mid-range projections to the ANC segments of adjoining suckers. The septa between the segments are employed as ANC nerve exits and as channels for ANC vasculature. Cellular analysis establishes that adjoining septa issue nerves with distinct fiber trajectories, which across two segments (or three septa) fully innervate the arm musculature. Sucker nerves also use the septa, setting up a nerve fiber "suckerotopy" in the sucker-side of the ANC. Comparative anatomy suggests a strong link between segmentation and flexible sucker-laden arms. In the squid Doryteuthis pealeii, the arms and the sucker-rich club of the tentacles have segments, but the sucker-poor stalk of the tentacles does not. The neural modules described here provide a new template for understanding the motor control of octopus soft tissues. In addition, this finding represents the first demonstration of nervous system segmentation in a mollusc.

Ancestral state reconstructions trace mitochondria but not phagocytosis to the last eukaryotic common ancestor

Two main theories have been put forward to explain the origin of mitochondria in eukaryotes: phagotrophic engulfment (undigested food) and microbial symbiosis (physiological interactions). The two theories generate mutually exclusive predictions about the order in which mitochondria and phagocytosis arose. To discriminate the alternatives, we have employed ancestral state reconstructions (ASR) for phagocytosis as a trait, phagotrophy as a feeding habit, the presence of mitochondria, the presence of plastids, and the multinucleated organization across major eukaryotic lineages. To mitigate the bias introduced by assuming a particular eukaryotic phylogeny, we reconstructed the appearance of these traits across 1789 different rooted gene trees, each having species from opisthokonts, mycetozoa, hacrobia, excavate, archeplastida, and Stramenopiles, Alveolates and Rhizaria. The trees reflect conflicting relationships and different positions of the root. We employed a novel phylogenomic test that summarizes ASR across trees which reconstructs a last eukaryotic common ancestor that possessed mitochondria, was multinucleated, lacked plastids, and was non-phagotrophic as well as non-phagocytic. This indicates that both phagocytosis and phagotrophy arose subsequent to the origin of mitochondria, consistent with findings from comparative physiology. Furthermore, our ASRs uncovered multiple origins of phagocytosis and of phagotrophy across eukaryotes, indicating that, like wings in animals, these traits are useful but neither ancestral nor homologous across groups. The data indicate that mitochondria preceded the origin of phagocytosis, such that phagocytosis cannot have been the mechanism by which mitochondria were acquired.

How many single-copy orthologous genes from whole genomes reveal deep gastropod relationships?

The Gastropoda contains 80% of existing mollusks and is the most diverse animal class second only to the Insecta. However, the deep phylogeny of gastropods has been controversial for a long time. Especially the position of Patellogastropoda is a major uncertainty. Morphology and some mitochondria studies concluded that Patellogastropoda is likely to be sister to all other gastropods (Orthogastropoda hypothesis), while transcriptomic and other mitogenomic studies indicated that Patellogastropoda and Vetigastropoda are sister taxa (Psilogastropoda). With the release of high-quality genomes, orthologous genes can be better identified and serve as powerful candidates for phylogenetic analysis. The question is, given the current limitations on the taxon sampling side, how many markers are needed to provide robust results. Here, we identified single-copy orthologous genes (SOGs) from 14 gastropods species with whole genomes available which cover five main gastropod subclasses. We generated different datasets from 395 to 1610 SOGs by allowing species missing in different levels. We constructed gene trees of each SOG, and inferred species trees from different collections of gene trees. We found as the number of SOGs increased, the inferred topology changed from Patellogastropoda being sister to all other gastropods to Patellogastropoda being sister to Vetigastropoda + Neomphalina (Psilogastropoda s.l.), with considerable support. Our study thus rejects the Orthogastropoda concept showing that the selection of the representative species and use of sufficient informative sites greatly influence the analysis of deep gastropod phylogeny.

Conflict over the eukaryote root resides in strong outliers, mosaics and missing data sensitivity of site-specific (CAT) mixture models

Phylogenetic reconstruction using concatenated loci ("phylogenomics" or "supermatrix phylogeny") is a powerful tool for solving evolutionary splits that are poorly resolved in single gene/protein trees (SGTs). However, recent phylogenomic attempts to resolve the eukaryote root have yielded conflicting results, along with claims of various artefacts hidden in the data. We have investigated these conflicts using two new methods for assessing phylogenetic conflict. ConJak uses whole marker (gene or protein) jackknifing to assess deviation from a central mean for each individual sequence, while ConWin uses a sliding window to screen for incongruent protein fragments (mosaics). Both methods allow selective masking of individual sequences or sequence fragments in order to minimize missing data, an important consideration for resolving deep splits with limited data. Analyses focused on a set of 76 eukaryotic proteins of bacterial-ancestry previously used in various combinations to assess the branching order among the three major divisions of eukaryotes: Amorphea (mainly animals, fungi and Amoebozoa), Diaphoretickes (most other well-known eukaryotes and nearly all algae) and Excavata, represented here by Discoba (Jakobida, Heterolobosea, and Euglenozoa). ConJak analyses found strong outliers to be concentrated in under-sampled lineages, while ConWin analyses of Discoba, the most under-sampled of the major lineages, detected potentially incongruent fragments scattered throughout. Phylogenetic analyses of the full data using an LG-gamma model support a Discoba sister scenario (neozoan-excavate root), which rises to 99-100% bootstrap support with data masked according to either protocol. However, analyses with two site-specific (CAT) mixture models yielded widely inconsistent results and a striking sensitivity to missing data. The neozoan-excavate root places Amorphea and Diaphoretickes as more closely related to each other than either is to Discoba, a fundamental relationship that should remain unaffected by additional taxa.

Transcriptomics provides a robust framework for the relationships of the major clades of cladobranch sea slugs (Mollusca, Gastropoda, Heterobranchia), but fails to resolve the position of the enigmatic genus Embletonia

Background

The soft-bodied cladobranch sea slugs represent roughly half of the biodiversity of marine nudibranch molluscs on the planet. Despite their global distribution from shallow waters to the deep sea, from tropical into polar seas, and their important role in marine ecosystems and for humans (as targets for drug discovery), the evolutionary history of cladobranch sea slugs is not yet fully understood.

Results

To enlarge the current knowledge on the phylogenetic relationships, we generated new transcriptome data for 19 species of cladobranch sea slugs and two additional outgroup taxa (Berthella plumula and Polycera quadrilineata). We complemented our taxon sampling with previously published transcriptome data, resulting in a final data set covering 56 species from all but one accepted cladobranch superfamilies. We assembled all transcriptomes using six different assemblers, selecting those assemblies that provided the largest amount of potentially phylogenetically informative sites. Quality-driven compilation of data sets resulted in four different supermatrices: two with full coverage of genes per species (446 and 335 single-copy protein-coding genes, respectively) and two with a less stringent coverage (667 genes with 98.9% partition coverage and 1767 genes with 86% partition coverage, respectively). We used these supermatrices to infer statistically robust maximum-likelihood trees. All analyses, irrespective of the data set, indicate maximal statistical support for all major splits and phylogenetic relationships at the family level. Besides the questionable position of Noumeaella rubrofasciata, rendering the Facelinidae as polyphyletic, the only notable discordance between the inferred trees is the position of Embletonia pulchra. Extensive testing using Four-cluster Likelihood Mapping, Approximately Unbiased tests, and Quartet Scores revealed that its position is not due to any informative phylogenetic signal, but caused by confounding signal.

Conclusions

Our data matrices and the inferred trees can serve as a solid foundation for future work on the taxonomy and evolutionary history of Cladobranchia. The placement of E. pulchra, however, proves challenging, even with large data sets and various optimization strategies. Moreover, quartet mapping results show that confounding signal present in the data is sufficient to explain the inferred position of E. pulchra, again leaving its phylogenetic position as an enigma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01944-0.

Comprehensive phenotypic phylogenetic analysis supports the monophyly of stromateiform fishes (Teleostei: Percomorphacea)

More than half the ray-finned fishes and about one-quarter of all living vertebrates belong to Percomorphacea. Among its 30 orders, Stromateiformes encompass 77 species in 16 genera and six families. Stromateiform monophyly has never been tested using morphology, and it has been rejected by molecular analyses. This comprehensive revision of Stromateiformes includes all its valid genera of all percomorph families previously aligned with the order. We sampled 207 phenotypic characters in 66 terminal taxa representing 14 orders and 46 acanthopterygian families. This dataset significantly surpasses all previous phenotype-based phylogenies of Stromateiformes, which analysed only a fraction of these characters. Stromateiformes is recovered as monophyletic, supported by eight unequivocal synapomorphies. Amarsipidae is the sister group of all other Stromateiformes (= Stromateoidei). Centrolophidae is paraphyletic, with three of its genera allocated into an early-diverging clade and the other four appearing as successive sister groups to a lineage containing the remaining stromateiforms. All other stromateoid families are monophyletic, with the following cladistic arrangement: (Nomeidae (Stromateidae (Tetragonuridae, Ariommatidae))). Our analysis convincingly refutes recent molecular phylogenetic interpretations that fail to recover a monophyletic Stromateiformes. These findings call into question large-scale conclusions of percomorph relationships and trait evolution based solely on molecular data.

Geographic isolation and human-assisted dispersal in land snails: a Mediterranean story of Helix borealis and its relatives (Gastropoda: Stylommatophora: Helicidae)

The Mediterranean basin is a major centre for land-snail diversity, with many localized endemics, but there are also species widely spread by humans. Both endemics and introductions can be found in the snail genus Helix, which comprises many large-bodied species used for human consumption in the past and present. The Mediterranean clade of Helix is currently distributed throughout this region, but the phylogenetic and biogeographic relationships among its forms from different parts of the basin remain enigmatic. The reasons include insufficient sampling, taxa with unclear taxonomy and a significant impact of human-assisted transport obscuring the natural distribution of phylogenetic lineages. We provide evidence that European and Anatolian populations of H. cincta and its relatives are not native to those regions, but originate from the northern Levant. These results have implications for taxonomy of the genus, but also for the understanding of its evolutionary history. We posit that the Mediterranean clade consists of four geographically separated groups, which diversified in Northern Africa, the Apennine Peninsula and Corsica, the Aegean and Greece, and the northern Levant. This geographic pattern has been subsequently blurred by multiple instances of human-assisted dispersal. However, revealing the founding populations with certainty requires thorough sampling in currently inaccessible countries.

An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola)

BACKGROUND

The latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations.

RESULTS

Our results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera.

CONCLUSION

Our extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Topological support and data quality can only be assessed through multiple tests in reviewing Blattodea phylogeny

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness.

Evaluation of Systematic Position of Helicoprorodontids and Chaeneids (Ciliophora, Litostomatea): An Attempt to Break Long Branches in 18S rRNA Gene Phylogenies

Phylogenetic position of some free-living litostomatean taxa has not been correctly determined because of long-branch artifacts in 18S rRNA gene trees. The main aim of this study was to test the effectiveness of various masking algorithms, tree-building techniques, binarization of DNA data as well as combining morphological and molecular data to eliminate long-branch attraction of two problematic groups, helicoprorodontids and chaeneids. Guidance and SlowFaster masking in a combination with PhyloBayesian tree construction erased the artifactual positions of helicoprorodontids and chaeneids. On the other hand, binarization of DNA sequences and the strategy of combining morphological and molecular data eliminated only the artifactual position of chaeneids but not that of helicoprorodontids which were still being attracted by out-group taxa. According to statistical tree topology tests and comparative morphological studies, helicoprorodontids are classified as a distinct order while chaeneids are considered to be fast evolving members of the order Lacrymariida. The high body contractility, "cephalization" of the anterior body end, and helicalization of the anterior portion of some or all somatic ciliary rows indicate relatedness of helicoprorodontids, chaeneids, and lacrymariids. On the other hand, the dorsal brush separated from the circumoral kinety by dense ciliary files supports kinships of chaeneids, lacrymariids, and didiniids.

Consensus and confusion in molluscan trees: evaluating morphological and molecular phylogenies

Mollusks are the most morphologically disparate living animal phylum, they have diversified into all habitats, and have a deep fossil record. Monophyly and identity of their eight living classes is undisputed, but relationships between these groups and patterns of their early radiation have remained elusive. Arguments about traditional morphological phylogeny focus on a small number of topological concepts but often without regard to proximity of the individual classes. In contrast, molecular studies have proposed a number of radically different, inherently contradictory, and controversial sister relationships. Here, we assembled a data set of 42 unique published trees describing molluscan interrelationships. We used these data to ask several questions about the state of resolution of molluscan phylogeny compared with a null model of the variation possible in random trees constructed from a monophyletic assemblage of eight terminals. Although 27 different unique trees have been proposed from morphological inference, the majority of these are not statistically different from each other. Within the available molecular topologies, only four studies to date have included the deep sea class Monoplacophora; but 36.4% of all trees are not significantly different. We also present supertrees derived from two data partitions and three methods, including all available molecular molluscan phylogenies, which will form the basis for future hypothesis testing. The supertrees presented here were not constructed to provide yet another hypothesis of molluscan relationships, but rather to algorithmically evaluate the relationships present in the disparate published topologies. Based on the totality of available evidence, certain patterns of relatedness among constituent taxa become clear. The internodal distance is consistently short between a few taxon pairs, particularly supporting the relatedness of Monoplacophora and the chitons, Polyplacophora. Other taxon pairs are rarely or never found in close proximity, such as the vermiform Caudofoveata and Bivalvia. Our results have specific utility for guiding constructive research planning to better test relationships in Mollusca as well as other problematic groups. Taxa with consistently proximate relationships should be the focus of a combined approach in a concerted assessment of potential genetic and anatomical homology, whereas unequivocally distant taxa will make the most constructive choices for exemplar selection in higher level phylogenomic analyses.

The complete mitochondrial genome of Scutopus ventrolineatus (Mollusca: Chaetodermomorpha) supports the Aculifera hypothesis

Background

With more than 100000 living species, mollusks are the second most diverse metazoan phylum. The current taxonomic classification of mollusks recognizes eight classes (Neomeniomorpha, Chaetodermomorpha, Polyplacophora, Monoplacophora, Cephalopoda, Gastropoda, Bivalvia, and Scaphopoda) that exhibit very distinct body plans. In the past, phylogenetic relationships among mollusk classes have been contentious due to the lack of indisputable morphological synapomorphies. Fortunately, recent phylogenetic analyses based on multi-gene data sets are rendering promising results. In this regard, mitochondrial genomes have been widely used to reconstruct deep phylogenies. For mollusks, complete mitochondrial genomes are mostly available for gastropods, bivalves, and cephalopods, whereas other less-diverse lineages have few or none reported.

Results

The complete DNA sequence (14662 bp) of the mitochondrial genome of the chaetodermomorph Scutopus ventrolineatus Salvini-Plawen, 1968 was determined. Compared with other mollusks, the relative position of protein-coding genes in the mitochondrial genome of S. ventrolineatus is very similar to those reported for Polyplacophora, Cephalopoda and early-diverging lineages of Bivalvia and Gastropoda (Vetigastropoda and Neritimorpha; but not Patellogastropoda). The reconstructed phylogenetic tree based on combined mitochondrial and nuclear sequence data recovered monophyletic Aplacophora, Aculifera, and Conchifera. Within the latter, Cephalopoda was the sister group of Gastropoda and Bivalvia + Scaphopoda.

Conclusions

Phylogenetic analyses of mitochondrial sequences showed strong among-lineage rate heterogeneity that produced long-branch attraction biases. Removal of long branches (namely those of bivalves and patellogastropods) ameliorated but not fully resolved the problem. Best results in terms of statistical support were achieved when mitochondrial and nuclear sequence data were concatenated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0197-9) contains supplementary material, which is available to authorized users.

A protocol for species delineation of public DNA databases, applied to the Insecta

Public DNA databases are composed of data from many different taxa, although the taxonomic annotation on sequences is not always complete, which impedes the utilization of mined data for species-level applications. There is much ongoing work on species identification and delineation based on the molecular data itself, although applying species clustering to whole databases requires consolidation of results from numerous undefined gene regions, and introduces significant obstacles in data organization and computational load. In the current paper, we demonstrate an approach for species delineation of a sequence database. All DNA sequences for the insects were obtained and processed. After filtration of duplicated data, delineation of the database into species or molecular operational taxonomic units (MOTUs) followed a three-step process in which (i) the genetic loci L are partitioned, (ii) the species S are delineated within each locus, then (iii) species units are matched across loci to form the matrix L × S, a set of global (multilocus) species units. Partitioning the database into a set of homologous gene fragments was achieved by Markov clustering using edge weights calculated from the amount of overlap between pairs of sequences, then delineation of species units and assignment of species names were performed for the set of genes necessary to capture most of the species diversity. The complexity of computing pairwise similarities for species clustering was substantial at the cytochrome oxidase subunit I locus in particular, but made feasible through the development of software that performs pairwise alignments within the taxonomic framework, while accounting for the different ranks at which sequences are labeled with taxonomic information. Over 24 different homologs, the unidentified sequences numbered approximately 194,000, containing 41,525 species IDs (98.7% of all found in the insect database), and were grouped into 59,173 single-locus MOTUs by hierarchical clustering under parameters optimized independently for each locus. Species units from different loci were matched using a multipartite matching algorithm to form multilocus species units with minimal incongruence between loci. After matching, the insect database as represented by these 24 loci was found to be composed of 78,091 species units in total. 38,574 of these units contained only species labeled data, 34,891 contained only unlabeled data, leaving 4,626 units composed both of labeled and unlabeled sequences. In addition to giving estimates of species diversity of sequence repositories, the protocol developed here will facilitate species-level applications of modern-day sequence data sets. In particular, the L × S matrix represents a post-taxonomic framework that can be used for species-level organization of metagenomic data, and incorporation of these methods into phylogenetic pipelines will yield matrices more representative of species diversity.

The continuing debate on deep molluscan phylogeny: evidence for Serialia (Mollusca, Monoplacophora + Polyplacophora)

Molluscs are a diverse animal phylum with a formidable fossil record. Although there is little doubt about the monophyly of the eight extant classes, relationships between these groups are controversial. We analysed a comprehensive multilocus molecular data set for molluscs, the first to include multiple species from all classes, including five monoplacophorans in both extant families. Our analyses of five markers resolve two major clades: the first includes gastropods and bivalves sister to Serialia (monoplacophorans and chitons), and the second comprises scaphopods sister to aplacophorans and cephalopods. Traditional groupings such as Testaria, Aculifera, and Conchifera are rejected by our data with significant Approximately Unbiased (AU) test values. A new molecular clock indicates that molluscs had a terminal Precambrian origin with rapid divergence of all eight extant classes in the Cambrian. The recovery of Serialia as a derived, Late Cambrian clade is potentially in line with the stratigraphic chronology of morphologically heterogeneous early mollusc fossils. Serialia is in conflict with traditional molluscan classifications and recent phylogenomic data. Yet our hypothesis, as others from molecular data, implies frequent molluscan shell and body transformations by heterochronic shifts in development and multiple convergent adaptations, leading to the variable shells and body plans in extant lineages.

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.

Dollo's law and the irreversibility of digit loss in Bachia

Several recent studies conclude that exceptions to Dollo's law are more common than used to be thought. If the claims are true this would change our view on the role of developmental constraints in the evolution of body plans. One study claims the reevolution of lost digits in the lizard genus Bachia (Kohlsdorf and Wagner 2006). We evaluate this claim. We conclude that the proposed molecular phylogenetic tree is in conflict with evolutionary mechanisms concerning the biogeography of lizards and with morphology-based phylogenies. A reanalysis of the molecular data does not support the topology of the published tree. Furthermore, two implicit assumptions, that digit numbers are fixed and that polydactyly evolves independently from other characters, are incorrect. We conclude that there is no convincing support for reevolution of digits in Bachia. We discuss our findings in the light of the current evidence for the reversal of losses of complex traits. We conclude that in metazoans, exceptions to Dollo's law are mainly found among meristic traits that originate relatively late during embryogenesis, when developmental systems are more compartmentalized. Finally, our study shows that phylogenetic analyses should incorporate evolutionary mechanisms including constraints, variation, and selection, not only for correct phylogenetic reconstruction, but also for correct evolutionary inference.

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.