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

Coevolution of the Tlx homeobox gene with medusa development (Cnidaria: Medusozoa)

Cnidarians display a wide diversity of life cycles. Among the main cnidarian clades, only Medusozoa possesses a swimming life cycle stage called the medusa, alternating with a benthic polyp stage. The medusa stage was repeatedly lost during medusozoan evolution, notably in the most diverse medusozoan class, Hydrozoa. Here, we show that the presence of the homeobox gene Tlx in Cnidaria is correlated with the presence of the medusa stage, the gene having been lost in clades that ancestrally lack a medusa (anthozoans, endocnidozoans) and in medusozoans that secondarily lost the medusa stage. Our characterization of Tlx expression indicate an upregulation of Tlx during medusa development in three distantly related medusozoans, and spatially restricted expression patterns in developing medusae in two distantly related species, the hydrozoan Podocoryna carnea and the scyphozoan Pelagia noctiluca. These results suggest that Tlx plays a key role in medusa development and that the loss of this gene is likely linked to the repeated loss of the medusa life cycle stage in the evolution of Hydrozoa.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus. Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities similar to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far particular to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus . Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities analogous to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far unique to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Venom system variation and the division of labor in the colonial hydrozoan Hydractinia symbiolongicarpus

Cnidarians (jellyfish, hydroids, sea anemones, and corals) possess a unique method for venom production, maintenance, and deployment through a decentralized system composed of different types of venom-filled stinging structures called nematocysts. In many species, nematocyst types are distributed heterogeneously across functionally distinct tissues. This has led to a prediction that different nematocyst types contain specific venom components. The colonial hydrozoan, Hydractinia symbiolongicarpus, is an ideal system to study the functional distribution of nematocyst types and their venoms, given that they display a division of labor through functionally distinct polyps within the colony. Here, we characterized the composition and distribution of nematocysts (cnidome) in the different polyp types and show that the feeding polyp (gastrozooid) has a distinct cnidome compared to the reproductive (gonozooid) and predatory polyp (dactylozooid). We generated a nematocyst-specific reporter line to track nematocyst development (nematogenesis) in H. symbiolongicarpus, and were able to confirm that nematogenesis primarily occurs in the mid-region of the gastrozooid and throughout stolons (tubes of epithelia that connect the polyps in the colony). This reporter line enabled us to isolate a nematocyst-specific lineage of cells for de novo transcriptome assembly, annotate venom-like genes (VLGs) and determine differential expression (DE) across polyp types. We show that a majority of VLGs are upregulated in gastrozooids, consistent with it being the primary site of active nematogenesis. However, despite gastrozooids producing more nematocysts, we found a number of VLGs significantly upregulated in dactylozooids, suggesting that these VLGs may be important for prey-capture. Our transgenic Hydractinia reporter line provides an opportunity to explore the complex interplay between venom composition, nematocyst diversity, and ecological partitioning in a colonial hydrozoan that displays a division of labor.

•

Functionally specific polyp types in Hydractinia symbiolongicarpus have distinct cnidomes.

•

We present a nematocyst-targeted transgenic line for H. symbiolongicarpus, showcasing active areas of nematogenesis.

•

105 venom-like genes (VLGs) were annotated from an assembled nematocyst-enriched transcriptome.

•

Several VLGs were significantly upregulated in feeding polyps, consistent with being a site of active nematogenesis.

•

Differential expression analysis suggests that different polyp types express distinct combinations of VLGs.

Phylogenetic and Selection Analysis of an Expanded Family of Putatively Pore-Forming Jellyfish Toxins (Cnidaria: Medusozoa)

Many jellyfish species are known to cause a painful sting, but box jellyfish (class Cubozoa) are a well-known danger to humans due to exceptionally potent venoms. Cubozoan toxicity has been attributed to the presence and abundance of cnidarian-specific pore-forming toxins called jellyfish toxins (JFTs), which are highly hemolytic and cardiotoxic. However, JFTs have also been found in other cnidarians outside of Cubozoa, and no comprehensive analysis of their phylogenetic distribution has been conducted to date. Here, we present a thorough annotation of JFTs from 147 cnidarian transcriptomes and document 111 novel putative JFTs from over 20 species within Medusozoa. Phylogenetic analyses show that JFTs form two distinct clades, which we call JFT-1 and JFT-2. JFT-1 includes all known potent cubozoan toxins, as well as hydrozoan and scyphozoan representatives, some of which were derived from medically relevant species. JFT-2 contains primarily uncharacterized JFTs. Although our analyses detected broad purifying selection across JFTs, we found that a subset of cubozoan JFT-1 sequences are influenced by gene-wide episodic positive selection compared with homologous toxins from other taxonomic groups. This suggests that duplication followed by neofunctionalization or subfunctionalization as a potential mechanism for the highly potent venom in cubozoans. Additionally, published RNA-seq data from several medusozoan species indicate that JFTs are differentially expressed, spatially and temporally, between functionally distinct tissues. Overall, our findings suggest a complex evolutionary history of JFTs involving duplication and selection that may have led to functional diversification, including variability in toxin potency and specificity.

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

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

The evolution and development of coloniality in hydrozoans

Hydrozoan colonies display a variety of shapes and sizes including encrusting, upright, and pelagic forms. Phylogenetic patterns reveal a complex evolutionary history of these distinct colony forms, as well as colony loss. Within a species, phenotypic variation in colonies as a response to changing environmental cues and resources has been documented. The patterns of branching of colony specific tissue, called stolons in encrusting colonies and stalks in upright colonies, are likely under the control of signaling mechanisms whose changing expression in evolution and development are responsible for the diversity of hydrozoan colony forms. Although mechanisms of polyp development have been well studied, little research has focused on colony development and patterning. In the few studies that investigated mechanisms governing colony patterning, the Wnt signaling pathway has been implicated. The diversity of colony form, evolutionary patterns, and mechanisms of colony variation in Hydrozoa are reviewed here.

Frizzled3 expression and colony development in hydractiniid hydrozoans

Hydractiniid hydrozoan colonies are comprised of individual polyps connected by tube-like stolons or a sheet-like mat. Mat and stolons function to integrate the colony through continuous epithelia and shared gastrovascular cavity. Although mechanisms of hydrozoan polyp development have been well studied, little is known about the signaling processes governing the patterning of colonies. Here we investigate the Wnt receptor family Frizzled. Phylogenetic analysis reveals that hydrozoans possess four Frizzled orthologs. We find that one of these genes, Frizzled3, shows a spatially restricted expression pattern in colony-specific tissue in two hydractiniid hydrozoans, Hydractinia symbiolongicarpus and Podocoryna carnea, in a manner that corresponds to their distinct colony forms (stolonal mat in Hydractinia and free stolons in Podocoryna). Interestingly, Frizzled3 was lost in the genome of Hydra, which is a solitary polyp and thus lacks colony-specific tissue. Current evidence suggests that the Wnt signaling pathway plays a key role in the evolution of colony diversity and colony loss in Hydrozoa.

Transcriptome Characterization of Reverse Development in Turritopsis dohrnii (Hydrozoa, Cnidaria)

Medusae of Turritopsis dohrnii undergo reverse development in response to physical damage, adverse environmental conditions, or aging. Senescent, weakened or damaged medusae transform into a cluster of poorly differentiated cells (known as the cyst stage), which metamorphose back into a preceding life cycle stage, the polyp. During the metamorphosis, cell transdifferentiation occurs. The cyst represents the intermediate stage between a reverting medusa and a healthy polyp, during which cell transdifferentiation and tissue reorganization take place. Here we characterize and compare the transcriptomes of the polyp and newborn medusa stages of T. dohrnii with that of the cyst, to identify biological networks potentially involved in the reverse development and transdifferentiation processes. The polyp, medusa and cyst of T. dohrnii were sequenced through Illumina RNA-sequencing and assembled using a de novo approach, resulting in 92,569, 74,639 and 86,373 contigs, respectively. The transcriptomes were annotated and comparative analyses among the stages identified biological networks that were significantly over-and under-expressed in the cyst as compared to the polyp and medusa stages. Biological processes that occur at the cyst stage such as telomerase activity, regulation of transposable elements and DNA repair systems, and suppression of cell signaling pathways, mitotic cell division and cellular differentiation and development may be involved in T. dohrnii's reverse development and transdifferentiation. Our results are the first attempt to understand T. dohrnii's life-cycle reversal at the genetic level, and indicate possible avenues of future research on developmental strategies, cell transdifferentiation, and aging using T. dohrnii as a non-traditional in vivo system.

A case of nascent speciation: unique polymorphism of gonophores within hydrozoan Sarsia lovenii

Revealing the mechanisms of life cycle changes is critical for understanding the processes driving hydrozoan evolution. Our analysis of mitochondrial (COI, 16S) and nuclear (ITS1 and ITS2) gene fragments resulted in the discovery of unique polymorphism in the life cycle of Sarsia lovenii from the White Sea. This polymorphic species exhibits two types of gonophores: hydroids produce both free-swimming medusae and attached medusoids (phenotypic polymorphism). Our phylogenetic analysis revealed the intrinsic genetic structure of S. lovenii (genetic polymorphism). Two haplogroups inhabiting the White Sea differ in their reproductive modes. Haplogroup 1 produces attached medusoids, and haplogroup 2 produces free-swimming medusae. Our experiments indicated the possibility of free interbreeding between haplogroups that likely is a rare event in the sea. We propose that inter-haplogroup crossing of S. lovenii in the White Sea may be limited by discordance in periods of spawning or by spatial differences in habitat of spawning specimens. Our finding can be interpreted as a case of nascent speciation that illustrates the patterns of repeated medusa loss in hydrozoan evolution. Life cycle traits of S. lovenii may be useful for elucidating the molecular mechanisms of medusa reduction in hydrozoans.

The genome of the jellyfish Clytia hemisphaerica and the evolution of the cnidarian life-cycle

Jellyfish (medusae) are a distinctive life-cycle stage of medusozoan cnidarians. They are major marine predators, with integrated neurosensory, muscular and organ systems. The genetic foundations of this complex form are largely unknown. We report the draft genome of the hydrozoan jellyfish Clytia hemisphaerica and use multiple transcriptomes to determine gene use across life-cycle stages. Medusa, planula larva and polyp are each characterized by distinct transcriptome signatures reflecting abrupt life-cycle transitions and all deploy a mixture of phylogenetically old and new genes. Medusa-specific transcription factors, including many with bilaterian orthologues, associate with diverse neurosensory structures. Compared to Clytia, the polyp-only hydrozoan Hydra has lost many of the medusa-expressed transcription factors, despite similar overall rates of gene content evolution and sequence evolution. Absence of expression and gene loss among Clytia orthologues of genes patterning the anthozoan aboral pole, secondary axis and endomesoderm support simplification of planulae and polyps in Hydrozoa, including loss of bilateral symmetry. Consequently, although the polyp and planula are generally considered the ancestral cnidarian forms, in Clytia the medusa maximally deploys the ancestral cnidarian-bilaterian transcription factor gene complement.

Comparative transcriptomics reveals potential genes involved in the vegetative growth of Morchella importuna

True morels (Morchella spp.) are edible, medicinal mushrooms which have recently been artificially cultivated in China but stable production remains a problem. Here, we describe complete and comprehensive transcriptome of Morchella importuna at the stages of vegetative mycelium (VM), initial sclerotium (IS) and mature sclerotium (MS) by deep transcriptional sequencing and de novo assembly for the first time and which will potentially provide useful information for improving its cultivation. A total of 26,496 genes were identified with a contig N50 length of 1763 bp and an average length of over 1064 bp. Additionally, 11,957 open reading frames (ORFs) were predicted and 9676 of them (80.9%) were annotated. The 2605 differentially expressed genes (DEGs) identified by gene expression clustering were mainly involved with energy metabolism and could be divided into three broad clusters, of which genes in cluster_1 and cluster_2 were involved in the metabolic process of carbohydrate, polysaccharide, hydrolase, caprolactam, beta-galactosidase, and disaccharide, respectively. Genes in cluster_3 were the largest category, mainly identified with the catalytic activity and transporter activity. Overall, the enzymes involved in the carbohydrate metabolism were highly expressed, and the CAZyme (carbohydrate-active enzyme) genes were significantly expressed within cluster_3. For expression verification, 16 CAZYme genes were selected for qRT-PCR, and the results suggested that the catabolism of carbohydrates occurs mainly in the vegetative mycelium stage, and the anabolism of the energy-rich substances is the main event of mycelial growth and sclerotial morphogenesis of M. importuna.

The genome of the jellyfish Aurelia and the evolution of animal complexity

We present the genome of the moon jellyfish Aurelia, a genome from a cnidarian with a medusa life stage. Our analyses suggest that gene gain and loss in Aurelia is comparable to what has been found in its morphologically simpler relatives-the anthozoan corals and sea anemones. RNA sequencing analysis does not support the hypothesis that taxonomically restricted (orphan) genes play an oversized role in the development of the medusa stage. Instead, genes broadly conserved across animals and eukaryotes play comparable roles throughout the life cycle. All life stages of Aurelia are significantly enriched in the expression of genes that are hypothesized to interact in protein networks found in bilaterian animals. Collectively, our results suggest that increased life cycle complexity in Aurelia does not correlate with an increased number of genes. This leads to two possible evolutionary scenarios: either medusozoans evolved their complex medusa life stage (with concomitant shifts into new ecological niches) primarily by re-working genetic pathways already present in the last common ancestor of cnidarians, or the earliest cnidarians had a medusa life stage, which was subsequently lost in the anthozoans. While we favour the earlier hypothesis, the latter is consistent with growing evidence that many of the earliest animals were more physically complex than previously hypothesized.

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

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

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

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