Giants among Cnidaria: Large Nuclear Genomes and Rearranged Mitochondrial Genomes in Siphonophores

Siphonophores (Cnidaria: Hydrozoa) are abundant predators found throughout the ocean and are important constituents of the global zooplankton community. They range in length from a few centimeters to tens of meters. They are gelatinous, fragile, and difficult to collect, so many aspects of the biology of these roughly 200 species remain poorly understood. To survey siphonophore genome diversity, we performed Illumina sequencing of 32 species sampled broadly across the phylogeny. Sequencing depth was sufficient to estimate nuclear genome size from k-mer spectra in six specimens, ranging from 0.7 to 2.3 Gb, with heterozygosity estimates between 0.69% and 2.32%. Incremental k-mer counting indicates k-mer peaks can be absent with nearly 20× read coverage, suggesting minimum genome sizes range from 1.4 to 5.6 Gb in the 25 samples without peaks in the k-mer spectra. This work confirms most siphonophore nuclear genomes are large relative to the genomes of other cnidarians, but also identifies several with reduced size that are tractable targets for future siphonophore nuclear genome assembly projects. We also assembled complete mitochondrial genomes for 33 specimens from these new data, indicating a conserved gene order shared among nonsiphonophore hydrozoans, Cystonectae, and some Physonectae, revealing the ancestral mitochondrial gene order of siphonophores. Our results also suggest extensive rearrangement of mitochondrial genomes within other Physonectae and in Calycophorae. Though siphonophores comprise a small fraction of cnidarian species, this survey greatly expands our understanding of cnidarian genome diversity. This study further illustrates both the importance of deep phylogenetic sampling and the utility of k-mer-based genome skimming in understanding the genomic diversity of a clade.

Cell type evolution reconstruction across species through cell phylogenies of single-cell RNA sequencing data

The origin and evolution of cell types has emerged as a key topic in evolutionary biology. Driven by rapidly accumulating single-cell datasets, recent attempts to infer cell type evolution have largely been limited to pairwise comparisons because we lack approaches to build cell phylogenies using model-based approaches. Here we approach the challenges of applying explicit phylogenetic methods to single-cell data by using principal components as phylogenetic characters. We infer a cell phylogeny from a large, comparative single-cell dataset of eye cells from five distantly related mammals. Robust cell type clades enable us to provide a phylogenetic, rather than phenetic, definition of cell type, allowing us to forgo marker genes and phylogenetically classify cells by topology. We further observe evolutionary relationships between diverse vessel endothelia and identify the myelinating and non-myelinating Schwann cells as sister cell types. Finally, we examine principal component loadings and describe the gene expression dynamics underlying the function and identity of cell type clades that have been conserved across the five species. A cell phylogeny provides a rigorous framework towards investigating the evolutionary history of cells and will be critical to interpret comparative single-cell datasets that aim to ask fundamental evolutionary questions.

Normalizing need not be the norm: count-based math for analyzing single-cell data

Counting transcripts of mRNA are a key method of observation in modern biology. With advances in counting transcripts in single cells (single-cell RNA sequencing or scRNA-seq), these data are routinely used to identify cells by their transcriptional profile, and to identify genes with differential cellular expression. Because the total number of transcripts counted per cell can vary for technical reasons, the first step of many commonly used scRNA-seq workflows is to normalize by sequencing depth, transforming counts into proportional abundances. The primary objective of this step is to reshape the data such that cells with similar biological proportions of transcripts end up with similar transformed measurements. But there is growing concern that normalization and other transformations result in unintended distortions that hinder both analyses and the interpretation of results. This has led to an intense focus on optimizing methods for normalization and transformation of scRNA-seq data. Here, we take an alternative approach, by avoiding normalization and transformation altogether. We abandon the use of distances to compare cells, and instead use a restricted algebra, motivated by measurement theory and abstract algebra, that preserves the count nature of the data. We demonstrate that this restricted algebra is sufficient to draw meaningful and practical comparisons of gene expression through the use of the dot product and other elementary operations. This approach sidesteps many of the problems with common transformations, and has the added benefit of being simpler and more intuitive. We implement our approach in the package countland, available in python and R.

Not all clusters are equal: dynamics of molecular HIV-1 clusters in a statewide Rhode Island epidemic

OBJECTIVES

Molecular epidemiology is a powerful tool to characterize HIV epidemics and prioritize public health interventions. Typically, HIV clusters are assumed to have uniform patterns over time. We hypothesized that assessment of cluster evolution would reveal distinct cluster behavior, possibly improving molecular epidemic characterization, towards disrupting HIV transmission.

DESIGN

Retrospective cohort.

METHODS

Annual phylogenies were inferred by cumulative aggregation of all available HIV-1 pol sequences of individuals with HIV-1 in Rhode Island (RI) between 1990 and 2020, representing a statewide epidemic. Molecular clusters were detected in annual phylogenies by strict and relaxed cluster definition criteria, and the impact of annual newly-diagnosed HIV-1 cases to the structure of individual clusters was examined over time.

RESULTS

Of 2153 individuals, 31% (strict criteria) - 47% (relaxed criteria) clustered. Longitudinal tracking of individual clusters identified three cluster types: normal, semi-normal and abnormal. Normal clusters (83-87% of all identified clusters) showed predicted growing/plateauing dynamics, with approximately three-fold higher growth rates in large (15-18%) vs. small (∼5%) clusters. Semi-normal clusters (1-2% of all clusters) temporarily fluctuated in size and composition. Abnormal clusters (11-16% of all clusters) demonstrated collapses and re-arrangements over time. Borderline values of cluster-defining parameters explained dynamics of non-normal clusters.

CONCLUSIONS

Comprehensive tracing of molecular HIV clusters over time in a statewide epidemic identified distinct cluster types, likely missed in cross-sectional analyses, demonstrating that not all clusters are equal. This knowledge challenges current perceptions of consistent cluster behavior over time and could improve molecular surveillance of local HIV epidemics to better inform public health strategies.

The evolution of ovary-biased gene expression in Hawaiian Drosophila

With detailed data on gene expression accessible from an increasingly broad array of species, we can test the extent to which our developmental genetic knowledge from model organisms predicts expression patterns and variation across species. But to know when differences in gene expression across species are significant, we first need to know how much evolutionary variation in gene expression we expect to observe. Here we provide an answer by analyzing RNAseq data across twelve species of Hawaiian Drosophilidae flies, focusing on gene expression differences between the ovary and other tissues. We show that over evolutionary time, there exists a cohort of ovary specific genes that is stable and that largely corresponds to described expression patterns from laboratory model Drosophila species. Our results also provide a demonstration of the prediction that, as phylogenetic distance increases, variation between species overwhelms variation between tissue types. Using ancestral state reconstruction of expression, we describe the distribution of evolutionary changes in tissue-biased expression, and use this to identify gains and losses of ovary-biased expression across these twelve species. We then use this distribution to calculate the evolutionary correlation in expression changes between genes, and demonstrate that genes with known interactions in D. melanogaster are significantly more correlated in their evolution than genes with no or unknown interactions. Finally, we use this correlation matrix to infer new networks of genes that share evolutionary trajectories, and we present these results as a dataset of new testable hypotheses about genetic roles and interactions in the function and evolution of the Drosophila ovary.

Characterizing the secret diets of siphonophores (Cnidaria: Hydrozoa) using DNA metabarcoding

Siphonophores (Cnidaria: Hydrozoa) are abundant and diverse gelatinous predators in open-ocean ecosystems. Due to limited access to the midwater, little is known about the diets of most deep-dwelling gelatinous species, which constrains our understanding of food-web structure and nutrient flow in these vast ecosystems. Visual gut-content methods can rarely identify soft-bodied rapidly-digested prey, while observations from submersibles often overlook small prey items. These methods have been differentially applied to shallow and deep siphonophore taxa, confounding habitat and methodological biases. DNA metabarcoding can be used to assess both shallow and deep species' diets under a common methodological framework, since it can detect both small and gelatinous prey. We (1) further characterized the diets of open-ocean siphonophores using DNA metabarcoding, (2) compared the prey detected by visual and molecular methods to evaluate their technical biases, and (3) evaluated tentacle-based predictions of diet. To do this, we performed DNA metabarcoding analyses on the gut contents of 39 siphonophore species across depths to describe their diets, using six barcode regions along the 18S gene. Taxonomic identifications were assigned using public databases combined with local zooplankton sequences. We identified 55 unique prey items, including crustaceans, gelatinous animals, and fish across 47 siphonophore specimens in 24 species. We reported 29 novel predator-prey interactions, among them the first insights into the diets of nine siphonophore species, many of which were congruent with the dietary predictions based on tentilla morphology. Our analyses detected both small and gelatinous prey taxa underrepresented by visual methods in species from both shallow and deep habitats, indicating that siphonophores play similar trophic roles across depth habitats. We also reveal hidden links between siphonophores and filter-feeders near the base of the food web. This study expands our understanding of the ecological roles of siphonophores in the open ocean, their trophic roles within the 'jelly-web', and the importance of their diversity for nutrient flow and ecosystem functioning. Understanding these inconspicuous yet ubiquitous predator-prey interactions is critical to predict the impacts of climate change, overfishing, and conservation policies on oceanic ecosystems.

Incorporating Within-Host Diversity in Phylogenetic Analyses for Detecting Clusters of New HIV Diagnoses

Background

Phylogenetic analyses of HIV sequences are used to detect clusters and inform public health interventions. Conventional approaches summarize within-host HIV diversity with a single consensus sequence per host of the pol gene, obtained from Sanger or next-generation sequencing (NGS). There is growing recognition that this approach discards potentially important information about within-host sequence variation, which can impact phylogenetic inference. However, whether alternative summary methods that incorporate intra-host variation impact phylogenetic inference of transmission network features is unknown.

Methods

We introduce profile sampling, a method to incorporate within-host NGS sequence diversity into phylogenetic HIV cluster inference. We compare this approach to Sanger- and NGS-derived pol and near-whole-genome consensus sequences and evaluate its potential benefits in identifying molecular clusters among all newly-HIV-diagnosed individuals over six months at the largest HIV center in Rhode Island.

Results

Profile sampling cluster inference demonstrated that within-host viral diversity impacts phylogenetic inference across individuals, and that consensus sequence approaches can obscure both magnitude and effect of these impacts. Clustering differed between Sanger- and NGS-derived consensus and profile sampling sequences, and across gene regions.

Discussion

Profile sampling can incorporate within-host HIV diversity captured by NGS into phylogenetic analyses. This additional information can improve robustness of cluster detection.

Evolution of gene expression across species and specialized zooids in Siphonophora

Siphonophores are complex colonial animals, consisting of asexually produced bodies (zooids) that are functionally specialized for specific tasks, including feeding, swimming, and sexual reproduction. Though this extreme functional specialization has captivated biologists for generations, its genomic underpinnings remain unknown. We use RNA-seq to investigate gene expression patterns in five zooids and one specialized tissue across seven siphonophore species. Analyses of gene expression across species present several challenges, including identification of comparable expression changes on gene trees with complex histories of speciation, duplication, and loss. We examine gene expression within species, conduct classical analyses examining expression patterns between species, and introduce species branch filtering, which allows us to examine the evolution of expression across species in a phylogenetic framework. Within and across species, we identified hundreds of zooid-specific and species-specific genes, as well as a number of putative transcription factors showing differential expression in particular zooids and developmental stages. We found that gene expression patterns tended to be largely consistent in zooids with the same function across species, but also some large lineage-specific shifts in gene expression. Our findings show that patterns of gene expression have the potential to define zooids in colonial organisms. Traditional analyses of the evolution of gene expression focus on the tips of gene phylogenies, identifying large-scale expression patterns that are zooid or species variable. The new explicit phylogenetic approach we propose here focuses on branches (not tips) offering a deeper evolutionary perspective into specific changes in gene expression within zooids along all branches of the gene (and species) trees.

Rooting the Animal Tree of Life

Identifying our most distant animal relatives has emerged as one of the most challenging problems in phylogenetics. This debate has major implications for our understanding of the origin of multicellular animals and of the earliest events in animal evolution, including the origin of the nervous system. Some analyses identify sponges as our most distant animal relatives (Porifera-sister hypothesis), and others identify comb jellies (Ctenophora-sister hypothesis). These analyses vary in many respects, making it difficult to interpret previous tests of these hypotheses. To gain insight into why different studies yield different results, an important next step in the ongoing debate, we systematically test these hypotheses by synthesizing 15 previous phylogenomic studies and performing new standardized analyses under consistent conditions with additional models. We find that Ctenophora-sister is recovered across the full range of examined conditions, and Porifera-sister is recovered in some analyses under narrow conditions when most outgroups are excluded and site-heterogeneous CAT models are used. We additionally find that the number of categories in site-heterogeneous models is sufficient to explain the Porifera-sister results. Furthermore, our cross-validation analyses show CAT models that recover Porifera-sister have hundreds of additional categories and fail to fit significantly better than site-heterogenuous models with far fewer categories. Systematic and standardized testing of diverse phylogenetic models suggests that we should be skeptical of Porifera-sister results both because they are recovered under such narrow conditions and because the models in these conditions fit the data no better than other models that recover Ctenophora-sister.

The Evolutionary History of Siphonophore Tentilla: Novelties, Convergence, and Integration

Synopsis Siphonophores are free-living predatory colonial hydrozoan cnidarians found in every region of the ocean. Siphonophore tentilla (tentacle side branches) are unique biological structures for prey capture, composed of a complex arrangement of cnidocytes (stinging cells) bearing different types of nematocysts (stinging capsules) and auxiliary structures. Tentilla present an extensive morphological and functional diversity across species. While associations between tentillum form and diet have been reported, the evolutionary history giving rise to this morphological diversity is largely unexplored. Here we examine the evolutionary gains and losses of novel tentillum substructures and nematocyst types on the most recent siphonophore phylogeny. Tentilla have a precisely coordinated high-speed strike mechanism of synchronous unwinding and nematocyst discharge. Here we characterize the kinematic diversity of this prey capture reaction using high-speed video and find relationships with morphological characters. Since tentillum discharge occurs in synchrony across a broad morphological diversity, we evaluate how phenotypic integration is maintaining character correlations across evolutionary time. We found that the tentillum morphospace has low dimensionality, identified instances of heterochrony and morphological convergence, and generated hypotheses on the diets of understudied siphonophore species. Our findings indicate that siphonophore tentilla are phenotypically integrated structures with a complex evolutionary history leading to a phylogenetically-structured diversity of forms that are predictive of kinematic performance and feeding habits.

A genome-scale phylogeny of the kingdom Fungi

Phylogenomic studies using genome-scale amounts of data have greatly improved understanding of the tree of life. Despite the diversity, ecological significance, and biomedical and industrial importance of fungi, evolutionary relationships among several major lineages remain poorly resolved, especially those near the base of the fungal phylogeny. To examine poorly resolved relationships and assess progress toward a genome-scale phylogeny of the fungal kingdom, we compiled a phylogenomic data matrix of 290 genes from the genomes of 1,644 species that includes representatives from most major fungal lineages. We also compiled 11 data matrices by subsampling genes or taxa from the full data matrix based on filtering criteria previously shown to improve phylogenomic inference. Analyses of these 12 data matrices using concatenation- and coalescent-based approaches yielded a robust phylogeny of the fungal kingdom, in which ∼85% of internal branches were congruent across data matrices and approaches used. We found support for several historically poorly resolved relationships as well as evidence for polytomies likely stemming from episodes of ancient diversification. By examining the relative evolutionary divergence of taxonomic groups of equivalent rank, we found that fungal taxonomy is broadly aligned with both genome sequence divergence and divergence time but also identified lineages where current taxonomic circumscription does not reflect their levels of evolutionary divergence. Our results provide a robust phylogenomic framework to explore the tempo and mode of fungal evolution and offer directions for future fungal phylogenetic and taxonomic studies.

Revising transcriptome assemblies with phylogenetic information

A common transcriptome assembly error is to mistake different transcripts of the same gene as transcripts from multiple closely related genes. This error is difficult to identify during assembly, but in a phylogenetic analysis such errors can be diagnosed from gene phylogenies where they appear as clades of tips from the same species with improbably short branch lengths. treeinform is a method that uses phylogenetic information across species to refine transcriptome assemblies within species. It identifies transcripts of the same gene that were incorrectly assigned to multiple genes and reassign them as transcripts of the same gene. The treeinform method is implemented in Agalma, available at https://bitbucket.org/caseywdunn/agalma, and the general approach is relevant in a variety of other contexts.

Empirical comparison of analytical approaches for identifying molecular HIV-1 clusters

Public health interventions guided by clustering of HIV-1 molecular sequences may be impacted by choices of analytical approaches. We identified commonly-used clustering analytical approaches, applied them to 1886 HIV-1 Rhode Island sequences from 2004-2018, and compared concordance in identifying molecular HIV-1 clusters within and between approaches. We used strict (topological support ≥ 0.95; distance 0.015 substitutions/site) and relaxed (topological support 0.80-0.95; distance 0.030-0.045 substitutions/site) thresholds to reflect different epidemiological scenarios. We found that clustering differed by method and threshold and depended more on distance than topological support thresholds. Clustering concordance analyses demonstrated some differences across analytical approaches, with RAxML having the highest (91%) mean summary percent concordance when strict thresholds were applied, and three (RAxML-, FastTree regular bootstrap- and IQ-Tree regular bootstrap-based) analytical approaches having the highest (86%) mean summary percent concordance when relaxed thresholds were applied. We conclude that different analytical approaches can yield diverse HIV-1 clustering outcomes and may need to be differentially used in diverse public health scenarios. Recognizing the variability and limitations of commonly-used methods in cluster identification is important for guiding clustering-triggered interventions to disrupt new transmissions and end the HIV epidemic.

Treeio: An R Package for Phylogenetic Tree Input and Output with Richly Annotated and Associated Data

Phylogenetic trees and data are often stored in incompatible and inconsistent formats. The outputs of software tools that contain trees with analysis findings are often not compatible with each other, making it hard to integrate the results of different analyses in a comparative study. The treeio package is designed to connect phylogenetic tree input and output. It supports extracting phylogenetic trees as well as the outputs of commonly used analytical software. It can link external data to phylogenies and merge tree data obtained from different sources, enabling analyses of phylogeny-associated data from different disciplines in an evolutionary context. Treeio also supports export of a phylogenetic tree with heterogeneous-associated data to a single tree file, including BEAST compatible NEXUS and jtree formats; these facilitate data sharing as well as file format conversion for downstream analysis. The treeio package is designed to work with the tidytree and ggtree packages. Tree data can be processed using the tidy interface with tidytree and visualized by ggtree. The treeio package is released within the Bioconductor and rOpenSci projects. It is available at https://www.bioconductor.org/packages/treeio/.

Morphology and development of the Portuguese man of war, Physalia physalis

The Portuguese man of war, Physalia physalis, is one of the most conspicuous, but poorly understood members of the pleuston, a community of organisms that occupy a habitat at the sea-air interface. Physalia physalis is a siphonophore that uses a gas-filled float as a sail to catch the wind. The development, morphology, and colony organization of P. physalis is very different from all other siphonophores. Here, we look at live and fixed larval and juvenile specimens, and use optical projection tomography to build on existing knowledge about the morphology and development of this species. We also propose a framework for homologizing the axes with other siphonophores, and also suggest that the tentacle bearing zooids should be called tentacular palpons. Previous descriptions of P. physalis larvae, especially descriptions of budding order, were often framed with the mature colony in mind. However, we use the simpler organization of larvae and the juvenile specimens to inform our understanding of the morphology, budding order, and colony organization in the mature specimen. Finally, we review what is known about the ecology and lifecycle of P. physalis.

We are not so special

New sequence data from choanoflagellates improves our understanding of the genetic changes that occurred along the branch of the evolutionary tree that gave rise to animals.

Improved phylogenetic resolution within Siphonophora (Cnidaria) with implications for trait evolution

Siphonophores are a diverse group of hydrozoans (Cnidaria) that are found at most depths of the ocean - from the surface, like the familiar Portuguese man of war, to the deep sea. They play important roles in ocean ecosystems, and are among the most abundant gelatinous predators. A previous phylogenetic study based on two ribosomal RNA genes provided insight into the internal relationships between major siphonophore groups. There was, however, little support for many deep relationships within the clade Codonophora. Here, we present a new siphonophore phylogeny based on new transcriptome data from 29 siphonophore species analyzed in combination with 14 publicly available genomic and transcriptomic datasets. We use this new phylogeny to reconstruct several traits that are central to siphonophore biology, including sexual system (monoecy vs. dioecy), gain and loss of zooid types, life history traits, and habitat. The phylogenetic relationships in this study are largely consistent with the previous phylogeny, but we find strong support for new clades within Codonophora that were previously unresolved. These results have important implications for trait evolution within Siphonophora, including favoring the hypothesis that monoecy arose at least twice.

Description of Tottonophyes enigmatica gen. nov., sp. nov. (Hydrozoa, Siphonophora, Calycophorae), with a reappraisal of the function and homology of nectophoral canals

A new species of calycophoran siphonophore, Tottonophyes enigmatica gen. nov, sp. nov., is described. It has a unique combination of traits, some shared with prayomorphs (including two rounded nectophores) and some with clausophyid diphyomorphs (the nectophores are dissimilar, with one slightly larger and slightly to the anterior of the other, and both possess a somatocyst). Molecular phylogenetic analyses indicate that the new species is the sister group to all other diphyomorphs. A new family, Tottonophyidae, is established for it. Its phylogenetic position and distinct morphology help clarify diphyomorph evolution. The function and homology of the nectophoral canals and somatocyst is also re-examined and further clarification is given to their nomenclature.

Animal Evolution: Are Phyla Real?

A recent study compares developmental gene expression among very distantly related animal species and interprets the results as providing biological justification for phyla. Several problems in the design and interpretation of the study call this conclusion into question.

Pairwise comparisons across species are problematic when analyzing functional genomic data

There is considerable interest in comparing functional genomic data across species. One goal of such work is to provide an integrated understanding of genome and phenotype evolution. Most comparative functional genomic studies have relied on multiple pairwise comparisons between species, an approach that does not incorporate information about the evolutionary relationships among species. The statistical problems that arise from not considering these relationships can lead pairwise approaches to the wrong conclusions and are a missed opportunity to learn about biology that can only be understood in an explicit phylogenetic context. Here, we examine two recently published studies that compare gene expression across species with pairwise methods, and find reason to question the original conclusions of both. One study interpreted pairwise comparisons of gene expression as support for the ortholog conjecture, the hypothesis that orthologs tend to have more similar attributes (expression in this case) than paralogs. The other study interpreted pairwise comparisons of embryonic gene expression across distantly related animals as evidence for a distinct evolutionary process that gave rise to phyla. In each study, distinct patterns of pairwise similarity among species were originally interpreted as evidence of particular evolutionary processes, but instead, we find that they reflect species relationships. These reanalyses concretely show the inadequacy of pairwise comparisons for analyzing functional genomic data across species. It will be critical to adopt phylogenetic comparative methods in future functional genomic work. Fortunately, phylogenetic comparative biology is also a rapidly advancing field with many methods that can be directly applied to functional genomic data.

Indoles induce metamorphosis in a broad diversity of jellyfish, but not in a crown jelly (Coronatae)

Many animals go through one or more metamorphoses during their lives, however, the molecular underpinnings of metamorphosis across diverse species are not well understood. Medusozoa (Cnidaria) is a clade of animals with complex life cycles, these life cycles can include a polyp stage that metamorphoses into a medusa (jellyfish). Medusae are produced through a variety of different developmental mechanisms-in some species polyps bud medusae (Hydrozoa), in others medusae are formed through polyp fission (Scyphozoa), while in others medusae are formed through direct transformation of the polyp (Cubozoa). To better understand the molecular mechanisms that may coordinate these different forms of metamorphosis, we tested two compounds first identified to induce metamorphosis in the moon jellyfish Aurelia aurita (indomethacin and 5-methoxy-2-methylindole) on a broad diversity of medusozoan polyps. We discovered that indole-containing compounds trigger metamorphosis across a broad diversity of species. All tested discomedusan polyps metamorphosed in the presence of both compounds, including species representatives of several major lineages within the clade (Pelagiidae, Cyaneidae, both clades of Rhizostomeae). In a cubozoan, low levels of 5-methoxy-2-methylindole reliably induced complete and healthy metamorphosis. In contrast, neither compound induced medusa metamorphosis in a coronate scyphozoan, or medusa production in either hydrozoan tested. Our results support the hypothesis that metamorphosis is mediated by a conserved induction pathway within discomedusan scyphozoans, and possibly cubozoans. However, failure of these compounds to induce metamorphosis in a coronate suggests this induction mechanism may have been lost in this clade, or is convergent between Scyphozoa and Cubozoa.

Phylogenomic Analyses Support Traditional Relationships within Cnidaria

Cnidaria, the sister group to Bilateria, is a highly diverse group of animals in terms of morphology, lifecycles, ecology, and development. How this diversity originated and evolved is not well understood because phylogenetic relationships among major cnidarian lineages are unclear, and recent studies present contrasting phylogenetic hypotheses. Here, we use transcriptome data from 15 newly-sequenced species in combination with 26 publicly available genomes and transcriptomes to assess phylogenetic relationships among major cnidarian lineages. Phylogenetic analyses using different partition schemes and models of molecular evolution, as well as topology tests for alternative phylogenetic relationships, support the monophyly of Medusozoa, Anthozoa, Octocorallia, Hydrozoa, and a clade consisting of Staurozoa, Cubozoa, and Scyphozoa. Support for the monophyly of Hexacorallia is weak due to the equivocal position of Ceriantharia. Taken together, these results further resolve deep cnidarian relationships, largely support traditional phylogenetic views on relationships, and provide a historical framework for studying the evolutionary processes involved in one of the most ancient animal radiations.

A phylogenetic backbone for Bivalvia: an RNA-seq approach

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

The evolution of animal genomes

Genome sequences are now available for hundreds of species sampled across the animal phylogeny, bringing key features of animal genome evolution into sharper focus. The field of animal evolutionary genomics has focused on identifying and classifying the diversity genomic features, reconstructing the history of evolutionary changes in animal genomes, and testing hypotheses about the evolutionary relationships of animals. The grand challenges moving forward are to connect evolutionary changes in genomes with particular evolutionary changes in phenotypes, and to determine which changes are driven by selection. This will require far greater genome sampling both across and within species, extensive phenotype data, a well resolved animal phylogeny, and advances in comparative methods.

Re-evaluation of characters in Apolemiidae (Siphonophora), with description of two new species from Monterey Bay, California

Siphonophores are polymorphic planktonic marine Cnidarians. The family Apolemiidae is sister to all other species of physonect and calycophoran siphonophores. Although this enigmatic group arguably includes the longest animal species on the planet, their colony-level organization and growth patterns are not well understood. Here we describe two new apolemiid species: Apolemia lanosa sp. nov. and A. rubriversa sp. nov. We provide detailed descriptions of zooid budding and the organization of mature zooids within the siphosome. Our findings reveal that at least two distinct general patterns of siphosomal organization are found in different Apolemia species. In the first pattern, dispersed organization, zooids independently attach directly to the siphosomal stem. In the second pattern, pedunculate organization, only the gastrozooid is attached directly to the stem, and the other zooids of the cormidium branch from its peduncle. This diversity within Apolemia indicates that fundamental aspects of zooid budding and organization are homoplastic within Siphonophora, as both patterns are also found in other siphonophores. The observations presented here greatly clarify the interpretation of diagnostic characters within Apolemiidae, bear on the status of the three previously described species, provide critical detail for understanding the diversity of colony-level organization in siphonophores, and establish a foundation for the description of additional apolemiid species.

Fluorescent proteins function as a prey attractant: experimental evidence from the hydromedusa Olindias formosus and other marine organisms

Although proteins in the green fluorescent protein family (GFPs) have been discovered in a wide array of taxa, their ecological functions in these organisms remain unclear. Many hypothesized roles are related to modifying bioluminescence spectra or modulating the light regime for algal symbionts, but these do not explain the presence of GFPs in animals that are non-luminous and non-symbiotic. Other hypothesized functions are unrelated to the visual signals themselves, including stress responses and antioxidant roles, but these cannot explain the localization of fluorescence in particular structures on the animals. Here we tested the hypothesis that fluorescence might serve to attract prey. In laboratory experiments, the predator was the hydromedusa Olindias formosus (previously known as O. formosa), which has fluorescent and pigmented patches on the tips of its tentacles. The prey, juvenile rockfishes in the genus Sebastes, were significantly more attracted (P<1×10⁻⁵) to the medusa's tentacles under lighting conditions where fluorescence was excited and tentacle tips were visible above the background. The fish did not respond significantly when treatments did not include fluorescent structures or took place under yellow or white lights, which did not generate fluorescence visible above the ambient light. Furthermore, underwater observations of the behavior of fishes when presented with a brightly illuminated point showed a strong attraction to this visual stimulus. In situ observations also provided evidence for fluorescent lures as supernormal stimuli in several other marine animals, including the siphonophore Rhizophysa eysenhardti. Our results support the idea that fluorescent structures can serve as prey attractants, thus providing a potential function for GFPs and other fluorescent proteins in a diverse range of organisms.

Automation and Evaluation of the SOWH Test with SOWHAT

The Swofford–Olsen–Waddell–Hillis (SOWH) test evaluates statistical support for incongruent phylogenetic topologies. It is commonly applied to determine if the maximum likelihood tree in a phylogenetic analysis is significantly different than an alternative hypothesis. The SOWH test compares the observed difference in log-likelihood between two topologies to a null distribution of differences in log-likelihood generated by parametric resampling. The test is a well-established phylogenetic method for topology testing, but it is sensitive to model misspecification, it is computationally burdensome to perform, and its implementation requires the investigator to make several decisions that each have the potential to affect the outcome of the test. We analyzed the effects of multiple factors using seven data sets to which the SOWH test was previously applied. These factors include a number of sample replicates, likelihood software, the introduction of gaps to simulated data, the use of distinct models of evolution for data simulation and likelihood inference, and a suggested test correction wherein an unresolved “zero-constrained” tree is used to simulate sequence data. To facilitate these analyses and future applications of the SOWH test, we wrote SOWHAT, a program that automates the SOWH test. We find that inadequate bootstrap sampling can change the outcome of the SOWH test. The results also show that using a zero-constrained tree for data simulation can result in a wider null distribution and higher p-values, but does not change the outcome of the SOWH test for most of the data sets tested here. These results will help others implement and evaluate the SOWH test and allow us to provide recommendations for future applications of the SOWH test. SOWHAT is available for download from https://github.com/josephryan/SOWHAT.

Phylogenomic analyses of deep gastropod relationships reject Orthogastropoda

Gastropods are a highly diverse clade of molluscs that includes many familiar animals, such as limpets, snails, slugs and sea slugs. It is one of the most abundant groups of animals in the sea and the only molluscan lineage that has successfully colonized land. Yet the relationships among and within its constituent clades have remained in flux for over a century of morphological, anatomical and molecular study. Here, we re-evaluate gastropod phylogenetic relationships by collecting new transcriptome data for 40 species and analysing them in combination with publicly available genomes and transcriptomes. Our datasets include all five main gastropod clades: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia. We use two different methods to assign orthology, subsample each of these matrices into three increasingly dense subsets, and analyse all six of these supermatrices with two different models of molecular evolution. All 12 analyses yield the same unrooted network connecting the five major gastropod lineages. This reduces deep gastropod phylogeny to three alternative rooting hypotheses. These results reject the prevalent hypothesis of gastropod phylogeny, Orthogastropoda. Our dated tree is congruent with a possible end-Permian recovery of some gastropod clades, namely Caenogastropoda and some Heterobranchia subclades.

The histology of Nanomia bijuga (Hydrozoa: Siphonophora)

The siphonophore Nanomia bijuga is a pelagic hydrozoan (Cnidaria) with complex morphological organization. Each siphonophore is made up of many asexually produced, genetically identical zooids that are functionally specialized and morphologically distinct. These zooids predominantly arise by budding in two growth zones, and are arranged in precise patterns. This study describes the cellular anatomy of several zooid types, the stem, and the gas‐filled float, called the pneumatophore. The distribution of cellular morphologies across zooid types enhances our understanding of zooid function. The unique absorptive cells in the palpon, for example, indicate specialized intracellular digestive processing in this zooid type. Though cnidarians are usually thought of as mono‐epithelial, we characterize at least two cellular populations in this species which are not connected to a basement membrane. This work provides a greater understanding of epithelial diversity within the cnidarians, and will be a foundation for future studies on N. bijuga, including functional assays and gene expression analyses. J. Exp. Zool. (Mol. Dev. Evol.) 324B:435–449, 2015. © 2015 The Authors. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution Published by Wiley Periodicals, Inc.

Stem cells in Nanomia bijuga (Siphonophora), a colonial animal with localized growth zones

BACKGROUND

Siphonophores (Hydrozoa) have unparalleled colony-level complexity, precision of colony organization, and functional specialization between zooids (i.e., the units that make up colonies). Previous work has shown that, unlike other colonial animals, most growth in siphonophores is restricted to one or two well-defined growth zones that are the sites of both elongation and zooid budding. It remained unknown, however, how this unique colony growth and development is realized at the cellular level.

RESULTS

To understand the colony-level growth and development of siphonophores at the cellular level, we characterize the distribution of proliferating cells and interstitial stem cells (i-cells) in the siphonophore Nanomia bijuga. Within the colony, we find evidence that i-cells are present at the tip of the horn, the structure within the growth zone that gives rise to new zooids. Co-localized gene expression of vasa-1, pl10, piwi, nanos-1, and nanos-2 suggests that i-cells persist in the youngest zooid buds and that i-cells become progressively restricted to specific regions within the zooids until they are mostly absent from the oldest zooids. The examined genes remain expressed in gametogenic regions. No evidence for i-cells is found in the stem between maturing zooids. Domains of high cell proliferation include regions where the examined genes are expressed, but also include some areas in which the examined genes were not expressed such as the stem within the growth zones. Cell proliferation in regions devoid of vasa-1, pl10, piwi, nanos-1, and nanos-2 expression indicates the presence of mitotically active epithelial cell lineages and, potentially, progenitor cell populations.

CONCLUSIONS

We provide the first evidence for i-cells in a siphonophore. Our findings suggest maintenance of i-cell populations at the sites of growth zones and that these sites are the main source of i-cells. This restriction of stem cells to particular regions in the colony, in combination with localized budding and spatial patterning during pro-bud subdivision, may play a major role in facilitating the precision of siphonophore growth. Spatially restricted maintenance of i-cells in mature zooids and absence of i-cells along the stem may explain the reduced developmental plasticity in older parts of the colony.

Par system components are asymmetrically localized in ectodermal epithelia, but not during early development in the sea anemone Nematostella vectensis

Background

The evolutionary origins of cell polarity in metazoan embryos are unclear. In most bilaterian animals, embryonic and cell polarity are set up during embryogenesis with the same molecules being utilized to regulate tissue polarity at different life stages. Atypical protein kinase C (aPKC), lethal giant larvae (Lgl), and Partitioning-defective (Par) proteins are conserved components of cellular polarization, and their role in establishing embryonic asymmetry and tissue polarity have been widely studied in model bilaterian groups. However, the deployment and role of these proteins in animals outside Bilateria has not been studied. We address this by characterizing the localization of different components of the Par system during early development of the sea anemone Nematostella vectensis, a member of the clade Cnidaria, the sister group to bilaterian animals.

Results

Immunostaining using specific N. vectensis antibodies and the overexpression of mRNA-reporter constructs show that components of the N. vectensis Par system (NvPar-1, NvPar-3, NvPar-6, NvaPKC, and NvLgl) distribute throughout the microtubule cytoskeleton of eggs and early embryos without clear polarization along any embryonic axis. However, they become asymmetrically distributed at later stages, when the embryo forms an ectodermal epithelial layer. NvLgl and NvPar-1 localize in the basolateral cortex, and NvaPKC, NvPar-6, and NvPar-3 at the apical zone of the cell in a manner seen in bilaterian animals.

Conclusions

The cnidarian N. vectensis exhibits clear polarity at all stages of early embryonic development, which appears to be established independent of the Par system reported in many bilaterian embryos. However, in N. vectensis, using multiple immunohistochemical and fluorescently labeled markers in vivo, components of this system are deployed to organize epithelial cell polarity at later stages of development. This suggests that Par system proteins were co-opted to organize early embryonic cell polarity at the base of the Bilateria and that, therefore, different molecular mechanisms operate in early cnidarian embryogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0014-6) contains supplementary material, which is available to authorized users.

Comparative muscle development of scyphozoan jellyfish with simple and complex life cycles

Background

Simple life cycles arise from complex life cycles when one or more developmental stages are lost. This raises a fundamental question - how can an intermediate stage, such as a larva, be removed, and development still produce a normal adult? To address this question, we examined the development in several species of pelagiid jellyfish. Most members of Pelagiidae have a complex life cycle with a sessile polyp that gives rise to ephyrae (juvenile medusae); but one species within Pelagiidae, Pelagia noctiluca, spends its whole life in the water column, developing from a larva directly into an ephyra. In many complex life cycles, adult features develop from cell populations that remain quiescent in larvae, and this is known as life cycle compartmentalization and may facilitate the evolution of direct life cycles. A second type of metamorphic processes, known as remodeling, occurs when adult features are formed through modification of already differentiated larval structures. We examined muscle morphology to determine which of these alternatives may be present in Pelagiidae.

Results

We first examined the structure and development of polyp and ephyra musculature in Chrysaora quinquecirrha, a close relative of P. noctiluca with a complex life cycle. Using phallotoxin staining and confocal microscopy, we verified that polyps have four to six cord muscles that persist in strobilae and discovered that cord muscles is physically separated from ephyra muscle. When cord muscle is removed from ephyra segments, normal ephyra muscle still develops. This suggests that polyp cord muscle is not necessary for ephyra muscle formation. We also found no evidence of polyp-like muscle in P. noctiluca. In both species, we discovered that ephyra muscle arises de novo in a similar manner, regardless of the life cycle.

Conclusions

The separate origins of polyp and ephyra muscle in C. quinquecirrha and the absence of polyp-like muscle in P. noctiluca suggest that polyp muscle is not remodeled to form ephyra muscle in Pelagiidae. Life cycle stages in Scyphozoa may instead be compartmentalized. Because polyp muscle is not directly remodeled, this may have facilitated the loss of the polyp stage in the evolution of P. noctiluca.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0005-7) contains supplementary material, which is available to authorized users.

The genome of the ctenophore Mnemiopsis leidyi and its implications for cell type evolution

An understanding of ctenophore biology is critical for reconstructing events that occurred early in animal evolution. Toward this goal, we have sequenced, assembled, and annotated the genome of the ctenophore Mnemiopsis leidyi. Our phylogenomic analyses of both amino acid positions and gene content suggest that ctenophores rather than sponges are the sister lineage to all other animals. Mnemiopsis lacks many of the genes found in bilaterian mesodermal cell types, suggesting that these cell types evolved independently. The set of neural genes in Mnemiopsis is similar to that of sponges, indicating that sponges may have lost a nervous system. These results present a newly supported view of early animal evolution that accounts for major losses and/or gains of sophisticated cell types, including nerve and muscle cells.

Agalma: an automated phylogenomics workflow

BACKGROUND

In the past decade, transcriptome data have become an important component of many phylogenetic studies. They are a cost-effective source of protein-coding gene sequences, and have helped projects grow from a few genes to hundreds or thousands of genes. Phylogenetic studies now regularly include genes from newly sequenced transcriptomes, as well as publicly available transcriptomes and genomes. Implementing such a phylogenomic study, however, is computationally intensive, requires the coordinated use of many complex software tools, and includes multiple steps for which no published tools exist. Phylogenomic studies have therefore been manual or semiautomated. In addition to taking considerable user time, this makes phylogenomic analyses difficult to reproduce, compare, and extend. In addition, methodological improvements made in the context of one study often cannot be easily applied and evaluated in the context of other studies.

RESULTS

We present Agalma, an automated tool that constructs matrices for phylogenomic analyses. The user provides raw Illumina transcriptome data, and Agalma produces annotated assemblies, aligned gene sequence matrices, a preliminary phylogeny, and detailed diagnostics that allow the investigator to make extensive assessments of intermediate analysis steps and the final results. Sequences from other sources, such as externally assembled genomes and transcriptomes, can also be incorporated in the analyses. Agalma is built on the BioLite bioinformatics framework, which tracks provenance, profiles processor and memory use, records diagnostics, manages metadata, installs dependencies, logs version numbers and calls to external programs, and enables rich HTML reports for all stages of the analysis. Agalma includes a small test data set and a built-in test analysis of these data. In addition to describing Agalma, we here present a sample analysis of a larger seven-taxon data set. Agalma is available for download at https://bitbucket.org/caseywdunn/agalma.

CONCLUSIONS

Agalma allows complex phylogenomic analyses to be implemented and described unambiguously as a series of high-level commands. This will enable phylogenomic studies to be readily reproduced, modified, and extended. Agalma also facilitates methods development by providing a complete modular workflow, bundled with test data, that will allow further optimization of each step in the context of a full phylogenomic analysis.

Toward a statistically explicit understanding of de novo sequence assembly

MOTIVATION

Draft de novo genome assemblies are now available for many organisms. These assemblies are point estimates of the true genome sequences. Each is a specific hypothesis, drawn from among many alternative hypotheses, of the sequence of a genome. Assembly uncertainty, the inability to distinguish between multiple alternative assembly hypotheses, can be due to real variation between copies of the genome in the sample, errors and ambiguities in the sequenced data and assumptions and heuristics of the assemblers. Most assemblers select a single assembly according to ad hoc criteria, and do not yet report and quantify the uncertainty of their outputs. Those assemblers that do report uncertainty take different approaches to describing multiple assembly hypotheses and the support for each.

RESULTS

Here we review and examine the problem of representing and measuring uncertainty in assemblies. A promising recent development is the implementation of assemblers that are built according to explicit statistical models. Some new assembly methods, for example, estimate and maximize assembly likelihood. These advances, combined with technical advances in the representation of alternative assembly hypotheses, will lead to a more complete and biologically relevant understanding of assembly uncertainty. This will in turn facilitate the interpretation of downstream analyses and tests of specific biological hypotheses.

Phylogenetic analysis of gene expression

Phylogenetic analyses of gene expression have great potential for addressing a wide range of questions. These analyses will, for example, identify genes that have evolutionary shifts in expression that are correlated with evolutionary changes in morphological, physiological, and developmental characters of interest. This will provide entirely new opportunities to identify genes related to particular phenotypes. There are, however, 3 key challenges that must be addressed for such studies to realize their potential. First, data on gene expression must be measured from multiple species, some of which may be field-collected, and parameterized in such a way that they can be compared across species. Second, it will be necessary to develop comparative phylogenetic methods suitable for large multidimensional datasets. In most phylogenetic comparative studies to date, the number n of independent observations (independent contrasts) has been greater than the number p of variables (characters). The behavior of comparative methods for these classic problems is now well understood under a wide variety of conditions. In studies of gene expression, and in studies based on other high-throughput tools, the number n of samples is dwarfed by the number p of variables. The estimated covariance matrices will be singular, complicating their analysis and interpretation, and prone to spurious results. Third, new approaches are needed to investigate the expression of the many genes whose phylogenies are not congruent with species phylogenies due to gene loss, gene duplication, and incomplete lineage sorting. Here we outline general considerations of project design for phylogenetic analyses of gene expression and suggest solutions to these three categories of challenges. These topics are relevant to high-throughput phenotypic data well beyond gene expression.

Differential gene expression in the siphonophore Nanomia bijuga (Cnidaria) assessed with multiple next-generation sequencing workflows

We investigated differential gene expression between functionally specialized feeding polyps and swimming medusae in the siphonophore Nanomia bijuga (Cnidaria) with a hybrid long-read/short-read sequencing strategy. We assembled a set of partial gene reference sequences from long-read data (Roche 454), and generated short-read sequences from replicated tissue samples that were mapped to the references to quantify expression. We collected and compared expression data with three short-read expression workflows that differ in sample preparation, sequencing technology, and mapping tools. These workflows were Illumina mRNA-Seq, which generates sequence reads from random locations along each transcript, and two tag-based approaches, SOLiD SAGE and Helicos DGE, which generate reads from particular tag sites. Differences in expression results across workflows were mostly due to the differential impact of missing data in the partial reference sequences. When all 454-derived gene reference sequences were considered, Illumina mRNA-Seq detected more than twice as many differentially expressed (DE) reference sequences as the tag-based workflows. This discrepancy was largely due to missing tag sites in the partial reference that led to false negatives in the tag-based workflows. When only the subset of reference sequences that unambiguously have tag sites was considered, we found broad congruence across workflows, and they all identified a similar set of DE sequences. Our results are promising in several regards for gene expression studies in non-model organisms. First, we demonstrate that a hybrid long-read/short-read sequencing strategy is an effective way to collect gene expression data when an annotated genome sequence is not available. Second, our replicated sampling indicates that expression profiles are highly consistent across field-collected animals in this case. Third, the impacts of partial reference sequences on the ability to detect DE can be mitigated through workflow choice and deeper reference sequencing.

Assessing the root of bilaterian animals with scalable phylogenomic methods

A clear picture of animal relationships is a prerequisite to understand how the morphological and ecological diversity of animals evolved over time. Among others, the placement of the acoelomorph flatworms, Acoela and Nemertodermatida, has fundamental implications for the origin and evolution of various animal organ systems. Their position, however, has been inconsistent in phylogenetic studies using one or several genes. Furthermore, Acoela has been among the least stable taxa in recent animal phylogenomic analyses, which simultaneously examine many genes from many species, while Nemertodermatida has not been sampled in any phylogenomic study. New sequence data are presented here from organisms targeted for their instability or lack of representation in prior analyses, and are analysed in combination with other publicly available data. We also designed new automated explicit methods for identifying and selecting common genes across different species, and developed highly optimized supercomputing tools to reconstruct relationships from gene sequences. The results of the work corroborate several recently established findings about animal relationships and provide new support for the placement of other groups. These new data and methods strongly uphold previous suggestions that Acoelomorpha is sister clade to all other bilaterian animals, find diminishing evidence for the placement of the enigmatic Xenoturbella within Deuterostomia, and place Cycliophora with Entoprocta and Ectoprocta. The work highlights the implications that these arrangements have for metazoan evolution and permits a clearer picture of ancestral morphologies and life histories in the deep past.

Broad phylogenomic sampling improves resolution of the animal tree of life

Long-held ideas regarding the evolutionary relationships among animals have recently been upended by sometimes controversial hypotheses based largely on insights from molecular data. These new hypotheses include a clade of moulting animals (Ecdysozoa) and the close relationship of the lophophorates to molluscs and annelids (Lophotrochozoa). Many relationships remain disputed, including those that are required to polarize key features of character evolution, and support for deep nodes is often low. Phylogenomic approaches, which use data from many genes, have shown promise for resolving deep animal relationships, but are hindered by a lack of data from many important groups. Here we report a total of 39.9 Mb of expressed sequence tags from 29 animals belonging to 21 phyla, including 11 phyla previously lacking genomic or expressed-sequence-tag data. Analysed in combination with existing sequences, our data reinforce several previously identified clades that split deeply in the animal tree (including Protostomia, Ecdysozoa and Lophotrochozoa), unambiguously resolve multiple long-standing issues for which there was strong conflicting support in earlier studies with less data (such as velvet worms rather than tardigrades as the sister group of arthropods), and provide molecular support for the monophyly of molluscs, a group long recognized by morphologists. In addition, we find strong support for several new hypotheses. These include a clade that unites annelids (including sipunculans and echiurans) with nemerteans, phoronids and brachiopods, molluscs as sister to that assemblage, and the placement of ctenophores as the earliest diverging extant multicellular animals. A single origin of spiral cleavage (with subsequent losses) is inferred from well-supported nodes. Many relationships between a stable subset of taxa find strong support, and a diminishing number of lineages remain recalcitrant to placement on the tree.

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.

The evolution of colony-level development in the Siphonophora (Cnidaria:Hydrozoa)

Evolutionary developmental biology has focused almost exclusively on multicellular organisms, but there are other relevant levels of biological organization that have remained largely neglected. Animal colonies are made up of multiple physiologically integrated and genetically identical units called zooids that are each homologous to solitary, free-living animals. Siphonophores, a group of pelagic hydrozoans (Cnidaria), have the most complex colony-level organization of all animals. Here the colony-level development of five siphonophore species, strategically sampled across the siphonophore phylogeny, is described from specimens collected using deep-sea submersibles and by self-contained underwater breathing apparatus diving. These species include three cystonects, Bathyphysa sibogae, Rhizophysa filiformis, and Rhizophysa eysenhardti, and two "physonects", Agalma elegans and Nanomia bijuga. These data, together with previous findings, are analyzed in a phylogenetic framework to reconstruct key features of the history of colony-level organization and development in the Siphonophora. It is shown that gonodendra and gastrozooids of the examined cystonects arise as independent buds directly on the stem, whereas probud subdivision (the origin of feeding, reproductive, and other zooids from a single bud) is a synapomorphy of the Codonophora. The origin of probud subdivision is associated with the origin of cormidia as integrated units of colony organization, and may have allowed for greater morphological and ecological diversification in the Codonophora relative to the Cystonectae. It is also found that symmetry is labile in siphonophores, with multiple gains and/or losses of directional asymmetry in the group. This descriptive work will enable future mechanistic and molecular studies of colony-level development in the siphonophores.