Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration

Annelid adult cell type diversity and their pluripotent cellular origins

Many annelids can regenerate missing body parts or reproduce asexually, generating all cell types in adult stages. However, the putative adult stem cell populations involved in these processes, and the diversity of cell types generated by them, are still unknown. To address this, we recover 75,218 single cell transcriptomes of the highly regenerative and asexually-reproducing annelid Pristina leidyi. Our results uncover a rich cell type diversity including annelid specific types as well as novel types. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. We find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors, in piwi+ cells. Finally, lineage reconstruction analyses reveal computational differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids by single cell transcriptomics and suggest that a piwi+ cell population with a pluripotent stem cell signature is associated with adult cell type differentiation.

Cell Proliferation Indices in Regenerating Alitta virens (Annelida, Errantia)

In recent years, interest in the possible molecular regulators of cell proliferation and differentiation in a wide range of regeneration models has grown significantly, but the cell kinetics of this process remain largely a mystery. Here we try to elucidate the cellular aspects of regeneration by EdU incorporation in intact and posteriorly amputated annelid Alitta virens using quantitative analysis. We found that the main mechanism of blastema formation in A. virens is local dedifferentiation; mitotically active cells of intact segments do not significantly contribute to the blastemal cellular sources. Amputation-induced proliferation occurred predominantly within the epidermal and intestinal epithelium, as well as wound-adjacent muscle fibers, where clusters of cells at the same stage of the cell cycle were found. The resulting regenerative bud had zones of high proliferative activity and consisted of a heterogeneous population of cells that differed in their anterior-posterior positions and in their cell cycle parameters. The data presented allowed for the quantification of cell proliferation in the context of annelid regeneration for the first time. Regenerative cells showed an unprecedentedly high cycle rate and an exceptionally large growth fraction, making this regeneration model especially valuable for studying coordinated cell cycle entry in vivo in response to injury.

Annelid adult cell type diversity and their pluripotent cellular origins

Annelids are a broadly distributed, highly diverse, economically and environmentally important group of animals. Most species can regenerate missing body parts, and many are able to reproduce asexually. Therefore, many annelids can generate all adult cell types in adult stages. However, the putative adult stem cell populations involved in these processes, as well as the diversity of adult cell types generated by them, are still unknown. Here, we recover 75,218 single cell transcriptomes of Pristina leidyi, a highly regenerative and asexually-reproducing freshwater annelid. We characterise all major annelid adult cell types, and validate many of our observations by HCR in situ hybridisation. Our results uncover complex patterns of regionally expressed genes in the annelid gut, as well as neuronal, muscle and epidermal specific genes. We also characterise annelid-specific cell types such as the chaetal sacs and globin+ cells, and novel cell types of enigmatic affinity, including a vigilin+ cell type, a lumbrokinase+ cell type, and a diverse set of metabolic cells. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. In these piwi+ cells, we also find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors. Finally, lineage reconstruction analyses reveal the existence of differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids for the first time and serve as a resource for studying annelid cell types and their evolution. On the other hand, our characterisation of a piwi+ cell population with a pluripotent stem cell signature will serve as a platform for the study of annelid stem cells and their role in regeneration.

Distinct patterns of gene expression during regeneration and asexual reproduction in the annelid Pristina leidyi

Regeneration, the ability to replace lost body parts, is a widespread phenomenon in the animal kingdom often connected to asexual reproduction or fission, since the only difference between the two appears to be the stimulus that triggers them. Both developmental processes have largely been characterized; however, the molecular toolkit and genetic mechanisms underlying these events remain poorly unexplored. Annelids, in particular the oligochaete Pristina leidyi, provide a good model system to investigate these processes as they show diverse ways to regenerate, and can reproduce asexually through fission under laboratory conditions. Here, we used a comparative transcriptomics approach based on RNA-sequencing and differential gene expression analyses to understand the molecular mechanisms involved in anterior regeneration and asexual reproduction. We found 291 genes upregulated during anterior regeneration, including several regeneration-related genes previously reported in other annelids such as frizzled, paics, and vdra. On the other hand, during asexual reproduction, 130 genes were found upregulated, and unexpectedly, many of them were related to germline development during sexual reproduction. We also found important differences between anterior regeneration and asexual reproduction, with the latter showing a gene expression profile more similar to that of control individuals. Nevertheless, we identified 35 genes that were upregulated in both conditions, many of them related to cell pluripotency, stem cells, and cell proliferation. Overall, our results shed light on the molecular mechanisms that control anterior regeneration and asexual reproduction in annelids and reveal similarities with other animals, suggesting that the genetic machinery controlling these processes is conserved across metazoans.

Journey beyond the embryo: The beauty of Pristina and naidine annelids for studying regeneration and agametic reproduction

Regeneration, asexual agametic reproduction, and other forms of postembryonic development are ecologically important and widely variable across animals, yet our understanding of this important aspect of animal diversity remains limited. A significant limitation has been the dearth of adequate study systems for exploring the mechanisms and evolution of these processes. Here I describe key parts of our journey in developing naid annelids as a study system for investigating the evolution and development of regeneration and fission. Naids are small freshwater annelids that provide numerous advantages for studying postembryonic development: they are small and reproduce readily by fission, they include species with diverse regenerative abilities, and many species are easy to culture. Among the naids, Pristina leidyi is a particularly useful study species, being largely transparent and quite robust to a variety of experimental manipulations. Building on a sparse but long history of past research on these animals, we have developed this system by establishing methods and generating resources for working with them. Naids are yielding novel insights into the evolution of regeneration and fission, providing one of many examples of the value of developing new study species to enable the exploration of fundamental and understudied questions in biology. Establishing new study systems comes with challenges but is exciting and rewarding, and I provide perspectives from my own experiences with the hope of encouraging the further expansion of study systems in biology.

Morphological and molecular features of early regeneration in the marine annelid Ophryotrocha xiamen

Regeneration capability varies in the phylum Annelida making them an excellent group to investigate the differences between closely related organisms. Several studies have described the process of regeneration, while the underlying molecular mechanism remains unclear, especially during the early stage (wound healing and blastema formation). In this study, the newly identified Ophryotrocha xiamen was used to explore the early regeneration. The detailed morphological and molecular analyses positioned O. xiamen within 'labronica' clade. We analyzed the morphological changes during regeneration process (0-3 days post amputation) and molecular changes during the early regeneration stage (1 day post amputation). Wound healing was achieved within one day and a blastema formed one day later. A total of 243 DEGs were mainly involved in metabolism and signal transduction. Currently known regeneration-related genes were identified in O. xiamen which could help with exploring the functions of genes involved in regeneration processes. According to their conserved motif, we identified 8 different Hox gene fragments and Hox5 and Lox2 were found to be absent in early regeneration and during regular growth. Our data can promote further use of O. xiamen which can be used as an experimental model for resolving crucial problems of developmental biology in marine invertebrates.

Effects of GSK3β inhibition in the regeneration of Syllis malaquini (Syllidae, Annelida)

The Wnt/β-catenin signaling pathway has been widely associated to the reestablishment of anteroposterior body polarities in the embryonic development and regeneration in animals. For instance, in annelids, cellular proliferation, wound healing, and blastema development can be affected when this pathway is disrupted. However, very little is known about the genetic regulatory processes involved in these anomalies. Here, we investigate the morphological effects of 1-azakenpaullone, a pharmacological inhibitor of GSK3β that is supposed to over-activate the Wnt/β-catenin pathway, during the anterior and posterior regeneration of the annelid Syllis malaquini. The results showed that high concentrations of 1-azakenpaullone affect the stages of blastema differentiation and resegmentation. Therefore, GSK3β-associated gene regulatory networks are candidate to investigate the genetic mechanisms involved in the regular course of S. malaquini regeneration.

Comparative Aspects of Annelid Regeneration: Towards Understanding the Mechanisms of Regeneration

The question of why animals vary in their ability to regenerate remains one of the most intriguing questions in biology. Annelids are a large and diverse phylum, many members of which are capable of extensive regeneration such as regrowth of a complete head or tail and whole-body regeneration, even from few segments. On the other hand, some representatives of both of the two major annelid clades show very limited tissue regeneration and are completely incapable of segmental regeneration. Here we review experimental and descriptive data on annelid regeneration, obtained at different levels of organization, from data on organs and tissues to intracellular and transcriptomic data. Understanding the variety of the cellular and molecular basis of regeneration in annelids can help one to address important questions about the role of stem/dedifferentiated cells and "molecular morphallaxis" in annelid regeneration as well as the evolution of regeneration in general.

Injury-Induced Innate Immune Response During Segment Regeneration of the Earthworm, Eisenia andrei

Regeneration of body parts and their interaction with the immune response is a poorly understood aspect of earthworm biology. Consequently, we aimed to study the mechanisms of innate immunity during regeneration in Eisenia andrei earthworms. In the course of anterior and posterior regeneration, we documented the kinetical aspects of segment restoration by histochemistry. Cell proliferation peaked at two weeks and remitted by four weeks in regenerating earthworms. Apoptotic cells were present throughout the cell renewal period. Distinct immune cell (e.g., coelomocyte) subsets were accumulated in the newly-formed blastema in the close proximity of the apoptotic area. Regenerating earthworms have decreased pattern recognition receptors (PRRs) (e.g., TLR, except for scavenger receptor) and antimicrobial peptides (AMPs) (e.g., lysenin) mRNA patterns compared to intact earthworms. In contrast, at the protein level, mirroring regulation of lysenins became evident. Experimental coelomocyte depletion caused significantly impaired cell divisions and blastema formation during anterior and posterior regeneration. These obtained novel data allow us to gain insight into the intricate interactions of regeneration and invertebrate innate immunity.

The transcriptome of anterior regeneration in earthworm Eudrilus eugeniae

The oligochaete earthworm, Eudrilus eugeniae is capable of regenerating both anterior and posterior segments. The present study focuses on the transcriptome analysis of earthworm E. eugeniae to identify and functionally annotate the key genes supporting the anterior blastema formation and regulating the anterior regeneration of the worm. The Illumina sequencing generated a total of 91,593,182 raw reads which were assembled into 105,193 contigs using CLC genomics workbench. In total, 40,946 contigs were annotated against the NCBI nr and SwissProt database and among them, 15,702 contigs were assigned to 14,575 GO terms. Besides a total of 9389 contigs were mapped to 416 KEGG biological pathways. The RNA-Seq comparison study identified 10,868 differentially expressed genes (DEGs) and of them, 3986 genes were significantly upregulated in the anterior regenerated blastema tissue samples of the worm. The GO enrichment analysis showed angiogenesis and unfolded protein binding as the top enriched functions and the pathway enrichment analysis denoted TCA cycle as the most significantly enriched pathway associated with the upregulated gene dataset of the worm. The identified DEGs and their function and pathway information can be effectively utilized further to interpret the key cellular, genetic and molecular events associated with the regeneration of the worm.

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.

Large scale changes in the transcriptome of Eisenia fetida during regeneration

Earthworms show a wide spectrum of regenerative potential with certain species like Eisenia fetida capable of regenerating more than two-thirds of their body while other closely related species, such as Paranais litoralis seem to have lost this ability. Earthworms belong to the phylum Annelida, in which the genomes of the marine oligochaete Capitella telata and the freshwater leech Helobdella robusta have been sequenced and studied. Herein, we report the transcriptomic changes in Eisenia fetida (Indian isolate) during regeneration. Following injury, E. fetida regenerates the posterior segments in a time spanning several weeks. We analyzed gene expression changes both in the newly regenerating cells and in the adjacent tissue, at early (15days post amputation), intermediate (20days post amputation) and late (30 days post amputation) by RNAseq based de novo assembly and comparison of transcriptomes. We also generated a draft genome sequence of this terrestrial red worm using short reads and mate-pair reads. An in-depth analysis of the miRNome of the worm showed that many miRNA gene families have undergone extensive duplications. Sox4, a master regulator of TGF-beta mediated epithelial-mesenchymal transition was induced in the newly regenerated tissue. Genes for several proteins such as sialidases and neurotrophins were identified amongst the differentially expressed transcripts. The regeneration of the ventral nerve cord was also accompanied by the induction of nerve growth factor and neurofilament genes. We identified 315 novel differentially expressed transcripts in the transcriptome, that have no homolog in any other species. Surprisingly, 82% of these novel differentially expressed transcripts showed poor potential for coding proteins, suggesting that novel ncRNAs may play a critical role in regeneration of earthworm.

Molecular mechanisms of fission in echinoderms: Transcriptome analysis

Echinoderms are capable of asexual reproduction by fission. An individual divides into parts due to changes in the strength of connective tissue of the body wall. The structure of connective tissue and the mechanisms of variations in its strength in echinoderms remain poorly studied. An analysis of transcriptomes of individuals during the process of fission provides a new opportunity to understand the mechanisms of connective tissue mutability. In the holothurian Cladolabes schmeltzii, we have found a rather complex organization of connective tissue. Transcripts of genes encoding a wide range of structural proteins of extracellular matrix, as well as various proteases and their inhibitors, have been discovered. All these molecules may constitute a part of the mechanism of connective tissue mutability. According to our data, the extracellular matrix of echinoderms is substantially distinguished from that of vertebrates by the lack of elastin, fibronectins, and tenascins. In case of fission, a large number of genes of transcription factors and components of different signaling pathways are expressed. Products of these genes are probably involved in regulation of asexual reproduction, connective tissue mutability, and preparation of tissues for subsequent regeneration. It has been shown that holothurian tensilins are a special group of tissue inhibitors of metalloproteinases, which has formed within the class Holothuroidea and is absent from other echinoderms. Our data can serve a basis for the further study of the mechanisms of extracellular matrix mutability, as well as the mechanisms responsible for asexual reproduction in echinoderms.

Regeneration mechanisms in Syllidae (Annelida)

Syllidae is one of the most species-rich groups within Annelida, with a wide variety of reproductive modes and different regenerative processes. Syllids have striking ability to regenerate their body anteriorly and posteriorly, which in many species is redeployed during sexual (schizogamy) and asexual (fission) reproduction. This review summarizes the available data on regeneration in syllids, covering descriptions of regenerative mechanisms in different species as well as regeneration in relation to reproductive modes. Our survey shows that posterior regeneration is widely distributed in syllids, whereas anterior regeneration is limited in most of the species, excepting those reproducing by fission. The latter reproductive mode is well known for a few species belonging to Autolytinae, Eusyllinae, and Syllinae. Patterns of fission areas have been studied in these animals. Deviations of the regular regeneration pattern or aberrant forms such as bifurcated animals or individuals with multiple heads have been reported for several species. Some of these aberrations show a deviation of the bilateral symmetry and antero-posterior axis, which, interestingly, can also be observed in the regular branching body pattern of some species of syllids.

Developmental and molecular biology of annelid regeneration: a comparative review of recent studies

Studies of annelid regeneration have greatly increased in frequency in recent years, providing new insights into the developmental basis and evolution of regeneration. In this review, we summarize recent findings related to regeneration in annelids, focusing on molecular and developmental studies of epimorphic (blastema-based) regeneration, morphallactic (tissue-remodeling based) regeneration, and development and regeneration of putative stem cells of the posterior growth zone and germline. Regeneration is being investigated in a broad range of annelids spanning the phylum, and comparing findings among species reveals both widely conserved features that may be ancestral for the phylum as well as features that are variable across the group.

A transcriptional blueprint for a spiral-cleaving embryo

Background

The spiral cleavage mode of early development is utilized in over one-third of all animal phyla and generates embryonic cells of different size, position, and fate through a conserved set of stereotypic and invariant asymmetric cell divisions. Despite the widespread use of spiral cleavage, regulatory and molecular features for any spiral-cleaving embryo are largely uncharted. To address this gap we use RNA-sequencing on the spiralian model Platynereis dumerilii to capture and quantify the first complete genome-wide transcriptional landscape of early spiral cleavage.

Results

RNA-sequencing datasets from seven stages in early Platynereis development, from the zygote to the protrochophore, are described here including the de novo assembly and annotation of ~17,200 Platynereis genes. Depth and quality of the RNA-sequencing datasets allow the identification of the temporal onset and level of transcription for each annotated gene, even if the expression is restricted to a single cell. Over 4000 transcripts are maternally contributed and cleared by the end of the early spiral cleavage phase. Small early waves of zygotic expression are followed by major waves of thousands of genes, demarcating the maternal to zygotic transition shortly after the completion of spiral cleavages in this annelid species.

Conclusions

Our comprehensive stage-specific transcriptional analysis of early embryonic stages in Platynereis elucidates the regulatory genome during early spiral embryogenesis and defines the maternal to zygotic transition in Platynereis embryos. This transcriptome assembly provides the first systems-level view of the transcriptional and regulatory landscape for a spiral-cleaving embryo.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2860-6) contains supplementary material, which is available to authorized users.

Gonad establishment during asexual reproduction in the annelid Pristina leidyi

Animals that can reproduce by both asexual agametic reproduction and sexual reproduction must transmit or re-establish their germ line post-embryonically. Although such a dual reproductive mode has evolved repeatedly among animals, how asexually produced individuals establish their germ line remains poorly understood in most groups. We investigated germ line development in the annelid Pristina leidyi, a species that typically reproduces asexually by paratomic fission, intercalating a new tail and head in the middle of the body followed by splitting. We found that in fissioning individuals, gonads occur in anterior segments in the anterior-most individual as well as in new heads forming within fission zones. Homologs of the germ line/multipotency genes piwi, vasa, and nanos are expressed in the gonads, as well as in proliferative tissues including the posterior growth zone, fission zone, and regeneration blastema. In fissioning animals, certain cells on the ventral nerve cord express a homolog of piwi, are abundant near fission zones, and sometimes make contact with gonads. Such cells are typically undetectable near the blastema and posterior growth zone. Time-lapse imaging provides direct evidence that cells on the ventral nerve cord migrate preferentially towards fission zones. Our findings indicate that gonads form routinely in fissioning individuals, that a population of piwi-positive cells on the ventral nerve cord is associated with fission and gonads, and that cells resembling these piwi-positive cells migrate along the ventral nerve cord. We suggest that the piwi-positive ventral cells are germ cells that transmit the germ line across asexually produced individuals via migration along the ventral nerve cord.

De novo assembly and characterization of the fruit transcriptome of Chinese jujube (Ziziphus jujuba Mill.) Using 454 pyrosequencing and the development of novel tri-nucleotide SSR markers

Chinese jujube (Ziziphus jujuba Mill.) is an economically important deciduous tree that has high therapeutic value and health benefits. However, a lack of sequence data and molecular markers have constrained genetic and breeding studies for better fruit quality and other traits in Chinese jujube. In this study, two combined cDNA libraries of ‘Dongzao’ fruit representing the early and late stages of fruit development were constructed and sequenced on the 454 GS FLX Titanium platform. In total, 1,124,197 reads were generated and then de novo assembled into 97,479 unigenes. A total of 52,938 unigenes were homologous to genes in the NCBI non-redundant sequence database. A total of 33,123 unigenes were assigned to one or more Gene Ontology terms, and 16,693 unigenes were classified into 319 Kyoto Encyclopedia of Genes and Genomes pathways. The results showed that the Smirnoff-Wheeler pathway was the main pathway for the biosynthesis of ascorbic acid in Chinese jujube. The number of differentially expressed genes between the two stages of fruit development was 1,764, among which 974 and 790 genes were up-regulated and down-regulated, respectively. Furthermore, 9,893 sequences were identified containing SSRs. 93 primer pairs designed from the sequences with a tri-nucleotide repeat showed successful PCR amplification and could be validated in Chinese jujube accessions and Z. mauritiana Lam and Z. acidojujuba as well, of which 71 primer pairs were polymorphic. The obtained transcriptome provides a most comprehensive resource currently available for gene discovery and the development of functional markers in Z. jujuba. The newly developed microsatellite markers could be used in applications such as genetic linkage analysis and association studies, diversity analysis, and marker-assisted selection in Chinese jujube and related species.

Characterization of Liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis

Background

Liaoning cashmere goat is a famous goat breed for cashmere wool. In order to increase the transcriptome data and accelerate genetic improvement for this breed, we performed denovo transcriptome sequencing to generate the first expressed sequence tag dataset for the Liaoning cashmere goat, using next-generation sequencing technology.

Results

Transcriptome sequencing of Liaoning cashmere goat on a Roche 454 platform yielded 804,601 high-quality reads. Clustering and assembly of these reads produced a non-redundant set of 117,854 unigenes, comprising 13,194 isotigs and 104,660 singletons. Based on similarity searches with known proteins, 17,356 unigenes were assigned to 6,700 GO categories, and the terms were summarized into three main GO categories and 59 sub-categories. 3,548 and 46,778 unigenes had significant similarity to existing sequences in the KEGG and COG databases, respectively. Comparative analysis revealed that 42,254 unigenes were aligned to 17,532 different sequences in NCBI non-redundant nucleotide databases. 97,236 (82.51%) unigenes were mapped to the 30 goat chromosomes. 35,551 (30.17%) unigenes were matched to 11,438 reported goat protein-coding genes. The remaining non-matched unigenes were further compared with cattle and human reference genes, 67 putative new goat genes were discovered. Additionally, 2,781 potential simple sequence repeats were initially identified from all unigenes.

Conclusion

The transcriptome of Liaoning cashmere goat was deep sequenced, denovo assembled, and annotated, providing abundant data to better understand the Liaoning cashmere goat transcriptome. The potential simple sequence repeats provide a material basis for future genetic linkage and quantitative trait loci analyses.