Differential gene expression between functionally specialized polyps of the colonial hydrozoan Hydractinia symbiolongicarpus (Phylum Cnidaria)
BMC Genomics volume 15, Article number: 406 (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.
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].
The main labor divisions are between feeding, reproduction, and defense, where specialized polyp types are morphologically distinct, reflecting their particular functions. Hydractinia symbiolongicarpus has four different polyp types (Figure 1). The feeding polyp (called a gastrozooid) has a mouth and tentacles, which are absent in the reproductive polyp (gonozooid), defensive and food gathering polyp (dactylozooid), and the less common defensive polyp (tentaculozooid, not shown). The dactylozooid has an elongated body column with an abundance of epithelial muscular cells and nematocytes (stinging cells). The gonozooid bears gonophores, which house the gametes. The gonozooid and dactylozooid are similar in their distal ends, with clusters of nematocysts and lacking a functional mouth and elongate tentacles. The tentaculozooid resembles a tentacle of the gastrozooid, but is the size of an individual polyp.
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
From the three normalized libraries, 49,759, 43,776, and 142,408 contigs were assembled for the gastrozooid, gonozooid, and dactylozooid, respectively. Individual transcriptomes were merged into a single assembly of 101,518 unique transcripts using cuffmerge (Figure 2). Cuffmerge merges novel and common transcripts into a single assembly and removes artifact constructions, improving the overall quality of the assembly. This step allows for easy annotation and differential expression analyses of a single assembly, without concerns regarding orthology assignments between multiple assemblies. After filtering for transcripts less than 200 bp in length, our final assembly consisted of 66,508 transcripts, with an N50 of 1,451 bp (Additional file 1).
Approximately 54% of the transcriptome (35,636 transcripts) was annotated using Blast2GO, CEGMA, orthoMCL, and HMMscan (Additional file 2), with these transcripts showing significant similarity to sequences in at least one database in our annotation pipeline (Figures 2 and 3). These include 416 (91%) of the “core” and 238 (96%) of the “ultra-conserved” eukaryotic genes identified using CEGMA (Additional file 3). Figure 3 shows the number of transcripts annotated by one or more of the annotation methods.
Differential expression analyses
Statistically significant differences in expression between different polyp types were detected using two DE packages, DESeq and edgeR. Figure 4C and D shows the effect of heterogametic adjustments on the Euclidean distances (sum of the pairwise distance across all transcripts) between libraries. Both DESeq and edgeR reveal that dactylozooids and gastrozooids share the fewest number of DE transcripts and the smallest change in the number of DEs recovered after the heterogametic adjustments (Table 1), while DE analyses including gonozooids show a much larger increase in the number of DE transcripts after those adjustments. This large increase can be explained by the huge amount of variability found when ignoring the sexual differences between gametic tissues in gonozooid samples. The DE analysis between the male and female gonozooid libraries identified 11,798 (DESeq) and 12,886 (edgeR) transcripts significantly up- or down-regulated (Table 1, Figure 4B). Removal of all male/female DE transcripts clusters gonozooid samples by polyp type rather than sex and increases the distance between gastrozooids and dactylozooids (Figure 4D), while treating male and female gonozooid libraries as different conditions reduces the average dispersion estimate for each transcript, essentially increasing the power of the DE analyses (Figure 4E and F). Yet, even after the heterogametic adjustments, gonozooids still have the largest number of polyp-specific transcripts (transcripts that are strictly up- or down-regulated in a particular polyp type when compared to other polyps; Figure 5, Table 2, Additional file 4).
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].
Figure 6 shows whole mount in situ hybridization (ISH) results of several polyp-specific transcripts identified through the DE analyses (listed in Table 3). DE analyses reported several different toxin transcripts to be differentially expressed between the different polyps. Polyp specificity of one of the three toxins identified as gastrozooid-specific by DE analyses, referred to here as toxin_5320, was confirmed by ISH. This transcript was expressed solely in specific gastrodermal cells around the base of the hypostome/tentacle margin of the gastrozooids. Three distinct cell types populate the gastrodermis of the hypostome in Hydractinia: gastrodermal epithelia (including digestive cells) and two glandular cell types (spumeous and spherulous cells) [40–42]. Toxin_5320 expression appears to be limited to the spherulous cells of the hypostome (Figure 6, Additional file 10). DE analyses found toxin_3875 to be dactylozooid-specific and ISH found expression to be limited to nematocytes primarily found in the proximal portion of the body column of the dactylozooid (Figure 6, Additional file 10).
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.
In an effort to distinguish between gametogenic-specific expression and expression specific to gonozooid polyp identity, several preliminary DE analyses were conducted to adjust for gametogenic expression (Figure 7). First, a DE analysis was conducted between male and female gonozooid libraries (step 1, Figure 7), identifying significantly (padj < 0.05) up- or down-regulated transcripts (step 2, Figure 7). Second, transcripts found to be significant were excluded from the template pool; a DE analysis was then performed on non-significant transcripts that included counts from both male and female libraries (step 3 shown in black, Figure 7). However, this analysis excluded maternal transcripts that could also play a role in somatic morphogenesis. To include these maternal transcripts, a second DE analysis was conducted on transcripts up-regulated in female gonozooids (putative maternal transcripts). For these comparisons, only the expression counts from the male gonozooid libraries were used (step 3 shown in red, Figure 7). In this approach, expression patterns consistent with developmental patterning and functional specialization were less likely obscured by the expression of genes specific to gametogenesis in the DE analyses. Results from both analyses (step 3, Figure 7) were combined.
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.
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.
Berrill NJ: Growth and form in gymnoblastic hydroids; polymorphism within the Hydractiniidae. J Morph. 1953, 92: 241-272. 10.1002/jmor.1050920203.
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.
Cartwright P, Nawrocki AM: Character evolution in Hydrozoa (phylum Cnidaria). J Integ Comp Biol. 2010, 50: 456-472. 10.1093/icb/icq089.
Agassiz A: North American Acalephae. 1865, Welsh, Bigelow, Cambridge, U.K: University Press
Haeckel E: Report on the Siphonophorae Collected by H. M. S. Challenger During the Years 1873–76. 1888, London: The Challenger Reports
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 Society
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.
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.
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.
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.
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.
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.
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.127233
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.
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.
Bebenek HG, Gates RD, Morris J, Hartenstein V, Jacobs DK: sine oculis in basal Metazoa. Dev Genes Evol. 2004, 214: 342-351.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.x
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.
Boero F, Sará M: Motile sexual stages and evolution of Leptomedusae (Cnidaria). Boll Zool. 1987, 54: 131-139. 10.1080/11250008709355572.
Frazão B, Vasconcelos V, Antunes A: Sea anemone (Cnidaria, Anthozoa, Actiniaria) toxins: an overview. Mar Drugs. 2012, 10: 1812-1851. 10.3390/md10081812.
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.
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.
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.
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.
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.
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.
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.
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.
Hazen AP: Regeneration in Hydractinia and Podocoryne. Amer Nat. 1902, 36: 193-200. 10.1086/278099.
Campbell RD: Cell proliferation and morphological patterns in the hydroids Tubularia and Hydractinia. J Embryol Exp Morph. 1967, 17: 607-616.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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/gkr1065
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,
Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of orthologous groups for Eukaryotic genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.
Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9: 35-359.
Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.
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.
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.
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.
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.
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.
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.
Stamatakis A, Hoover P, Rougemont J: A fast bootstrapping alogorithm for the RAxML web-servers. Sys Biol. 2008, 57: 758-771. 10.1080/10635150802429642.
CIPRES Science Gateway. http://www.phylo.org/index.php/portal/,
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.
The authors declare that they have no competing interests.
SMS performed tissue collections for the 12 non-normalized libraries, assembled and annotated transcriptome, designed the heterogametic adjustment algorithm, performed DE analyses and ISH, and drafted the manuscript. MS performed tissue collection, RNA isolation, normalizations, and library construction of the 3 normalized libraries. MS also wrote the python script that extracts counts from the .sam files. PC oversaw the design and implementation of the study and helped draft the manuscript. All authors have read and approved of the final manuscript.
Electronic supplementary material
Additional file 1: Histogram of the size distribution of assembled transcripts. This does not include transcripts that were removed because they were < 200 bp in length. Inset table displays assembly numbers and size statistics before and after filtering out the <200 bp transcripts. (PDF 145 KB)
Additional file 4: All polyp-specific transcripts. This list includes the assembly sequence ID, top BLASTX hit, number of gene ontology IDs, top HMM protein domain, polyp specificity, top significance threshold, and transcript sequence. Polyp specificity is defined in two separate columns, ‘Polyp’ and ‘Direction’ (example: ‘Gono, DOWN, .05,’ would be a transcript that is significantly down regulated in the gonozooid when compared to the other two polyps at a significance level between .01 and .05). (CSV 9 MB)
Additional file 5: List of polyp-specific DEs discussed in Results and discussion section. This table is a subset of Additional file 4 and contains all the same information, but with BS support and the accession number of the closest cnidarian ortholog, if molecular phylogenetic analyses performed. In cases where phylogenetic analyses recovered well supported sister relationships between two monophyloetic clades with no clear one to one orthologous relationships, a single accession number was selected from the non-H. symbiolongicarpus monophyly. List is divided up by color according to the subheading in the Results and discussion section they are discussed in: Blue – ‘Gametogenic expression’; Yellow – ‘Homeodomains’; Pink – ‘Myosins’; Green – ‘Toxins’; Orange – ‘Astacins’. (*) Marks transcripts with whole mount in situ hybridizations data in this study. (XLSX 152 KB)
Additional file 6: Cnidarian homeodomain gene tree. Molecular phylogeny of cnidarian homeodomains sampled from GenBank’s nr database. Accession numbers are appended to the ends of the tip labels. Only polyp-specific homeodomains from H. symbiolongicarpus (highlighted in red) were included in the analysis. Fasta and alignment file available upon request. (PDF 12 KB)
Additional file 7: Cnidarian myosin gene tree. Molecular phylogeny of cnidarian myosins sampled from GenBank’s nr database. Accession numbers are appended to the ends of the tip labels. Only polyp-specific homeodomains from H. symbiolongicarpus (highlighted in red) were included in the analysis. Fasta and alignment file available upon request. (PDF 4 KB)
Additional file 8: Cnidarian toxin gene tree. Molecular phylogeny of cnidarian toxins sampled from GenBank’s nr database. Accession numbers are appended to the ends of the tip labels. Only polyp-specific homeodomains from H. symbiolongicarpus (highlighted in red) were included in the analysis. Fasta and alignment file available upon request. (PDF 4 KB)
Additional file 9: Cnidarian astacin gene tree. Molecular phylogeny of cnidarian astacins sampled from GenBank’s nr database. Accession numbers are appended to the ends of the tip labels. Only polyp-specific astacins from H. symbiolongicarpus (highlighted in red) were included in the analysis. Fasta and alignment file available upon request. (PDF 5 KB)
Additional file 10: In situ hybridization (higher magnification) of toxins. A. toxin_5320. B. toxin_3875. sp = spumeous cells; ns = nematocyst; nc = nematocyte. (PDF 5 MB)
Additional file 11: Python script for extracting count data from .sam files. Save as counts-paired.py. Use: python counts-paired.py infile.sam.(CSV 2 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Sanders, S.M., Shcheglovitova, M. & Cartwright, P. Differential gene expression between functionally specialized polyps of the colonial hydrozoan Hydractinia symbiolongicarpus (Phylum Cnidaria). BMC Genomics 15, 406 (2014). https://doi.org/10.1186/1471-2164-15-406