De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans

De novo transcriptome sequencing and annotation of the Antarctic polychaete Microspio moorei (Spionidae) with its characterization of the heat stress-related proteins (HSP, SOD & CAT)

We present a de novo transcriptome assembly for the non-model Antarctic polychaete worm Microspio moorei (Spionidae) collected during Antarctic field expedition in Fildes Bay, King George Island, Antarctic Peninsula, in 2017. Here, we report the first transcriptome reference array for Microspio spp. The gene sequences of the spionid worm were annotated from a wide range of functions (i.e., biological, and metabolic processes, catalytic processes, and catalytic activity). HSP70, HSP90 SOD and CAT families were compared to reported annelid transcriptomes and proteomes. The phylogenetic analysis using COI, 16S, and 18S markers effectively clusters the species within the family. However, it also casts uncertainty on the monophyletic nature of the Microspio genera, indicating the necessity for additional data and potentially requiring a reevaluation of its grouping. Within these protein families, 3D model software was used to create one representative of their protein structures. Structural predictions were compared with related reported annelids living at different temperatures and a human X-ray reference. We found structural differences (RMSE >1.8) between the human HSP proteins but no significant differences between the polychaete-predicted proteins (RMSE <1.2). These results encourage further research of heat stress-related proteins, the development of genetic markers for climate change-induced temperature stress, and the study of the underlying mechanisms of the heat response. Moreover, these results motivate the extension of these findings to congeneric species.

Droplet digital PCR as a tool for investigating dynamics of cryptic symbionts

Interactions among symbiotic organisms and their hosts are major drivers of ecological and evolutionary processes. Monitoring the infection patterns among natural populations and identifying factors affecting these interactions are critical for understanding symbiont-host relationships. However, many of these interactions remain understudied since the knowledge about the symbiont species is lacking, which hinders the development of appropriate tools. In this study, we developed a digital droplet PCR (ddPCR) assay based on apicomplexan COX1 gene to detect an undescribed agamococcidian symbiont. We show that the method gives precise and reproducible results and enables detecting cryptic symbionts in low target concentration. We further exemplify the assay's use to survey seasonally sampled natural host (Pygospio elegans) populations for symbiont infection dynamics. We found that symbiont prevalence differs spatially but does not show seasonal changes. Infection load differed between populations and was low in spring and significantly increased towards fall in all populations. We also found that the symbiont prevalence is affected by host length and population density. Larger hosts were more likely to be infected, and high host densities were found to have a lower probability of infection. The observed variations could be due to characteristics of both symbiont and host biology, especially the seasonal variation in encounter rates. Our findings show that the developed ddPCR assay is a robust tool for detecting undescribed symbionts that are otherwise difficult to quantify, enabling further insight into the impact cryptic symbionts have on their hosts.

RNA-Seq reveals divergent gene expression between larvae with contrasting trophic modes in the poecilogonous polychaete Boccardia wellingtonensis

The polychaete Boccardia wellingtonensis is a poecilogonous species that produces different larval types. Females may lay Type I capsules, in which only planktotrophic larvae are present, or Type III capsules that contain planktotrophic and adelphophagic larvae as well as nurse eggs. While planktotrophic larvae do not feed during encapsulation, adelphophagic larvae develop by feeding on nurse eggs and on other larvae inside the capsules and hatch at the juvenile stage. Previous works have not found differences in the morphology between the two larval types; thus, the factors explaining contrasting feeding abilities in larvae of this species are still unknown. In this paper, we use a transcriptomic approach to study the cellular and genetic mechanisms underlying the different larval trophic modes of B. wellingtonensis. By using approximately 624 million high-quality reads, we assemble the de novo transcriptome with 133,314 contigs, coding 32,390 putative proteins. We identify 5221 genes that are up-regulated in larval stages compared to their expression in adult individuals. The genetic expression profile differed between larval trophic modes, with genes involved in lipid metabolism and chaetogenesis over expressed in planktotrophic larvae. In contrast, up-regulated genes in adelphophagic larvae were associated with DNA replication and mRNA synthesis.

Comparative transcriptomics in Syllidae (Annelida) indicates that posterior regeneration and regular growth are comparable, while anterior regeneration is a distinct process

Background

Annelids exhibit remarkable postembryonic developmental abilities. Most annelids grow during their whole life by adding segments through the action of a segment addition zone (SAZ) located in front of the pygidium. In addition, they show an outstanding ability to regenerate their bodies. Experimental evidence and field observations show that many annelids are able to regenerate their posterior bodies, while anterior regeneration is often limited or absent. Syllidae, for instance, usually show high abilities of posterior regeneration, although anterior regeneration varies across species. Some syllids are able to partially restore the anterior end, while others regenerate all lost anterior body after bisection. Here, we used comparative transcriptomics to detect changes in the gene expression profiles during anterior regeneration, posterior regeneration and regular growth of two syllid species: Sphaerosyllis hystrix and Syllis gracilis; which exhibit limited and complete anterior regeneration, respectively.

Results

We detected a high number of genes with differential expression: 4771 genes in S. hystrix (limited anterior regeneration) and 1997 genes in S. gracilis (complete anterior regeneration). For both species, the comparative transcriptomic analysis showed that gene expression during posterior regeneration and regular growth was very similar, whereas anterior regeneration was characterized by up-regulation of several genes. Among the up-regulated genes, we identified putative homologs of regeneration-related genes associated to cellular proliferation, nervous system development, establishment of body axis, and stem-cellness; such as rup and JNK (in S. hystrix); and glutamine synthetase, elav, slit, Hox genes, β-catenin and PL10 (in S. gracilis).

Conclusions

Posterior regeneration and regular growth show no significant differences in gene expression in the herein investigated syllids. However, anterior regeneration is associated with a clear change in terms of gene expression in both species. Our comparative transcriptomic analysis was able to detect differential expression of some regeneration-related genes, suggesting that syllids share some features of the regenerative mechanisms already known for other annelids and invertebrates.

Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

BACKGROUND

De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or "pipelines," on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research.

RESULTS

New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla.

CONCLUSIONS

Given current bioinformatics approaches, there is no single "best" reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

Delegating Sex: Differential Gene Expression in Stolonizing Syllids Uncovers the Hormonal Control of Reproduction

Stolonization in syllid annelids is a unique mode of reproduction among animals. During the breeding season, a structure resembling the adult but containing only gametes, called stolon, is formed generally at the posterior end of the animal. When stolons mature, they detach from the adult and gametes are released into the water column. The process is synchronized within each species, and it has been reported to be under environmental and endogenous control, probably via endocrine regulation. To further understand reproduction in syllids and to elucidate the molecular toolkit underlying stolonization, we generated Illumina RNA-seq data from different tissues of reproductive and nonreproductive individuals of Syllis magdalena and characterized gene expression during the stolonization process. Several genes involved in gametogenesis (ovochymase, vitellogenin, testis-specific serine/threonine-kinase), immune response (complement receptor 2), neuronal development (tyrosine-protein kinase Src42A), cell proliferation (alpha-1D adrenergic receptor), and steroid metabolism (hydroxysteroid dehydrogenase 2) were found differentially expressed in the different tissues and conditions analyzed. In addition, our findings suggest that several neurohormones, such as methyl farnesoate, dopamine, and serotonin, might trigger stolon formation, the correct maturation of gametes and the detachment of stolons when gametogenesis ends. The process seems to be under circadian control, as indicated by the expression patterns of r-opsins. Overall, our results shed light into the genes that orchestrate the onset of gamete formation and improve our understanding of how some hormones, previously reported to be involved in reproduction and metamorphosis processes in other invertebrates, seem to also regulate reproduction via stolonization.

The transcriptome of the Bermuda fireworm Odontosyllis enopla (Annelida: Syllidae): A unique luciferase gene family and putative epitoky-related genes

The Bermuda fireworm Odontosyllis enopla exhibits an extremely tight circalunar circadian behavior that results in an impressive bioluminescent mating swarm, thought to be due to a conventional luciferase-mediated oxidation of a light-emitting luciferin. In addition, the four eyes become hypertrophied and heavily pigmented, and the nephridial system is modified to store and release gametes and associated secretions. In an effort to elucidate transcripts related to bioluminescence, circadian or circalunar periodicity, as well as epitoky-related changes of the eyes and nephridial system, we examined the transcriptomic profile of three female O. enopla during a bioluminescent swarm in Ferry Reach, Bermuda. Using the well-characterized luciferase gene of the Japanese syllid Odontosyllis undecimdonta as a reference, a complete best-matching luciferase open reading frame (329 amino acids in length) was found in all three individuals analyzed in addition to numerous other paralogous sequences in this new gene family. No photoproteins were detected. We also recovered a predicted homolog of 4-coumarate-CoA ligase (268 amino acids in length) that best matched luciferase of the firefly Luciola with the best predicted template being the crystal structure of luciferase for Photinus pyralis, the common eastern firefly. A wide variety of genes associated with periodicity were recovered including predicted homologs of clock, bmal1, period, and timeless. Several genes corresponding to putative epitoky-related changes of the eyes were recovered including predicted homologs of a phototransduction gene, a retinol dehydrogenase and carotenoid isomerooxygenase as well as a visual perception related retinal rod rhodopsin-sensitive cGMP 3',5'-cyclic phosphodiesterase. Genes correlating to putative epitoky-related changes of the nephridia included predicted homologs of nephrocystin-3 and an egg-release sex peptide receptor.

Seasonal genetic variation associated with population dynamics of a poecilogonous polychaete worm

Poecilogonous species show variation in developmental mode, with larvae that differ both morphologically and ecologically. The spionid polychaete Pygospio elegans shows variation in developmental mode not only between populations, but also seasonally within populations. We investigated the consequences of this developmental polymorphism on the spatial and seasonal genetic structure of P. elegans at four sites in the Danish Isefjord-Roskilde-Fjord estuary at six time points, from March 2014 until February 2015. We found genetic differentiation between our sampling sites as well as seasonal differentiation at two of the sites. The seasonal genetic shift correlated with the appearance of new size cohorts in the populations. Additionally, we found that the genetic composition of reproductive individuals did not always reflect the genetic composition of the entire sample, indicating that variance in reproductive success among individuals is a likely explanation for the patterns of chaotic genetic patchiness observed during this and previous studies. The heterogeneous, unpredictable character of the estuary might maintain poecilogony in P. elegans as a bet-hedging strategy in the Isefjord-Roskilde-Fjord complex in comparison with other sites where P. elegans are expected to be fixed to a certain mode of development.