Kevin M. Kocot

The ctenophore genome and the evolutionary origins of neural systems

The origins of neural systems remain unresolved. In contrast to other basal metazoans, ctenophores (comb jellies) have both complex nervous and mesoderm-derived muscular systems. These holoplanktonic predators also have sophisticated ciliated locomotion, behaviour and distinct development. Here we present the draft genome of Pleurobrachia bachei, Pacific sea gooseberry, together with ten other ctenophore transcriptomes, and show that they are remarkably distinct from other animal genomes in their content of neurogenic, immune and developmental genes. Our integrative analyses place Ctenophora as the earliest lineage within Metazoa. This hypothesis is supported by comparative analysis of multiple gene families, including the apparent absence of HOX genes, canonical microRNA machinery, and reduced immune complement in ctenophores. Although two distinct nervous systems are well recognized in ctenophores, many bilaterian neuron-specific genes and genes of ‘classical’ neurotransmitter pathways either are absent or, if present, are not expressed in neurons. Our metabolomic and physiological data are consistent with the hypothesis that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals.

Phylogenomics reveals deep molluscan relationships

Evolutionary relationships among the eight major lineages of Mollusca have remained unresolved despite their diversity and importance. Previous investigations of molluscan phylogeny, based primarily on nuclear ribosomal gene sequences or morphological data, have been unsuccessful at elucidating these relationships. Recently, phylogenomic studies using dozens to hundreds of genes have greatly improved our understanding of deep animal relationships. However, limited genomic resources spanning molluscan diversity has prevented use of a phylogenomic approach. Here we use transcriptome and genome data from all major lineages (except Monoplacophora) and recover a well-supported topology for Mollusca. Our results strongly support the Aculifera hypothesis placing Polyplacophora (chitons) in a clade with a monophyletic Aplacophora (worm-like molluscs). Additionally, within Conchifera, a sister-taxon relationship between Gastropoda and Bivalvia is supported. This grouping has received little consideration and contains most (>95%) molluscan species. Thus we propose the node-based name Pleistomollusca. In light of these results, we examined the evolution of morphological characters and found support for advanced cephalization and shells as possibly having multiple origins within Mollusca.

Error, signal, and the placement of Ctenophora sister to all other animals

Significance Traditional interpretation of animal phylogeny suggests traits, such as mesoderm, muscles, and neurons, evolved only once given the assumed placement of sponges as sister to all other animals. In contrast, placement of ctenophores as the first branching animal lineage raises the possibility of multiple origins of many complex traits considered important for animal diversification and success. We consider sources of potential error and increase taxon sampling to find a single, statistically robust placement of ctenophores as our most distant animal relatives, contrary to the traditional understanding of animal phylogeny. Furthermore, ribosomal protein genes are identified as creating conflict in signal that caused some past studies to recover a sister relationship between ctenophores and cnidarians. Elucidating relationships among early animal lineages has been difficult, and recent phylogenomic analyses place Ctenophora sister to all other extant animals, contrary to the traditional view of Porifera as the earliest-branching animal lineage. To date, phylogenetic support for either ctenophores or sponges as sister to other animals has been limited and inconsistent among studies. Lack of agreement among phylogenomic analyses using different data and methods obscures how complex traits, such as epithelia, neurons, and muscles evolved. A consensus view of animal evolution will not be accepted until datasets and methods converge on a single hypothesis of early metazoan relationships and putative sources of systematic error (e.g., long-branch attraction, compositional bias, poor model choice) are assessed. Here, we investigate possible causes of systematic error by expanding taxon sampling with eight novel transcriptomes, strictly enforcing orthology inference criteria, and progressively examining potential causes of systematic error while using both maximum-likelihood with robust data partitioning and Bayesian inference with a site-heterogeneous model. We identified ribosomal protein genes as possessing a conflicting signal compared with other genes, which caused some past studies to infer ctenophores and cnidarians as sister. Importantly, biases resulting from elevated compositional heterogeneity or elevated substitution rates are ruled out. Placement of ctenophores as sister to all other animals, and sponge monophyly, are strongly supported under multiple analyses, herein.

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.

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.

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.

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

Phylogenomics supports Panpulmonata: Opisthobranch paraphyly and key evolutionary steps in a major radiation of gastropod molluscs

Pulmonates, with over 30,000 described species, represent the largest radiation of non-marine animals outside of Arthropoda. The pulmonate lung was a key evolutionary innovation enabling diversification of terrestrial and freshwater snails and slugs. However, recent studies drew conflicting conclusions about pulmonate monophyly, and support for a sister group is lacking, hindering our understanding of this major animal radiation. Analyses of mitochondrial protein-coding genes recovered a paraphyletic Pulmonata grading into a monophyletic Opisthobranchia, a traditional group of sea slugs long considered sister to pulmonates. Conversely, analyses of datsets dominated by nuclear rDNA indicated Opisthobranchia is paraphyletic with respect to Pulmonata. No study resolved the placement of two key taxa: Sacoglossa, an opisthobranch group including photosynthetic sea slugs, and Siphonarioidea, intertidal limpet-like snails traditionally in Pulmonata. To examine evolutionary relationships at the base of the pulmonate radiation, we performed a phylogenomic analysis of 102 nuclear protein-coding gene regions for 19 gastropods. Opisthobranchia was recovered as paraphyletic with respect to Panpulmonata, a clade in which Sacoglossa was sister to Pulmonata, with Siphonarioidea as the basal pulmonate lineage. Siphonarioideans share a similar gill structure with shelled sacoglossans but lack the contractile pneumostome of pulmonates, suggesting descent from an evolutionary intermediate that facilitated the pulmonate radiation into non-marine habitats.

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

Phylogenomics of Lophotrochozoa with consideration of systematic error

&NA; 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.]