Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration
© Nyberg et al.; licensee BioMed Central Ltd. 2012
Received: 23 February 2012
Accepted: 29 June 2012
Published: 29 June 2012
The naid annelids contain a number of species that vary in their ability to regenerate lost body parts, making them excellent candidates for evolution of regeneration studies. However, scant sequence data exists to facilitate such studies. We constructed a cDNA library from the naid Pristina leidyi, a species that is highly regenerative and also reproduces asexually by fission, using material from a range of regeneration and fission stages for our library. We then sequenced the transcriptome of P. leidyi using 454 technology.
454 sequencing produced 1,550,174 reads with an average read length of 376 nucleotides. Assembly of 454 sequence reads resulted in 64,522 isogroups and 46,679 singletons for a total of 111,201 unigenes in this transcriptome. We estimate that over 95% of the transcripts in our library are present in our transcriptome. 17.7% of isogroups had significant BLAST hits to the UniProt database and these include putative homologs of a number of genes relevant to regeneration research. Although many sequences are incomplete, the mean sequence length of transcripts (isotigs) is 707 nucleotides. Thus, many sequences are large enough to be immediately useful for downstream applications such as gene expression analyses. Using in situ hybridization, we show that two Wnt/β-catenin pathway genes (homologs of frizzled and β-catenin) present in our transcriptome are expressed in the regeneration blastema of P. leidyi, demonstrating the usefulness of this resource for regeneration research.
454 sequencing is a rapid and efficient approach for identifying large numbers of genes in an organism that lacks a sequenced genome. This transcriptome dataset will be a valuable resource for molecular analyses of regeneration in P. leidyi and will serve as a starting point for comparisons to non-regenerating naids. It also contributes significantly to the still limited genomic resources available for annelids and lophotrochozoans more generally.
KeywordsTranscriptome 454 pyrosequencing Regeneration Evolution Asexual reproduction Fission Annelid Pristina leidyi Cell signaling Wnt signaling
The process of regeneration, or the replacement of lost body parts, has long captured the interest of biologists. While early experiments on crayfish  and Hydra demonstrated the remarkable abilities of some animals to develop lost parts anew, it is also clear that many animals, including humans, do not possess such abilities. The ability to regenerate is thought to have been lost over the course of evolution in many animal lineages [3–5]. Despite recent advances in knowledge of the molecular and developmental basis of regeneration in a variety of animal systems [6–8], little is currently known about the developmental and evolutionary mechanisms that drive loss of regeneration ability . Understanding this phenomenon requires a comparative approach and the identification and development of animal systems that show variation in regeneration ability among closely related species.
The naid annelids are among a small number of documented groups in which regeneration ability varies among close relatives [5, 9–14], making them a good model for studying the loss of regeneration. Naids (the minimal clade including both the Naidinae and Pristininae) are a group of small aquatic oligochaete worms, many of which can reproduce asexually by fission . Many naids, including Pristina leidyi, possess excellent regeneration abilities, being able to regrow both their heads and tails after amputation. Following amputation, tissues at the wound site actively proliferate and form a regeneration blastema (a mass of undifferentiated cells) which ultimately differentiates to give rise to regenerated structures . The ability to regenerate anteriorly and posteriorly is thought to be ancestral for the clade. However, recent experiments indicate that head regeneration ability has been lost at least three times within the naids, allowing multiple independent comparisons between regenerating and non-regenerating species [9, 10, 17]. The degree of loss of the regeneration machinery can vary between lineages, suggesting that different developmental mechanisms may underlie independent evolutionary losses of regeneration . Thus, in the naids, evolution has crafted an ideal experiment for investigating loss of regeneration.
Much recent work on the developmental basis of regeneration has focused on the role of signaling pathways, such as the Wnt pathway, in recruiting stem cells and promoting morphogenesis in regenerated tissues [18–29]. In order to investigate the role of signaling pathways and other molecules in variation of regeneration ability, genomic resources are needed for the naids. Recent advances in high-throughput sequencing and bioinformatic analyses have made transcriptome sequencing feasible for discovering novel genes in non-model systems.
454 pyrosequencing, with sequence reads now approaching the length of traditional Sanger sequences, is ideal for transcriptome sequencing in a model that lacks a sequenced genome [30, 31]. While the sequencing depth of 454 is modest compared to that of other deep sequencing technologies, 454 does offer depth orders of magnitude above what can be obtained via Sanger sequencing . In addition, recent versions of the Newbler assembler from 454 allow for assembling sequences from cDNA, grouping presumptive gene isoforms together. Here, we describe the sequencing and assembly of a full run of 454 GS FLX sequencing with Titanium reagents from the regenerating annelid P. leidyi.
Results and discussion
Genome size estimation
We estimated the genome size of P. leidyi and four other naid species currently used in comparative regeneration studies . Using Feulgen densitometry analysis, we estimated a C-value of 1.37 pg for P. leidyi and C-values ranging from 0.54 to 1.09 pg for the four other naid species (Additional file 1). Previously published estimates for two naid species are 1.53 and 3.23 pg, and the mean of values reported for oligochaetes is ~1.6 pg (range: 0.43 to 7.64) [33, 34]. Thus, the genome size of P. leidyi is typical for this group.
Construction of a partially normalized cDNA library
RNA was extracted from whole worms at multiple time points between the initiation of regeneration and its completion and from unamputated worms that were actively undergoing fission. cDNA was synthesized with an oligo-dT primer and a MINT full-length reverse-transcription kit. PCR assays indicated that the dried Spirulina powder used as food was not metabolically active and could not be detected by RT-PCR in the cDNA sample (Additional file 2).
The overall amount of cDNA is greatly reduced during normalization, making PCR amplification necessary to produce a sufficient quantity of cDNA for 454 sequencing. Because PCR has its own biases, particularly against large amplicons, we pooled an unamplified, non-normalized sample with a normalized sample in a 1:2 ratio to increase the representation of longer transcripts (Figure 1). This pooled cDNA library was used for 454 pyrosequencing.
454 pyrosequencing and transcriptome assembly
Sequence and assembly output
454 sequencing results
mean read length
Newbler Assembler v2.3 assembly results
mean isotig length
mean max isotig length
fraction captured transcripts
gene discovery rate
33.21 reads/new gene
454 sequence reads were assembled using the Newbler Assembler v2.3 . Sequence output using the cDNA option of Newbler v2.3 differs from that of traditional genomic assemblers (e.g. SeqMan NGen 2.0, CAP3, Newbler 2.2) by taking into account the possibility that multiple isoforms (e.g. alternative splice variants) of a gene may be present. Overlapping sequence reads are assembled into contigs, much like traditional assemblers. However, if multiple isoforms are present, a sequence read may contain a portion that aligns perfectly with the previously constructed contig and a portion that does not (with the point of divergence being, for example, an exon-exon junction). When this occurs, Newbler v2.3 breaks up the aligned sequences into multiple contigs. Sequences shared between multiple isoforms are retained as unique contigs, and any adjacent variant sequences are split off as their own unique contigs. Thus, a single gene isoform might be assembled into multiple contigs, and the same contig might be shared across multiple isoforms. Each putative isoform identified by Newbler v2.3 is termed an isotig, and the multiple isoforms for each gene are organized into isogroups, representing putative gene loci.
BLAST analysis of 454 isotigs
Annotation of isotigs
EST2uni BLAST results
64,522 (max isotig per isogroup)
11,388 (17.7%, E-value < e-10)
Taxonomic identity of top UniProt hit
Total BLAST hits
Isotig matches to UniProt database
UniProt ref. size
5’ + 3’
BLAST results suggest that our efforts to minimize environmental contamination were successful. BLAST searches against the P. leidyi isotig dataset using either a cnidarian or human 16 S sequence [Hydra magnipapillata: GenBank:NC_011220|:307-2044; Homo sapiens: GenBank: FJ794693.1|:1673-3230] retrieve isotigs matching to Pristina 16S as the only hits with any reasonable significance (E-value < 0.1), suggesting no metazoan contamination. Furthermore, only a small number of isotigs matched prokaryotic genes (1.4%) (Table 2) and BLAST searches of the P. leidyi isotigs using bacterial 16S RNA from either the proteobacterium Escherichia coli [GenBank:4924485] or the cyanobacterium Arthrospira platensis (Spirulina) [GenBank:FJ798612.1] return a very limited number of isotigs (only nine, and the same nine, isotigs for both searches; E-values < e-2). Interestingly, four of these isotigs match 16S from bacterial genera known to be common endosymbionts in animal intestines (gammaproteobacteria Edwardsiella/Xenorhabdus/Photorhabdus; bacteroidetes Paenicardinium/Cardinium). Thus, some bacterial sequences present in the transcriptome may represent the endemic gut flora of P. leidyi.
To assess how well 5’ and 3’ ends were captured in our dataset, isotigs with significant BLAST hits were compared to their counterparts in UniProt (Table 3). The proportion of isotigs with captured ends (within 10 amino acids of the corresponding end of the UniProt sequence) varies with the length of the coding sequence in UniProt. Isotigs matching shorter UniProt sequences are more likely to be complete on both the 5’ and 3’ ends than isotigs matching longer UniProt sequences. 31.7% of isotigs matching UniProt sequences of less than 250 amino acids are complete on both ends, while only a single isotig matching UniProt sequences greater than 750 amino acids is complete. In total, 6.5% of isotigs can be considered complete, 13.0% have captured the 5’ end, 16.1% have captured the 3’ end, and 64.3% have no matches against either end. While we estimate that we have captured the vast majority of transcripts in the original cDNA library (Table 1) and that there is no strong bias towards either 5’ or 3’ ends, it is clear that most of our unigenes consist of only partial transcript sequences. Further sequencing, either in a high-throughput manner or on a targeted basis with genes of interest, will be necessary to fill in these gaps.
Gene ontology analysis of 454 isotigs
Identification of candidate regeneration genes
BLAST results for candidate regeneration genes
07867, 10773, 21383, 44961, 47307
03233, 13291, 15038, 40749, 49382
03101, 07423, 11748, 64463
00856, 52001, 54688
Independent confirmation of assembled transcripts
The utility of this transcriptome will ultimately be determined by whether these assembled sequences can be independently confirmed and manipulated for further studies. Because 454 sequencing has a nebulization step and is not performed using intact transcripts, Newbler v2.3 is unable to reconstruct with complete accuracy the actual gene isoforms present in vivo. Therefore, isotigs should be treated as predicted gene isoforms that require independent confirmation, such as by PCR assay. Contigs, on the other hand, already are well supported via the original sequencing and are expected to be true contiguous sequence and thus amplifiable by PCR.
Expression of wnt/β-catenin pathway genes during regeneration
To further demonstrate the utility of our sequencing effort for regeneration studies, we examined the expression patterns of two genes present in the transcriptome, homologs of frizzled (fz) and β-catenin (β-cat). These genes were chosen because they are components of the Wnt/β-catenin pathway, an important cell signaling pathway implicated in numerous developmental processes, including regeneration [18–29]. We identified from our transcriptome several homologs of fz (a Wnt ligand receptor) and a single homolog of β-cat (a multifunctional protein that acts as a transcription factor when Wnt signaling is activated). We examined expression of one of these fz homologs, Pl-fzA, and the homolog of β-cat Pl-β-cat, by whole mount in situ hybridization of regenerating and fissioning P. leidyi.
Our results provide the first expression data for frizzled or β-catenin genes during annelid regeneration and strongly implicate Wnt signaling in P. leidyi regeneration. They also add to the accumulating data showing a close developmental relationship between regeneration and fission in these animals. More broadly, these findings demonstrate that sequences from this transcriptome can provide new insights into annelid development, setting the stage for future comparative studies of annelid regeneration.
Applications for further regeneration research
This transcriptome dataset provides a valuable new resource for regeneration research in annelids. The work described here already demonstrates the utility of this dataset: our GO analysis suggests that a large number of genes relevant for regeneration research are represented, our sequence confirmation assays show that putative transcripts from our assembly can be independently validated, and our expression studies demonstrate that genes expressed during regeneration are indeed present in this transcriptome and can provide new insight into annelid regeneration.
This transcriptome resource promises to accelerate regeneration research in P. leidyi and provides a stepping-stone to studies of regeneration failure in closely related, non-regenerating naid species. This dataset greatly facilitates gene discovery, allowing genes of interest to be quickly identified and characterized by RT-PCR or in situ hybridization. Importantly, this resource can also provide a reference transcriptome for larger, genome-scale studies, such as high-throughput analyses of gene expression by microarrays or RNA-Seq [78–80]. Extending these approaches to closely related regenerating and non-regenerating naid species holds great promise for elucidating the genetic basis of both regeneration success and failure.
This transcriptome sequencing project has produced the first genomic-type sequence data for any of the naid annelids, a promising group for understanding regeneration loss. This dataset also represents the first regeneration-based, large-scale sequence database for annelids as a whole and thus provides a valuable resource for regeneration research more broadly. Our approach of using mixed-stage starting material and combining normalized/non-normalized cDNA pools was successful in producing a transcriptome with high gene representation. Based on BLAST searches for known regeneration genes and relevant GO analyses matches, we conclude that our methods captured a significant number of genes that may be involved in regeneration. This transcriptome resource enabled gene expression studies that have provided novel insight into annelid regeneration, yielding the first evidence suggesting that a cell-signaling pathway important in other regenerating systems, Wnt/β-catenin signaling, is initiated during annelid regeneration. Thus, this dataset promises to be instrumental in determining which genes are involved in regeneration processes in P. leidyi and will subsequently inform evolution of regeneration studies in the naids as a whole.
The development of genomic resources for the Lophotrochozoa (the large clade of bilaterians including platyhelminths, annelids, and molluscs, among others) has lagged considerably behind that of other major groups of animals. With the advent of less expensive sequencing technologies and an increased appreciation of the value of non-traditional model systems, genomic resources for this group are finally becoming available. Transcriptomes have recently been generated for a range of lophotrochozoan taxa [81–87], including the highly regenerative planarians [88, 89]. The growing number of genomic resources for lophotrochozoans promises to help fuel research on a broad range of questions in this large and diverse clade.
Genome size measurement
The genome sizes of five naid species were estimated using the Feulgen image analysis densitometry method . Individuals were obtained from laboratory cultures of the following species: P. leidyi (Carolina Biological Supply Company), Allonais paraguayensis (Wards Natural Science), Dero digitata (originally collected from Edwards Lake, University of Maryland at College Park, USA), Dero furcata (Connecticut Valley Biological Supply), and Paranais litoralis (originally collected from Herrington Bay, Deale, MD, USA). Fifty or more nuclei were measured from each sample. The Integrated Optical Density of the sample was converted to a genome size value (in picograms) using Gallus gallus domesticus (1.25 pg) as a standard.
Worm culture, sampling, and RNA extraction
To generate material for this sequencing effort, we established twelve replicate lab cultures of a single clonal line of Pristina leidyi (PRIle(cbs)cloneA). Each culture was initially started with 100 worms and was maintained at room temperature in 20 cm glass bowls filled with ~1 liter of commercially purchased Poland Spring Water (PSW). To ensure purity of the samples, worm cultures were rinsed frequently to remove algae and debris and cultures were routinely inspected visually for the presence of small metazoans (e.g. rotifers). Worms were fed dried Spirulina powder twice weekly, and water was changed at least 1-2 times per week.
Possible contamination by the dried Spirulina food source was assayed via PCR, using cDNA samples derived from live Arthrospira platensis (Spirulina) as a reference. RNA was extracted using TRIReagent (Applied Biosystems), and cDNA was constructed using random oligos and Superscript III reverse transcriptase (Invitrogen). Primers were constructed against the large subunit of rubisco (rbcL) [GenBank:AY147205.1] and c-phycocyanin (cpc) [GenBank:AF164139.1] genes of A. platensis (Additional file 5).
Worms were collected at a range of stages of regeneration and fission (Figure 1). For the fission material, 1,000 worms that were actively growing and fissioning were collected and starved for 24 hours. To generate regenerating material, 3,485 worms were amputated anteriorly and posteriorly and allowed to regenerate for various lengths of time before collection. Most worms were actively fissioning and consisted of chains of linked zooids at time of amputation. A cut was made 2 body segments anterior to the most anterior fission zone to elicit posterior regeneration. A second cut was made after the 6th body segment of the most posterior zooid to elicit anterior regeneration. If a worm did not consist of at least two nearly formed zooids, a single cut after the 6th body segment was made to elicit anterior regeneration. Because the initiation of regeneration processes holds particular significance for future studies, 1,985 of the regenerating worms in the sample were allowed to regenerate between 0 and 24 hours, which is roughly coincident with the start of blastema formation. Batches of 250 worms were also collected at 1.25 days post-amputation (dpa), 1.75 dpa, 2 dpa, 2.5 dpa, 3 dpa, and 3.5 dpa, when differentiation of adult morphology is nearly complete.
Fissioning and regenerating worms were washed 5x in PSW prior to RNA extraction. RNA was extracted using TRIReagent (Applied Biosystems), and RNA from all samples was then pooled together.
cDNA library construction
We constructed a pooled cDNA library consisting of a normalized fraction to capture lowly expressed transcripts and a non-normalized fraction to capture large transcripts that might be lost during the PCR amplification steps of the normalization process.
First-strand cDNA (F.S. cDNA) was made using a MINT full-length cDNA synthesis kit (Evrogen) and manufacturer’s instructions. A modified oligo-dT primer with breaks in the homopolymer-T run was used to minimize the negative effects of an extensive homopolymer run on 454 sequence quality (Additional file 5). A portion of the F.S. cDNA was incubated for 2 hours at 15°C with NEB Buffer 2, DNA Polymerase I (New England Biolabs), and RNase H in order to make full-length double-strand cDNA. A fraction of F.S. cDNA was then normalized with Evrogen’s duplex-specific nuclease (DSN). F.S. cDNA-RNA duplexes in hybridization buffer were denatured at 98° for 3 minutes and then allowed to hybridize at 70°C for 5 hours. Preheated DSN at 1/4x concentration was then added and incubated for 20 minutes at 70°C. DSN stop solution was then added, and the sample was incubated for 5 minutes at 70°C. Normalized cDNA was then PCR amplified using an Encyclo PCR kit (Evrogen). PCR conditions were: 1 cycle × 95°C-1 min.; 17 cycles × 95°C-15 sec., 66°C-20 sec., 72°C-3 min.; 1 cycle × 66°C-15 sec., 72°C-3 min. Normalization efficiency was assayed via gel smear and qPCR of genes with known relative abundance (Additional file 5). qPCR analysis was performed using LinRegPCR [91, 92].
The non-normalized and normalized cDNA libraries were then pooled in a 1:2 ratio (Figure 1).
Five μg of pooled cDNA library was sent to the Roy J. Carver Biotechnology Center at the University of Illinois for sequencing. The cDNA library was sheared to 500-800 bp in length, 454 sequencing adaptors were ligated onto ends (Additional file 5), and the library was then converted to a single-stranded template library. Three titration runs (each of 1/16 lane) were performed to optimize sequencing conditions. A full plate was then sequenced on a Roche/454 GS FLX Sequencer using Titanium reagents.
Assembly of 454 sequence reads
Reads from the full plate and three titration runs were assembled using the Newbler Assembler v2.3 (Roche) using default parameters under the cDNA option. Prior to assembly, specified primers and adaptors were trimmed, namely the oligo-dT primer, the MINT PlugOligo adapter and PCR primer (Evrogen), and the 454 sequencing adaptors (Additional file 5).
Determination of fraction of captured transcripts
The coverage statistic developed by Susko and Roger (2004) estimates the proportion of genes from a cDNA library that is actually represented in the sequence data . Using this method, the unbiased estimate of coverage was calculated for our transcriptome with the equation = 1 – n 1 /n, where n 1 is the number of singletons in the assembly and n is the total number of reads [39, 40]. The new gene discovery rate was estimated using the term 1/(1 – ).
Annotation and analysis of BLAST hits
The set of 95,644 isotigs was input into the EST2uni annotation pipeline using default parameters, but with PCR marker integration, microarray printing, reciprocal BLAST for orthologues, Gene Ontology, and RFLP integration options turned off . Within EST2uni, the CAP3 assembly parameters were adjusted to –f 2 –g 100 –p 100 –d 110 to produce an assembly of all singletons, thereby preserving the isotigs produced by Newbler. Isotigs were annotated if they produced BLASTX matches against UniProtKB Release 2010_04 (23-Mar-2010) with E-values less than e-10. Parsing of BLAST data for Table 2 was done with custom Perl scripts (available upon request).
Completeness of annotated isotigs in Table 3 was performed using only the largest isotig from each isogroup as a representative for its putative gene locus. A Perl script utilizing BioPerl modules (available upon request) was used for completeness analysis. An isotig was considered complete on either end if it matched within ten amino acids of the corresponding end of the UniProt sequence .
The set of max isotigs was also used to identify Gene Ontology (GO) designations using the program Blast2GO [44, 45]. Results from BLAST searches against the UniProt database were imported from EST2uni, and GO annotation in Blast2GO was performed with default parameters (E-value threshold of e-6).
Identification of candidate regeneration genes
TBLASTX was used to search the P. leidyi isotig dataset for homologs of genes implicated in animal regeneration in the literature. A reciprocal TBLASTX search was then performed against UniProt or the nr database via NCBI to verify the putative identity of candidate regeneration genes in P. leidyi.
Validation of transcript assembly
Transcript validation assays were performed for two isogroups, isogroup08478 (a putative piwi-like homolog) and isogroup03233 (a putative frizzled homolog). Isotigs of each isogroup were aligned together using ClustalX v2.1 with manual editing by Seaview v4.0 [93, 94]. PCR was then performed to verify contigs and isotigs of each isogroup (Additional file 5). PCR amplicons of the expected size were either sequenced directly or cloned into the pGEM-T Easy vector (Promega) prior to sequencing. Sequencing was performed using an Applied Biosystems 3730 × l DNA Analyzer.
Transcript assembly was also verified for four previously known P. leidyi genes. BLAST searches against the 454 dataset were performed using the published gene sequences of Pl-en [GenBank: AF336055.1], Pl-otx1 [GenBank: AF336056.1], Pl-otx2 [GenBank: AF336057.1], and Pl-nos [GenBank: GQ369728.1]. GenBank sequences were aligned to transcriptome sequences using Sequencher v.4.7 (Gene Codes Corporation).
Analysis of gene expression by whole mount in situ hybridization
A ~1250 bp fragment of Pl-fzA (isogroup23343) and a ~1300 bp fragment of Pl-β-cat (isogroup01340) were amplified by PCR (Additional file 5). Synthesis of sense and antisense riboprobes and in situ hybridization were performed as previously described .
We thank Yinon Bentor for help with PERL and other programming issues, the University of Illinois’ Roy J. Carver Biotechnology Center for performing 454 sequencing and helping with subsequent analyses, Karen L. Carleton for assistance with DSN normalization and qPCR, Duygu Özpolat for assistance with in situ hybridizations, the Delwiche Lab for live cultures of Arthrospira platensis, and Thomas D. Kocher for providing bioinformatic resources. This work was funded by NSF grant IOS-0920502 to AEB.
- Réaumur RAF: Sur les diverses reproductions qui se font dans les Ecrevisse, les Omars, les Crabes, etc. et entr'autres sur celles de leurs Jambes et de leurs Ecailles. Mem Acad Royal Sci. 1712, : 223-245.
- Trembley A: Memoires pour servir a l'Histoire d'un Genre de Polypes d'Eau douce, a Bras en Forme de Cornes. 1744
- Bely AE, Nyberg KG: Evolution of animal regeneration: re-emergence of a field. Trends Ecol Evol. 2010, 25: 161-170. 10.1016/j.tree.2009.08.005.View ArticlePubMed
- Goss RJ: The evolution of regeneration: adaptive or inherent?. J Theor Biol. 1992, 159: 241-260. 10.1016/S0022-5193(05)80704-0.View ArticlePubMed
- Bely AE: Evolutionary loss of animal regeneration: pattern and process. Integr Comp Biol. 2010, 50 (4): 515-527. 10.1093/icb/icq118.View ArticlePubMed
- Sánchez Alvarado A, Tsonis PA: Bridging the regeneration gap: genetic insights from diverse animal models. Nat Rev Genet. 2006, 7: 873-884.View ArticlePubMed
- Brockes JP, Kumar A: Comparative aspects of animal regeneration. Annu Rev Cell Dev Biol. 2008, 24: 525-549. 10.1146/annurev.cellbio.24.110707.175336.View ArticlePubMed
- Tanaka EM, Reddien PW: The cellular basis for animal regeneration. Dev Cell. 2011, 21: 172-185. 10.1016/j.devcel.2011.06.016.PubMed CentralView ArticlePubMed
- Bely AE: Distribution of segment regeneration ability in the Annelida. Integr Comp Biol. 2006, 46: 508-518. 10.1093/icb/icj051.View ArticlePubMed
- Bely AE, Sikes JM: Latent regeneration abilities persist following recent evolutionary loss in asexual annelids. Proc Natl Acad Sci USA. 2010, 107: 1464-1469. 10.1073/pnas.0907931107.PubMed CentralView ArticlePubMed
- Scadding SR: Phylogenic distribution of limb regeneration capacity in adult Amphibia. J Exp Zool. 1977, 202: 57-68. 10.1002/jez.1402020108.View Article
- Scadding SR: Limb regeneration in adult amphibia. Can J Zool. 1981, 59: 34-46. 10.1139/z81-007.View Article
- Wagner GP, Misof BY: Evolutionary modification of regenerative capability in vertebrates: a comparative study on teleost pectoral fin regeneration. J Exp Zool. 1992, 261: 62-78. 10.1002/jez.1402610108.View ArticlePubMed
- Vollrath F: Leg regeneration in web spiders and its implications for orb weaver phylogeny. Bull Br Arachnol Soc. 1990, 8: 177-184.
- Brinkhurst RO, Jamieson BGM: Aquatic oligochaeta of the world. 1971, Oliver and Boyd, Edinburgh
- Zattara EE, Bely AE: Evolution of a novel developmental trajectory: fission is distinct from regeneration in the annelid Pristina leidyi. Evol Dev. 2011, 13 (1): 80-95. 10.1111/j.1525-142X.2010.00458.x.View ArticlePubMed
- Bely AE: Decoupling of fission and regenerative capabilities in an asexual oligochaete. Hydrobiologia. 1999, 406: 243-251.View Article
- 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.View ArticlePubMed
- Gurley KA, Rink JC, Sánchez Alvarado A: β-catenin defines head versus tail identity during planarian regeneration and homeostasis. Science. 2008, 319: 323-327. 10.1126/science.1150029.PubMed CentralView ArticlePubMed
- Peterson CP, Reddien PW: Smed-Bcatenin-1 is required for anteroposterior blastema polarity in planarian regeneration. Science. 2008, 319: 327-330. 10.1126/science.1149943.View Article
- Adell T, Saló E, Boutros M, Bartscherer K: Smed-Evi/Wntless is required for beta-catenin-dependent and -independent processes during planarian regeneration. Development. 2009, 136: 905-910. 10.1242/dev.033761.View ArticlePubMed
- Almuedo-Castillo M, Saló E, Adell T: Disheveled is essential for neural connectivity and planar cell polarity in planarians. Proc Natl Acad Sci USA. 2011, 108: 2813-2818. 10.1073/pnas.1012090108.PubMed CentralView ArticlePubMed
- Iglesias M, Gomez-Skarmeta JL, Saló E, Adell T: Silencing of Smed-beta-catenin1 generates radial-like hypercephalized planarians. Development. 2008, 135: 1215-1221. 10.1242/dev.020289.View ArticlePubMed
- Kawakami Y, Rodríguez Esteban C, Raya M, Kawakami H, Marti M, Dubova I, Izpisua Belmonte JC: Wnt/beta-catenin signaling regulates vertebrate limb regeneration. Genes Dev. 2006, 20: 3232-3237. 10.1101/gad.1475106.PubMed CentralView ArticlePubMed
- Lin G, Slack JM: Requirement for Wnt and FGF signaling in Xenopus tadpole tail regeneration. Dev Biol. 2008, 316: 323-335. 10.1016/j.ydbio.2008.01.032.View ArticlePubMed
- McClure KD, Schubiger G: Transdetermination: drosophila imaginal disc cells exhibit stem cell-like potency. Int J Biochem Cell Biol. 2007, 39: 1105-1118. 10.1016/j.biocel.2007.01.007.PubMed CentralView ArticlePubMed
- Peterson CP, Reddien PW: A wound-induced Wnt expression program controls planarian regeneration polarity. Proc Natl Acad Sci USA. 2009, 106: 17061-17066. 10.1073/pnas.0906823106.View Article
- Schubiger M, Sustar A, Schubiger G: Regeneration and transdetermination: the role of wingless and its regulation. Dev Biol. 2010, 347: 315-324. 10.1016/j.ydbio.2010.08.034.PubMed CentralView ArticlePubMed
- Stoick-Cooper CL, Moon RT, Weidinger G: Advances in signaling in vertebrate regeneration as a prelude to regenerative medicine. Genes Dev. 2007, 21 (11): 1292-1315. 10.1101/gad.1540507.View ArticlePubMed
- Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang H, Landherr L, Tomsho LP, Hu Y, Carlson JE, Ma H, Schuster SC, Soltis DE, Soltis PS, Altman N, dePamphilis CW: Comparison of next generation sequencing technologies for transcriptome characterization. BMC Genomics. 2009, 10: 347-10.1186/1471-2164-10-347.PubMed CentralView ArticlePubMed
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer MLI, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL, Lu H, Makhijani VB, McDade KE, McKenna MP, Myers EW, Nickerson E, Nobile JR, Plant R, Puc BP, Ronan MT, Roth GT, Sarkis GJ, Simons JF, Simpson JW, Srinivasan M, Tartaro KR, Tomasz A, Vogt KA, Volkmer GA, Wang SH, Wang Y, Weiner MP, Yu P, Begley RF, Rothberg JM: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.PubMed CentralPubMed
- Rothberg JM, Leamon JH: The development and impact of 454 sequencing. Nat Biotechnol. 2008, 26: 1117-1124. 10.1038/nbt1485.View ArticlePubMed
- Gregory TR, Hebert PDN: Genome size estimates for some oligochaete annelids. Can J Zool. 2002, 80 (8): 1485-1489. 10.1139/z02-145.View Article
- Gregory TR: Animal genome size database. 2012,http://www.genomesize.com,
- Bely AE, Wray GA: Evolution of regeneration and fission in annelids: insights from engrailed- and orthodenticle-class gene expression. Development. 2001, 128: 2781-2791.PubMed
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, Shagin DA: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res. 2004, 32 (3): e37-10.1093/nar/gnh031.PubMed CentralView ArticlePubMed
- Bogdanova EA, Shagin DA, Lukyanov SA: Normalization of full-length enriched cDNA. Mol Biosyst. 2008, 4: 205-212. 10.1039/b715110c.View ArticlePubMed
- NCBI Short Read Archive.http://www.ncbi.nlm.nih.gov/Traces/sra/,
- Susko E, Roger AJ: Estimating and comparing the rates of gene discovery and expressed sequence tag (EST) frequencies in EST surveys. Bioinformatics. 2004, 20: 2279-2287. 10.1093/bioinformatics/bth239.View ArticlePubMed
- Lee BY, Howe AE, Conte MA, D'Cotta H, Pepey E, Baroiller JF, di Palma F, Carleton KL, Kocher TD: An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags. BMC Genomics. 2010, 11: 278-10.1186/1471-2164-11-278.PubMed CentralView ArticlePubMed
- BouillaBase EST database.http://cichlid.umd.edu/est2uni//login.php,
- Bely lab resources.http://www.life.umd.edu/biology/bely/web/belylab/resources.html,
- Forment J, Gilabert F, Robles A, Conejero V, Nuez F, Blanca JM: EST2uni: an open, parallel tool for automated EST analysis and database creation, with a data mining web interface and microarray expression data integration. BMC Bioinformatics. 2008, 9: 5-10.1186/1471-2105-9-5.PubMed CentralView ArticlePubMed
- Conesa A, Gotz 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 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMed
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralView ArticlePubMed
- Campbell LJ, Crews CM: Wound epidermis formation and function in urodele amphibian limb regeneration. Cell Mol Life Sci. 2008, 65: 73-79. 10.1007/s00018-007-7433-z.View ArticlePubMed
- Bely AE, Sikes JM: Acoel and platyhelminth models for stem-cell research. J Biol. 2010, 9 (2): 14-10.1186/jbiol223.PubMed CentralView ArticlePubMed
- Cho S-J, Valles Y, Giani VC, Seaver EC, Weisblat DA: Evolutionary dynamics of the wnt gene family: a lophotrochozoan perspective. Mol Biol Evol. 2010, 27 (7): 1645-1658. 10.1093/molbev/msq052.PubMed CentralView ArticlePubMed
- Riddiford N, Olson PD: Wnt gene loss in flatworms. Dev Genes Evol. 2011, 221 (4): 187-197. 10.1007/s00427-011-0370-8.View ArticlePubMed
- Kato T, Miyazaki K, Shimizu-Nishikawa K, Koshiba K, Obara M, Mishima HK, Yoshizato K: Unique expression patterns of matrix metalloproteinases in regenerating newt limbs. Dev Dyn. 2003, 226 (2): 366-376. 10.1002/dvdy.10247.View ArticlePubMed
- Leontovich AA, Zhang JS, Shimokawa K, Nagase H, Sarras MP: A novel hydra matrix metalloproteinase (HMMP) functions in extracellular matrix degradation, morphogenesis and the maintenance of differentiated cells in the foot process. Development. 2000, 127 (4): 907-920.PubMed
- Altincicek B, Vilcinskas A: Comparative analysis of septic injury-inducible genes in phylogenetically distant model organisms of regeneration and stem cell research, the planarian Schmidtea mediterranea and the cnidarian Hydra vulgaris. Frontiers Zool. 2008, 5: 6-10.1186/1742-9994-5-6.View Article
- Reddien PW, Oviedo NJ, Jennings JR, Jenkin JC, Sánchez Alvarado A: SMEDWI-2 is a PIWI-like protein that regulates planarian stem cells. Science. 2005, 310: 1327-1330. 10.1126/science.1116110.View ArticlePubMed
- De Mulder K, Pfister D, Kuales G, Egger B, Salvenmoser W, Willems M, Steger J, Fauster K, Micura R, Borgonie G, Ladurner P: Stem cells are differentially regulated during development, regeneration, and homeostasis in flatworms. Dev Biol. 2009, 334: 198-212. 10.1016/j.ydbio.2009.07.019.View ArticlePubMed
- De Mulder K, Kuales G, Pfister D, Willems M, Egger B, Salvenmoser W, Thaler M, Gorny A-K, Hrouda M, Borgonie G, Ladurner P: Characterization of the stem cell system of the acoel Isodiametra pulchra. BMC Dev Biol. 2009, 9: 69-10.1186/1471-213X-9-69.PubMed CentralView ArticlePubMed
- Seipel K, Yanze N, Schmid V: The germ line and somatic stem cell gene Cniwi in the jellyfish Podocoryne carnea. Int J Dev Biol. 2004, 48 (1): 1-7. 10.1387/ijdb.15005568.View ArticlePubMed
- Handberg-Thorsager M, Saló E: The planarian nanos-like gene Smednos is expressed in germline and eye precursor cells during development and regeneration. Dev Genes Evol. 2007, 217 (5): 403-411. 10.1007/s00427-007-0146-3.View ArticlePubMed
- Wang Y, Zayas RM, Guo T, Newmark PA: nanos function is essential for development and regeneration of planarian germ cells. Proc Natl Acad Sci USA. 2007, 104: 5901-5906. 10.1073/pnas.0609708104.PubMed CentralView ArticlePubMed
- Hayashi T, Mizuno N, Takada R, Takada S, Kondoh H: Determinative role of Wnt signals in dorsal iris-derived lens regeneration in newt eye. Mech Dev. 2006, 123 (11): 793-800. 10.1016/j.mod.2006.08.009.View ArticlePubMed
- Adell T, Marsal M, Saló E: Planarian GSK3s are involved in neural regeneration. Dev Genes Evol. 2008, 218 (2): 89-103. 10.1007/s00427-007-0199-3.View ArticlePubMed
- D'Jamoos CA, McMahon G, Tsonis PA: Fibroblast growth factor receptors regulate the ability for hindlimb regeneration in Xenopus laevis. Wound Repair Regen. 1998, 6 (4): 388-397.View ArticlePubMed
- Poss FD, Shen JX, Nechiporuk A, McMahon G, Thisse B, Thisse C, Keating MT: Roles for Fgf signaling during zebrafish fin regeneration. Dev Biol. 2000, 222 (2): 347-358. 10.1006/dbio.2000.9722.View ArticlePubMed
- Martin P, Parkhurst SM: Parallels between tissue repair and embryo morphogenesis. Development. 2004, 131 (13): 3021-3034. 10.1242/dev.01253.View ArticlePubMed
- Tasaki J, Shibata N, Sakurai T, Agata K, Umesono Y: Role of c-Jun N-terminal kinase activation in blastema formation during planarian regeneration. Dev Growth Differ. 2011, 53 (3): 389-400. 10.1111/j.1440-169X.2011.01254.x.View ArticlePubMed
- Reddien PW, Bermange AL, Kicza AM, Sánchez Alvarado A: BMP signaling regulates the dorsal planarian midline and is needed for asymmetric regeneration. Development. 2007, 134 (22): 4043-4051. 10.1242/dev.007138.View ArticlePubMed
- Molina MD, Saló E, Cebrià F: The BMP pathway is essential for re-specification and maintenance of the dorsoventral axis in regenerating and intact planarians. Dev Biol. 2007, 311 (1): 79-94. 10.1016/j.ydbio.2007.08.019.View ArticlePubMed
- Pearl EJ, Barker D, Day RC, Beck CW: Identification of genes associated with regenerative success of Xenopus laevis hindlimbs. BMC Dev Biol. 2008, 8: 66-10.1186/1471-213X-8-66.PubMed CentralView ArticlePubMed
- Molina MD, Neto A, Maeso I, Luis Gomez-Skarmeta J, Saló E, Cebrià F: Noggin and noggin-like genes control dorsoventral axis regeneration in planarians. Curr Biol. 2011, 21 (4): 300-305. 10.1016/j.cub.2011.01.016.View ArticlePubMed
- Rink JC, Gurley KA, Elliot SA, Sánchez Alvarado A: Planarian Hh signaling regulates regeneration polarity and links Hh pathway evolution to cilia. Science. 2009, 326: 1406-1410. 10.1126/science.1178712.PubMed CentralView ArticlePubMed
- Yazawa S, Umesono Y, Hayashi T, Tarui H, Agata K: Planarian Hedgehog/Patched establishes anterior-posterior polarity by regulating Wnt signaling. Proc Natl Acad Sci USA. 2009, 106 (52): 22329-22334. 10.1073/pnas.0907464106.PubMed CentralView ArticlePubMed
- Schnapp E, Kragl M, Rubin L, Tanaka EM: Hedgehog signaling controls dorsoventral patterning, blastema cell proliferation and cartilage induction during axolotl tail regeneration. Development. 2005, 132 (14): 3243-3253. 10.1242/dev.01906.View ArticlePubMed
- Tsonis PA, Vergara MN, Spence JR, Madhavan M, Kramer EL, Call MK, Santiago WG, Vallance JE, Robbins DJ, Del Rio-Tsonis K: A novel role of the hedgehog pathway in lens regeneration. Dev Biol. 2004, 267 (2): 450-461. 10.1016/j.ydbio.2003.12.005.View ArticlePubMed
- Mannini L, Deri P, Gremigni V, Rossi L, Salvetti A, Batistoni R: Two msh/msx-related genes, Djmsh1 and Djmsh2, contribute to the early blastema growth during planarian head regeneration. Int J Dev Biol. 2008, 52 (7): 943-952. 10.1387/ijdb.072476lm.View ArticlePubMed
- Carlson MRJ, Bryant SV, Gardiner DM: Expression of Msx-2 during development, regeneration, and wound healing in axolotl limbs. J Exp Zool. 1998, 282 (6): 715-723. 10.1002/(SICI)1097-010X(19981215)282:6<715::AID-JEZ7>3.0.CO;2-F.View ArticlePubMed
- Cho S-J, Lee MS, Tak ES, Lee E, Koh KS, Ahn CH, Park SC: Gene expression profile in the anterior regeneration of the earthworm using expressed sequence tags. Biosci Biotechnol Biochem. 2009, 73 (1): 29-34. 10.1271/bbb.80391.View ArticlePubMed
- Mullen LM, Bryant SV, Torok MA, Blumberg B, Gardiner DM: Nerve dependency of regeneration: the role of Distal-less and FGF signaling in amphibian limb regeneration. Development. 1996, 122 (11): 3487-3497.PubMed
- Juliano CE, Swartz SZ, Wessel GM: A conserved germline multipotency program. Development. 2010, 137 (24): 4113-4126. 10.1242/dev.047969.PubMed CentralView ArticlePubMed
- Neira-Oviedo M, Tsyganov-Bodounov A, Lycett GJ, Kokoza V, Raikhel AS, Krzywinski J: The RNA-Seq approach to studying the expression of mosquito mitochondrial genes. Insect Mol Biol. 2011, 20 (2): 141-152. 10.1111/j.1365-2583.2010.01053.x.View ArticlePubMed
- Wilhelm BT, Landry J-R: RNA-Seq-quantitative measurement of expression through massively parallel RNA-sequencing. Methods. 2009, 48 (3): 249-257. 10.1016/j.ymeth.2009.03.016.View ArticlePubMed
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMed
- Henry JQ, Perry KJ, Fukui L, Alvi N: Differential localization of mRNAs during early development in the mollusc, Crepidula fornicata. Integr Comp Biol. 2010, 50 (5): 720-733. 10.1093/icb/icq088.View ArticlePubMed
- Lambert JD, Chan XY, Spiecker B, Sweet HC: Characterizing the embryonic transcriptome of the snail Ilyanassa. Integr Comp Biol. 2010, 50 (5): 768-777. 10.1093/icb/icq121.View ArticlePubMed
- Clark MS, Thorne MAS, Vieira FA, Cardoso JCR, Power DM, Peck LS: Insights into shell deposition in the Antarctic bivalve Laternula elliptica: gene discovery in the mantle transcriptome using 454 pyrosequencing. BMC Genomics. 2010, 11: 362-10.1186/1471-2164-11-362.PubMed CentralView ArticlePubMed
- Hou R, Bao Z, Wang S, Su H, Li Y, Du H, Hu J, Wang S, Hu X: Transcriptome sequencing and de novo analysis for Yesso scallop (Patinopecten yessoensis) using 454 GS FLX. PLoS One. 2011, 6 (6): e21560-10.1371/journal.pone.0021560.PubMed CentralView ArticlePubMed
- Gong P, Pirooznia M, Guan X, Perkins EJ: Design, validation and annotation of transcriptome-wide oligonucleotide probes for the oligochaete annelid Eisenia fetida. PLoS One. 2010, 5 (12): e14266-10.1371/journal.pone.0014266.PubMed CentralView ArticlePubMed
- Heyland A, Vue Z, Voolstra CR, Medina M, Moroz LL: Developmental transcriptome of Aplysia californica. J Exp Zool B. 2011, 316B (2): 113-134. 10.1002/jez.b.21383.View Article
- Milan M, Coppe A, Reinhardt R, Cancela LM, Leite RB, Saavedra C, Ciofi C, Chelazzi G, Patarnello T, Bortoluzzi S, Bargelloni L: Transcriptome sequencing and microarray development for the Manila clam, Ruditapes philippinarum: genomic tools for environmental monitoring. BMC Genomics. 2011, 12: 234-10.1186/1471-2164-12-234.PubMed CentralView ArticlePubMed
- Abril JF, Cebrià F, Rodríguez Esteban G, Horn T, Fraguas S, Calvo B, Bartscherer K, Saló E: Smed454 dataset: unraveling the transcriptome of Schmidtea mediterranea. BMC Genomics. 2010, 11: 731-10.1186/1471-2164-11-731.PubMed CentralView ArticlePubMed
- Adamidi C, Wang Y, Gruen D, Mastrobuoni G, You X, Tolle D, Dodt M, Mackowiak SD, Gogol-Doering A, Oenal P, Rybak A, Ross E, Sánchez Alvarado A, Kempa S, Dieterich C, Rajewsky N, Chen W: De novo assembly and validation of planaria transcriptome by massive parallel sequencing and shotgun proteomics. Genome Res. 2011, 21: 1193-1200. 10.1101/gr.113779.110.PubMed CentralView ArticlePubMed
- Hardie DC, Gregory TR, Hebert PDN: From pixels to picograms: a beginners' guide to genome quantification by Feulgen image analysis densitometry. J Histochem Cytochem. 2002, 50 (6): 735-749. 10.1177/002215540205000601.View ArticlePubMed
- Ramakers C, Ruijter JM, Deprez RHL, Moorman AFM: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003, 339 (1): 62-66. 10.1016/S0304-3940(02)01423-4.View ArticlePubMed
- Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, Moorman AFM: Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009, 37 (6): e45-10.1093/nar/gkp045.PubMed CentralView ArticlePubMed
- Gouy M, Guindon S, Gascuel O: SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010, 27 (2): 221-224. 10.1093/molbev/msp259.View ArticlePubMed
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.