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

Genetic evidence against monophyly of Oniscidea implies a need to revise scenarios for the origin of terrestrial isopods

Among the few crustacean taxa that managed to inhabit terrestrial environments, Oniscidea includes the most successful colonizers in terms of species richness and abundance. However, neither morphological traits nor molecular markers have definitively resolved phylogenetic relationships among major Oniscidea clades or established the monophyly of the taxon. Herein, we employed the highly conserved, nuclear protein-coding genes Sodium-Potassium Pump (NAK) and Phosphoenolpyruvate Carboxykinase (PEPCK), along with the traditionally used 18 s and 28 s ribosomal RNA genes, in an attempt to clarify these questions. Our dataset included sequences representing all major Oniscidea clades and closely related aquatic taxa, as suggested by previous studies. We applied Bayesian Inference and Maximum Likelihood methods and produced a robust and fully resolved phylogenetic tree that offers strong evidence against the monophyly of Oniscidea. The amphibious genus Ligia appears to be more closely related to representatives of marine suborders, while the phylogenetic pattern of the remaining Oniscidea implies a complex history of the transition from the marine environment to land. With the exception of the basal clade, all other established major clades have been recovered as monophyletic, even though relationships within these clades call for a revised interpretation of morphological characters used in terrestrial isopod taxonomy.

New ways to acquire resistance: imperfect convergence in insect adaptations to a potent plant toxin

Evolution of insensitivity to the toxic effects of cardiac glycosides has become a model in the study of convergent evolution, as five taxonomic orders of insects use the same few similar amino acid substitutions in the otherwise highly conserved Na,K-ATPase α. We show here that insensitivity in pyrgomorphid grasshoppers evolved along a slightly divergent path. As in other lineages, duplication of the Na,K-ATPase α gene paved the way for subfunctionalization: one copy maintains the ancestral, sensitive state, while the other copy is resistant. Nonetheless, in contrast with all other investigated insects, the grasshoppers' resistant copy shows length variation by two amino acids in the first extracellular loop, the main part of the cardiac glycoside-binding pocket. RT-qPCR analyses confirmed that this copy is predominantly expressed in tissues exposed to the toxins, while the ancestral copy predominates in the nervous tissue. Functional tests with genetically engineered Drosophila Na,K-ATPases bearing the first extracellular loop of the pyrgomorphid genes showed the derived form to be highly resistant, while the ancestral state is sensitive. Thus, we report convergence in gene duplication and in the gene targets for toxin insensitivity; however, the means to the phenotypic end have been novel in pyrgomorphid grasshoppers.

Phylogenomics of Lophotrochozoa with Consideration of Systematic Error

Phylogenomic studies have improved understanding of deep metazoan phylogeny and show promise for resolving incongruences among analyses based on limited numbers of loci. One region of the animal tree that has been especially difficult to resolve, even with phylogenomic approaches, is relationships within Lophotrochozoa (the animal clade that includes molluscs, annelids, and flatworms among others). Lack of resolution in phylogenomic analyses could be due to insufficient phylogenetic signal, limitations in taxon and/or gene sampling, or systematic error. Here, we investigated why lophotrochozoan phylogeny has been such a difficult question to answer by identifying and reducing sources of systematic error. We supplemented existing data with 32 new transcriptomes spanning the diversity of Lophotrochozoa and constructed a new set of Lophotrochozoa-specific core orthologs. Of these, 638 orthologous groups (OGs) passed strict screening for paralogy using a tree-based approach. In order to reduce possible sources of systematic error, we calculated branch-length heterogeneity, evolutionary rate, percent missing data, compositional bias, and saturation for each OG and analyzed increasingly stricter subsets of only the most stringent (best) OGs for these five variables. Principal component analysis of the values for each factor examined for each OG revealed that compositional heterogeneity and average patristic distance contributed most to the variance observed along the first principal component while branch-length heterogeneity and, to a lesser extent, saturation contributed most to the variance observed along the second. Missing data did not strongly contribute to either. Additional sensitivity analyses examined effects of removing taxa with heterogeneous branch lengths, large amounts of missing data, and compositional heterogeneity. Although our analyses do not unambiguously resolve lophotrochozoan phylogeny, we advance the field by reducing the list of viable hypotheses. Moreover, our systematic approach for dissection of phylogenomic data can be applied to explore sources of incongruence and poor support in any phylogenomic data set. [Annelida; Brachiopoda; Bryozoa; Entoprocta; Mollusca; Nemertea; Phoronida; Platyzoa; Polyzoa; Spiralia; Trochozoa.].

RNA sequencing analysis of transcriptional change in the freshwater mussel Elliptio complanata after environmentally relevant sodium chloride exposure

To identify potential biomarkers of salt stress in a freshwater sentinel species, we examined transcriptional responses of the common mussel Elliptio complanata to controlled sodium chloride (NaCl) exposures. Ribonucleic acid sequencing (RNA-Seq) of mantle tissue identified 481 transcripts differentially expressed in adult mussels exposed to 2 ppt NaCl (1.2 ppt chloride) for 7 d, of which 290 had nonoverlapping intervals. Differentially expressed gene categories included ion and transmembrane transport, oxidoreductase activity, maintenance of protein folding, and amino acid metabolism. The rate-limiting enzyme for synthesis of taurine, an amino acid frequently linked to osmotic stress in aquatic species, was upregulated, as was the transmembrane ion pump sodium/potassium adenosine 5'-triphosphatase. These patterns confirm a primary transcriptional response to the experimental dose, albeit likely overlapping with nonspecific secondary stress responses. Substantial involvement of the heat shock protein 70 chaperone family and the water-transporting aquaporin family was not detected, however, in contrast to some studies in other bivalves. A subset of the most significantly regulated genes was confirmed by quantitative polymerase chain reaction in an independent sample. Cluster analysis showed separation of mussels exposed to 2 ppt NaCl from control mussels in multivariate space, but mussels exposed to 1 ppt NaCl were largely indistinguishable from controls. Transcriptome-scale analysis of salt exposure under laboratory conditions efficiently identified candidate biomarkers for further functional analysis and field validation. Environ Toxicol Chem 2017;36:2352-2366. © Published 2017 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Structural analysis of the α subunit of Na(+)/K(+) ATPase genes in invertebrates

The Na(+)/K(+) ATPase is a ubiquitous pump coordinating the transport of Na(+) and K(+) across the membrane of cells and its role is fundamental to cellular functions. It is heteromer in eukaryotes including two or three subunits (α, β and γ which is specific to the vertebrates). The catalytic functions of the enzyme have been attributed to the α subunit. Several complete α protein sequences are available, but only few gene structures were characterized. We identified the genomic sequences coding the α-subunit of the Na(+)/K(+) ATPase, from the whole-genome shotgun contigs (WGS), NCBI Genomes (chromosome), Genomic Survey Sequences (GSS) and High Throughput Genomic Sequences (HTGS) databases across distinct phyla. One copy of the α subunit gene was found in Annelida, Arthropoda, Cnidaria, Echinodermata, Hemichordata, Mollusca, Placozoa, Porifera, Platyhelminthes, Urochordata, but the nematodes seem to possess 2 to 4 copies. The number of introns varied from 0 (Platyhelminthes) to 26 (Porifera); and their localization and length are also highly variable. Molecular phylogenies (Maximum Likelihood and Maximum Parsimony methods) showed some clusters constituted by (Chordata/(Echinodermata/Hemichordata)) or (Plathelminthes/(Annelida/Mollusca)) and a basal position for Porifera. These structural analyses increase our knowledge about the evolutionary events of the α subunit genes in the invertebrates.

Sternites and spiracles - the unclear homology of ventral sclerites in the basal millipede order Glomeridesmida (Myriapoda, Diplopoda)

We report the discovery of a ventral plate in the basal and little-known chilognath millipede order Glomeridesmida. This ventral plate, interpreted here as a 'true sternite', is clearly separate from both the coxa and the more lateral stigma-carrying plates commonly referred to as 'diplopod sternites'. Therefore, the lateral, stigma-carrying plates of the Diplopoda, previously referred to as sternites, are not sternal elements, but subcoxal elements associated with the limb base. This discovery changes the nomenclature used for the ventral plates in Diplopoda, with the formerly named 'sternite' better referred to as 'stigma-carrying plate'. In helminthomorph Diplopoda, the stigma-carrying plates are apparently secondarily fused with the sternite. The main argument for the independent evolution of tracheae in insects and myriapods, the different location of their respective spiracles, no longer holds true. In all Myriapoda and Hexapoda the spiracles associated with subcoxal elements are located lateral to the limb base. This discovery shows that the arguments for an independent origin of tracheae in insects and myriapods are not uncontestable.

New phylogenomic data support the monophyly of Lophophorata and an Ectoproct-Phoronid clade and indicate that Polyzoa and Kryptrochozoa are caused by systematic bias

BACKGROUND

Within the complex metazoan phylogeny, the relationships of the three lophophorate lineages, ectoprocts, brachiopods and phoronids, are particularly elusive. To shed further light on this issue, we present phylogenomic analyses of 196 genes from 58 bilaterian taxa, paying particular attention to the influence of compositional heterogeneity.

RESULTS

The phylogenetic analyses strongly support the monophyly of Lophophorata and a sister-group relationship between Ectoprocta and Phoronida. Our results contrast previous findings based on rDNA sequences and phylogenomic datasets which supported monophyletic Polyzoa (= Bryozoa sensu lato) including Ectoprocta, Entoprocta and Cycliophora, Brachiozoa including Brachiopoda and Phoronida as well as Kryptrochozoa including Brachiopoda, Phoronida and Nemertea, thus rendering Lophophorata polyphyletic. Our attempts to identify the causes for the conflicting results revealed that Polyzoa, Brachiozoa and Kryptrochozoa are supported by character subsets with deviating amino acid compositions, whereas there is no indication for compositional heterogeneity in the character subsets supporting the monophyly of Lophophorata.

CONCLUSION

Our results indicate that the support for Polyzoa, Brachiozoa and Kryptrochozoa gathered so far is likely an artifact caused by compositional bias. The monophyly of Lophophorata implies that the horseshoe-shaped mesosomal lophophore, the tentacular feeding apparatus of ectoprocts, phoronids and brachiopods is, indeed, a synapomorphy of the lophophorate lineages. The same may apply to radial cleavage. However, among phoronids also spiral cleavage is known. This suggests that the cleavage pattern is highly plastic and has changed several times within lophophorates. The sister group relationship of ectoprocts and phoronids is in accordance with the interpretation of the eversion of a ventral invagination at the beginning of metamorphosis as a common derived feature of these taxa.

Molecular systematics of the deep-sea hydrothermal vent endemic Brachyuran family Bythograeidae: a comparison of three Bayesian species tree methods

Brachyuran crabs of the family Bythograeidae are endemic to deep-sea hydrothermal vents and represent one of the most successful groups of macroinvertebrates that have colonized this extreme environment. Occurring worldwide, the family includes six genera (Allograea, Austinograea, Bythograea, Cyanagraea, Gandalfus, and Segonzacia) and fourteen formally described species. To investigate their evolutionary relationships, we conducted Maximum Likelihood and Bayesian molecular phylogenetic analyses, based on DNA sequences from fragments of three mitochondrial genes (16S rDNA, Cytochrome oxidase I, and Cytochrome b) and three nuclear genes (28S rDNA, the sodium-potassium ATPase a-subunit 'NaK', and Histone H3A). We employed traditional concatenated (i.e., supermatrix) phylogenetic methods, as well as three recently developed Bayesian multilocus methods aimed at inferring species trees from potentially discordant gene trees. We found strong support for two main clades within Bythograeidae: one comprising the members of the genus Bythograea; and the other comprising the remaining genera. Relationships within each of these two clades were partially resolved. We compare our results with an earlier hypothesis on the phylogenetic relationships among bythograeid genera based on morphology. We also discuss the biogeography of the family in the light of our results. Our species tree analyses reveal differences in how each of the three methods weighs conflicting phylogenetic signal from different gene partitions and how limits on the number of outgroup taxa may affect the results.

The complete mitochondrial genome of Flustra foliacea (Ectoprocta, Cheilostomata) - compositional bias affects phylogenetic analyses of lophotrochozoan relationships

Background

The phylogenetic relationships of the lophophorate lineages, ectoprocts, brachiopods and phoronids, within Lophotrochozoa are still controversial. We sequenced an additional mitochondrial genome of the most species-rich lophophorate lineage, the ectoprocts. Although it is known that there are large differences in the nucleotide composition of mitochondrial sequences of different lineages as well as in the amino acid composition of the encoded proteins, this bias is often not considered in phylogenetic analyses. We applied several approaches for reducing compositional bias and saturation in the phylogenetic analyses of the mitochondrial sequences.

Results

The complete mitochondrial genome (16,089 bp) of Flustra foliacea (Ectoprocta, Gymnolaemata, Cheilostomata) was sequenced. All protein-encoding, rRNA and tRNA genes are transcribed from the same strand. Flustra shares long intergenic sequences with the cheilostomate ectoproct Bugula, which might be a synapomorphy of these taxa. Further synapomorphies might be the loss of the DHU arm of the tRNA L(UUR), the loss of the DHU arm of the tRNA S(UCN) and the unique anticodon sequence GAG of the tRNA L(CUN). The gene order of the mitochondrial genome of Flustra differs strongly from that of the other known ectoprocts. Phylogenetic analyses of mitochondrial nucleotide and amino acid data sets show that the lophophorate lineages are more closely related to trochozoan phyla than to deuterostomes or ecdysozoans confirming the Lophotrochozoa hypothesis. Furthermore, they support the monophyly of Cheilostomata and Ectoprocta. However, the relationships of the lophophorate lineages within Lophotrochozoa differ strongly depending on the data set and the used method. Different approaches for reducing heterogeneity in nucleotide and amino acid data sets and saturation did not result in a more robust resolution of lophotrochozoan relationships.

Conclusion

The contradictory and usually weakly supported phylogenetic reconstructions of the relationships among lophotrochozoan phyla based on mitochondrial sequences indicate that these alone do not contain enough information for a robust resolution of the relations of the lophotrochozoan phyla. The mitochondrial gene order is also not useful for inferring their phylogenetic relationships, because it is highly variable in ectoprocts, brachiopods and some other lophotrochozoan phyla. However, our study revealed several rare genomic changes like the evolution of long intergenic sequences and changes in the structure of tRNAs, which may be helpful for reconstructing ectoproct phylogeny.

FPGA Acceleration of the phylogenetic likelihood function for Bayesian MCMC inference methods

Background

Likelihood (ML)-based phylogenetic inference has become a popular method for estimating the evolutionary relationships among species based on genomic sequence data. This method is used in applications such as RAxML, GARLI, MrBayes, PAML, and PAUP. The Phylogenetic Likelihood Function (PLF) is an important kernel computation for this method. The PLF consists of a loop with no conditional behavior or dependencies between iterations. As such it contains a high potential for exploiting parallelism using micro-architectural techniques. In this paper, we describe a technique for mapping the PLF and supporting logic onto a Field Programmable Gate Array (FPGA)-based co-processor. By leveraging the FPGA's on-chip DSP modules and the high-bandwidth local memory attached to the FPGA, the resultant co-processor can accelerate ML-based methods and outperform state-of-the-art multi-core processors.

Results

We use the MrBayes 3 tool as a framework for designing our co-processor. For large datasets, we estimate that our accelerated MrBayes, if run on a current-generation FPGA, achieves a 10× speedup relative to software running on a state-of-the-art server-class microprocessor. The FPGA-based implementation achieves its performance by deeply pipelining the likelihood computations, performing multiple floating-point operations in parallel, and through a natural log approximation that is chosen specifically to leverage a deeply pipelined custom architecture.

Conclusions

Heterogeneous computing, which combines general-purpose processors with special-purpose co-processors such as FPGAs and GPUs, is a promising approach for high-performance phylogeny inference as shown by the growing body of literature in this field. FPGAs in particular are well-suited for this task because of their low power consumption as compared to many-core processors and Graphics Processor Units (GPUs) [1].

Tentaculate fossils from the Cambrian of Canada (British Columbia) and China (Yunnan) interpreted as primitive deuterostomes

Molecular and morphological evidence unite the hemichordates and echinoderms as the Ambulacraria, but their earliest history remains almost entirely conjectural. This is on account of the morphological disparity of the ambulacrarians and a paucity of obvious stem-groups. We describe here a new taxon Herpetogaster collinsi gen. et sp. nov. from the Burgess Shale (Middle Cambrian) Lagerstätte. This soft-bodied vermiform animal has a pair of elongate dendritic oral tentacles, a flexible stolon with an attachment disc, and a re-curved trunk with at least 13 segments that is directed dextrally. A differentiated but un-looped gut is enclosed in a sac suspended by mesenteries. It consists of a short pharynx, a conspicuous lenticular stomach, followed by a narrow intestine sub-equal in length. This new taxon, together with the Lower Cambrian Phlogites and more intriguingly the hitherto enigmatic discoidal eldoniids (Cambrian-Devonian), form a distinctive clade (herein the cambroernids). Although one hypothesis of their relationships would look to the lophotrochozoans (specifically the entoprocts), we suggest that the evidence is more consistent with their being primitive deuterostomes, with specific comparisons being made to the pterobranch hemichordates and pre-radial echinoderms. On this basis some of the earliest ambulacrarians are interpreted as soft-bodied animals with a muscular stalk, and possessing prominent tentacles.

Bilaterian phylogeny: a broad sampling of 13 nuclear genes provides a new Lophotrochozoa phylogeny and supports a paraphyletic basal acoelomorpha

During the past decade, great progress has been made in clarifying the relationships among bilaterian animals. Studies based on a limited number of markers established new hypotheses such as the existence of three superclades (Deuterostomia, Ecdysozoa, and Lophotrochozoa) but left major questions unresolved. The data sets used to the present either bear few characters for many taxa (i.e., the ribosomal genes) or present many characters but lack many phyla (such as recent phylogenomic approaches) failing to provide definitive answers for all the regions of the bilaterian tree. We performed phylogenetic analyses using a molecular matrix with a high number of characters and bilaterian phyla. This data set is built from 13 genes (8,880 bp) belonging to 90 taxa from 27 bilaterian phyla. Probabilistic analyses robustly support the three superclades, the monophyly of Chordata, a spiralian clade including Brachiozoa, the basal position of a paraphyletic Acoelomorpha, and point to an ecdysozoan affiliation for Chaetognatha. This new phylogeny not only agrees with most classical molecular results but also provides new insights into the relationships between lophotrochozoans and challenges the results obtained using high-throughput strategies, highlighting the problems associated with the current trend to increase gene number rather than taxa.

Complete mitochondrial genome of Bugula neritina (Bryozoa, Gymnolaemata, Cheilostomata): phylogenetic position of Bryozoa and phylogeny of lophophorates within the Lophotrochozoa

BACKGROUND

The phylogenetic position of Bryozoa is one of the most controversial issues in metazoan phylogeny. In an attempt to address this issue, the first bryozoan mitochondrial genome from Flustrellidra hispida (Gymnolaemata, Ctenostomata) was recently sequenced and characterized. Unfortunately, it has extensive gene translocation and extremely reduced size. In addition, the phylogenies obtained from the result were conflicting, so they failed to assign a reliable phylogenetic position to Bryozoa or to clarify lophophorate phylogeny. Thus, it is necessary to characterize further mitochondrial genomes from slowly-evolving bryozoans to obtain a more credible lophophorate phylogeny.

RESULTS

The complete mitochondrial genome (15,433 bp) of Bugula neritina (Bryozoa, Gymnolaemata, Cheilostomata), one of the most widely distributed cheliostome bryozoans, is sequenced. This second bryozoan mitochondrial genome contains the set of 37 components generally observed in other metazoans, differing from that of F. hispida (Bryozoa, Gymnolaemata, Ctenostomata), which has only 36 components with loss of tRNAser(ucn) genes. The B. neritina mitochondrial genome possesses 27 multiple noncoding regions. The gene order is more similar to those of the two remaining lophophorate phyla (Brachiopoda and Phoronida) and a chiton Katharina tunicate than to that of F. hispida. Phylogenetic analyses based on the nucleotide sequences or amino acid residues of 12 protein-coding genes showed consistently that, within the Lophotrochozoa, the monophyly of the bryozoan class Gymnolaemata (B. neritina and F. hispida) was strongly supported and the bryozoan clade was grouped with brachiopods. Echiura appeared as a subtaxon of Annelida, and Entoprocta as a sister taxon of Phoronida. The clade of Bryozoa + Brachiopoda was clustered with either the clade of Annelida-Echiura or that of Phoronida + Entoprocta.

CONCLUSION

This study presents the complete mitochondrial genome of a cheliostome bryozoan, B. neritina. The phylogenetic analyses suggest a close relationship between Bryozoa and Brachiopoda within the Lophotrochozoa. However, the sister group of Bryozoa + Brachiopoda is still ambiguous, although it has some attractions with Annelida-Echiura or Phoronida + Entoprocta. If the latter is a true phylogeny, lophophorate monophyly including Entoprocta is supported. Consequently, the present results imply that Brachiozoa (= Brachiopoda + Phoronida) and the recently-resurrected Bryozoa concept comprising Ectoprocta and Entoprocta may be refuted.

Evolutionary history of Na,K-ATPases and their osmoregulatory role

The Na/K pump, or Na,K-ATPase, is a key enzyme to the homeostasis of osmotic pressure, cell volume, and the maintenance of electrochemical gradients. Its alpha subunit, which holds most of its functions, belongs to a large family of ATPases known as P-type, and to the subfamily IIC, which also includes H,K-ATPases. In this study, we attempt to describe the evolutionary history of IIC ATPases by doing phylogenetic analysis with most of the currently available protein sequences (over 200), and pay special attention to the relationship between their diversity and their osmoregulatory role. We include proteins derived from many completed or ongoing genome projects, many of whose IIC ATPases have not been phylogenetically analyzed previously. We show that the most likely origin of IIC proteins is prokaryotic, and that many of them are present in non-metazoans, such as algae, protozoans or fungi. We also suggest that the pre-metazoan ancestor, represented by the choanoflagellate Monosiga brevicollis, whose genome has been sequenced, presented at least two IIC-type proteins. One of these proteins would have given rise to most current animal IIC ATPases, whereas the other apparently evolved into a lineage that, so far, has only been found in nematodes. We also propose that early deuterostomes presented a single IIC gene, from which all the extant diversity of vertebrate IIC proteins originated by gene and genome duplications.

Phylogenomic analyses of lophophorates (brachiopods, phoronids and bryozoans) confirm the Lophotrochozoa concept

Based on embryological and morphological evidence, Lophophorata was long considered to be the sister or paraphyletic stem group of Deuterostomia. By contrast, molecular data have consistently indicated that the three lophophorate lineages, Ectoprocta, Brachiopoda and Phoronida, are more closely related to trochozoans (annelids, molluscs and related groups) than to deuterostomes. For this reason, the lophophorate groups and Trochozoa were united to Lophotrochozoa. However, the relationships of the lophophorate lineages within Lophotrochozoa are still largely unresolved. Maximum-likelihood and Bayesian analyses were performed based on a dataset comprising 11,445 amino acid positions derived from 79 ribosomal proteins of 39 metazoan taxa including new sequences obtained from a brachiopod and a phoronid. These analyses show that the three lophophorate lineages are affiliated with trochozoan rather than deuterostome phyla. All hypotheses claiming that they are more closely related to Deuterostomia than to Protostomia can be rejected by topology testing. Monophyly of lophophorates was not recovered but that of Bryozoa including Ectoprocta and Entoprocta and monophyly of Brachiozoa including Brachiopoda and Phoronida were strongly supported. Alternative hypotheses that are refuted include (i) Brachiozoa as the sister group of Mollusca, (ii) ectoprocts as sister to all other Lophotrochozoa including Platyzoa, and (iii) ectoprocts as sister or to all other protostomes except chaetognaths.

Complete nucleotide sequences of mitochondrial genomes of two solitary entoprocts, Loxocorone allax and Loxosomella aloxiata: implications for lophotrochozoan phylogeny

The complete nucleotide sequences of the mitochondrial (mt) genomes of the entoprocts Loxocorone allax and Loxosomella aloxiata were determined. Both species carry the typical gene set of metazoan mt genomes and have similar organizations of their mt genes. However, they show differences in the positions of two tRNA(Leu) genes. Additionally, the tRNA(Val) gene, and half of the long non-coding region, is duplicated and inverted in the Loxos. aloxiata mt genome. The initiation codon of the Loxos. aloxiata cytochrome oxidase subunit I gene is expected to be ACG rather than AUG. The mt gene organizations in these two entoproct species most closely resemble those of mollusks such as Katharina tunicata and Octopus vulgaris, which have the most evolutionarily conserved mt gene organization reported to date in mollusks. Analyses of the mt gene organization in the lophotrochozoan phyla (Annelida, Brachiopoda, Echiura, Entoprocta, Mollusca, Nemertea, and Phoronida) suggested a close phylogenetic relationship between Brachiopoda, Annelida, and Echiura. However, Phoronida was excluded from this grouping. Molecular phylogenetic analyses based on the sequences of mt protein-coding genes suggested a possible close relationship between Entoprocta and Phoronida, and a close relationship among Brachiopoda, Annelida, and Echiura.

EST sequencing of Onychophora and phylogenomic analysis of Metazoa

Onychophora (velvet worms) represent a small animal taxon considered to be related to Euarthropoda. We have obtained 1873 5' cDNA sequences (expressed sequence tags, ESTs) from the velvet worm Epiperipatus sp., which were assembled into 833 contigs. BLAST similarity searches revealed that 51.9% of the contigs had matches in the protein databases with expectation values lower than 10(-4). Most ESTs had the best hit with proteins from either Chordata or Arthropoda (approximately 40% respectively). The ESTs included sequences of 27 ribosomal proteins. The orthologous sequences from 28 other species of a broad range of phyla were obtained from the databases, including other EST projects. A concatenated amino acid alignment comprising 5021 positions was constructed, which covers 4259 positions when problematic regions were removed. Bayesian and maximum likelihood methods place Epiperipatus within the monophyletic Ecdysozoa (Onychophora, Arthropoda, Tardigrada and Nematoda), but its exact relation to the Euarthropoda remained unresolved. The "Articulata" concept was not supported. Tardigrada and Nematoda formed a well-supported monophylum, suggesting that Tardigrada are actually Cycloneuralia. In agreement with previous studies, we have demonstrated that random sequencing of cDNAs results in sequence information suitable for phylogenomic approaches to resolve metazoan relationships.

Comparative effects of aluminum and ouabain on synaptosomal choline uptake, acetylcholine release and (Na+/K+)ATPase

Closing the gap between adverse health effects of aluminum and its mechanisms of action still represents a huge challenge. Cholinergic dysfunction has been implicated in neuronal injury induced by aluminum. Previously reported data also indicate that in vivo and in vitro exposure to aluminum inhibits the mammalian (Na(+)/K(+))ATPase, an ubiquitous plasma membrane pump. This study was undertaken with the specific aim of determining whether in vitro exposure to AlCl(3) and ouabain, the foremost utilized selective inhibitor of (Na(+)/K(+))ATPase, induce similar functional modifications of cholinergic presynaptic nerve terminals, by comparing their effects on choline uptake, acetylcholine release and (Na(+)/K(+))ATPase activity, on subcellular fractions enriched in synaptic nerve endings isolated from rat brain, cuttlefish optic lobe and torpedo electric organ. Results obtained show that choline uptake by rat synaptosomes was inhibited by submillimolar AlCl(3), whereas the amount of choline taken up by synaptosomes isolated from cuttlefish and torpedo remained unchanged. Conversely, choline uptake was reduced by ouabain to a large extent in all synaptosomal preparations analyzed. In contrast to ouabain, which modified the K(+) depolarization evoked release of acetylcholine by rat, cuttlefish and torpedo synaptosomal fractions, AlCl(3) induced reduction of stimulated acetylcholine release was only observed when rat synaptosomes were challenged. Finally, it was observed that the aluminum effect on cuttlefish and torpedo synaptosomal (Na(+)/K(+))ATPase activity was slight when compared to its inhibitory action on mammalian (Na(+)/K(+))ATPase. In conclusion, inhibition of (Na(+)/K(+))ATPase by AlCl(3) and ouabain jeopardized the high-affinity (Na(+)-dependent, hemicholinium-3 sensitive) uptake of choline and the Ca(2+)-dependent, K(+) depolarization evoked release of acetylcholine by rat, cuttlefish and torpedo synaptosomal fractions. The effects of submillimolar AlCl(3) on choline uptake and acetylcholine release only resembled those of ouabain when rat synaptosomes were assayed. Therefore, important differences were found between the species regarding the cholinotoxic action of aluminum. The variability of (Na(+)/K(+))ATPase sensitivity to aluminum of cholinergic neurons might contribute to their differential susceptibility to this neurotoxic agent.

Local decoding of sequences and alignment-free comparison

Subword composition plays an important role in a lot of analyses of sequences. Here we define and study the "local decoding of order N of sequences," an alternative that avoids some drawbacks of "subwords of length N" approaches while keeping informations about environments of length N in the sequences ("decoding" is taken here in the sense of hidden Markov modeling, i.e., associating some state to all positions of the sequence). We present an algorithm for computing the local decoding of order N of a given set of sequences. Its complexity is linear in the total length of the set (whatever the order N) both in time and memory space. In order to show a use of local decoding, we propose a very basic dissimilarity measure between sequences which can be computed both from local decoding of order N and composition in subwords of length N. The accuracies of these two dissimilarities are evaluated, over several datasets, by computing their linear correlations with a reference alignment-based distance. These accuracies are also compared to the one obtained from another recent alignment-free comparison.

The short-lived fish Nothobranchius furzeri as a new model system for aging studies

Genetic and pharmacological research on aging is hampered by the lifespan of available vertebrate models. We recently initiated studies on Nothobranchius furzeri, a species with a maximum life expectancy in captivity of just three months which represents the shortest documented captive lifespan for a vertebrate. Further research on N. furzeri has demonstrated that 1. Short lifespan is tied with explosive growth and accelerated sexual maturation. 2. Short lifespan is correlated with expression of age-related behavioral and histological changes. 3. Lifespan and expression of age-related markers can be modulated by water temperature. 4. Resveratrol, a drug characterized for its life-extending action in Caenorhabditis elegans and Drosophila, increases lifespan and retards expression of age-related markers. 5. Aging-related genes can be easily isolated by homology cloning. Finally, different populations or species of Nothobranchius show large-scale differences in captive lifespan. In the last three years, N. furzeri has moved from biological curiosity to a promising model system for drug validation. Furthermore, this species occupies a favorable position in the Teleost's "tree of life". It is very close to the Japanese Medaka, and close to the pufferfishes and stickleback and might represent a very useful model for comparative genomics of aging.

The partial mitochondrial genome of the Cephalothrix rufifrons (Nemertea, Palaeonemertea): characterization and implications for the phylogenetic position of Nemertea

A continuous 10.1kb fragment of the Cephalothrix rufifrons (Nemertea, Palaeonemertea) mitochondrial genome was sequenced and characterized to further assess organization of protostome mitochondrial genomes and evaluate the phylogenetic potential of gene arrangement and amino acid characters. The genome is A-T rich (72%), and this biased base composition is partly reflected in codon usage. Inferred tRNA secondary structures are typical of those reported for other metazoan mitochondrial DNAs. The arrangement of the 26 genes contained in the fragment exhibits marked similarity to those of many protostome taxa, most notably molluscs with highly conserved arrangements and a phoronid. Separate and simultaneous phylogenetic analyses of inferred amino acid sequences and gene adjacencies place the nemertean within the protostomes among coelomate lophotrochozoan taxa, but do not find a well-supported sister taxon link.

Darwin's dilemma: the realities of the Cambrian 'explosion'

The Cambrian 'explosion' is widely regarded as one of the fulcrum points in the history of life, yet its origins and causes remain deeply controversial. New data from the fossil record, especially of Burgess Shale-type Lagerstätten, indicate, however, that the assembly of bodyplans is not only largely a Cambrian phenomenon, but can already be documented in fair detail. This speaks against a much more ancient origin of the metazoans, and current work is doing much to reconcile the apparent discrepancies between the fossil record, including the Ediacaran assemblages of latest Neoproterozoic age and molecular 'clocks'. Hypotheses to explain the Cambrian 'explosion' continue to be generated, but the recurrent confusion of cause and effect suggests that the wrong sort of question is being asked. Here I propose that despite its step-like function this evolutionary event is the inevitable consequence of Earth and biospheric change.

Lophotrochozoan phylogeny assessed with LSU and SSU data: Evidence of lophophorate polyphyly

Of the three major bilaterian clades, Lophotrochozoa has the greatest diversity and disparity of body forms and is the least understood in terms of phylogenetic history. Within this clade, small nuclear ribosomal subunit (SSU or 18S) studies have failed to provide resolution and other molecular markers have insufficient taxon sampling. To examine relationships within Lophotrochozoa, we collected and complied complete SSU data and nearly complete (>90%) large nuclear ribosomal subunit (LSU or 28S) data totaling approximately 5kb per taxon, for 36 lophotrochozoans. Results of LSU and combined SSU+LSU likelihood analyses provide topologies more consistent with morphological data than analyses of SSU data alone. Namely, most phyla recognized on morphological grounds are recovered as monophyletic entities when the LSU data is considered (contra SSU data alone). These new data show with significant support that "Lophophorata" (traditionally recognized to include Brachiopoda, Phoronida, and Bryozoa) is not a monophyletic entity. Further, the data suggest that Platyzoa is real and may be derived within lophotrochozans rather than a basal or sister taxon. The recently discovered Cycliophora are allied to entoprocts, consistent with their initial placement based on morphology. Additional evidence for Syndermata (i.e., Rotifera+Acanthocephala) is also found. Although relationships among groups with trochophore-like larvae could not be resolved and nodal support values are generally low, the addition of LSU data is a considerable advance in our understanding of lophotrochozoan phylogeny from the molecular perspective.

The complete mitochondrial genome of Flustrellidra hispida and the phylogenetic position of Bryozoa among the Metazoa

The complete mitochondrial genome of Flustrellidra hispida (Bryozoa, Ctenostomata, Flustrellidridae) was sequenced using a transposon-mediated approach. All but one of the 36 genes were identified (trnS2). The genome is 13,026 bp long, being one of the smallest metazoan mitochondrial genomes sequenced to date with a unique gene order when compared to other Metazoa. The genome has an overall AT richness of 59.4%. We found seven regions of overlaps between tRNAs and protein-coding genes ranging from 2 to 11 nt, and seven regions of overlap between tRNAs, ranging from 1 to 8 nt, resulting in a total number of 46 overlapping nucleotides. Genes nad4, cox2, atp8, and nad3 are terminated by the abbreviated stop codon T and cytb is suggested to terminate on (ACT)AA; we postulate that mRNA editing is required to remove AC for TAA to be functional in terminating translation. Phylogenetic analysis of nucleotide and amino acid data place Flustrellidra in the Lophotrochozoa. DNA for this study originated from two populations resulting in a contig consisting of multiple haplotypes. Twenty-seven SNP sites were detected, the majority occurring in cox1 and nad5. With cox1 already established as a marker in bryozoan studies, we advocate the further testing of nad5.

Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements

Phylogenetic relationships among the metazoan phyla are the subject of an ongoing controversy. Analysis of mitochondrial gene arrangements is a powerful tool to investigate these relationships; however, its previous application outside of individual animal phyla has been hampered by the lack of informative out-group data. To address this shortcoming, we determined complete mitochondrial DNA sequences for the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. With sponges as an outgroup, we investigated phylogenetic relationships of nine bilaterian phyla using both breakpoint analysis of global mitochondrial gene arrangements and maximum parsimony analysis of mitochondrial gene adjacencies. Our results provide strong support for a group that includes protostome (but not deuterostome) coelomate, pseudocoelomate, and acoelomate animals, thus clearly rejecting the Coelomata hypothesis. Two other groups of bilaterian animals, Lophotrochozoa and Ambulacraria, are also supported by our analyses. However, due to the remarkable stability of mitochondrial gene arrangements in Deuterostomia and the Ecdysozoa, conclusions on their evolutionary history cannot be drawn.

Genome-scale evidence of the nematode-arthropod clade

The most extensive phylogenetic analysis carried out to date, including 11 complete genomes, is shown to support the Ecdysozoa hypothesis in the open-ended debate of the Coelomata-Ecdysozoa evolutionary problem.

Background

The issue of whether coelomates form a single clade, the Coelomata, or whether all animals that moult an exoskeleton (such as the coelomate arthropods and the pseudocoelomate nematodes) form a distinct clade, the Ecdysozoa, is the most puzzling issue in animal systematics and a major open-ended subject in evolutionary biology. Previous single-gene and genome-scale analyses designed to resolve the issue have produced contradictory results. Here we present the first genome-scale phylogenetic evidence that strongly supports the Ecdysozoa hypothesis.

Results

Through the most extensive phylogenetic analysis carried out to date, the complete genomes of 11 eukaryotic species have been analyzed in order to find homologous sequences derived from 18 human chromosomes. Phylogenetic analysis of datasets showing an increased adjustment to equal evolutionary rates between nematode and arthropod sequences produced a gradual change from support for Coelomata to support for Ecdysozoa. Transition between topologies occurred when fast-evolving sequences of Caenorhabditis elegans were removed. When chordate, nematode and arthropod sequences were constrained to fit equal evolutionary rates, the Ecdysozoa topology was statistically accepted whereas Coelomata was rejected.

Conclusions

The reliability of a monophyletic group clustering arthropods and nematodes was unequivocally accepted in datasets where traces of the long-branch attraction effect were removed. This is the first phylogenomic evidence to strongly support the 'moulting clade' hypothesis.

Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia

Almost a decade ago, a new phylogeny of bilaterian animals was inferred from small-subunit ribosomal RNA (rRNA) that claimed the monophyly of two major groups of protostome animals: Ecdysozoa (e.g., arthropods, nematodes, onychophorans, and tardigrades) and Lophotrochozoa (e.g., annelids, molluscs, platyhelminths, brachiopods, and rotifers). However, it received little additional support. In fact, several multigene analyses strongly argued against this new phylogeny. These latter studies were based on a large amount of sequence data and therefore showed an apparently strong statistical support. Yet, they covered only a few taxa (those for which complete genomes were available), making systematic artifacts of tree reconstruction more probable. Here we expand this sparse taxonomic sampling and analyze a large data set (146 genes, 35,371 positions) from a diverse sample of animals (35 species). Our study demonstrates that the incongruences observed between rRNA and multigene analyses were indeed due to long-branch attraction artifacts, illustrating the enormous impact of systematic biases on phylogenomic studies. A refined analysis of our data set excluding the most biased genes provides strong support in favor of the new animal phylogeny and in addition suggests that urochordates are more closely related to vertebrates than are cephalochordates. These findings have important implications for the interpretation of morphological and genomic data.