Differential gene expression between functionally specialized polyps of the colonial hydrozoan Hydractinia symbiolongicarpus (Phylum Cnidaria)
© Sanders et al.; licensee BioMed Central Ltd. 2014
Received: 24 February 2014
Accepted: 20 May 2014
Published: 28 May 2014
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.
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.
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.
KeywordsHydractinia symbiolongicarpus RNA-Seq Polymorphism Differential expression Transcriptome assembly Annotation
Colonial hydrozoans are composed of individual polyps connected through continuous epithelia and a shared gastrovascular cavity. Hydrozoans are members of the phylum Cnidaria, which are characterized by their diploblastic construction, comprising only two epithelial layers, the epidermis and gastrodermis. Despite their simple epithelial construction, many hydrozoan species evolved complex colonies through functional specialization of genetically identical yet morphologically distinct polyp types, conferring a division of labor within the colony [1–3]. This division of labor is known as polyp polymorphism [1–3].
It has long been hypothesized that these specialized polyps arose through evolutionary alterations in oral-aboral patterning in the ancestral gastrozooid [5–7]. Previous studies using candidate gene approaches have identified patterning genes specific to different polyp types. Cartwright et al.  focused on the involvement of Cnox-2, a parahox gene, in patterning these different polyp types of H. symbiolongicarpus. Immunolocalization of the Cnox-2 protein showed expression in body column tissue and down-regulation in oral structures of the gastrozooid. Mokady et al.  compared expression of Cn-ems (empty spiracles homolog) between gastrozooids and gonozooids of H. symbiolongicarpus. Whole mount in situ hybridization revealed no expression of Cn-ems in the gonozooid, while mRNAs were detected in the gastrodermal epithelia (“digestive cells”) of the gastrozooid.
More recently, Siebert et al.  used an RNA-Seq approach to examine differential expression between several polyp types of another hydrozoan, the siphonophore Nanomia bijuga. Although the focus of their paper was to evaluate next generation sequencing (NGS) platforms for differential expression (DE), they confirmed, through whole mount in situ hybridization, that at least one gene identified through their DE analyses (isogroup03256) was expressed in a polyp specific manner.
With the advent of NGS technologies, an unbiased approach to identify genes involved in the differentiation of different tissues (e.g. ) and developmental stages (e.g. ), or as well as those that are differentially expressed between species (e.g. ) can be made without reference to particular candidate genes. We report a transcriptome assembly, annotations, and DE analyses between three different polyp types in H. symbiolongicarpus. Our results, confirmed, with whole mount in situ hybridization, that DE analyses using RNA-Seq is a powerful approach for identifying genes and pathways involved in conferring a division of labor within this colonial organism.
Results and discussion
Transcriptome assembly and annotation
Differential expression analyses
Number of DE transcripts in different pairwise comparisons of the libraries
Dact vs Gast
Dact vs Gono
Gast vs Gono
Male vs Female
Number of transcripts identified as always up- or down-regulated in a specific polyp
padj < .05
Our DE analyses revealed several polyp-specific genes that are consistent with previous studies using candidate gene approaches in cnidarians as discussed below. Furthermore, our analyses revealed additional genes that were not previously considered to play specific developmental, functional and/or structural roles in cnidarians. Below we summarize of few of these results and suggest areas of interest for further study.
While adjusting for differences in gene expression between males and females greatly reduced the effect of heterogametic expression on the DE analyses, genes likely involved in non-sex specific gametogenesis were found up-regulated in the gonozooids. Of the 76 polyp-specific DE transcripts annotated with functional terms that include mitosis, cell cycle, and germline maintenance (Additional file 5), 69 are up-regulated in the gonozooids, including four DE transcripts annotated as known hydrozoan stem cell markers nanos, vasa, and piwi. This is consistent with expression studies of nanos and vasa genes in a closely related species, Hydractinia echinata[14, 15].
Several homeobox transcripts are differentially expressed between different polyp types (Additional file 5). Homeobox genes up-regulated in the gastrozooid include members of the LIM (lhx), sine oculus (six), empty spiracles (ems), and PRD classes, confirmed by molecular phylogenetic analysis of cnidarian homeodomains (Additional file 5; Additional file 6). The up-regulation of the empty spiracles homolog (100% bootstrap [BS] support; Additional file 5; Additional file 6), Cn-ems, is consistent with the findings of Mokady et al.  discussed previously. Up-regulated gastrozooid expression of two lhx-like transcripts, one six-like and one orthopedia (PRD class) transcript (100%, 99%, and 98% BS support, respectively; Additional file 5; Additional file 6) is also consistent with expression studies in other cnidarians, including Aurelia, Nematostella[18–20], Craspedacusta, Cladonema, and Podocoryna, where their expression was found in regions specific to feeding and/or digestion, including tentacles and gastric tissue.
One of the homeodomain-containing transcripts up-regulated in the gonozooid belongs to the POU class (Additional file 5; Additional file 6). Expression of POU homeodomain transcription factors has also been categorized in other cnidarians, including Aurelia and H. echinata. In H. echinata, the POU gene, pln, is expressed around interstitial stem cell (i-cells) . The H. symbiolongicarpus ortholog to pln (100% BS support; Additional file 5; Additional file 6) is up-regulated in the gonozooid, which is consistent with that of the other stem cell markers mentioned previously.
Myosin genes are a superfamily of molecular motor proteins, primarily associated with muscular contraction and cell movement. Here we find a complex pattern of differential expression of several different myosin transcripts up-regulated in each polyp type (four, six, and four transcripts in the gastrozooids, gonozooids, and dactylozooids, respectively), spanning several myosin classes (Additional file 5, Additional file 7). Of particular note is the up-regulation of a tropomyosin transcript in the gonozooids. In the hydrozoan Podocoryna carnea, a tropomyosin, tpm2, is expressed solely in the striated muscle of the developing and adult medusa life cycle stages and not in the polyp , as opposed to tpm1, which is ubiquitously expressed in both polyp and medusae stages . In Hydractinia, gonophore development is greatly truncated and never reaches the medusae stage. Instead Hydractinia forms sporosacs, which are believed to lack all medusae like features, including striated muscle necessary for medusae to swim [27–29]. Phylogenetic analysis of cnidarian myosins did not recover any well-supported orthologous relationship between this polyp-specific tropomyosin and other known cnidarian tropomyosins, although orthology assignments of several other polyp-specific myosins were revealed (Additional file 7). Further discovery of tropomyosin genes in additional cnidarian taxa are necessary to determine if different tropomyosin orthologs are specific to certain medusae features and/or reduced developmental forms.
While research into the characterization and properties of cnidarian toxins is on the rise, very little is known of their function and location of endogenous expression . We identified 13 DE transcripts annotated as some type of toxin (three up-regulated in the gastrozooids, seven in dactylozooids, and one in gonozooids; two down-regulated in gonozooids; Additional file 5). Phylogenetic analysis of cnidarian toxins recovered a monophyletic cluster of six H. symbiolongicarpus ‘echotoxin’ transcripts as sister to a group of anthozoan toxin genes (60% BS support; Additional file 8), and a strongly supported (92% BS; Additional file 8) sister relationship between a four H. symbiolongicarpus toxins and two scyphozoan toxins from Aurelia (TX1 and TX2; Additional file 8). The remaining three polyp-specific toxins were not placed in any well-supported orthologous groups. Further study is warranted to determine if these toxins each play a unique role in different functions, such as prey capture, defense, and/or digestion.
A large number of transcripts belonging to the astacin subfamily are up-regulated in the gastrozooid (44 total, Additional file 5; Additional file 9), consistent with one of their roles as digestive enzymes in other metazoans [31–34]. Expression studies of several astacin genes in hydrozoans also suggest a role in digestion. In P. carnea, pmp1 is expressed in both the mouth of the polyp and the manubrium of the medusa stage . Immunolocalization of the HMP1 protein found it expressed in the head and tentacle regions of Hydra, while Kumpfmüller et al.  found farm1 expressed in both the epi and gastrodermal layers of gastric region of Hydra.
Its important to note that digestion is just one function of the astacin subfamily. Another function is in regeneration, as shown in H. echinata, where Möhrlen et al.  found astacins hea1 and hea2 expressed throughout development and soon after the gastrozooid is subjected to tissue injury (expression in other polyp types not mentioned). HMP1 was also up-regulated during head regeneration in Hydra. Orthologs of hea1 and hea2 were among the 44 gastrozooid-specific astacins in our study (94% and 99% BS support, respectively; see Additional file 5; Additional file 9). Up-regulation of these transcripts may be a result of tissue damage response during dissections prior to RNA extractions. However, it is interesting that they are specific to the gastrozooid, suggesting that gonozooids and dactylozooids may have different regenerative properties than gastrozooids [1, 3, 39].
Polyp-specific DE transcripts analyzed with whole mount in situ hybridization
Top blast hit
Sea anemone cytotoxic protein
myosin heavy chain isoform a
Sea anemone cytotoxic protein
transcription factor capicua
HMG (high mobility group)
indian hedgehog b
Protein kinase domain
A myosin gene, referred to here as myosin_48857, was identified as a dactylozooid-specific gene by DE analyses. ISH confirmed this, recovering expression limited to the ectoderm of the body column of the extended side when the dactylozooid is curled in on itself (Figure 6). Minor expression is also detected around the base of the gonophores and on the body column of some gonozooids (not shown).
The gene cerberus is also found to be dactylozooid-specific by DE analysis. This gene is only expressed in the gastrodermis beneath the clusters of nematocysts at the distal end of the dactylozooid (Figure 6). Expression studies of cerberus in other metazoans have shown it to act as an antagonist of TGF-ß and Wnt signaling [43, 44]. Here, expression in the dactylozooids is consistent with its antagonist role in Wnt signaling. H. symbiolongicarpus’ canonical wnt, HsWnt3 [GenBank:KF745052], is expressed at the distal tip of the dactylozooid (not shown, unpublished). Cerberus is expressed at the proximal boundary of HsWnt3 expression, potentially acting to maintain HsWnt3’s expression boundary. This however appears to be specific to the dactylozooids, as H. symbiolongicarpus feeding polyps express Wnt3 (not shown, unpublished), similar to other hydrozoan feeding polyps including H. echinata[45–47], P. carnea (unpublished), and Hydra[48–52]), but do not express cerberus (Figure 6).
ISH also confirmed the specificity of several gonozooid-specific DE transcripts. Expression of the hedgehog homolog is restricted to the gastrodermis of both male and female gonophores (Figure 6). Expression of a bmp receptor gene, tentatively called bmpR_54452, and capicua are primarily limited to developing oocytes in females and the gastrodermis of male gonophores (Figure 6). ISH expression patterns of these transcripts suggest their involvement in some stage of meiotic/mitotic division during gametogenesis. In Hydractinia, oogenesis begins in the germinal zone (body column) of female gonozooids and oocyte differentiation continues after moving into the gonophores [1, 2, 53], while spermatogenesis takes place entirely in the gastrodermis of the male gonophores [2, 53].
For several of these transcripts, expression in the females might not only be associated with germline proliferation, but with maternal transcript generation as well. Maternal expression of capicua and BMP receptors in early embryonic development has been reported in other metazoans [54, 55]. Expression corresponding to maternal transcript generation is consistent with strong expression around developing oocytes in the germinal zone. By contrast, ISH of hedgehog in Hydractinia recovered no expression in the germinal zone of female gonophores (Figure 6). Instead, its expression was limited to the gastrodermal tissues surrounding maturing oocytes in female gonophores. This is consistent with hedgehog genes implicated in germline proliferation and differentiation in other metazoans [56, 57] and in Nematostella, where one hedgehog appears to be involved in germline proliferation, but lacks maternal expression .
Our non-biased approach of characterizing differential expression in different polyp types enabled us to identify key genes potentially involved in the morphological and functional differences between these different polyps. However, in interpreting results from a DE analysis, it is important to understand the distinction between biological relevance and statistical significance. We do not propose that every transcript in our list of putative polyp-specific genes is involved in the patterning or function of these different polyps, nor do we report to have captured all polyp-specific genes. One type of information not captured in this method would be those genes whose spatial or temporal expression (but not abundance) confers differences between polyp types. For example, the parahox gene Cnox-2, which was shown to be expressed in all polyps uniformly except for the oral region of the gastrozooid  was not recovered in the DE analysis. This is likely due to the fact that Cnox-2 has different patterns of expression but not distinct differences in abundance between polyp types.
Even given the potential limitations, this unbiased approach of RNA-Seq DE analysis, validated through in situ hybridization, identified many potential patterning and functional/structural genes without limiting our investigations to particular candidate genes. While potentially originating through simple changes in patterning, polymorphic polyps in Hydractinia are the result of differentially expressed functional, histological, and patterning genes. The DE 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.
Colonies of H. symbiolongicarpus encrusting gastropod shells occupied by the hermit crab Pagurus longicarpus were purchased from Marine Biological Laboratories (Woods Hole MA). Some colonies of H. symbiolongicarpus were surgically explanted onto microscope slides, placed in slide racks kept in seawater (REEF CRYSTALS, Aquarium Systems) aquaria, maintained at 21°C, and fed 2-3-day-old nauplii of Artemia three times a week. P. longicarpus were maintained in similar conditions and fed frozen shrimp three times a week.
Tissue collection and RNA isolation
Tissue and RNA preps were divided into two categories based on the ultimate use of the samples (transcriptome assembly or DE analyses). Gastrozooids, gonozooids, and dactylozooids were individually dissected and collected from colonies encrusting the gastropod shells inhabited by P. longicarpus. The fourth polyp type (tentaculozooid) was not collected due to its rare occurrence in a colony. Excised polyps were immediately flash-frozen and stored at -80°C until RNA extractions were performed. Care was taken to only include polyp tissue and to exclude tissue from the stolons and stolonal mat of the colony. In order to obtain sufficient quantities of tissue, polyps from multiple colonies were often pooled together.
RNA extractions were carried out on pooled samples of approximately 100 individuals of a single polyp type. Total RNA was isolated using the TriReagent isolation protocol (Invitrogen) followed by a DNase treatment using the TURBO DNase kit (Ambion) or performed at the University of Kansas Medical Genome Sequencing Facility (KUMC-GSF) according to standard Illumina protocols. In samples collected for transcriptome assembly, gonozooid samples were from both male and female colonies and were pooled together during RNA extraction whereas, for the gonozooid samples collected for the downstream DE analyses, males and females were kept separate from tissue collection through sequencing.
Library construction and sequencing
RNA libraries were constructed according to the TruSeq RNA Sample Preparation Guide (Illumina) using the TruSeq RNA Sample Preparation Kit (Box A). To increase transcript discovery, libraries used for transcriptome assembly were normalized using the Evrogen duplex-specific thermostable nuclease (DSN) kit following the Illumina DSN Normalization protocol. DNA fragments with adapters ligated on both ends were PCR-enriched after DSN normalization. Three normalized libraries were constructed with an average insert size of 160 bp and subsequently barcoded, pooled, and multiplexed across three lanes of an Illumina HiSeq2000 flowcell.
For DE analyses, a total of twelve other libraries (four for each polyp type, including two male and two female gonozooid libraries) were constructed similarly, but without DSN normalization at KUMC-GSF. These samples were barcoded, pooled, and multiplexed on a single lane of an Illumina HiSeq2500 flowcell. All libraries were 100 bp paired-end and sequenced at KUMC-GSF.
Transcriptome assembly and annotation
The workflow from sequencing through transcript annotation and differential expression analyses is shown in Figure 2. Raw reads from all three normalized libraries were filtered based on quality score and separately mapped to a set of unpublished genomic scaffolds of H. symbiolongicarpus using TopHat 2.0.6 . TopHat alignments were assembled into transcripts using Cufflinks 2.1.1 , generating three separate assemblies, one for each library. These assemblies were then merged into a single assembly using the cuffmerge function from Cufflinks . This assembly was then filtered by transcript size, removing assembled transcripts less than 200 bp in length. This assembly has been submitted to the NCBI Transcriptome Shotgun Assembly (TSA) database (Accession Number GAWH00000000 ). The raw reads have been submitted to the NCBI Sequence Read Archive (SRA; Project Number: SRX474462).
Transcripts were annotated using several different methods. Gene Ontology (GO) terms were added with Blast2GO [62, 63], using the BLASTX algorithm and a significance threshold of 1 × 10-03 to search against NCBI’s non-redundant (NR) protein database. Annotation names from the GO analysis represent the top BLAST hit (Additional file 2). A set of conserved eukaryotic genes was identified with CEGMA v2.4  (Additional file 3). HMM (hidden markov model) protein families from the PFAM  and TIGR  databases were assigned to the amino acid translation of the most likely reading frame (identified using an open reading frame prediction tool ) of each transcript using HMMscan  under default settings. HMMscan annotations were constrained to a significance threshold of 0.01 (Additional file 2). Orthogroups were assigned to the same amino acid translations using the orthoMCL web server  (Additional file 2).
Differential expression analyses
Reads from the 12 non-normalized RNASeq libraries were mapped to the transcriptome assembly using Bowtie2 2.0.2 . The raw reads from these libraries have been submitted to the NCBI Sequence Read Archive (SRA; Project Number: SRX474878). Counts for transcripts for each library were extracted from the bowtie output (.sam files) using a python script that only counts reads in which both paired reads mapped to the same transcript (Additional file 11). The count data for each library was then fed through the DESeq  and edgeR  packages to assess statistically significant DE between all pairwise combinations of polyp types, including a comparison between male and female gonozooids. Both methods were used because they often give distinctly different results, with DESeq generally being more conservative in its assessment [73–75].
Given that the goal of this study was to identify differential gene expression between somatic tissues in the different polyp types, it was necessary to reduce the effect of gametogenic expression for the DE analyses. In Hydractinia, there are no discernable morphological differences between male and female gonozooids aside from the type of gametes present. Thus it can be assumed that any differences in expression between male and female gonozooids can be attributed to differences in gametogenesis (heterogametic expression) and need to be accounted for prior to DE analyses between polyp types.
These DE analyses produced a list of DE transcripts specific to one or more of the three pairwise comparisons made between the different libraries, but not ones truly specific to a certain polyp type. In order to identify these polyp-specific transcripts, only transcripts significantly up- or down-regulated (padj < 0.05) in a particular polyp when compared to either of the other two polyp types (must be significant in only two of the three pairwise comparisons) from both edgeR and DESeq were considered polyp-specific. Figure 5 is a Venn diagram that lists the number of transcripts that meet these requirements (Additional file 4).
Probe Synthesis and in situhybridization
Several polyp-specific transcripts identified during the DE analyses were selected for confirmation and further investigation with whole mount in situ hybridization (ISH) experiments (Table 3). Sequences for these transcripts were identified in the assembly, amplified from cDNA, cloned using the Invitrogen TOPO-TA Cloning Kit, and anti-DIG labeled riboprobes were synthesized from clones using the Invitrogen T7/T3 Megascript kit. ISH of these transcripts were performed following methods from Nawrocki & Cartwright .
Molecular phylogenetic analyses
Several gene trees were constructed of select gene families, including homeodomains, myosins, toxins, and astacins (Additional file 6; Additional file 7; Additional file 8; Additional file 9). Cnidarian sequences belonging to families of interest were mined from the nr NCBI database and aligned using Mafft . Depending on the family, either the L-insi or E-insi alignment algorithm was used. Only polyp-specific H. symbiolongicarpus sequences annotated with these families were included in the alignments. Maximum likelihood estimates of the molecular phylogenies of these gene families were then produced using RAxML  on the CIRPES portal  using the rapid bootstrapping (-f a) algorithm with 1000 bootstrap replicates under the PROTGAMMA + WAG model.
Differential expression or differentially expressed
Next generation sequencing
Whole mount in situ hybridization
Duplex-Specific thermostable nuclease
Hidden markov model
Open reading frame.
We would like to thank Uri Frank for help with the ISH protocol. We thank Clark Bloomer at KUMC-GSF for help in RNA isolation, library construction, and sequencing. We are grateful for the draft genomic scaffolds provided by Christine Schnitzler and Andy Baxevanis and for their comments on an earlier draft of this manuscript. This works was supported by NSF grant DEB‒095357 and MDIBL Visiting Scientist Fellowship to PC, and the MDIBL Towle Fellowship to SMS.
- Muller W: Experimentelle Untersuchungen Über Stockentwicklung, Polypendifferenzierung, und sexualchimären bei Hydractinia echinata. Wilhelm Roux’ Arch Entwicklungsmech Org. 1964, 155: 181-268. 10.1007/BF00573905.View ArticleGoogle Scholar
- Berrill NJ: Growth and form in gymnoblastic hydroids; polymorphism within the Hydractiniidae. J Morph. 1953, 92: 241-272. 10.1002/jmor.1050920203.View ArticleGoogle Scholar
- Burnett AL, Sindelar W, Diehl N: An examination of polymorphism in the hydroid Hydractinia echinata. J Mar Biol Assoc UK. 1967, 47: 645-658. 10.1017/S0025315400035256.View ArticleGoogle Scholar
- Cartwright P, Nawrocki AM: Character evolution in Hydrozoa (phylum Cnidaria). J Integ Comp Biol. 2010, 50: 456-472. 10.1093/icb/icq089.View ArticleGoogle Scholar
- Agassiz A: North American Acalephae. 1865, Welsh, Bigelow, Cambridge, U.K: University PressView ArticleGoogle Scholar
- Haeckel E: Report on the Siphonophorae Collected by H. M. S. Challenger During the Years 1873–76. 1888, London: The Challenger ReportsView ArticleGoogle Scholar
- Huxley TH: The Oceanic Hydrozoa: A Description of the Calycophoridae and Physophoridae Observed During the Voyage of the H. M. S. “Rattlesnake” in the years 1846–1850: With a General Introduction. 1859, London: Ray SocietyGoogle Scholar
- Cartwright P, Bowsher J, Buss LW: Expression of the HOX type gene, Cnox-2, and the division of labor in a colonial hydroid. Proc Natl Acad Sci USA. 1999, 96: 2183-2186. 10.1073/pnas.96.5.2183.PubMed CentralPubMedView ArticleGoogle Scholar
- Mokady O, Dick MH, Lackschewitz D, Schierwater B, Buss LW: Over one-half billion years of head conservation? Expression of an ems class gene in Hydractinia symbiolongicarpus (Cnidaria: Hydrozoa). Proc Natl Acad Sci USA. 1998, 95: 3673-3678. 10.1073/pnas.95.7.3673.PubMed CentralPubMedView ArticleGoogle Scholar
- Siebert S, Robinson MD, Tintori SC, Goetz F, Helm RR, Smith SA, Shaner N, Haddock SHD, Dunn CW: Differential gene expression in the siphonophore Nanomia bijuga (Cnidaria) assessed with multiple next-generation sequencing workflows. PLoS One. 2011, 6: e22953-10.1371/journal.pone.0022953.PubMed CentralPubMedView ArticleGoogle Scholar
- Hao DC, Ge G, Xiao P, Zhang Y, Yang L: The first insight into the tissue specific Taxus transcriptome via Illumina second generation sequencing. PLoS One. 2011, 6: e21220-10.1371/journal.pone.0021220.PubMed CentralView ArticleGoogle Scholar
- Helm RR, Siebert S, Tulin S, Smith J, Dunn CW: Characterization of differential transcript abundance through time during Nematostella vectensis development. BMC Genomics. 2013, 14: 266-10.1186/1471-2164-14-266.PubMed CentralPubMedView ArticleGoogle Scholar
- Lui S, Lin L, Jiang P, Wang D, Xing Y: A comparison of RNA-Seq and high-density exon array for detecting differential expression between closely related species. Nucl Acids Res. 2010, 39: 578-588.Google Scholar
- Kanska J, Frank U: Novels roles for Nanos in neural cell fate determination revealed by studies in a cnidarian. J Cell Science. 2013, doi:10.1242/ jcs.127233Google Scholar
- Rebscher N, Volk C, Teo R, Plickert G: The germ plasm component vasa allows tracing of the interstitial stem cells in the Cnidarian Hydractinia echinata. Dev Dynamics. 2008, 237: 1736-1745. 10.1002/dvdy.21562.View ArticleGoogle Scholar
- Seipel K, Yanze N, Schmid V: The germ line and somatic stem cell gene Cniwi in the jellyfish Podocoryne carnea. J Dev Biol. 2004, 48: 1-7. 10.1387/ijdb.15005568.View ArticleGoogle Scholar
- Bebenek HG, Gates RD, Morris J, Hartenstein V, Jacobs DK: sine oculis in basal Metazoa. Dev Genes Evol. 2004, 214: 342-351.PubMedView ArticleGoogle Scholar
- Martindale MQ, Pang K, Finnerty JR: Investigating the origins of triploblasty: ‘mesodermal’ gene expression in a diploblastic animal, the sea anemone Nematostella vectensis (phylum, Cnidaria; class, Anthozoa). Dev. 2004, 131: 2463-2474. 10.1242/dev.01119.View ArticleGoogle Scholar
- Srivastava M, Larroux C, Lu DR, Mohanty K, Chapman J, Degnan BM, Rokhsar DS: Early evolution of the LIM homeobox gene family. BMC Biol. 2010, 8: 4-10.1186/1741-7007-8-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Mazza ME, Pang K, Reitzel AM, Martindale MQ, Finnerty JR: A conserved cluster of three PRD-class homeobox genes (homeobrain, rx and orthopedia) in the Cnidaria and Protostomia. EvoDevo. 2010, 1: 3-10.1186/2041-9139-1-3.PubMed CentralPubMedView ArticleGoogle Scholar
- Hroudova M, Vojta P, Strnad H, Krejcik Z, Ridl J, Paces J, Vlcek C, Paces V: Diversity, phylogeny and expression patterns of pou and six homeodomain transcription factors in hydrozoan jellyfish Craspedacusta sowerbyi. PLoS One. 2012, 7: e36420-10.1371/journal.pone.0036420.PubMed CentralPubMedView ArticleGoogle Scholar
- Stierwald M, Yanze N, Bamert RP, Kammermeier L, Schmid V: The Sine oculis/Six class family of homeobox genes in jellyfish with and without eyes: development and eye regeneration. Dev Biol. 2004, 274: 70-81. 10.1016/j.ydbio.2004.06.018.PubMedView ArticleGoogle Scholar
- Nakanishi N, Yuan D, Hartenstein V, Jacobs DK: Evolutionary origin of rhopalia: insights from cellular-level analyses of Otx and POU expression patterns in the developing rhopalial nervous system. Evol Dev. 2010, 12: 404-415. 10.1111/j.1525-142X.2010.00427.x.PubMedView ArticleGoogle Scholar
- Millane RC, Kanska J, Duffy DJ, Seoighe C, Cunningham S, Plickert G, Frank U: Induced stem cell neoplasia in a cnidarian by ectopic expression of a POU domain transcription factor. Dev. 2011, 138: 2429-2439. 10.1242/dev.064931.View ArticleGoogle Scholar
- Gröger H, Callaerts P, Gehring WJ, Schmid V: Gene duplication and recruitment of a specific tropomyosin into striated muscle cells in the jellyfish Podocoryne carnea. J Exp Zoo. 1999, 285: 378-386. 10.1002/(SICI)1097-010X(19991215)285:4<378::AID-JEZ5>3.0.CO;2-C.View ArticleGoogle Scholar
- Baader CD, Schmid V, Schuchert P: Characterization of a tropomyosin cDNA from the hydrozoan Podocoryne carnea. Fed Eur Biochem Soc. 1993, 328: 63-66. 10.1016/0014-5793(93)80966-X.View ArticleGoogle Scholar
- Miglietta MP, Cunningham CW: Evolution of life cycle, colony morphology, and host specificity in the family Hydractiniidae (Hydrozoa, Cnidaria). Evol. 2012, doi:10.1111/j.1558-5646.2012.01717.xGoogle Scholar
- Bouillon J, Medel D, Cantero ALP: The taxonomic status of the genus Stylactaria Stechow, 1921 (Hydroidomedusae, Anthomedusae, Hydractiniidae), with the description of a new species. Sci. 1997, 61: 471-486.Google Scholar
- Boero F, Sará M: Motile sexual stages and evolution of Leptomedusae (Cnidaria). Boll Zool. 1987, 54: 131-139. 10.1080/11250008709355572.View ArticleGoogle Scholar
- Frazão B, Vasconcelos V, Antunes A: Sea anemone (Cnidaria, Anthozoa, Actiniaria) toxins: an overview. Mar Drugs. 2012, 10: 1812-1851. 10.3390/md10081812.PubMed CentralPubMedView ArticleGoogle Scholar
- Yasumasu S, Yamada K, Akasaka K, Mitsunaga K, Iuchi I, Shimada H, Yamagami K: Isolation of cDNAs for LCE and HCE, two constituent proteases of the hatching enzyme of Oryzias latipes, and concurrent expression of their mRNAs during development. Dev Biol. 1992, 153: 250-258. 10.1016/0012-1606(92)90110-3.PubMedView ArticleGoogle Scholar
- Hishida R, Ishihara T, Kondo K, Katsura I: hch-1, a gene required for normal hatching and normal migration of a neuroblast in C. elegans, encodes a protein related to TOLLOID and BMP-1. EMBO J. 1996, 15: 411-4122.Google Scholar
- Geier G, Zwilling R: Cloning and characterization of a cDNA coding for Astacus embryonic astacin, a member of the astacin family of metalloproteases from the crayfish Astacus astacus. Eur J Biochem. 1998, 253: 796-803. 10.1046/j.1432-1327.1998.2530796.x.PubMedView ArticleGoogle Scholar
- Fan TJ, Katagiri C: Properties of the hatching enzyme from Xenopus laevis. Eur J Biochem. 2001, 268: 4892-4898. 10.1046/j.0014-2956.2001.02420.x.PubMedView ArticleGoogle Scholar
- Pan T, Gröger H, Schmid V, Spring J: A toxin homology domain in an astacing-like metallopteinase of the jellyfish of Podocoryne carnea with a dual role in digestion and development. Dev Genes Evol. 1998, 208: 259-266. 10.1007/s004270050180.PubMedView ArticleGoogle Scholar
- Yan L, Pollock GH, Nagase H, Sarras MP: A 25.7×103M r hydra metalloproteinase (HMP1), a member of the astacin family, localizes to the extracellular matrix of Hydra vulgaris in a head-specific manner and has a developmental function. Dev. 1995, 121: 1591-1602.Google Scholar
- Kumpfmüller G, Rybakine V, Takahashi T, Fujisawa T, Bosch TCG: Identification of an astacin matrix metalloprotease as target gene for Hydra foot activator peptides. Dev Genes Evol. 1995, 209: 601-607.Google Scholar
- Möhrlen F, Maniura M, Plickert G, Frohme M, Frank U: Evolution of astacin-like metalloproteases in animals and their function in development. Evol Dev. 2006, 8: 223-231. 10.1111/j.1525-142X.2006.00092.x.PubMedView ArticleGoogle Scholar
- Hazen AP: Regeneration in Hydractinia and Podocoryne. Amer Nat. 1902, 36: 193-200. 10.1086/278099.View ArticleGoogle Scholar
- Campbell RD: Cell proliferation and morphological patterns in the hydroids Tubularia and Hydractinia. J Embryol Exp Morph. 1967, 17: 607-616.PubMedGoogle Scholar
- Braverman M: Studies on hydroid differentiation; the replacement of hypostomal gland cells of Podocoryne carnea. J Morph. 1968, 126: 95-106. 10.1002/jmor.1051260106.PubMedView ArticleGoogle Scholar
- Bouillon J: Les cellules glandulaires des hydroïdes et des hydroméduses; leur structure et la mature de leurs sécrétions. Cahiers Biologie Marine. 1966, 7: 157-205.Google Scholar
- Hsu DR, Economides AN, Wang X, Eimon PM, Harland RM: The Xenpus dorsalizing factor Gremlin indentifies a novel family of secreted proteins that antagonize BMP activities. Mol Cell. 1998, 1: 673-683. 10.1016/S1097-2765(00)80067-2.PubMedView ArticleGoogle Scholar
- Piccolo S, Agius E, Leyns L, Bhattacharyya B, Grunz H, Bouwmeester T, De Robertis EM: The head inducer Cerberus is a multifunctional antagonist of Nodal, BMP and Wnt signals. Nature. 1999, 397: 707-710. 10.1038/17820.PubMed CentralPubMedView ArticleGoogle Scholar
- Duffy DJ, Plickert G, Kuenzel T, Tilmann W, Frank U: Wnt signaling promotes oral but suppresses aboral structures in Hydractinia metamorphosis and regeneration. Dev. 2010, 137: 3057-3066. 10.1242/dev.046631.View ArticleGoogle Scholar
- Müller W, Frank U, Teo R, Mokady O, Guette C, Plickert G: Wnt signaling in hydroid development: ectopic heads and giant buds induced by GSK-3ß inhibitors. Int J Dev Biol. 2007, 51: 211-220. 10.1387/ijdb.062247wm.PubMedView ArticleGoogle Scholar
- Plickert G, Jacoby V, Frank U, Müller WA, Mokady O: Wnt signaling in hydroid development: formation of the primary body axis in embryogenesis and its subsequent patterning. Dev Biol. 2006, 298: 368-378. 10.1016/j.ydbio.2006.06.043.PubMedView ArticleGoogle Scholar
- Broun M, Gee L, Reinhardt B, Bode HR: Formation of the head organizer in hydra involves the canonical Wnt pathway. Dev. 2005, 132: 2907-2916. 10.1242/dev.01848.View ArticleGoogle Scholar
- Gee L, Hartig J, Law L, Wittlieb J, Khalturin K, Bosch TCG, Bode HR: ß-catenin plays a central role in setting up the head organizer in Hydra. Dev Biol. 2010, 340: 116-124. 10.1016/j.ydbio.2009.12.036.PubMedView ArticleGoogle Scholar
- Guder C, Pinho S, Nacak TG, Schmidt HA, Hobmayer B, Niehrs C, Holstein TW: An ancient Wnt-Dickkopf anatagonism in Hydra. Dev. 2005, 133: 901-911.View ArticleGoogle Scholar
- Lengfeld T, Watanabe H, Simakov O, Lindgens D, Gee L, Law L, Schmidt HA, Özbek S, Bode H, Holstein TW: Multiple Wnts are involved in Hydra organizer formation and regeneration. Dev Biol. 2009, 330: 186-199. 10.1016/j.ydbio.2009.02.004.PubMedView ArticleGoogle Scholar
- Nakamura Y, Tsiairis CD, Özbek S, Holstein TW: Autoregulatory and repressive inputs localize Hydra Wnt3 to the head organizer. Proc Natl Acad Sci USA. 2011, 108: 9137-9142. 10.1073/pnas.1018109108.PubMed CentralPubMedView ArticleGoogle Scholar
- Bunting M: The origin of the sex-cells in Hydractinia and Podocoryne; and the development of Hydractinia. J Morph. 1894, 9: 203-236. 10.1002/jmor.1050090204.View ArticleGoogle Scholar
- Jiménez G, Guichet A, Ephrussi A, Casanova J: Relief of gene repression by Torso RTK signaling: role of capicua in Drosophila terminal and dorsoventral patterning. Genes Dev. 1894, 14: 224-231.Google Scholar
- Monteiro R, van Dinther M, Bakkers J, Wilkinson R, Patient R, ten Diijke P, Mummery C: Two novel type II receptors mediate BMP signaling and are required to establish left-right asymmetry in zebrafish. Dev Biol. 2008, 315: 55-71. 10.1016/j.ydbio.2007.11.038.PubMedView ArticleGoogle Scholar
- Deshpande G, Swanhart L, Chiang P, Schedl P: Hedgehog signaling in germ cell migration. Cell. 2001, 106: 759-769. 10.1016/S0092-8674(01)00488-3.PubMedView ArticleGoogle Scholar
- Bitgood MJ, Shen L, McMahon AP: Sertoli cell signaling by Desert hedgehog regulates the male germline. Current Biol. 1996, 6: 298-304. 10.1016/S0960-9822(02)00480-3.View ArticleGoogle Scholar
- Matus DQ, Magie CR, Pang K, Martindale MQ, Thomsen GH: The Hedgehog gene family of the cnidarian, Nematostella vectensis, and implications for understanding metazoan Hedgehog pathway evolution. Dev Biol. 2008, 313: 501-518. 10.1016/j.ydbio.2007.09.032.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertation, deletions and gene fusions. Genome Biol. 2013, 14: R36-10.1186/gb-2013-14-4-r36.PubMed CentralPubMedView ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotech. 2010, 28: 511-515. 10.1038/nbt.1621.View ArticleGoogle Scholar
- This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GAWH00000000. The version described in this paper is the first version, GAWH01000000
- Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucl Acids Res. 2008, 36: 3420-3435. 10.1093/nar/gkn176.PubMed CentralPubMedView ArticleGoogle Scholar
- Parra G, Bradnam K, Korf I: CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007, 23: 1061-1067. 10.1093/bioinformatics/btm071.PubMedView ArticleGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer ELL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucl Acid Res. 2012, Database Issue: doi:10.1093/nar/gkr1065Google Scholar
- TIGR Database. http://blast.jcvi.org/web-hmm/,
- Min XJ, Butler G, Storms R, Tsang A: OrfPredictor: predicting protein-coding regions in EST-derived sequences. Nucleic Acids Res. 2005, W677-W680. http://proteomics.ysu.edu/tools/OrfPredictor.html, Web Server Issue,
- HMMER. http://hmmer.org/,
- Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of orthologous groups for Eukaryotic genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.PubMed CentralPubMedView ArticleGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9: 35-359.View ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMed CentralPubMedView ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductorpackage for differential expression analysis of digital gene expression data. Bioinformatics. 2010, 26: 139-140. 10.1093/bioinformatics/btp616.PubMed CentralPubMedView ArticleGoogle Scholar
- Rapaport F, Khanin R, Liang Y, Pirun M, Krek A, Zumbo P, Mason CE, Socci ND, Betel D: Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data. Genome Biol. 2013, 14: R95-10.1186/gb-2013-14-9-r95.PubMed CentralPubMedView ArticleGoogle Scholar
- Anders S, McCarthy DJ, Chen Y, Okoniewski M, Smyth GK, Huber W, Robinson MD: Count-based differential expression analysis of RNA sequencing data using R and Bioconductor. Nature Prot. 2013, 8: 1765-1786. 10.1038/nprot.2013.099.View ArticleGoogle Scholar
- Soneson C, Delorenzi M: A comparison of methods for differential expression analysis of RNA-seq data. BMC Bioinformatics. 2013, 14: 91-10.1186/1471-2105-14-91.PubMed CentralPubMedView ArticleGoogle Scholar
- Nawrocki AM, Cartwright P: Expression of Wnt pathway genes in polyps and medusa-like structures of Ectopleura larynx (Cnidaria: Hydrozoa). Evol Dev. 2013, 15: 373-384.PubMedView ArticleGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in the accuracy of multiple sequence alignment. Nucl Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.PubMed CentralPubMedView ArticleGoogle Scholar
- Stamatakis A, Hoover P, Rougemont J: A fast bootstrapping alogorithm for the RAxML web-servers. Sys Biol. 2008, 57: 758-771. 10.1080/10635150802429642.View ArticleGoogle Scholar
- CIPRES Science Gateway. http://www.phylo.org/index.php/portal/,
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.