Phylogenomic reconstruction of Solenogastres (Mollusca, Aplacophora) informs hypotheses on body size evolution

Body size is a fundamental characteristic of animals that impacts every aspect of their biology from anatomical complexity to ecology. In Mollusca, Solenogastres has been considered important to understanding the group's early evolution as most morphology-based phylogenetic reconstructions placed it as an early branching molluscan lineage. Under this scenario, molluscs were thought to have evolved from a small, turbellarian-like ancestor and small (i.e., macrofaunal) body size was inferred to be plesiomorphic for Solenogastres. More recently, phylogenomic studies have shown that aplacophorans (Solenogastres + Caudofoveata) form a clade with chitons (Polyplacophora), which is sister to all other molluscs, suggesting a relatively large-bodied (i.e., megafaunal) ancestor for Mollusca. Meanwhile, recent investigations into aplacophoran phylogeny have called the assumption that the last common ancestor of Solenogastres was small-bodied into question, but sampling of meiofaunal species was limited, biasing these studies towards large-bodied taxa and leaving fundamental questions about solenogaster body size evolution unanswered. Here, we supplemented available data with transcriptomes from eight diverse meiofaunal species of Solenogastres and conducted phylogenomic analyses on datasets of up to 949 genes. Maximum likelihood analyses support the meiofaunal family Meiomeniidae as the sister group to all other solenogasters, congruent with earlier ideas of a small-bodied ancestor of Solenogastres. In contrast, Bayesian Inference analyses support the large-bodied family Amphimeniidae as the sister group to all other solenogasters. Investigation of phylogenetic signal by comparing site-wise likelihood scores for the two competing hypotheses support the Meiomeniidae-first topology. In light of these results, we performed ancestral character state reconstruction to explore the implications of both hypotheses on understanding of Solenogaster evolution and review previous hypotheses about body size evolution and its potential consequences for solenogaster biology. Both hypotheses imply that body size evolution has been highly dynamic over the course of solenogaster evolution and that their relatively static body plan has successfully allowed for evolutionary transitions between meio-, macro- and megafaunal size ranges.

Genomic and transcriptomic survey of bryozoan Hox and ParaHox genes with emphasis on phylactolaemate bryozoans

BACKGROUND

Bryozoans are mostly sessile aquatic colonial invertebrates belonging to the clade Lophotrochozoa, which unites many protostome bilaterian phyla such as molluscs, annelids and brachiopods. While Hox and ParaHox genes have been extensively studied in various lophotrochozoan lineages, investigations on Hox and ParaHox gene complements in bryozoans are scarce.

RESULTS

Herein, we present the most comprehensive survey of Hox and ParaHox gene complements in bryozoans using four genomes and 35 transcriptomes representing all bryozoan clades: Cheilostomata, Ctenostomata, Cyclostomata and Phylactolaemata. Using similarity searches, phylogenetic analyses and detailed manual curation, we have identified five Hox genes in bryozoans (pb, Dfd, Lox5, Lox4 and Post2) and one ParaHox gene (Cdx). Interestingly, we observed lineage-specific duplication of certain Hox and ParaHox genes (Dfd, Lox5 and Cdx) in some bryozoan lineages.

CONCLUSIONS

The bryozoan Hox cluster does not retain the ancestral lophotrochozoan condition but appears relatively simple (includes only five genes) and broken into two genomic regions, characterized by the loss and duplication of serval genes. Importantly, bryozoans share the lack of two Hox genes (Post1 and Scr) with their proposed sister-taxon, Phoronida, which suggests that those genes were missing in the most common ancestor of bryozoans and phoronids.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

The embryology, metamorphosis, and muscle development of Schizocardium karankawa sp. nov. (Enteropneusta) from the Gulf of Mexico

Schizocardium karankawa sp. nov. has been collected from subtidal muds of the Laguna Madre, Texas, and the Mississippi coast, Gulf of Mexico. The Texas population is reproductive from early February to mid-April. Gametes are liberated by a small incision in a gonad. Oocyte germinal vesicle breakdown is increased in the presence of sperm, and the highest fertilization success was in the artificial seawater Jamarin U. Manually dechorionated embryos develop normally. Development was asynchronous via a tornaria larva, metamorphosis and maintained to the juvenile worm 6 gill-pore stage. Phalloidin-labeled late-stage tornaria revealed retractor muscles that connect the pericardial sac with the apical tuft anteriorly, the oesophagus ventrally, and muscle cells of the early mesocoels. The muscle development of early juvenile worms began with dorso-lateral trunk muscles, lateral trunk bands, and sphincters around the gill pores and anus. Adult worms are characterized by a stomochord that bifurcates anteriorly into paired vermiform processes, gill bars that extend almost the entire dorsal to ventral branchial region resulting in a narrow ventral hypobranchial ridge, and an elaborate epibranchial organ with six zones of discrete cell types. The trunk has up to three rows of liver sacs, and lateral gonads. The acorn worm evo-devo model species Saccoglossus kowalevskii, Ptychodera flava, and Schizocardium californicum are phylogenetically distant with disparate life histories. S. karnakawa from S. californicum are phylogenetically close, and differences between them that become apparent as adult worms include the number of gill pores and hepatic sacs, and elaborations of the heart-kidney-stomochord complex. An important challenge for evolutionary developmental biology is to form links from phylogenetically distant and large-scale differences to phylogenetically close and small-scale differences. This description of the embryology, development, and adult morphology of S. karankawa permits investigations into how acorn worm development evolves at fine scales.

Phylogenomics reveals deep relationships and diversification within phylactolaemate bryozoans

Bryozoans are mostly sessile colonial invertebrates that inhabit all kinds of aquatic ecosystems. Extant bryozoan species fall into two clades with one of them, Phylactolaemata, being the only exclusively freshwater clade. Phylogenetic relationships within the class Phylactolaemata have long been controversial owing to their limited distinguishable characteristics that reflect evolutionary relationships. Here, we present the first phylogenomic analysis of Phylactolaemata using transcriptomic data combined with dense taxon sampling of six families to better resolve the interrelationships and to estimate divergence time. Using maximum-likelihood and Bayesian inference approaches, we recovered a robust phylogeny for Phylactolaemata in which the interfamilial relationships are fully resolved. We show Stephanellidae is the sister taxon of all other phylactolaemates and confirm that Lophopodidae represents the second offshoot within the phylactolaemate tree. Plumatella fruticosa clearly falls outside Plumatellidae as previous investigations have suggested, and instead clusters with Pectinatellidae and Cristatellidae as the sister taxon of Fredericellidae. Our results demonstrate that cryptic speciation is very likely in F. sultana and in two species of Plumatella (P. repens and P. casmiana). Divergence time estimates show that Phylactolaemata appeared at the end of the Ediacaran and started to diverge in the Silurian, although confidence intervals were large for most nodes. The radiation of most extant phylactolaemate families occurred mainly in the Palaeogene and Neogene highlighting post-extinction diversification.

Different phylogenomic methods support monophyly of enigmatic 'Mesozoa' (Dicyemida + Orthonectida, Lophotrochozoa)

Dicyemids and orthonectids were traditionally classified in a group called Mesozoa, but their placement in a single clade has been contested and their position(s) within Metazoa is uncertain. Here, we assembled a comprehensive matrix of Lophotrochozoa (Metazoa) and investigated the position of Dicyemida (= Rhombozoa) and Orthonectida, employing multiple phylogenomic approaches. We sequenced seven new transcriptomes and one draft genome from dicyemids (Dicyema, Dicyemennea) and two transcriptomes from orthonectids (Rhopalura). Using these and published data, we assembled and analysed contamination-filtered datasets with up to 987 genes. Our results recover Mesozoa monophyletic and as a close relative of Platyhelminthes or Gnathifera. Because of the tendency of the long-branch mesozoans to group with other long-branch taxa in our analyses, we explored the impact of approaches purported to help alleviate long-branch attraction (e.g. taxon removal, coalescent inference, gene targeting). None of these were able to break the association of Orthonectida with Dicyemida in the maximum-likelihood trees. Contrastingly, the Bayesian analysis and site-specific frequency model in maximum-likelihood did not recover a monophyletic Mesozoa (but only when using a specific 50 gene matrix). The classic hypothesis on monophyletic Mesozoa is possibly reborn and should be further tested.

Genome of the lepidopleurid chiton Hanleya hanleyi (Mollusca, Polyplacophora)

Mollusca is the second most species-rich phylum and includes animals as disparate as octopuses, clams, and chitons. Dozens of molluscan genomes are available, but only one representative of the subphylum Aculifera, the sister taxon to all other molluscs, has been sequenced to date, hindering comparative and evolutionary studies. To facilitate evolutionary studies across Mollusca, we sequenced the genome of a second aculiferan mollusc, the lepidopleurid chiton Hanleya hanleyi (Bean 1844), using a hybrid approach combining Oxford Nanopore and Illumina reads. After purging redundant haplotigs and removing contamination from this 1.3% heterozygous genome, we produced a 2.5 Gbp haploid assembly (>4X the size of the other chiton genome sequenced to date) with an N50 of 65.0 Kbp. Despite a fragmented assembly, the genome is rather complete (92.0% of BUSCOs detected; 79.4% complete plus 12.6% fragmented). Remarkably, the genome has the highest repeat content of any molluscan genome reported to date (>66%). Our gene annotation pipeline predicted 69,284 gene models (92.9% of BUSCOs detected; 81.8% complete plus 11.1% fragmented) of which 35,362 were supported by transcriptome and/or protein evidence. Phylogenomic analysis recovered Polyplacophora sister to all other sampled molluscs with maximal support. The Hanleya genome will be a valuable resource for studies of molluscan biology with diverse potential applications ranging from evolutionary and comparative genomics to molecular ecology.

Phylogenomic resolution of the root of Panpulmonata, a hyperdiverse radiation of gastropods: new insight into the evolution of air breathing

Transitions to terrestriality have been associated with major animal radiations including land snails and slugs in Stylommatophora (>20 000 described species), the most successful lineage of 'pulmonates' (a non-monophyletic assemblage of air-breathing gastropods). However, phylogenomic studies have failed to robustly resolve relationships among traditional pulmonates and affiliated marine lineages that comprise clade Panpulmonata (Mollusca, Gastropoda), especially two key taxa: Sacoglossa, a group including photosynthetic sea slugs, and Siphonarioidea, intertidal limpet-like snails with a non-contractile pneumostome (narrow opening to a vascularized pallial cavity). To clarify the evolutionary history of the panpulmonate radiation, we performed phylogenomic analyses on datasets of up to 1160 nuclear protein-coding genes for 110 gastropods, including 40 new transcriptomes for Sacoglossa and Siphonarioidea. All 18 analyses recovered Sacoglossa as the sister group to a clade we named Pneumopulmonata, within which Siphonarioidea was sister to the remaining lineages in most analyses. Comparative modelling indicated shifts to marginal habitat (estuarine, mangrove and intertidal zones) preceded and accelerated the evolution of a pneumostome, present in the pneumopulmonate ancestor along with a one-sided plicate gill. These findings highlight key intermediate stages in the evolution of air-breathing snails, supporting the hypothesis that adaptation to marginal zones played an important role in major sea-to-land transitions.

Phylogenomic Analysis of the Phylum Nematoda: Conflicts and Congruences With Morphology, 18S rRNA, and Mitogenomes

Phylogenetic relationships within many lineages of the phylum Nematoda remain unresolved, despite numerous morphology-based and molecular analyses. We performed several phylogenomic analyses using 286 published genomes and transcriptomes and 19 new transcriptomes by focusing on Trichinellida, Spirurina, Rhabditina, and Tylenchina separately, and by analyzing a selection of species from the whole phylum Nematoda. The phylogeny of Trichinellida supported the division of Trichinella into encapsulated and non-encapsulated species and placed them as sister to Trichuris. The Spirurina subtree supported the clades formed by species from Ascaridomorpha and Spiruromorpha respectively, but did not support Dracunculoidea. The analysis of Tylenchina supported a clade that included all sampled species from Tylenchomorpha and placed it as sister to clades that included sampled species from Cephalobomorpha and Panagrolaimomorpha, supporting the hypothesis that postulates the single origin of the stomatostylet. The Rhabditina subtree placed a clade composed of all sampled species from Diplogastridae as sister to a lineage consisting of paraphyletic Rhabditidae, a single representative of Heterorhabditidae and a clade composed of sampled species belonging to Strongylida. It also strongly supported all suborders within Strongylida. In the phylum-wide analysis, a clade composed of all sampled species belonging to Enoplia were consistently placed as sister to Dorylaimia + Chromadoria. The topology of the Nematoda backbone was consistent with previous studies, including polyphyletic placement of sampled representatives of Monhysterida and Araeolaimida.

Benchmarking Oxford Nanopore read assemblers for high-quality molluscan genomes

Choosing the optimum assembly approach is essential to achieving a high-quality genome assembly suitable for comparative and evolutionary genomic investigations. Significant recent progress in long-read sequencing technologies such as PacBio and Oxford Nanopore Technologies (ONT) has also brought about a large variety of assemblers. Although these have been extensively tested on model species such as Homo sapiens and Drosophila melanogaster, such benchmarking has not been done in Mollusca, which lacks widely adopted model species. Molluscan genomes are notoriously rich in repeats and are often highly heterozygous, making their assembly challenging. Here, we benchmarked 10 assemblers based on ONT raw reads from two published molluscan genomes of differing properties, the gastropod Chrysomallon squamiferum (356.6 Mb, 1.59% heterozygosity) and the bivalve Mytilus coruscus (1593 Mb, 1.94% heterozygosity). By optimizing the assembly pipeline, we greatly improved both genomes from previously published versions. Our results suggested that 40-50X of ONT reads are sufficient for high-quality genomes, with Flye being the recommended assembler for compact and less heterozygous genomes exemplified by C. squamiferum, while NextDenovo excelled for more repetitive and heterozygous molluscan genomes exemplified by M. coruscus. A phylogenomic analysis using the two updated genomes with 32 other published high-quality lophotrochozoan genomes resulted in maximum support across all nodes, and we show that improved genome quality also leads to more complete matrices for phylogenomic inferences. Our benchmarking will ensure efficiency in future assemblies for molluscs and perhaps also for other marine phyla with few genomes available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

On the diversity of abyssal Dondersiidae (Mollusca: Aplacophora) with the description of a new genus, six new species, and a review of the family

So far, of the 292 known species of solenogasters (Mollusca, Aplacophora), 62 belong to the clade Pholidoskepia Salvini-Plawen, 1978. Of these, only two have an abyssal distribution (3500-6000 m depth). Among Pholidoskepia, Dondersiidae Simroth, 1893 is the most diverse family. This study contributes to the knowledge of this family with the description of one new genus and six new species from the abyssal South Atlantic Ocean: Dondersia ? foraminosa sp. n., Nematomenia divae sp. n., Nematomenia brasiliensis sp. n., Nematomenia ? guineana sp. n., Helluoherpia vieiralaneroi sp. n. and Inopinatamenia (gen. n.) calamitosa sp. n. Specimens were collected during DIVA (Latitudinal Gradients of Deep-Sea BioDIVersity in the Atlantic Ocean) expeditions in the Guinea (DIVA 2 Me 63/2, 2005) and Brazil (DIVA 3 Me 79/1, 2008) Basins. Specimens were characterized based primarily on the sclerites and internal anatomy, which was studied using histology. The importance of the radula and mantle sclerites for taxonomy is emphasized. Amended diagnoses for the family and some genera within this family are provided. This contribution increases the described diversity of Dondersiidae to ten genera and 38 species and highlights the need for more study of solenogasters in the deep sea.

Assessment of mitochondrial genomes for heterobranch gastropod phylogenetics

Background

Heterobranchia is a diverse clade of marine, freshwater, and terrestrial gastropod molluscs. It includes such disparate taxa as nudibranchs, sea hares, bubble snails, pulmonate land snails and slugs, and a number of (mostly small-bodied) poorly known snails and slugs collectively referred to as the “lower heterobranchs”. Evolutionary relationships within Heterobranchia have been challenging to resolve and the group has been subject to frequent and significant taxonomic revision. Mitochondrial (mt) genomes can be a useful molecular marker for phylogenetics but, to date, sequences have been available for only a relatively small subset of Heterobranchia.

Results

To assess the utility of mitochondrial genomes for resolving evolutionary relationships within this clade, eleven new mt genomes were sequenced including representatives of several groups of “lower heterobranchs”. Maximum likelihood analyses of concatenated matrices of the thirteen protein coding genes found weak support for most higher-level relationships even after several taxa with extremely high rates of evolution were excluded. Bayesian inference with the CAT + GTR model resulted in a reconstruction that is much more consistent with the current understanding of heterobranch phylogeny. Notably, this analysis recovered Valvatoidea and Orbitestelloidea in a polytomy with a clade including all other heterobranchs, highlighting these taxa as important to understanding early heterobranch evolution. Also, dramatic gene rearrangements were detected within and between multiple clades. However, a single gene order is conserved across the majority of heterobranch clades.

Conclusions

Analysis of mitochondrial genomes in a Bayesian framework with the site heterogeneous CAT + GTR model resulted in a topology largely consistent with the current understanding of heterobranch phylogeny. However, mitochondrial genomes appear to be too variable to serve as good phylogenetic markers for robustly resolving a number of deeper splits within this clade.

The Iron-Responsive Genome of the Chiton Acanthopleura granulata

Molluscs biomineralize structures that vary in composition, form, and function, prompting questions about the genetic mechanisms responsible for their production and the evolution of these mechanisms. Chitons (Mollusca, Polyplacophora) are a promising system for studies of biomineralization because they build a range of calcified structures including shell plates and spine- or scale-like sclerites. Chitons also harden the calcified teeth of their rasp-like radula with a coat of iron (as magnetite). Here we present the genome of the West Indian fuzzy chiton Acanthopleura granulata, the first from any aculiferan mollusc. The A. granulata genome contains homologs of many genes associated with biomineralization in conchiferan molluscs. We expected chitons to lack genes previously identified from pathways conchiferans use to make biominerals like calcite and nacre because chitons do not use these materials in their shells. Surprisingly, the A. granulata genome has homologs of many of these genes, suggesting that the ancestral mollusc may have had a more diverse biomineralization toolkit than expected. The A. granulata genome has features that may be specialized for iron biomineralization, including a higher proportion of genes regulated directly by iron than other molluscs. A. granulata also produces two isoforms of soma-like ferritin: one is regulated by iron and similar in sequence to the soma-like ferritins of other molluscs, and the other is constitutively translated and is not found in other molluscs. The A. granulata genome is a resource for future studies of molluscan evolution and biomineralization.

Zooid morphology and molecular phylogeny of the graptolite Rhabdopleura annulata (Hemichordata, Pterobranchia) from Heron Island, Australia

Rhabdopleura Allman, 1869 is one of the longest surviving animal genera. The five-known species are the only living Graptolithina, a group well known from their diverse Paleozoic fossil record. Here we add information on the soft-bodied zooids and molecular phylogenetics of Rhabdopleura annulata Norman, 1921, which was previously only known from its tubes. Tubes and zooids were collected from Heron Island, Queensland, Australia. Zooids have a single pair of tentaculated arms. Dark pigment granules are found throughout the body, and particularly dense in the pair of arms and the anterior lip of the cephalic shield. Colonies grow encrusted in and on coral debris. The tubes are either creeping or erect, but no stolon has been found. Inside of the coral matrix lacunae, the tube cortex formed a parchment-like wallpaper. Phylogenetic analysis based on combined 18S+16S rRNA sequences placed R. annulata as sister to the remaining rhabdopleurids, albeit with weak support. The biogeographic range of R. annulata extends from Indonesia to Tasmania, and New Zealand. Its occurrence on Heron Island does not extend this range, but highlights that rhabdopleurids may be more common, and in shallower waters, than previously appreciated, permitting further studies that may shed light on graptolite paleobiology.

Phylogenomics of Aplacophora (Mollusca, Aculifera) and a solenogaster without a foot

Recent molecular phylogenetic investigations strongly supported the placement of the shell-less, worm-shaped aplacophoran molluscs (Solenogastres and Caudofoveata) and chitons (Polyplacophora) in a clade called Aculifera, which is the sister taxon of all other molluscs. Thus, understanding the evolutionary history of aculiferan molluscs is important for understanding early molluscan evolution. In particular, fundamental questions about evolutionary relationships within Aplacophora have long been unanswered. Here, we supplemented the paucity of available data with transcriptomes from 25 aculiferans and conducted phylogenomic analyses on datasets with up to 525 genes and 75 914 amino acid positions. Our results indicate that aplacophoran taxonomy requires revision as several traditionally recognized groups are non-monophyletic. Most notably, Cavibelonia, the solenogaster taxon defined by hollow sclerites, is polyphyletic, suggesting parallel evolution of hollow sclerites in multiple lineages. Moreover, we describe Apodomenia enigmatica sp. nov., a bizarre new species that appears to be a morphological intermediate between Solenogastres and Caudofoveata. This animal is not a missing link, however; molecular and morphological studies show that it is a derived solenogaster that lacks a foot, mantle cavity and radula. Taken together, these results shed light on the evolutionary history of Aplacophora and reveal a surprising degree of morphological plasticity within the group.

Complete mitochondrial genomes of two scaphopod molluscs

Complete mitochondrial genomes were determined for two scaphopod molluscs: the dentaliid Antalis entalis and an unidentified Antarctic gadilid. Both genomes are complete except, in Gadilida sp. indet., a short stretch of nad5 was undetermined and trnR could not be annotated. Organization of the Gadilida sp. genome is nearly identical to that previously reported for the gadilid Siphonodentalium whereas trnK, nad5, trnD, nad4, and nad4l are transposed to the opposite strand in the previously published Graptacme genome relative to that of Antalis. Phylogenetic analysis of the 13 protein-coding and 2 rRNA genes recovered Scaphopoda, Gadilida, and Dentaliida monophyletic with maximal support.

Revisiting metazoan phylogeny with genomic sampling of all phyla

Proper biological interpretation of a phylogeny can sometimes hinge on the placement of key taxa-or fail when such key taxa are not sampled. In this light, we here present the first attempt to investigate (though not conclusively resolve) animal relationships using genome-scale data from all phyla. Results from the site-heterogeneous CAT + GTR model recapitulate many established major clades, and strongly confirm some recent discoveries, such as a monophyletic Lophophorata, and a sister group relationship between Gnathifera and Chaetognatha, raising continued questions on the nature of the spiralian ancestor. We also explore matrix construction with an eye towards testing specific relationships; this approach uniquely recovers support for Panarthropoda, and shows that Lophotrochozoa (a subclade of Spiralia) can be constructed in strongly conflicting ways using different taxon- and/or orthologue sets. Dayhoff-6 recoding sacrifices information, but can also reveal surprising outcomes, e.g. full support for a clade of Lophophorata and Entoprocta + Cycliophora, a clade of Placozoa + Cnidaria, and raising support for Ctenophora as sister group to the remaining Metazoa, in a manner dependent on the gene and/or taxon sampling of the matrix in question. Future work should test the hypothesis that the few remaining uncertainties in animal phylogeny might reflect violations of the various stationarity assumptions used in contemporary inference methods.

Improved phylogenomic sampling of free-living nematodes enhances resolution of higher-level nematode phylogeny

Background

Nematodes are among the most diverse and abundant metazoans on Earth, but research on them has been biased toward parasitic taxa and model organisms. Free-living nematodes, particularly from the clades Enoplia and Dorylaimia, have been underrepresented in genome-scale phylogenetic analyses to date, leading to poor resolution of deep relationships within the phylum.

Results

We supplemented publicly available data by sequencing transcriptomes of nine free-living nematodes and two important outgroups and conducted a phylum-wide phylogenomic analysis including a total of 108 nematodes. Analysis of a dataset generated using a conservative orthology inference strategy resulted in a matrix with a high proportion of missing data and moderate to weak support for branching within and placement of Enoplia. A less conservative orthology inference approach recovered more genes and resulted in higher support for the deepest splits within Nematoda, recovering Enoplia as the sister taxon to the rest of Nematoda. Relationships within major clades were similar to those found in previously published studies based on 18S rDNA.

Conclusions

Expanded transcriptome sequencing of free-living nematodes has contributed to better resolution among deep nematode lineages, though the dataset is still strongly biased toward parasites. Inclusion of more free-living nematodes in future phylogenomic analyses will allow a clearer understanding of many interesting aspects of nematode evolution, such as morphological and molecular adaptations to parasitism and whether nematodes originated in a marine or terrestrial environment.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1444-x) contains supplementary material, which is available to authorized users.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

Conservation of mitochondrial genome arrangements in brittle stars (Echinodermata, Ophiuroidea)

Brittle stars are conspicuous members of benthic ecosystems, fill many ecological niches and are the most speciose of all classes of echinoderms. With high levels of biodiversity, elucidating the evolutionary history of this group is important. Understanding of higher-level relationships within Ophiuroidea has been aided by multilocus nuclear data and DNA barcoding. However, the degree of consistency between mitochondrial and nuclear data within ophiuroids remains unclear and deserves further assessment. In this study, 17 mitochondrial genomes spanning the taxonomic breadth of Ophiuroidea were utilized to explore evolutionary relationships through maximum likelihood analyses, Bayesian inference and comparative assessment of gene order. Our phylogenetic analyses, based on both nucleotide and amino acid residues, support recent findings based on multilocus nuclear data and morphology, in that the brittle star clades Ophintegrida and Euryophiurida were recovered as monophyletic with the latter comprising Euyalida, Ophiuridae and Ophiopyrgidae. Only three different arrangements of the 13 protein coding and 2 ribosomal RNA genes were observed. As expected, tRNA genes were more likely to have undergone rearrangement but the order of all 37 genes was found to be conserved in all sampled Euryalida and Ophiuridae. Both Euryalida and the clade comprised of Ophiuridae and Ophiopyrgidae, each had their own conserved rearrangement of protein coding genes and ribosomal genes, after divergence from their last common ancestor. Euryalida has a rearrangement of the two ribosomal RNA genes, rrnS and rrnL, in contrast to Ophiuridae and Ophiopyrgidae, which had an inversion of the genes nad1, nad2, and cob relative to Ophintegrida. Further, our data support the gene order found in all sampled Euryalida as the most likely ancestral order for all Ophiuroidea.

Mitogenomics reveals phylogenetic relationships of caudofoveate aplacophoran molluscs

The worm-shaped, shell-less aplacophoran molluscs Caudofoveata and Solenogastres have recently received attention as part of the clade Aculifera, but relationships within these two lineages are still largely unknown. Here, we use complete mitochondrial genomes to shed light on higher-level relationships within Caudofoveata. Mitochondrial genomes have been sequenced for many diverse molluscs, but only two mitochondrial genomes from aplacophoran molluscs (the caudofoveates Scutopus ventrolineatus and Chaetoderma nitidulum) are available to date. We sequenced and assembled complete or near complete mitochondrial genomes of five additional species of Caudofoveata (Falcidens acutargatus, Falcidens halanychi, Scutopus robustus, Psilodens balduri and Spathoderma clenchi) and one species of Solenogastres (Neomenia carinata) for comparison to available mitochondrial genomes of aculiferans. Comparison of mitochondrial gene order among different lineages revealed a highly conserved order of protein coding genes corresponding to the hypothesized ancestral gene order for Mollusca. Unique arrangements of tRNAs were found among lineages of aculiferan molluscs as well as among caudofoveate taxa. Phylogenetic analyses of amino acid sequences for the 13 protein-coding genes recovered a monophyletic Aplacophora. Within Caudofoveata, Chaetodermatidae, but not Limifossoridae, was recovered monophyletic. Taken together, our results suggest that mitochondrial genomes can serve as useful molecular markers for aculiferan phylogenetics, especially for more recent phylogenetic events.

Step-wise evolution of complex chemical defenses in millipedes: a phylogenomic approach

With fossil representatives from the Silurian capable of respiring atmospheric oxygen, millipedes are among the oldest terrestrial animals, and likely the first to acquire diverse and complex chemical defenses against predators. Exploring the origin of complex adaptive traits is critical for understanding the evolution of Earth's biological complexity, and chemical defense evolution serves as an ideal study system. The classic explanation for the evolution of complexity is by gradual increase from simple to complex, passing through intermediate "stepping stone" states. Here we present the first phylogenetic-based study of the evolution of complex chemical defenses in millipedes by generating the largest genomic-based phylogenetic dataset ever assembled for the group. Our phylogenomic results demonstrate that chemical complexity shows a clear pattern of escalation through time. New pathways are added in a stepwise pattern, leading to greater chemical complexity, independently in a number of derived lineages. This complexity gradually increased through time, leading to the advent of three distantly related chemically complex evolutionary lineages, each uniquely characteristic of each of the respective millipede groups.

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.].

Ctenophore relationships and their placement as the sister group to all other animals

Ctenophora, compromising approximately 200 described species, is an important lineage for understanding metazoan evolution and is of great ecological and economic importance. Ctenophore diversity includes species with unique colloblasts used for prey capture, smooth and striated muscles, benthic and pelagic lifestyles, and locomotion with ciliated paddles or muscular propulsion. However, ancestral states of traits are debated and relationships among many lineages are unresolved. Here, using 27 newly sequenced ctenophore transcriptomes, publicly available data, and methods to control systematic error we establish the placement of Ctenophora as the sister group to all other animals and refine phylogenetic relationships within ctenophores. Molecular clock analyses suggest modern ctenophore diversity originated approximately 350MYA ± 88 MY, conflicting with previous hypotheses of approximately 65 MYA. We recover Euplokamis dunlapae, a species with striated muscles, as the sister lineage to other sampled ctenophores. Ancestral state reconstruction shows the most recent common ancestor of extant ctenophores was pelagic, possessed tentacles, was bioluminescent, and did not have separate sexes. Our results imply at least two transitions from a pelagic to a benthic lifestyle within Ctenophora, suggesting such transitions were more common in animal diversification than appreciated.

The crown-of-thorns starfish genome as a guide for biocontrol of this coral reef pest

The crown-of-thorns starfish (COTS, the Acanthaster planci species group) is a highly fecund predator of reef-building corals throughout the Indo-Pacific region. COTS population outbreaks cause substantial loss of coral cover, diminishing the integrity and resilience of reef ecosystems. Here we sequenced genomes of COTS from the Great Barrier Reef, Australia and Okinawa, Japan to identify gene products that underlie species-specific communication and could potentially be used in biocontrol strategies. We focused on water-borne chemical plumes released from aggregating COTS, which make the normally sedentary starfish become highly active. Peptide sequences detected in these plumes by mass spectrometry are encoded in the COTS genome and expressed in external tissues. The exoproteome released by aggregating COTS consists largely of signalling factors and hydrolytic enzymes, and includes an expanded and rapidly evolving set of starfish-specific ependymin-related proteins. These secreted proteins may be detected by members of a large family of olfactory-receptor-like G-protein-coupled receptors that are expressed externally, sometimes in a sex-specific manner. This study provides insights into COTS-specific communication that may guide the generation of peptide mimetics for use on reefs with COTS outbreaks.

Molecular clocks indicate turnover and diversification of modern coleoid cephalopods during the Mesozoic Marine Revolution

Coleoid cephalopod molluscs comprise squid, cuttlefish and octopuses, and represent nearly the entire diversity of modern cephalopods. Sophisticated adaptations such as the use of colour for camouflage and communication, jet propulsion and the ink sac highlight the unique nature of the group. Despite these striking adaptations, there are clear parallels in ecology between coleoids and bony fishes. The coleoid fossil record is limited, however, hindering confident analysis of the tempo and pattern of their evolution. Here we use a molecular dataset (180 genes, approx. 36 000 amino acids) of 26 cephalopod species to explore the phylogeny and timing of cephalopod evolution. We show that crown cephalopods diverged in the Silurian-Devonian, while crown coleoids had origins in the latest Palaeozoic. While the deep-sea vampire squid and dumbo octopuses have ancient origins extending to the Early Mesozoic Era, 242 ± 38 Ma, incirrate octopuses and the decabrachian coleoids (10-armed squid) diversified in the Jurassic Period. These divergence estimates highlight the modern diversity of coleoid cephalopods emerging in the Mesozoic Marine Revolution, a period that also witnessed the radiation of most ray-finned fish groups in addition to several other marine vertebrates. This suggests that that the origin of modern cephalopod biodiversity was contingent on ecological competition with marine vertebrates.

Monoplacophoran mitochondrial genomes: convergent gene arrangements and little phylogenetic signal

BACKGROUND

Although recent studies have greatly advanced understanding of deep molluscan phylogeny, placement of some taxa remains uncertain as different datasets support competing class-relationships. Traditionally, morphologists have placed Monoplacophora, a group of morphologically simple, limpet-like molluscs as sister group to all other conchiferans (shelled molluscs other than Polyplacophora), a grouping that is supported by the latest large-scale phylogenomic study that includes Laevipilina. However, molecular datasets dominated by nuclear ribosomal genes support Monoplacophora + Polyplacophora (Serialia). Here, we evaluate the potential of mitochondrial genome data for resolving placement of Monoplacophora.

RESULTS

Two complete (Laevipilina antarctica and Vema ewingi) and one partial (Laevipilina hyalina) mitochondrial genomes were sequenced, assembled, and compared. All three genomes show a highly similar architecture including an unusually high number of non-coding regions. Comparison of monoplacophoran gene order shows a gene arrangement pattern not previously reported; there is an inversion of one large gene cluster. Our reanalyses of recently published polyplacophoran mitogenomes show, however, that this feature is also present in some chiton species. Maximum Likelihood and Bayesian Inference analyses of 13 mitochondrial protein-coding genes failed to robustly place Monoplacophora and hypothesis testing could not reject any of the evaluated placements of Monoplacophora.

CONCLUSIONS

Under both serialian or aculiferan-conchiferan scenarios, the observed gene cluster inversion appears to be a convergent evolution of gene arrangements in molluscs. Our phylogenetic results are inconclusive and sensitive to taxon sampling. Aculifera (Polyplacophora + Aplacophora) and Conchifera were never recovered. However, some analyses recovered Serialia (Monoplacophora + Polyplacophora), Diasoma (Bivalvia + Scaphopoda) or Pleistomollusca (Bivalvia + Gastropoda). Although we could not shed light on deep evolutionary traits of Mollusca we found unique patterns of gene arrangements that are common to monoplacophoran and chitonine polyplacophoran species but not to acanthochitonine Polyplacophora. Complete mitochondrial genome of Laevipilina antarctica.

Sea shell diversity and rapidly evolving secretomes: insights into the evolution of biomineralization

An external skeleton is an essential part of the body plan of many animals and is thought to be one of the key factors that enabled the great expansion in animal diversity and disparity during the Cambrian explosion. Molluscs are considered ideal to study the evolution of biomineralization because of their diversity of highly complex, robust and patterned shells. The molluscan shell forms externally at the interface of animal and environment, and involves controlled deposition of calcium carbonate within a framework of macromolecules that are secreted from the dorsal mantle epithelium. Despite its deep conservation within Mollusca, the mantle is capable of producing an incredible diversity of shell patterns, and macro- and micro-architectures. Here we review recent developments within the field of molluscan biomineralization, focusing on the genes expressed in the mantle that encode secreted proteins. The so-called mantle secretome appears to regulate shell deposition and patterning and in some cases becomes part of the shell matrix. Recent transcriptomic and proteomic studies have revealed marked differences in the mantle secretomes of even closely-related molluscs; these typically exceed expected differences based on characteristics of the external shell. All mantle secretomes surveyed to date include novel genes encoding lineage-restricted proteins and unique combinations of co-opted ancient genes. A surprisingly large proportion of both ancient and novel secreted proteins containing simple repetitive motifs or domains that are often modular in construction. These repetitive low complexity domains (RLCDs) appear to further promote the evolvability of the mantle secretome, resulting in domain shuffling, expansion and loss. RLCD families further evolve via slippage and other mechanisms associated with repetitive sequences. As analogous types of secreted proteins are expressed in biomineralizing tissues in other animals, insights into the evolution of the genes underlying molluscan shell formation may be applied more broadly to understanding the evolution of metazoan biomineralization.

Spider phylogenomics: untangling the Spider Tree of Life

Spiders (Order Araneae) are massively abundant generalist arthropod predators that are found in nearly every ecosystem on the planet and have persisted for over 380 million years. Spiders have long served as evolutionary models for studying complex mating and web spinning behaviors, key innovation and adaptive radiation hypotheses, and have been inspiration for important theories like sexual selection by female choice. Unfortunately, past major attempts to reconstruct spider phylogeny typically employing the "usual suspect" genes have been unable to produce a well-supported phylogenetic framework for the entire order. To further resolve spider evolutionary relationships we have assembled a transcriptome-based data set comprising 70 ingroup spider taxa. Using maximum likelihood and shortcut coalescence-based approaches, we analyze eight data sets, the largest of which contains 3,398 gene regions and 696,652 amino acid sites forming the largest phylogenomic analysis of spider relationships produced to date. Contrary to long held beliefs that the orb web is the crowning achievement of spider evolution, ancestral state reconstructions of web type support a phylogenetically ancient origin of the orb web, and diversification analyses show that the mostly ground-dwelling, web-less RTA clade diversified faster than orb weavers. Consistent with molecular dating estimates we report herein, this may reflect a major increase in biomass of non-flying insects during the Cretaceous Terrestrial Revolution 125-90 million years ago favoring diversification of spiders that feed on cursorial rather than flying prey. Our results also have major implications for our understanding of spider systematics. Phylogenomic analyses corroborate several well-accepted high level groupings: Opisthothele, Mygalomorphae, Atypoidina, Avicularoidea, Theraphosoidina, Araneomorphae, Entelegynae, Araneoidea, the RTA clade, Dionycha and the Lycosoidea. Alternatively, our results challenge the monophyly of Eresoidea, Orbiculariae, and Deinopoidea. The composition of the major paleocribellate and neocribellate clades, the basal divisions of Araneomorphae, appear to be falsified. Traditional Haplogynae is in need of revision, as our findings appear to support the newly conceived concept of Synspermiata. The sister pairing of filistatids with hypochilids implies that some peculiar features of each family may in fact be synapomorphic for the pair. Leptonetids now are seen as a possible sister group to the Entelegynae, illustrating possible intermediates in the evolution of the more complex entelegyne genitalic condition, spinning organs and respiratory organs.

Employing Phylogenomics to Resolve the Relationships among Cnidarians, Ctenophores, Sponges, Placozoans, and Bilaterians

Despite an explosion in the amount of sequence data, phylogenomics has failed to settle controversy regarding some critical nodes on the animal tree of life. Understanding relationships among Bilateria, Ctenophora, Cnidaria, Placozoa, and Porifera is essential for studying how complex traits such as neurons, muscles, and gastrulation have evolved. Recent studies have cast doubt on the historical viewpoint that sponges are sister to all other animal lineages with recent studies recovering ctenophores as sister. However, the ctenophore-sister hypothesis has been criticized as unrealistic and caused by systematic error. We review past phylogenomic studies and potential causes of systematic error in an effort to identify areas that can be improved in future studies. Increased sampling of taxa, less missing data, and a priori removal of sequences and taxa that may cause systematic error in phylogenomic inference will likely be the most fruitful areas of focus when assembling future datasets. Ultimately, we foresee metazoan relationships being resolved with higher support in the near future, and we caution against dismissing novel hypotheses merely because they conflict with historical viewpoints of animal evolution.

Reconstruction of cyclooxygenase evolution in animals suggests variable, lineage-specific duplications, and homologs with low sequence identity

Cyclooxygenase (COX) enzymatically converts arachidonic acid into prostaglandin G/H in animals and has importance during pregnancy, digestion, and other physiological functions in mammals. COX genes have mainly been described from vertebrates, where gene duplications are common, but few studies have examined COX in invertebrates. Given the increasing ease in generating genomic data, as well as recent, although incomplete descriptions of potential COX sequences in Mollusca, Crustacea, and Insecta, assessing COX evolution across Metazoa is now possible. Here, we recover 40 putative COX orthologs by searching publicly available genomic resources as well as ~250 novel invertebrate transcriptomic datasets. Results suggest the common ancestor of Cnidaria and Bilateria possessed a COX homolog similar to those of vertebrates, although such homologs were not found in poriferan and ctenophore genomes. COX was found in most crustaceans and the majority of molluscs examined, but only specific taxa/lineages within Cnidaria and Annelida. For example, all octocorallians appear to have COX, while no COX homologs were found in hexacorallian datasets. Most species examined had a single homolog, although species-specific COX duplications were found in members of Annelida, Mollusca, and Cnidaria. Additionally, COX genes were not found in Hemichordata, Echinodermata, or Platyhelminthes, and the few previously described COX genes in Insecta lacked appreciable sequence homology (although structural analyses suggest these may still be functional COX enzymes). This analysis provides a benchmark for identifying COX homologs in future genomic and transcriptomic datasets, and identifies lineages for future studies of COX.

Nemertean toxin genes revealed through transcriptome sequencing

Nemerteans are one of few animal groups that have evolved the ability to utilize toxins for both defense and subduing prey, but little is known about specific nemertean toxins. In particular, no study has identified specific toxin genes even though peptide toxins are known from some nemertean species. Information about toxin genes is needed to better understand evolution of toxins across animals and possibly provide novel targets for pharmaceutical and industrial applications. We sequenced and annotated transcriptomes of two free-living and one commensal nemertean and annotated an additional six publicly available nemertean transcriptomes to identify putative toxin genes. Approximately 63–74% of predicted open reading frames in each transcriptome were annotated with gene names, and all species had similar percentages of transcripts annotated with each higher-level GO term. Every nemertean analyzed possessed genes with high sequence similarities to known animal toxins including those from stonefish, cephalopods, and sea anemones. One toxin-like gene found in all nemerteans analyzed had high sequence similarity to Plancitoxin-1, a DNase II hepatotoxin that may function well at low pH, which suggests that the acidic body walls of some nemerteans could work to enhance the efficacy of protein toxins. The highest number of toxin-like genes found in any one species was seven and the lowest was three. The diversity of toxin-like nemertean genes found here is greater than previously documented, and these animals are likely an ideal system for exploring toxin evolution and industrial applications of toxins.

Phylogenomic resolution of the hemichordate and echinoderm clade

Ambulacraria, comprising Hemichordata and Echinodermata, is closely related to Chordata, making it integral to understanding chordate origins and polarizing chordate molecular and morphological characters. Unfortunately, relationships within Hemichordata and Echinodermata have remained unresolved, compromising our ability to extrapolate findings from the most closely related molecular and developmental models outside of Chordata (e.g., the acorn worms Saccoglossus kowalevskii and Ptychodera flava and the sea urchin Strongylocentrotus purpuratus). To resolve long-standing phylogenetic issues within Ambulacraria, we sequenced transcriptomes for 14 hemichordates as well as 8 echinoderms and complemented these with existing data for a total of 33 ambulacrarian operational taxonomic units (OTUs). Examination of leaf stability values revealed rhabdopleurid pterobranchs and the enteropneust Stereobalanus canadensis were unstable in placement; therefore, analyses were also run without these taxa. Analyses of 185 genes resulted in reciprocal monophyly of Enteropneusta and Pterobranchia, placed the deep-sea family Torquaratoridae within Ptychoderidae, and confirmed the position of ophiuroid brittle stars as sister to asteroid sea stars (the Asterozoa hypothesis). These results are consistent with earlier perspectives concerning plesiomorphies of Ambulacraria, including pharyngeal gill slits, a single axocoel, and paired hydrocoels and somatocoels. The resolved ambulacrarian phylogeny will help clarify the early evolution of chordate characteristics and has implications for our understanding of major fossil groups, including graptolites and somasteroideans.

Repurposed transcriptomic data facilitate discovery of innate immunity toll-like receptor (TLR) Genes across Lophotrochozoa

The growing volume of genomic data from across life represents opportunities for deriving valuable biological information from data that were initially collected for another purpose. Here, we use transcriptomes collected for phylogenomic studies to search for toll-like receptor (TLR) genes in poorly sampled lophotrochozoan clades (Annelida, Mollusca, Brachiopoda, Phoronida, and Entoprocta) and one ecdysozoan clade (Priapulida). TLR genes are involved in innate immunity across animals by recognizing potential microbial infection. They have an extracellular leucine-rich repeat (LRR) domain connected to a transmembrane domain and an intracellular toll/interleukin-1 receptor (TIR) domain. Consequently, these genes are important in initiating a signaling pathway to trigger defense. We found at least one TLR ortholog in all but two taxa examined, suggesting that a broad array of lophotrochozoans may have innate immune systems similar to those observed in vertebrates and arthropods. Comparison to the SMART database confirmed the presence of both the LRR and the TIR protein motifs characteristic of TLR genes. Because we looked at only one transcriptome per species, discovery of TLR genes was limited for most taxa. However, several TRL-like genes that vary in the number and placement of LRR domains were found in phoronids. Additionally, several contigs contained LRR domains but lacked TIR domains, suggesting they were not TLRs. Many of these LRR-containing contigs had other domains (e.g., immunoglobin) and are likely involved in innate immunity.

Immunofluorescence localization of the tubulin cytoskeleton during cell division and cell growth in members of the Coleochaetales (Streptophyta)

Study of charophycean green algae, including the Coleochaetales, may shed light on the evolutionary history of characters they share with their land plant relatives. We examined the tubulin cytoskeleton during mitosis, cytokinesis, and growth in members of the Coleochaetales with diverse morphologies to determine if phragmoplasts occurred throughout this order and to identify microtubular patterns associated with cell growth. Species representing three subgroups of Coleochaete and its sister genus Chaetosphaeridium were studied. Cytokinesis involving a phragmoplast was found in the four taxa examined. Differential interference contrast microscopy of living cells confirmed that polar cytokinesis like that described in the model flowering plant Arabidopsis occurred in all species when the forming cell plate traversed a vacuole. Calcofluor labeling of cell walls demonstrated directed growth from particular cell regions of all taxa. Electron microscopy confirmed directed growth in the unusual growth pattern of Chaetosphaeridium. All four species exhibited unordered microtubule patterns associated with diffuse growth in early cell expansion. In subsequent elongating cells, Coleochaete irregularis Pringsheim and Chaetosphaeridium globosum (Nordstedt) Klebahn exhibited tubulin cytoskeleton arrays corresponding to growth patterns associated with tip growth in plants, fungi, and other charophycean algae. Hoop-shaped microtubules frequently associated with diffuse growth of elongating cells in plants were not observed in any of these species. Presence of phragmoplasts in the diverse species studied supports the hypothesis that cytokinesis involving a phragmoplast originated in a common ancestor of the Coleochaetales, and possibly in a common ancestor of Charales, Coleochaetales, Zygnematales, and plants.

PhyloTreePruner: A Phylogenetic Tree-Based Approach for Selection of Orthologous Sequences for Phylogenomics

Molecular phylogenetics relies on accurate identification of orthologous sequences among the taxa of interest. Most orthology inference programs available for use in phylogenomics rely on small sets of pre-defined orthologs from model organisms or phenetic approaches such as all-versus-all sequence comparisons followed by Markov graph-based clustering. Such approaches have high sensitivity but may erroneously include paralogous sequences. We developed PhyloTreePruner, a software utility that uses a phylogenetic approach to refine orthology inferences made using phenetic methods. PhyloTreePruner checks single-gene trees for evidence of paralogy and generates a new alignment for each group containing only sequences inferred to be orthologs. Importantly, PhyloTreePruner takes into account support values on the tree and avoids unnecessarily deleting sequences in cases where a weakly supported tree topology incorrectly indicates paralogy. A test of PhyloTreePruner on a dataset generated from 11 completely sequenced arthropod genomes identified 2,027 orthologous groups sampled for all taxa. Phylogenetic analysis of the concatenated supermatrix yielded a generally well-supported topology that was consistent with the current understanding of arthropod phylogeny. PhyloTreePruner is freely available from http://sourceforge.net/projects/phylotreepruner/.

Hard and soft anatomy in two genera of Dondersiidae (Mollusca, Aplacophora, Solenogastres)

Phylogenetic relationships and identifications in the aplacophoran taxon Solenogastres (Neomeniomorpha) are in flux largely because descriptions of hard parts--sclerites, radulae, copulatory spicules--and body shape have often not been adequately illustrated or utilized. With easily recognizable and accessible hard parts, descriptions of Solenogastres are of greater use, not just to solenogaster taxonomists, but also to ecologists, paleontologists, and evolutionary biologists. Phylogenetic studies of Aplacophora, Mollusca, and the Lophotrochozoa as a whole, whether morphological or molecular, would be enhanced. As an example, morphologic characters, both isolated hard parts and internal anatomy, are provided for two genera in the Dondersiidae. Five species are described or redescribed and earlier descriptions corrected and enhanced. Three belong to Dondersia: D. festiva Hubrecht, D. incali (Scheltema), and D. namibiensis n. sp., the latter differentiated unambiguously from D. incali only by sclerites and copulatory spicules. Two species belong to Lyratoherpia: L. carinata Salvini-Plawen and L. californica (Heath). Notes are given for other species in Dondersiidae: L. bracteata Salvini-Plawen, Ichthyomenia ichthyodes (Pruvot), and Heathia porosa (Heath). D. indica Stork is synonymized with D. annulata. A cladistic morphological analysis was conducted to examine the utility of hard parts for reconstructing solenogaster phylogeny. Results indicate monophyly of Dondersia and Lyratoherpia as described here.

Transcriptomic evidence for the expression of horizontally transferred algal nuclear genes in the photosynthetic sea slug, Elysia chlorotica

Analysis of the transcriptome of the kleptoplastic sea slug, Elysia chlorotica, has revealed the presence of at least 101 chloroplast-encoded gene sequences and 111 transcripts matching 52 nuclear-encoded genes from the chloroplast donor, Vaucheria litorea. These data clearly show that the symbiotic chloroplasts are translationally active and, of even more interest, that a variety of functional algal genes have been transferred into the slug genome, as has been suggested by earlier indirect experiments. Both the chloroplast- and nuclear-encoded sequences were rare within the E. chlorotica transcriptome, suggesting that their copy numbers and synthesis rates are low, and required both a large amount of sequence data and native algal sequences to find. These results show that the symbiotic chloroplasts residing inside the host molluscan cell are maintained by an interaction of both organellar and host biochemistry directed by the presence of transferred genes.

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).