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

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.

Establishing Bilateral Symmetry in Hydrozoan Planula Larvae, a Review of Siphonophore Early Development

Siphonophores are colonial hydrozoans, characterized by complex colony organization and unparalleled zooid functional specialization. Recent genomic studies have offered an evolutionary perspective on how this morphological complexity arose, but a molecular characterization of symmetry breaking in siphonophore embryonic development is still largely missing. Here, bringing together historical data on early development with new immunohistochemical data, we review the diversity of developmental trajectories that lead to the formation of bilaterally symmetric planula larvae in siphonophores. Embryonic development, up to the planula stage, is remarkably similar across siphonophore phylogeny. Then, with the appearance of the lateral endodermal thickening (= ventral endoderm), larval development diverges between taxa, differing in the location and patterning of the primary buds, chronology of budding, establishment of growth zones, and retention of larval zooids. Our work also uncovers a number of open questions in siphonophore development, including homology of different zooids, mechanisms underlying formation and maintenance of spatially restricted growth zone(s), and molecular factors establishing a secondary dorsal-ventral axis in planulae. By discussing siphonophore development and body axes within the broader cnidarian context, we then set the framework for future work on siphonophores, which is finally achievable with the advent of culturing methods.

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.

Differential Gene Expression Between Polymorphic Zooids of the Marine Bryozoan Bugulina stolonifera

Bryozoans are a diverse phylum of marine and freshwater colonial invertebrates containing approximately 6,300 described living species. Bryozoans grow by budding new physiologically connected colony members (zooids) from a founding individual that forms from a metamorphosed larva. In some species these zooids come in different shapes and sizes and are specialized to serve different tasks within the colony. A complex interaction of genotype, environment, and developmental pathway shapes zooid fate, however, the specific mechanisms underlying the establishment of this division of labor remain unknown. Here, the first characterization of differential gene expression between polymorphic zooids of a bryozoan colony is presented. The development of different zooid types of lab-cultured Bugulina stolonifera colonies including feeding autozooids, avicularia (derived non-feeding zooids that are homologous to feeding autozooids but shaped like a bird’s beak), and rhizoids (a branching network of non-feeding anchoring zooids) was explored using RNA sequencing, de novo transcriptome assembly, and differential gene expression analyses. High throughput sequencing of cDNA libraries yielded an average of 14.9 ± 1.3 (SE) million high-quality paired-end reads per sample. Data for the first de novo transcriptome assemblies of B. stolonifera and the first characterization of genes involved in the formation and maintenance of zooid types within a bryozoan colony are presented. In a comparison between autozooid and avicularium tissues, 1,097 significant differentially expressed genes were uncovered. This work provides a much-needed foundation for understanding the mechanisms involved in the development of polymorphic zooids and the establishment of division of labor in bryozoans.

Coloniality, clonality, and modularity in animals: The elephant in the room

Nearly half of the animal phyla contain species that propagate asexually via agametic reproduction, often forming colonies of genetically identical modules, that is, ramets, zooids, or polyps. Clonal reproduction, colony formation, and modular organization have important consequences for many aspects of organismal biology. Theories in ecology, evolution, and development are often based on unitary and, mainly, strictly sexually reproducing organisms, and though colonial animals dominate many marine ecosystems and habitats, recognized concepts for the study of clonal species are often lacking. In this review, we present an overview of the study of colonial and clonal animals, from the historic interests in this subject to modern research in a range of topics, including immunology, stem cell biology, aging, biogeography, and ecology. We attempt to portray the fundamental questions lying behind the biology of colonial animals, focusing on how colonial animals challenge several dogmas in biology as well as the remaining puzzles still to be answered, of which there are many.

De Novo Assembly of the Transcriptome of Turritopsis, a Jellyfish that Repeatedly Rejuvenates

In most animals, aging is an irreversible process; however the species Turritopsis sp. has been observed to undergo a rejuvenation process as many as 14 times. In the present study, we used multiplexed RNA libraries to obtain the transcriptome from four developmental stages (St) of Turritopsis sp., including (I) immature medusa, (II) dumpling, (III) dumpling with a short stolon, and (IV) polyp, which had recently rejuvenated. A total of 4.02 billion paired-end reads were assembled de novo, yielding 90,327 contigs. Our analyses revealed that significant blast hits were recovered for 74% of the assembled contigs, and 19% were successfully annotated with gene ontology (GO) terms. A BLAST search demonstrated that 32% of the contigs were most similar to Hydra vulgarissequences. Raw reads from each sample were mapped against the contigs to find St-specific genes. This represents the first comprehensive set of de novo transcriptome data for this species, which may provide clues toward a better understanding of cyclical rejuvenation in multicellular animals.

PARRoT- a homology-based strategy to quantify and compare RNA-sequencing from non-model organisms

Background

Next-generation sequencing promises the de novo genomic and transcriptomic analysis of samples of interests. However, there are only a few organisms having reference genomic sequences and even fewer having well-defined or curated annotations. For transcriptome studies focusing on organisms lacking proper reference genomes, the common strategy is de novo assembly followed by functional annotation. However, things become even more complicated when multiple transcriptomes are compared.

Results

Here, we propose a new analysis strategy and quantification methods for quantifying expression level which not only generate a virtual reference from sequencing data, but also provide comparisons between transcriptomes. First, all reads from the transcriptome datasets are pooled together for de novo assembly. The assembled contigs are searched against NCBI NR databases to find potential homolog sequences. Based on the searched result, a set of virtual transcripts are generated and served as a reference transcriptome. By using the same reference, normalized quantification values including RC (read counts), eRPKM (estimated RPKM) and eTPM (estimated TPM) can be obtained that are comparable across transcriptome datasets. In order to demonstrate the feasibility of our strategy, we implement it in the web service PARRoT. PARRoT stands for Pipeline for Analyzing RNA Reads of Transcriptomes. It analyzes gene expression profiles for two transcriptome sequencing datasets. For better understanding of the biological meaning from the comparison among transcriptomes, PARRoT further provides linkage between these virtual transcripts and their potential function through showing best hits in SwissProt, NR database, assigning GO terms. Our demo datasets showed that PARRoT can analyze two paired-end transcriptomic datasets of approximately 100 million reads within just three hours.

Conclusions

In this study, we proposed and implemented a strategy to analyze transcriptomes from non-reference organisms which offers the opportunity to quantify and compare transcriptome profiles through a homolog based virtual transcriptome reference. By using the homolog based reference, our strategy effectively avoids the problems that may cause from inconsistencies among transcriptomes. This strategy will shed lights on the field of comparative genomics for non-model organism. We have implemented PARRoT as a web service which is freely available at http://parrot.cgu.edu.tw.

Tissue-Specific Venom Composition and Differential Gene Expression in Sea Anemones

Cnidarians represent one of the few groups of venomous animals that lack a centralized venom transmission system. Instead, they are equipped with stinging capsules collectively known as nematocysts. Nematocysts vary in abundance and type across different tissues; however, the venom composition in most species remains unknown. Depending on the tissue type, the venom composition in sea anemones may be vital for predation, defense, or digestion. Using a tissue-specific RNA-seq approach, we characterize the venom assemblage in the tentacles, mesenterial filaments, and column for three species of sea anemone (Anemonia sulcata, Heteractis crispa, and Megalactis griffithsi). These taxa vary with regard to inferred venom potency, symbiont abundance, and nematocyst diversity. We show that there is significant variation in abundance of toxin-like genes across tissues and species. Although the cumulative toxin abundance for the column was consistently the lowest, contributions to the overall toxin assemblage varied considerably among tissues for different toxin types. Our gene ontology (GO) analyses also show sharp contrasts between conserved GO groups emerging from whole transcriptome analysis and tissue-specific expression among GO groups in our differential expression analysis. This study provides a framework for future characterization of tissue-specific venom and other functionally important genes in this lineage of simple bodied animals.

Jellyfish Bioactive Compounds: Methods for Wet-Lab Work

The study of bioactive compounds from marine animals has provided, over time, an endless source of interesting molecules. Jellyfish are commonly targets of study due to their toxic proteins. However, there is a gap in reviewing successful wet-lab methods employed in these animals, which compromises the fast progress in the detection of related biomolecules. Here, we provide a compilation of the most effective wet-lab methodologies for jellyfish venom extraction prior to proteomic analysis-separation, identification and toxicity assays. This includes SDS-PAGE, 2DE, gel chromatography, HPLC, DEAE, LC-MS, MALDI, Western blot, hemolytic assay, antimicrobial assay and protease activity assay. For a more comprehensive approach, jellyfish toxicity studies should further consider transcriptome sequencing. We reviewed such methodologies and other genomic techniques used prior to the deep sequencing of transcripts, including RNA extraction, construction of cDNA libraries and RACE. Overall, we provide an overview of the most promising methods and their successful implementation for optimizing time and effort when studying jellyfish.

De Novo Assembly and Characterization of Four Anthozoan (Phylum Cnidaria) Transcriptomes

Many nonmodel species exemplify important biological questions but lack the sequence resources required to study the genes and genomic regions underlying traits of interest. Reef-building corals are famously sensitive to rising seawater temperatures, motivating ongoing research into their stress responses and long-term prospects in a changing climate. A comprehensive understanding of these processes will require extending beyond the sequenced coral genome (Acropora digitifera) to encompass diverse coral species and related anthozoans. Toward that end, we have assembled and annotated reference transcriptomes to develop catalogs of gene sequences for three scleractinian corals (Fungia scutaria, Montastraea cavernosa, Seriatopora hystrix) and a temperate anemone (Anthopleura elegantissima). High-throughput sequencing of cDNA libraries produced ~20-30 million reads per sample, and de novo assembly of these reads produced ~75,000-110,000 transcripts from each sample with size distributions (mean ~1.4 kb, N50 ~2 kb), comparable to the distribution of gene models from the coral genome (mean ~1.7 kb, N50 ~2.2 kb). Each assembly includes matches for more than half the gene models from A. digitifera (54-67%) and many reasonably complete transcripts (~5300-6700) spanning nearly the entire gene (ortholog hit ratios ≥0.75). The catalogs of gene sequences developed in this study made it possible to identify hundreds to thousands of orthologs across diverse scleractinian species and related taxa. We used these sequences for phylogenetic inference, recovering known relationships and demonstrating superior performance over phylogenetic trees constructed using single mitochondrial loci. The resources developed in this study provide gene sequences and genetic markers for several anthozoan species. To enhance the utility of these resources for the research community, we developed searchable databases enabling researchers to rapidly recover sequences for genes of interest. Our analysis of de novo assembly quality highlights metrics that we expect will be useful for evaluating the relative quality of other de novo transcriptome assemblies. The identification of orthologous sequences and phylogenetic reconstruction demonstrates the feasibility of these methods for clarifying the substantial uncertainties in the existing scleractinian phylogeny.

Interspecific Differential Expression Analysis of RNA-Seq Data Yields Insight into Life Cycle Variation in Hydractiniid Hydrozoans

Hydrozoans are known for their complex life cycles, which can alternate between an asexually reproducing polyp stage and a sexually reproducing medusa stage. Most hydrozoan species, however, lack a free-living medusa stage and instead display a developmentally truncated form, called a medusoid or sporosac, which generally remains attached to the polyp. Although evolutionary transitions in medusa truncation and loss have been investigated phylogenetically, little is known about the genes involved in the development and loss of this life cycle stage. Here, we present a new workflow for evaluating differential expression (DE) between two species using short read Illumina RNA-seq data. Through interspecific DE analyses between two hydractiniid hydrozoans, Hydractinia symbiolongicarpus and Podocoryna carnea, we identified genes potentially involved in the developmental, functional, and morphological differences between the fully developed medusa of P. carnea and reduced sporosac of H. symbiolongicarpus. A total of 10,909 putative orthologs of H. symbiolongicarpus and P. carnea were identified from de novo assemblies of short read Illumina data. DE analysis revealed 938 of these are differentially expressed between P. carnea developing and adult medusa, when compared with H. symbiolongicarpus sporosacs, the majority of which have not been previously characterized in cnidarians. In addition, several genes with no corresponding ortholog in H. symbiolongicarpus were expressed in developing medusa of P. carnea. Results presented here show interspecific DE analyses of RNA-seq data to be a sensitive and reliable method for identifying genes and gene pathways potentially involved in morphological and life cycle differences between species.

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.

Comparative transcriptome resources of two Dysosma species (Berberidaceae) and molecular evolution of the CYP719A gene in Podophylloideae

Dysosma species (Berberidaceae, Podophylloideae) are of great medicinal pharmacogenetic importance and used as model systems to study the drivers and mechanisms of species diversification of temperate plants in East Asia. Recently, we have sequenced the transcriptome of the low-elevation D. versipellis. In this study, we sequenced the transcriptome of the high-elevation D. aurantiocaulis and used comparative genomic approaches to investigate the transcriptome evolution of the two species. We retrieved 53,929 unigenes from D. aurantiocaulis by de novo transcriptome assemblies using the Illumina HiSeq 2000 platform. Comparing the transcriptomes of both species, we identified 4593 orthologs. Estimation of Ka/Ks ratios for 3126 orthologs revealed that none had a Ka/Ks significantly greater than 1, whereas 1273 (Ka/Ks < 0.5, P < 0.05) were inferred to be under purifying selection. A total of 51 primer pairs were successfully designed from 461 EST-SSRs contained in 4593 orthologs. Marker validation assay revealed that 26 (51%) and 41 (80.4%) produced clear fragments with the expected sizes in all Podophylloideae species. Specifically, 19 different sequences of CYP719A were identified from PCR-amplified genomic DNA of all 12 species of Podophylloideae using primers designed from the assembled transcripts. The data further indicated that CYP719A was likely subject to strong selective constraints maintaining only one copy per genome. In Dysosma, there was relaxed purifying selection or more positive selection for high-elevation species. Overall, this study has generated a wealth of molecular resources potentially useful for pharmacogenetic and evolutionary studies in Dysosma and allied taxa.

The genetics of colony form and function in Caribbean Acropora corals

Background

Colonial reef-building corals have evolved a broad spectrum of colony morphologies based on coordinated asexual reproduction of polyps on a secreted calcium carbonate skeleton. Though cnidarians have been shown to possess and use similar developmental genes to bilaterians during larval development and polyp formation, little is known about genetic regulation of colony morphology in hard corals. We used RNA-seq to evaluate transcriptomic differences between functionally distinct regions of the coral (apical branch tips and branch bases) in two species of Caribbean Acropora, the staghorn coral, A. cervicornis, and the elkhorn coral, A. palmata.

Results

Transcriptome-wide gene profiles differed significantly between different parts of the coral colony as well as between species. Genes showing differential expression between branch tips and bases were involved in developmental signaling pathways, such as Wnt, Notch, and BMP, as well as pH regulation, ion transport, extracellular matrix production and other processes. Differences both within colonies and between species identify a relatively small number of genes that may contribute to the distinct “staghorn” versus “elkhorn” morphologies of these two sister species.

Conclusions

The large number of differentially expressed genes supports a strong division of labor between coral branch tips and branch bases. Genes involved in growth of mature Acropora colonies include the classical signaling pathways associated with development of cnidarian larvae and polyps as well as morphological determination in higher metazoans.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1133) contains supplementary material, which is available to authorized users.

Gene co-expression modules underlying polymorphic and monomorphic zooids in the colonial hydrozoan, Hydractinia symbiolongicarpus

Advances in sequencing technology have forced a quantitative revolution in Evolutionary Biology. One important feature of this renaissance is that comprehensive genomic resources can be obtained quickly for almost any taxon, thus speeding the development of new model organisms. Here, we analyze 20 RNA-seq libraries from morphologically, sexually, and genetically distinct polyp types from the gonochoristic colonial hydrozoan, Hydractinia symbiolongicarpus (Cnidaria). Analyses of these data using weighted gene co-expression networks highlight deeply conserved genetic elements of animal spermatogenesis and demonstrate the utility of these methods in identifying modules of genes that correlate with different zooid types across various statistical contrasts. RNA-seq data and analytical scripts described here are deposited in publicly available databases.

Differential gene expression between functionally specialized polyps of the colonial hydrozoan Hydractinia symbiolongicarpus (Phylum Cnidaria)

Background

A colony of the hydrozoan Hydractinia symbiolongicarpus comprises genetically identical yet morphologically distinct and functionally specialized polyp types. The main labor divisions are between feeding, reproduction and defense. In H. symbiolongicarpus, the feeding polyp (called a gastrozooid) has elongated tentacles and a mouth, which are absent in the reproductive polyp (gonozooid) and defensive polyp (dactylozooid). Instead, the dactylozooid has an extended body column with an abundance of stinging cells (nematocysts) and the gonozooid bears gonophores on its body column. Morphological differences between polyp types can be attributed to simple changes in their axial patterning during development, and it has long been hypothesized that these specialized polyps arose through evolutionary alterations in oral-aboral patterning of the ancestral gastrozooid.

Results

An assembly of 66,508 transcripts (>200 bp) were generated using short-read Illumina RNA-Seq libraries constructed from feeding, reproductive, and defensive polyps of H. symbiolongicarpus. Using several different annotation methods, approximately 54% of the transcripts were annotated. Differential expression analyses were conducted between these three polyp types to isolate genes that may be involved in functional, histological, and pattering differences between polyp types. Nearly 7 K transcripts were differentially expressed in a polyp-specific manner, including members of the homeodomain, myosin, toxin and BMP gene families. We report the spatial expression of a subset of these polyp-specific transcripts to validate our differential expression analyses.

Conclusions

While potentially originating through simple changes in patterning, polymorphic polyps in Hydractinia are the result of differentially expressed functional, structural, and patterning genes.

The differentially expressed genes identified in our study provide a starting point for future investigations of the developmental patterning and functional differences that are displayed in the different polyp types that confer a division of labor within a colony of H. symbiolongicarpus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-406) contains supplementary material, which is available to authorized users.

Old cell, new trick? Cnidocytes as a model for the evolution of novelty

Understanding how new cell types arise is critical for understanding the evolution of organismal complexity. Questions of this nature, however, can be difficult to answer due to the challenge associated with defining the identity of a truly novel cell. Cnidarians (anemones, jellies, and their allies) provide a unique opportunity to investigate the molecular regulation and development of cell-novelty because they possess a cell that is unique to the cnidarian lineage and that also has a very well-characterized phenotype: the cnidocyte (stinging cell). Because cnidocytes are thought to differentiate from the cell lineage that also gives rise to neurons, cnidocytes can be expected to express many of the same genes expressed in their neural "sister" cells. Conversely, only cnidocytes posses a cnidocyst (the explosive organelle that gives cnidocytes their sting); therefore, those genes or gene-regulatory relationships required for the development of the cnidocyst can be expected to be expressed uniquely (or in unique combination) in cnidocytes. This system provides an important opportunity to: (1) construct the gene-regulatory network (GRN) underlying the differentiation of cnidocytes, (2) assess the relative contributions of both conserved and derived genes in the cnidocyte GRN, and (3) test hypotheses about the role of novel regulatory relationships in the generation of novel cell types. In this review, we summarize common challenges to studying the evolution of novelty, introduce the utility of cnidocyte differentiation in the model cnidarian, Nematostella vectensis, as a means of overcoming these challenges, and describe an experimental approach that leverages comparative tissue-specific transcriptomics to generate hypotheses about the GRNs underlying the acquisition of the cnidocyte identity.

Global diversity and review of Siphonophorae (Cnidaria: Hydrozoa)

In this review the history of discovery of siphonophores, from the first formal description by Carl Linnaeus in 1785 to the present, is summarized, and species richness together with a summary of world-wide distribution of this pelagic group within the clade Hydrozoa discussed. Siphonophores exhibit three basic body plans which are briefly explained and figured, whilst other atypical body plans are also noted. Currently, 175 valid siphonophore species are recognized in the latest WoRMS world list, including 16 families and 65 genera. Much new information since the last review in 1987 is revealed from the first molecular analysis of the group, enabling identification of some new morphological characters diagnostic for physonect siphonophores. Ten types of nematocysts (stinging cells) are identified in siphonophores, more than in any other cnidarian; these are incorporated into batteries in the side branches of the tentacles in most species (here termed tentilla), and tentilla are reviewed in the last section of this paper. Their discharge mechanisms are explained and also how the tentilla of several physonect siphonophores are modified into lures. Of particular interest is the recent discovery of a previously unknown red fluorescent lure in the tentilla of the deep sea physonect Erenna, the first described example of emission of red light by an invertebrate to attract prey.

Production of a reference transcriptome and transcriptomic database (EdwardsiellaBase) for the lined sea anemone, Edwardsiella lineata, a parasitic cnidarian

Background

The lined sea anemone Edwardsiella lineata is an informative model system for evolutionary-developmental studies of parasitism. In this species, it is possible to compare alternate developmental pathways leading from a larva to either a free-living polyp or a vermiform parasite that inhabits the mesoglea of a ctenophore host. Additionally, E. lineata is confamilial with the model cnidarian Nematostella vectensis, providing an opportunity for comparative genomic, molecular and organismal studies.

Description

We generated a reference transcriptome for E. lineata via high-throughput sequencing of RNA isolated from five developmental stages (parasite; parasite-to-larva transition; larva; larva-to-adult transition; adult). The transcriptome comprises 90,440 contigs assembled from >15 billion nucleotides of DNA sequence. Using a molecular clock approach, we estimated the divergence between E. lineata and N. vectensis at 215–364 million years ago. Based on gene ontology and metabolic pathway analyses and gene family surveys (bHLH-PAS, deiodinases, Fox genes, LIM homeodomains, minicollagens, nuclear receptors, Sox genes, and Wnts), the transcriptome of E. lineata is comparable in depth and completeness to N. vectensis. Analyses of protein motifs and revealed extensive conservation between the proteins of these two edwardsiid anemones, although we show the NF-κB protein of E. lineata reflects the ancestral structure, while the NF-κB protein of N. vectensis has undergone a split that separates the DNA-binding domain from the inhibitory domain. All contigs have been deposited in a public database (EdwardsiellaBase), where they may be searched according to contig ID, gene ontology, protein family motif (Pfam), enzyme commission number, and BLAST. The alignment of the raw reads to the contigs can also be visualized via JBrowse.

Conclusions

The transcriptomic data and database described here provide a platform for studying the evolutionary developmental genomics of a derived parasitic life cycle. In addition, these data from E. lineata will aid in the interpretation of evolutionary novelties in gene sequence or structure that have been reported for the model cnidarian N. vectensis (e.g., the split NF-κB locus). Finally, we include custom computational tools to facilitate the annotation of a transcriptome based on high-throughput sequencing data obtained from a “non-model system.”

Adhesion networks of cnidarians: a postgenomic view

Cell-extracellular matrix (ECM) and cell-cell adhesion systems are fundamental to the multicellularity of metazoans. Members of phylum Cnidaria were classified historically by their radial symmetry as an outgroup to bilaterian animals. Experimental study of Hydra and jellyfish has fascinated zoologists for many years. Laboratory studies, based on dissection, biochemical isolations, or perturbations of the living organism, have identified the ECM layer of cnidarians (mesoglea) and its components as important determinants of stem cell properties, cell migration and differentiation, tissue morphogenesis, repair, and regeneration. Studies of the ultrastructure and functions of intercellular gap and septate junctions identified parallel roles for these structures in intercellular communication and morphogenesis. More recently, the sequenced genomes of sea anemone Nematostella vectensis, Hydra magnipapillata, and coral Acropora digitifera have opened up a new frame of reference for analyzing the cell-ECM and cell-cell adhesion molecules of cnidarians and examining their conservation with bilaterians. This chapter integrates a review of literature on the structure and functions of cell-ECM and cell-cell adhesion systems in cnidarians with current analyses of genome-encoded repertoires of adhesion molecules. The postgenomic perspective provides a fresh view on fundamental similarities between cnidarian and bilaterian animals and is impelling wider adoption of species from phylum Cnidaria as model organisms.

Genome wide transcriptome analysis of dendritic cells identifies genes with altered expression in psoriasis

Activation of dendritic cells by different pathogens induces the secretion of proinflammatory mediators resulting in local inflammation. Importantly, innate immunity must be properly controlled, as its continuous activation leads to the development of chronic inflammatory diseases such as psoriasis. Lipopolysaccharide (LPS) or peptidoglycan (PGN) induced tolerance, a phenomenon of transient unresponsiveness of cells to repeated or prolonged stimulation, proved valuable model for the study of chronic inflammation. Thus, the aim of this study was the identification of the transcriptional diversity of primary human immature dendritic cells (iDCs) upon PGN induced tolerance. Using SAGE-Seq approach, a tag-based transcriptome sequencing method, we investigated gene expression changes of primary human iDCs upon stimulation or restimulation with Staphylococcus aureus derived PGN, a widely used TLR2 ligand. Based on the expression pattern of the altered genes, we identified non-tolerizeable and tolerizeable genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (Kegg) analysis showed marked enrichment of immune-, cell cycle- and apoptosis related genes. In parallel to the marked induction of proinflammatory mediators, negative feedback regulators of innate immunity, such as TNFAIP3, TNFAIP8, Tyro3 and Mer are markedly downregulated in tolerant cells. We also demonstrate, that the expression pattern of TNFAIP3 and TNFAIP8 is altered in both lesional, and non-lesional skin of psoriatic patients. Finally, we show that pretreatment of immature dendritic cells with anti-TNF-α inhibits the expression of IL-6 and CCL1 in tolerant iDCs and partially releases the suppression of TNFAIP8. Our findings suggest that after PGN stimulation/restimulation the host cell utilizes different mechanisms in order to maintain critical balance between inflammation and tolerance. Importantly, the transcriptome sequencing of stimulated/restimulated iDCs identified numerous genes with altered expression to date not associated with role in chronic inflammation, underlying the relevance of our in vitro model for further characterization of IFN-primed iDCs.

Genome wide transcriptome analysis of dendritic cells identifies genes with altered expression in psoriasis

Activation of dendritic cells by different pathogens induces the secretion of proinflammatory mediators resulting in local inflammation. Importantly, innate immunity must be properly controlled, as its continuous activation leads to the development of chronic inflammatory diseases such as psoriasis. Lipopolysaccharide (LPS) or peptidoglycan (PGN) induced tolerance, a phenomenon of transient unresponsiveness of cells to repeated or prolonged stimulation, proved valuable model for the study of chronic inflammation. Thus, the aim of this study was the identification of the transcriptional diversity of primary human immature dendritic cells (iDCs) upon PGN induced tolerance. Using SAGE-Seq approach, a tag-based transcriptome sequencing method, we investigated gene expression changes of primary human iDCs upon stimulation or restimulation with Staphylococcus aureus derived PGN, a widely used TLR2 ligand. Based on the expression pattern of the altered genes, we identified non-tolerizeable and tolerizeable genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (Kegg) analysis showed marked enrichment of immune-, cell cycle- and apoptosis related genes. In parallel to the marked induction of proinflammatory mediators, negative feedback regulators of innate immunity, such as TNFAIP3, TNFAIP8, Tyro3 and Mer are markedly downregulated in tolerant cells. We also demonstrate, that the expression pattern of TNFAIP3 and TNFAIP8 is altered in both lesional, and non-lesional skin of psoriatic patients. Finally, we show that pretreatment of immature dendritic cells with anti-TNF-α inhibits the expression of IL-6 and CCL1 in tolerant iDCs and partially releases the suppression of TNFAIP8. Our findings suggest that after PGN stimulation/restimulation the host cell utilizes different mechanisms in order to maintain critical balance between inflammation and tolerance. Importantly, the transcriptome sequencing of stimulated/restimulated iDCs identified numerous genes with altered expression to date not associated with role in chronic inflammation, underlying the relevance of our in vitro model for further characterization of IFN-primed iDCs.

Count-based differential expression analysis of RNA sequencing data using R and Bioconductor

RNA sequencing (RNA-seq) has been rapidly adopted for the profiling of transcriptomes in many areas of biology, including studies into gene regulation, development and disease. Of particular interest is the discovery of differentially expressed genes across different conditions (e.g., tissues, perturbations) while optionally adjusting for other systematic factors that affect the data-collection process. There are a number of subtle yet crucial aspects of these analyses, such as read counting, appropriate treatment of biological variability, quality control checks and appropriate setup of statistical modeling. Several variations have been presented in the literature, and there is a need for guidance on current best practices. This protocol presents a state-of-the-art computational and statistical RNA-seq differential expression analysis workflow largely based on the free open-source R language and Bioconductor software and, in particular, on two widely used tools, DESeq and edgeR. Hands-on time for typical small experiments (e.g., 4-10 samples) can be <1 h, with computation time <1 d using a standard desktop PC.

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.

Characterization of differential transcript abundance through time during Nematostella vectensis development

BACKGROUND

Nematostella vectensis, a burrowing sea anemone, has become a popular species for the study of cnidarian development. In previous studies, the expression of a variety of genes has been characterized during N. vectensis development with in situ mRNA hybridization. This has provided detailed spatial resolution and a qualitative perspective on changes in expression. However, little is known about broad transcriptome-level patterns of gene expression through time. Here we examine the expression of N. vectensis genes through the course of development with quantitative RNA-seq. We provide an overview of changes in the transcriptome through development, and examine the maternal to zygotic transition, which has been difficult to investigate with other tools.

RESULTS

We measured transcript abundance in N. vectensis with RNA-seq at six time points in development: zygote (2 hours post fertilization (HPF)), early blastula (7 HPF), mid-blastula (12 HPF), gastrula (24 HPF), planula (5 days post fertilization (DPF)) and young polyp (10 DPF). The major wave of zygotic expression appears between 7-12 HPF, though some changes occur between 2-7 HPF. The most dynamic changes in transcript abundance occur between the late blastula and early gastrula stages. More transcripts are upregulated between the gastrula and planula than downregulated, and a comparatively lower number of transcripts significantly change between planula and polyp. Within the maternal to zygotic transition, we identified a subset of maternal factors that decrease early in development, and likely play a role in suppressing zygotic gene expression. Among the first genes to be expressed zygotically are genes whose proteins may be involved in the degradation of maternal RNA.

CONCLUSIONS

The approach presented here is highly complementary to prior studies on spatial patterns of gene expression, as it provides a quantitative perspective on a broad set of genes through time but lacks spatial resolution. In addition to addressing the problems identified above, our work provides an annotated matrix that other investigators can use to examine genes and developmental events that we do not examine in detail here.

Pheromone evolution, reproductive genes, and comparative transcriptomics in mediterranean earthworms (annelida, oligochaeta, hormogastridae)

Animals inhabiting cryptic environments are often subjected to morphological stasis due to the lack of obvious agents driving selection, and hence chemical cues may be important drivers of sexual selection and individual recognition. Here, we provide a comparative analysis of de novo-assembled transcriptomes in two Mediterranean earthworm species with the objective to detect pheromone proteins and other reproductive genes that could be involved in cryptic speciation processes, as recently characterized in other earthworm species. cDNA libraries of unspecific tissue of Hormogaster samnitica and three different tissues of H. elisae were sequenced in an Illumina Genome Analyzer II or Hi-Seq. Two pheromones, Attractin and Temptin were detected in all tissue samples and both species. Attractin resulted in a reliable marker for phylogenetic inference. Temptin contained multiple paralogs and was slightly overexpressed in the digestive tissue, suggesting that these pheromones could be released with the casts. Genes involved in sexual determination and fertilization were highly expressed in reproductive tissue. This is thus the first detailed analysis of the molecular machinery of sexual reproduction in earthworms.

SOLiD-SAGE of endophyte-infected red fescue reveals numerous effects on host transcriptome and an abundance of highly expressed fungal secreted proteins

One of the most important plant-fungal symbiotic relationships is that of cool season grasses with endophytic fungi of the genera Epichloë and Neotyphodium. These associations often confer benefits, such as resistance to herbivores and improved drought tolerance, to the hosts. One benefit that appears to be unique to fine fescue grasses is disease resistance. As a first step towards understanding the basis of the endophyte-mediated disease resistance in Festuca rubra we carried out a SOLiD-SAGE quantitative transcriptome comparison of endophyte-free and Epichloë festucae-infected F. rubra. Over 200 plant genes involved in a wide variety of physiological processes were statistically significantly differentially expressed between the two samples. Many of the endophyte expressed genes were surprisingly abundant, with the most abundant fungal tag representing over 10% of the fungal mapped tags. Many of the abundant fungal tags were for secreted proteins. The second most abundantly expressed fungal gene was for a secreted antifungal protein and is of particular interest regarding the endophyte-mediated disease resistance. Similar genes in Penicillium and Aspergillus spp. have been demonstrated to have antifungal activity. Of the 10 epichloae whole genome sequences available, only one isolate of E. festucae and Neotyphodium gansuense var inebrians have an antifungal protein gene. The uniqueness of this gene in E. festucae from F. rubra, its transcript abundance, and the secreted nature of the protein, all suggest it may be involved in the disease resistance conferred to the host, which is a unique feature of the fine fescue-endophyte symbiosis.

Comparative description of ten transcriptomes of newly sequenced invertebrates and efficiency estimation of genomic sampling in non-model taxa

INTRODUCTION

Traditionally, genomic or transcriptomic data have been restricted to a few model or emerging model organisms, and to a handful of species of medical and/or environmental importance. Next-generation sequencing techniques have the capability of yielding massive amounts of gene sequence data for virtually any species at a modest cost. Here we provide a comparative analysis of de novo assembled transcriptomic data for ten non-model species of previously understudied animal taxa.

RESULTS

cDNA libraries of ten species belonging to five animal phyla (2 Annelida [including Sipuncula], 2 Arthropoda, 2 Mollusca, 2 Nemertea, and 2 Porifera) were sequenced in different batches with an Illumina Genome Analyzer II (read length 100 or 150 bp), rendering between ca. 25 and 52 million reads per species. Read thinning, trimming, and de novo assembly were performed under different parameters to optimize output. Between 67,423 and 207,559 contigs were obtained across the ten species, post-optimization. Of those, 9,069 to 25,681 contigs retrieved blast hits against the NCBI non-redundant database, and approximately 50% of these were assigned with Gene Ontology terms, covering all major categories, and with similar percentages in all species. Local blasts against our datasets, using selected genes from major signaling pathways and housekeeping genes, revealed high efficiency in gene recovery compared to available genomes of closely related species. Intriguingly, our transcriptomic datasets detected multiple paralogues in all phyla and in nearly all gene pathways, including housekeeping genes that are traditionally used in phylogenetic applications for their purported single-copy nature.

CONCLUSIONS

We generated the first study of comparative transcriptomics across multiple animal phyla (comparing two species per phylum in most cases), established the first Illumina-based transcriptomic datasets for sponge, nemertean, and sipunculan species, and generated a tractable catalogue of annotated genes (or gene fragments) and protein families for ten newly sequenced non-model organisms, some of commercial importance (i.e., Octopus vulgaris). These comprehensive sets of genes can be readily used for phylogenetic analysis, gene expression profiling, developmental analysis, and can also be a powerful resource for gene discovery. The characterization of the transcriptomes of such a diverse array of animal species permitted the comparison of sequencing depth, functional annotation, and efficiency of genomic sampling using the same pipelines, which proved to be similar for all considered species. In addition, the datasets revealed their potential as a resource for paralogue detection, a recurrent concern in various aspects of biological inquiry, including phylogenetics, molecular evolution, development, and cellular biochemistry.

RNA-Seq reveals complex genetic response to Deepwater Horizon oil release in Fundulus grandis

Background

The release of oil resulting from the blowout of the Deepwater Horizon (DH) drilling platform was one of the largest in history discharging more than 189 million gallons of oil and subject to widespread application of oil dispersants. This event impacted a wide range of ecological habitats with a complex mix of pollutants whose biological impact is still not yet fully understood. To better understand the effects on a vertebrate genome, we studied gene expression in the salt marsh minnow Fundulus grandis, which is local to the northern coast of the Gulf of Mexico and is a sister species of the ecotoxicological model Fundulus heteroclitus. To assess genomic changes, we quantified mRNA expression using high throughput sequencing technologies (RNA-Seq) in F. grandis populations in the marshes and estuaries impacted by DH oil release. This application of RNA-Seq to a non-model, wild, and ecologically significant organism is an important evaluation of the technology to quickly assess similar events in the future.

Results

Our de novo assembly of RNA-Seq data produced a large set of sequences which included many duplicates and fragments. In many cases several of these could be associated with a common reference sequence using blast to query a reference database. This reduced the set of significant genes to 1,070 down-regulated and 1,251 up-regulated genes. These genes indicate a broad and complex genomic response to DH oil exposure including the expected AHR-mediated response and CYP genes. In addition a response to hypoxic conditions and an immune response are also indicated. Several genes in the choriogenin family were down-regulated in the exposed group; a response that is consistent with AH exposure. These analyses are in agreement with oligonucleotide-based microarray analyses, and describe only a subset of significant genes with aberrant regulation in the exposed set.

Conclusion

RNA-Seq may be successfully applied to feral and extremely polymorphic organisms that do not have an underlying genome sequence assembly to address timely environmental problems. Additionally, the observed changes in a large set of transcript expression levels are indicative of a complex response to the varied petroleum components to which the fish were exposed.

RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics

Recent rapid advances in next generation RNA sequencing (RNA-Seq)-based provide researchers with unprecedentedly large data sets and open new perspectives in transcriptomics. Furthermore, RNA-Seq-based transcript profiling can be applied to non-model and newly discovered organisms because it does not require a predefined measuring platform (like e.g. microarrays). However, these novel technologies pose new challenges: the raw data need to be rigorously quality checked and filtered prior to analysis, and proper statistical methods have to be applied to extract biologically relevant information. Given the sheer volume of data, this is no trivial task and requires a combination of considerable technical resources along with bioinformatics expertise. To aid the individual researcher, we have developed RobiNA as an integrated solution that consolidates all steps of RNA-Seq-based differential gene-expression analysis in one user-friendly cross-platform application featuring a rich graphical user interface. RobiNA accepts raw FastQ files, SAM/BAM alignment files and counts tables as input. It supports quality checking, flexible filtering and statistical analysis of differential gene expression based on state-of-the art biostatistical methods developed in the R/Bioconductor projects. In-line help and a step-by-step manual guide users through the analysis. Installer packages for Mac OS X, Windows and Linux are available under the LGPL licence from http://mapman.gabipd.org/web/guest/robin.

RNA-Seq and molecular docking reveal multi-level pesticide resistance in the bed bug

BACKGROUND

Bed bugs (Cimex lectularius) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of C. lectularius has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance.

RESULTS

We performed a next-generation RNA sequencing (RNA-Seq) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from de novo assembled mRNA-Seq tags (30,404 ESTs) and our previous 454 pyrosequenced database (21,088 ESTs). The two-way GLMseq analysis revealed ~15,000 highly significant differentially expressed ESTs between the PR and PS strains. Among the top 5,000 differentially expressed ESTs, 109 putative defense genes (cuticular proteins, cytochrome P450s, antioxidant genes, ABC transporters, glutathione S-transferases, carboxylesterases and acetyl cholinesterase) involved in penetration resistance and metabolic resistance were identified. Tissue and development-specific expression of P450 CYP3 clan members showed high mRNA levels in the cuticle, Malpighian tubules, and midgut; and in early instar nymphs, respectively. Lastly, molecular modeling and docking of a candidate cytochrome P450 (CYP397A1V2) revealed the flexibility of the deduced protein to metabolize a broad range of insecticide substrates including DDT, deltamethrin, permethrin, and imidacloprid.

CONCLUSIONS

We developed significant molecular resources for C. lectularius putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-Seq profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in C. lectularius. Future research that is targeted towards RNA interference (RNAi) on the identified metabolic targets such as cytochrome P450s and cuticular proteins could lay the foundation for a better understanding of the genetic basis of insecticide resistance in C. lectularius.