Functional analysis of Girardia tigrina transcriptome seeds pipeline for anthelmintic target discovery
© Wheeler et al.; licensee BioMed Central. 2015
Received: 13 October 2014
Accepted: 23 December 2014
Published: 20 January 2015
Neglected diseases caused by helminth infections impose a massive hindrance to progress in the developing world. While basic research on parasitic flatworms (platyhelminths) continues to expand, researchers have yet to broadly adopt a free-living model to complement the study of these important parasites.
We report the high-coverage sequencing (RNA-Seq) and assembly of the transcriptome of the planarian Girardia tigrina across a set of dynamic conditions. The assembly was annotated and extensive orthology analysis was used to seed a pipeline for the rational prioritization and validation of putative anthelmintic targets. A small number of targets conserved between parasitic and free-living flatworms were comparatively interrogated.
240 million paired-end reads were assembled de novo to produce a strictly filtered predicted proteome consisting of over 22,000 proteins. Gene Ontology annotations were extended to 16,467 proteins. 2,693 sequences were identified in orthology groups spanning flukes, tapeworms and planaria, with 441 highlighted as belonging to druggable protein families. Chemical inhibitors were used on three targets in pharmacological screens using both planaria and schistosomula, revealing distinct motility phenotypes that were shown to correlate with planarian RNAi phenotypes.
This work provides the first comprehensive and annotated sequence resource for the model planarian G. tigrina, alongside a prioritized list of candidate drug targets conserved among parasitic and free-living flatworms. As proof of principle, we show that a simple RNAi and pharmacology pipeline in the more convenient planarian model system can inform parasite biology and serve as an efficient screening tool for the identification of lucrative anthelmintic targets.
KeywordsSchistosomiasis Girardia tigrina Planaria RNA-Seq Drug discovery
The Platyhelminthes (flatworms) comprise a diverse phylum of medically and economically important species. Trematodes (flukes) and cestodes (tapeworms) are the etiological agents of several Neglected Tropical Diseases (NTDs) that disproportionately devastate the health and economic prospects of the poor across much of the developing world. Schistosomes infect over 220 million in sub-Saharan Africa alone, and 600–800 million live at risk of infection worldwide . Echinococcosis and cysticercosis, while less prevalent than schistosomiasis, are zoonotic parasitic diseases of great public health importance. These neglected diseases inflict significant morbidity and mortality, accounting for upwards of 280,000 deaths  and an annual loss of between 3.5 - 70 million disability-adjusted life years (DALYs) .
The prioritization of flatworm-associated NTDs by the World Health Organization  underscores the urgency of efforts to control infection and to develop new anthelmintic treatments. The threat of drug resistance  further calls attention to the need for novel pipelines for drug target validation and drug discovery . Against this backdrop, free-living flatworms represent a new and potentially powerful screening model for parasite drug discovery efforts . Members of the free-living genus Planaria are widely interrogated in the realm of stem cell biology due to their remarkable regenerative abilities . In comparison to their parasitic counterparts, planarians are much more amenable to modern genetic protocols and their culture and maintenance within the laboratory is relatively inexpensive and simple. Many behavioral, biochemical, and morphological phenotypes have also been described for planaria [9-11], enabling straightforward inferences of function from combined reverse genetic, pharmacological, and phenotypic analyses.
In the case of schistosomes, it is necessary to maintain active populations of freshwater snails as intermediate hosts, manage periodic shedding of the infective cercariae, induce transformation to schistosomula or allow for penetrance into a definitive host (usually mice). The process is difficult, time consuming, moderately dangerous, and, for many labs, cost prohibitive. These concerns underpin efforts to extend the utility of planarian biology to the study of nearly-related parasites , mirroring the important role that Caenorhabditis elegans has played in furthering our understanding of the biology of parasitic nematodes .
A number of planarian species see use in the laboratory, with varying modes of reproduction, regenerative potential, and genome ploidy. Schmidtea mediterranea is among the most widely studied species. Clonal lines of S. mediterranea have been propagated to mitigate genetic heterogeneity, and both genomic  and transcriptomic  data have been published for this stable diploid [16-18]. Other notable planarian species include Girardia (formerly Dugesia) tigrina and Dugesia japonica . Genome assembly and analyses are partly complicated in these species due to their mixoploid genomes and the presence of large numbers of transposable elements . No significant sequence resources yet exist for G. tigrina, despite the convenient commercial availability of this species and its broad adoption in the fields of regeneration , pharmacology  and learning and memory .
The emergence of a comprehensive sequence resource for G. tigrina will open avenues for more precise biological manipulation of these planaria. RNA-Seq provides a powerful platform for producing a high-coverage transcriptome, without the complications of whole genome assembly. The selection of G. tigrina for this undertaking presents a reasonable trade-off, whereby some level of genetic heterogeneity is accepted for the greater ease of procuring, maintaining and scaling colonies, in comparison to clonally-derived laboratory strains. Genetic variation within this sexual strain is minor with respect to the accepted genetic distance between planaria and the flatworm parasites for which they are to serve as models. Although computationally intensive in the absence of a reference genome, a high-quality de novo transcriptome assembly would allow for closer examination of our overarching hypothesis that G. tigrina could provide a shortcut to identifying potential drug targets in the phylum.
Planarian culture and RNA isolation
G. tigrina (Ward’s Natural Science, Rochester, NY) colonies were maintained in aerated spring water. Planarians were starved for one week prior to RNA isolation. Five animals were randomly selected per experimental condition on day 7. Each group was washed repeatedly with spring water and tissue grinding was carried out using mortar and pestle in the presence of liquid nitrogen. A hybrid TRIzol (Invitrogen)/ RNeasy (Qiagen) protocol was used to isolate total RNA from ground tissue, whereby supernatants from the chloroform phase separation were combined with an equal volume of 100% ethanol and loaded into RNeasy columns for purification. Total RNA quality and concentration was assessed with an Agilent Bioanalyzer 2100. RNA integrity number (RIN) proved to be a poor benchmark of RNA quality, as the G. tigrina 28S rRNA subunit is evidently converted into fragments that co-migrate with 18S rRNA to produce a triple-peak, giving the misleading appearance of RNA degradation. All samples yielded at least 1 ug/ul of RNA when eluted in 40 ul of H20, with an OD A260/A280 of ∼ 2.1 and OD A260/A230 of ∼ 2.2.
Library preparation and Illumina paired-end RNA-seq
Illumina HiSeq 2000 paired-end (2x100 bp) library preparation and sequencing was carried out at the McGill University and Genome Quebec Innovation Center. The four RNA samples were multiplexed across two sequencing lanes with an average fragment size of ∼ 350 bp, corresponding to an average insert size of ∼ 225 bp. The sequencer run yielded ∼ 30 million paired-end reads per sample (241 million total paired-end reads) with an average Phred quality score of 37.
De novo transcriptome assembly
Adapter sequences were trimmed and reads were passed through a sliding window quality filter (window size = 4, minimum average quality score = 25) using Trimmomatic 0.22 . Paired-end reads and singletons ≥ 50 bp in length were retained. Overlapping paired-end reads were identified and merged using FLASH  with an expected insert size of 220 bp. Quality control and read collapsing led to a total filtered pool of 165 million paired-end reads and 55 million singletons. Surviving reads were combined and fed into the Trinity  pipeline for de novo assembly, performed on the GLUMEQ Guillimin supercomputer maintained by McGill University. Assembly optimization and runs were carried out on a 1 TB ScaleMP node that allows for a virtualized shared large memory environment required by the OpenMP standard. Final assembly was carried out with a minimum k-mer coverage of 2 and the default k-mer size of 25. Complex graphs that proved unresolvable within a 6 hour window were manually excised to allow the assembly to proceed. Separately, an available Python script was used to feed the same read pool into the Velvet  pipeline and to generate multiple k-mer assemblies (k = 21, 25, 29 and 33) for merging with Oases (k = 25) . The minimum contig or transcript length for both assembly pipelines was set to 200 nt. The statistical software R  was used to generate and evaluate assembly statistics. Further bioinformatic analysis was restricted to the Trinity transcriptome.
Transcriptome filtering and annotation
Filtered paired-end reads were mapped to the Trinity transcriptome with Bowtie . Abundance estimation with RSEM  was used to select for transcripts that accounted for at least 1% of the per-component (IsoPct) expression and that met a TPM cutoff of 1. Open reading frames (ORFs) with coding potential were predicted from the final transcriptome using log-likelihood scores based on codon usage with Transdecoder (http://transdecoder.sourceforge.net/). The resulting predicted proteome was further filtered with CD-HIT-EST  at a threshold of 0.95. BLAST2GO  was used to functionally annotate the G. tigrina predicted proteome and to assign GO terms to predicted proteins.
Differential expression analysis
Condition-specific abundance estimation was carried out with Bowtie and RSEM using the final filtered transcriptome. Existing Trinity scripts and the R/Bioconductor packages DESeq and edgeR of the statistical programming language R were used to identify differentially expressed transcripts. A threshold e-value of 10–3 and minimum four-fold expression changes were used to select and cluster transcripts as either up or down-regulated. Transcript sets were then mapped to previous annotations, where available.
Orthology analysis and drug target prioritization
Proteinortho  was used to detect ortholog groups between and among the G. tigrina, S. mansoni and E. multilocularis predicted proteomes, with default parameters and a coverage threshold of 0.5. The predicted parasite proteins were downloaded from GeneDB . The visualization tool Circos  was utilized to organize and display the orthologous relationships among these species in the context of Trematoda and Cestoda synteny. Provided that a transcriptome does not imply any spatial or chromosomal arrangement, the G. tigrina transcripts were arbitrarily arranged on a pseudo-chromosome - designated as an ideogram - to enable visualization. GFF files obtained from GeneDB containing gene coordinate data for the two parasitic species were parsed for necessary sequence features using a set of in-house Python scripts and used to draw orthologous “links” between ideograms. Heatmap data was created by running command-line BLAST . Similarity calculations were carried out with G. tigrina orthologs as queries against their corresponding parasitic orthologs, as well as with parasitic orthologs as queries against the RefSeq human proteome . Sequence homology was used to select orthologs that displayed high sequence similarity among the three examined flatworm species, as well those sufficiently diverged from their nearest-identifiable human homolog. These prioritized ortholog groups were mined with GO IDs, along with manually-selected ID’s, to extract and highlight annotated sequences belonging to notoriously druggable protein families .
Snails (Biomphalaria glabrata) infected with S. mansoni were provided by the Biomedical Research Institute (BRI) (Rockville, MD). Cercariae were shed from snails by light exposure and subsequently mechanically transformed to schistosomula in vitro per existing protocols . Somules were cultured in modified Basch medium (containing 10% Fetal Bovine Serum) at 37°C in 5% CO2 atmosphere.
Three small molecule inhibitors (CK-666, ARPC2; 3-NPA, SDH1; rotenone, NDUFV2) were purchased from Sigma-Aldrich (St. Louis, MO). Newly transformed schistosomula were incubated with various concentrations of each inhibitor for 30 minutes and recorded at 10× magnification for 5 minutes. Approximately 200 worms were assayed per experimental well, and 12 were randomly selected for video capture and analysis. The wrMTrck plugin of ImageJ  was used to track schistosomules and quantify motility in terms of contractile rate (body bends per second; BPS).
Planarian pharmacology and motility analysis
G. tigrina individuals were placed in 35 mm dishes filled with 4 mL of media supplemented with inhibitor at varying concentrations or an equal amount of solvent control, and the dishes were placed on a light box in a dark room. After 30 minutes of incubation, worms were recorded for 5 minutes by EthoVision , and motility was quantified by dividing the parameter DistanceTraveled (mm) by TimeInZone (s). Down-sampling was set to 5 to ensure that small bending and twisting motions were not factored. Tracking profiles were visually diagnosed for errors and manually edited where required. Errors were most often attributable to light reflections off of the surface of liquid media or imperfect arena definitions.
Planarian RNAi and regenerative assays
Total RNA was extracted from homogenized G. tigrina and converted to cDNA with Ambion’s RetroScript RT kit. 600 bp sequences were PCR amplified using primers designed with Primer3 . T7 promoter sites were added using a two-step PCR protocol, and dsRNA was transcribed with the Ambion MegaScript RNAi kit. dsRNA was added directly to homogenized liver paste according to prescribed methods  (10 umol/worm for one hour). Feedings were performed on Days 1, 3, and 5, and worms were bisected immediately above the pharynx on Day 6. A minimum of 10 planaria (20 worm halves) were used for each experiment. Two worms from each experimental group were set aside for semi-quantitative RT-PCR performed with Ambion’s QuantumRNA18S Internal Standards kit. In each experiment, 10 planaria (20 worm-halves) were observed for defects in regeneration over the full regeneration period (∼2-3 weeks).
Results and discussion
De novo transcriptome assembly
Protein prediction and transcriptome annotation
Predicted proteomes were created for each assembly using Transdecoder (http://transdecoder.sourceforge.net/) to evaluate the coding potential of open reading frames based on codon usage. After subtraction of redundant proteins, the Assembly-T proteome had a marginally larger unique protein count and was therefore used for all subsequent analysis. The 22,363 predicted proteins from this dataset were used in blastx  queries against the NCBI nr database. 16,467 sequences had at least one significant hit (e-value < 0.001) and the top 20 hits were retained for each sequence. We then applied Gene Ontology (GO) annotations to the de novo proteome using the Blast2GO  pipeline. The mapping of sequence-specific blast results to GO identifiers resulted in functional annotation of all 16,467 proteins with significant blast hits, accounting for 73% of the predicted proteome. This was complemented by InterProScan  domain mapping, resulting in the identification of at least one domain for 18,051 protein sequences.
Annotated predicted proteins were categorized according to their involvements in various biological processes at different hierarchy levels. Figure 2 depicts this categorization from more general level 2 categories through more specific level 6 categories, with many proteins binned into multiple categories. Separately, to further gauge coverage of core pathways, Kegg  pathway mapping was carried out. As a representative example, over 95% of core reference pathway enzymes were identified (data not shown). The very small numbers of unmapped enzymes in these and related metabolic pathways could result partly from fundamental biological differences, as opposed to gaps in our dataset.
Identification of differentially expressed transcripts
While the inclusion of different treatment conditions was primarily aimed at increasing transcript capture, it also presents an opportunity to identify transcripts that are significantly upregulated or downregulated with respect to these conditions. Previous investigators have carried out experiments to identify regeneration-associated genes in S. mediterranea  by performing RNAi screens to perturb normal regeneration. While many transcripts show greater than four-fold differences in expression between control and cut worms (e.g. RNA-dependent RNA polymerase, RNA helicases, reverse transcriptase), these data do not lend themselves to facile implications of molecular mechanisms of regeneration and are provided for further examination and investigation (Additional file 1).
Identification of planarian-parasite orthologs
Drug target prioritization
Ortholog groups with high sequence conservation through the phylum represent potential broad-spectrum therapeutic targets, and these can plausibly be interrogated using G. tigrina as a more tractable free-living model. Ideally, we want to prioritize planarian protein targets that share very high sequence similarity with both fluke and tapeworm homologs, and that exhibit lower levels of sequence similarity with any identifiable host (human) proteins. Figure 4C applies this selection logic towards exploitation of the available sequence data. Here, links between S. mansoni and E. multilocularis reveal synteny for those orthologous gene pairs that share a highly similar G. tigrina ortholog and which represent useful anthelmintic targets. The final links are restricted to a set of 441 putative drug targets, filtered from the initial set of 2,693 ortholog groups (Additional file 2).
To estimate the druggability of protein targets, we utilized GO annotations that are most often associated with established drug target classes (http://www.ebi.ac.uk/thorntonsrv/databases/cgi-bin/drugport/GetPage.pl?template=goanal.html). Highly-associated GO terms, manually supplemented, were collected into a list and used to extract specific sequences from the annotated ortholog dataset. Within this set, we looked to identify a handful of targets as proof of concept for our model paradigm. Specifically, three targets were chosen that showed high sequence similarity between free-living and parasitic flatworms and that we could potentially pharmacologically manipulate or inhibit with commercially available chemicals. These targets were actin-related protein complex 2/3 subunit 2 (ARPC2), succinate dehydrogenase (SDH1), and NADH dehydrogenase (ubiquinone) flavoprotein 2 (NDUFV2).
Comparative chemical screen of targets
To carry out a comparative first-pass phenotypic screen, chemical inhibitors for each target were used to treat newly transformed S. mansoni schistosomula and planaria across a range of concentrations. It should be noted that these chemical inhibitors have been shown to act on mammalian proteins, and the specificity of each interaction is therefore unknown. However, there is significant sequence conservation between protein domains in these proof of concept targets and their human homologs (Additional file 3), suggesting a high likelihood of a conserved mode of action. For example, G. tigrina ARPC2 shares 58% sequence identity with its human counterpart. While it can be hypothesized that CK-666, which locks this complex in an inactive state , performs this action in flatworms as well, it cannot be necessarily inferred.
RNAi in planaria is predictive of phenotypes in parasites
To assay regeneration, worms were bisected above the pharynx and each half was maintained in a separate well. Cephalic and caudal regeneration was observed over the course of 2–3 weeks. No developmental phenotype was observed for Dtig-NDUFV2(RNAi) worms, however, Dtig-ARPC2(RNAi) worms showed aberrant regeneration in comparison to control worms, consisting of a range of specific outcomes that included stalled or slowed regeneration, blastema malformation, an inability to seal the wound, and eventual death (Figures 6B and 6C). In these assays, readily identifiable phenotypes were observed in planaria that were predictive of visible phenotypes in the comparative chemical screen.
This proof of principle was limited to a handful of putative targets with commercially available inhibitors, but the pipeline can readily be scaled to larger numbers of targets. In the proposed scheme, highly conserved planarian-parasite orthologs are first interrogated with RNAi in planaria. The detection of phenotypes relating to planarian motility, morphology, or regeneration can serve as an efficient and high- throughput filter for target druggability in flatworm parasites. Parasite orthologs can then be investigated using available techniques such as RNAi and heterologous expression for drug target validation and functional characterization. This approach sidesteps the often prohibitive costs and technical challenges of carrying out large high-throughput screens in transient parasite life stages.
This work further promotes the adoption of planaria, and in particular Girardia tigrina, as a model screening organism for candidate drug targets in parasites. We provide a high-coverage annotated de novo transcriptome as a substrate for such efforts. The Assemblies, CDS predictions and annotations are available for download (http://goo.gl/BZU87d). The identification of ortholog groups that extend to planaria, blood flukes, and tapeworms, allows for the rational prioritization of likely broad-spectrum drug targets that can be readily screened in G. tigrina. We outline a pathway for the high-throughput evaluation of putative drug targets in planaria as a prelude to validation and more extensive characterization in parasitic flatworms. We further show how such screens can be predictive of biological phenotypes in parasites.
This study builds on other recent studies that have shown the utility of the planarian system in understanding parasite biology. For example, the antischistosomal praziquantel has been shown to lead to changes in planarian regenerative polarity through the action of voltage-operated calcium channel (VOCC) β subunits . In this conserved signaling pathway, regenerative polarity in planaria acts as a phenotypic correlate of drug efficacy and worm paralysis in schistosomes. This was followed by a more comprehensive investigation of the phenotypic correlates of manipulating signal transduction pathways in the planarian D. japonica and the parasite S. mansoni, as a predictive tool for the discovery of antischistosomal agents . The striking identification of adult stem cells in S. mansoni that resemble planarian neoblasts further strengthens the notion of fundamental biological conservation between free-living and parasitic flatworms . This annotated G. tigrina sequence resource, along with the orthology-based prioritization of putative drug targets, can help catalyze low cost and scalable in vivo pipelines for anthelmintic drug discovery.
This research was supported in part by a grant from the Fonds de Recherche du Quebec - Nature et Technologies (FRQNT) to PR (grant PR-147457). NW was supported by fellowship funds from the Office of Biotechnology, Iowa State University. B. glabrata snails exposed to Schistosoma mansoni, strain NMRI, were provided by the Schistosome Research Reagent Resource Center for distribution by BEI Resources, NIAID, NIH: NR-21962. We gratefully acknowledge Compute Canada/CLUMEQ for access to valuable computational resources.
- WHO. Schistosomiasis: population requiring preventive chemotherapy and number of people treated in 2010. Wkly Epidemiol Rec. 2012;87(4):37–44.Google Scholar
- van der Werf MJ, de Vlas SJ, Brooker S, Looman CWN, Nagelkerke NJD, Habbema JDF, et al. Quantification of clinical morbidity associated with schistosome infection in sub-Saharan Africa. Acta Trop. 2003;86(2–3):125–39. doi:10.1016/S0001-706X(03)00029-9.PubMedView ArticleGoogle Scholar
- Murray CJL, Vos T, Lozano R, Naghavi M, Flaxman AD, Michaud C, et al. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380(9859):2197–223. doi:10.1016/S0140-6736(12)61689-4.PubMedView ArticleGoogle Scholar
- Torreele E, Usdin M, Chirac P, Iyer BJK, Adams ER, Klatser PR: Priority medicines for Europe and the world: a public health approach to innovation. 2004Google Scholar
- Wang W, Wang L, Liang Y-S. Susceptibility or resistance of praziquantel in human schistosomiasis: a review. Parasitol Res. 2012;111(5):1871–7. doi:10.1007/s00436-012-3151-z.PubMedView ArticleGoogle Scholar
- Gilbert IH. Drug discovery for neglected diseases: molecular target-based and phenotypic approaches. J Med Chem. 2013;56(20):7719–26. doi:10.1021/jm400362b.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen B, Wen J-F. The adaptive evolution divergence of triosephosphate isomerases between parasitic and free-living flatworms and the discovery of a potential universal target against flatworm parasites. Parasitol Res. 2011;109(2):283–9. doi:10.1007/s00436-010-2249-4.PubMedView ArticleGoogle Scholar
- Scimone M, Kravarik K, Lapan S, Reddien P. Neoblast Specialization in Regeneration of the Planarian Schmidtea mediterranea. Stem Cell Rep. 2014;3(2):339–52. doi:10.1016/j.stemcr.2014.06.001.View ArticleGoogle Scholar
- Zhang J, Yuan Z, Zheng M, Sun Y, Wang Y, Yang S. Effects of N,N-dimethylformamide on behaviour and regeneration of planarian Dugesia japonica. Toxicol Ind Health. 2013;29(8):753–60. doi:10.1177/0748233712443148.PubMedView ArticleGoogle Scholar
- Evans DJ, Owlarn S, Tejada Romero B, Chen C, Aboobaker A. Combining classical and molecular approaches elaborates on the complexity of mechanisms underpinning anterior regeneration. PLoS One. 2011;6(11):27927. doi:10.1371/journal.pone.0027927.View ArticleGoogle Scholar
- Reddien PW, Bermange AL, Murfitt KJ, Jennings JR, Sanchez Alvarado A. Identification of genes needed for regeneration, stem cell function, and tissue homeostasis by systematic gene perturbation in planaria. Dev Cell. 2005;8(5):635–49. doi:10.1016/j.devcel.2005.02.014.PubMed CentralPubMedView ArticleGoogle Scholar
- Collins JJ, Hou X, Romanova EV, Lambrus BG, Miller CM, Saberi A, et al. Genome-wide analyses reveal a role for peptide hormones in planarian germline development. PLoS Biol. 2010;8(10):1000509. doi:10.1371/journal.pbio.1000509.View ArticleGoogle Scholar
- Gilleard JS. The use of Caenorhabditis elegans in parasitic nematode research. Parasitology. 2004;128(Suppl):49–70. doi:10.1017/S003118200400647X.Google Scholar
- Robb SMC, Ross E, Sa A. SmedGD : the Schmidtea mediterranea genome database. Nucleic Acids Res. 2008;36:599–606. doi:10.1093/nar/Gkm684.View ArticleGoogle Scholar
- Abril JF, Cebri’a F, Rodrıguez-Esteban G, Horn T, Fraguas S, Calvo B, et al. Smed454 dataset: unravelling the transcriptome of Schmidtea mediterranea. BMC Genomics. 2010;11(1):731. doi:10.1186/1471-2164-11-731.PubMed CentralPubMedView ArticleGoogle Scholar
- Galloni M. Global irradiation effects, stem cell genes and rare transcripts in the planarian transcriptome. Int J Dev Biol. 2012;56(1–3):103–16. doi:10.1387/ijdb.113455mg.PubMedView ArticleGoogle Scholar
- Resch AM, Palakodeti D, Lu Y-C, Horowitz M, Graveley BR. Transcriptome analysis reveals strain-specific and conserved stemness genes in Schmidtea mediterranea. PLoS One. 2012;7(4):34447. doi:10.1371/journal.pone.0034447.View ArticleGoogle Scholar
- Solana J, Kao D, Mihaylova Y, Jaber-Hijazi F, Malla S, Wilson R, et al. Defining the molecular profile of planarian pluripotent stem cells using a combinatorial RNAseq, RNA interference and irradiation approach. Genome Biol. 2012;13(3):19. doi:10.1186/gb-2012-13-3-r19.View ArticleGoogle Scholar
- Nishimura O, Hirao Y, Tarui H, Agata K. Comparative transcriptome analysis between planarian Dugesia japonica and other platyhelminth species. BMC Genomics. 2012;13(1):289. doi:10.1186/1471-2164-13-289.PubMed CentralPubMedView ArticleGoogle Scholar
- Garcia-Fernandez J, Ram J, Mar-any G, Mun AM. High Copy Number of Highly Similar mariner-like Transposons in Planarian ( Platyhelminthe ): Evidence for a Trans-Phyla Horizontal Transfer. Mol Biol Evol. 1995;12(3):421–31.PubMedGoogle Scholar
- Kreshchenko ND. Functions of flatworm neuropeptides NPF, GYIRF and FMRF in course of pharyngeal regeneration of anterior body fragments of planarian, Girardia tigrina. Acta Biol Hung. 2008;59(Suppl):199–207. doi:10.1556/ABiol.59.2008.Suppl.29.PubMedView ArticleGoogle Scholar
- Pag’an OR, Deats S, Baker D, Montgomery E, Wilk G, Tenaglia M, et al. Planarians require an intact brain to behaviorally react to cocaine, but not to react to nicotine. Neuroscience. 2013;246:265–70. doi:10.1016/j.neuroscience.2013.05.010.View ArticleGoogle Scholar
- Ramakrishnan L, Amatya C, DeSaer CJ, Dalhoff Z, Eggerichs MR. Galantamine reverses scopolamine-induced behavioral alterations in Dugesia tigrina. Invert Neurosci. 2014;14(2):91–101. doi:10.1007/s10158-013-0167-8.PubMedView ArticleGoogle Scholar
- Bolger AM, Lohse M, Usadel B.: Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 30(15):2114-20. doi:10.1093/bioinformatics/btu170Google Scholar
- Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63. doi:10.1093/bioinformatics/btr507.PubMed CentralPubMedView ArticleGoogle Scholar
- Grabherr M, Haas B, Yassour M. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nature. 2011;29(7):644–52. doi:10.1038/nbt.1883.Trinity.Google Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18(5):821–9. doi:10.1101/gr.074492.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012;28(8):1086–92. doi:10.1093/bioinformatics/bts094.PubMed CentralPubMedView ArticleGoogle Scholar
- Hornik, K.: The R FAQ (2014). http://cran.r-project.org/doc/FAQ/R-FAQ.html
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):25. doi:10.1186/gb-2009-10-3-r25.View ArticleGoogle Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. doi:10.1186/1471-2105-12-323.PubMed CentralPubMedView ArticleGoogle Scholar
- Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9. doi:10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- Conesa A, G’otz S, Garc’ıa-G’omez JM, Terol J, Tal’on M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6. doi:10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Lechner M, Findeiss S, Steiner L, Marz M, Stadler PF, Prohaska SJ. Proteinortho: detection of (co-)orthologs in large-scale analysis. BMC Bioinformatics. 2011;12(1):124. doi:10.1186/1471-2105-12-124.PubMed CentralPubMedView ArticleGoogle Scholar
- Logan-Klumpler FJ, De Silva N, Boehme U, Rogers MB, Velarde G, McQuillan J, et al. GeneDB–an annotation database for pathogens. Nucleic Acids Res. 2012;40(Database issue):98–108. doi:10.1093/nar/gkr1032.View ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45. doi:10.1101/gr.092759.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421. doi:10.1186/1471-2105-10-421.PubMed CentralPubMedView ArticleGoogle Scholar
- Pruitt KD, Brown GR, Hiatt SM, Thibaud-Nissen F, Astashyn A, Ermolaeva O, et al. RefSeq: an update on mammalian reference sequences. Nucleic Acids Res. 2014;42:756–63. doi:10.1093/nar/gkt1114.View ArticleGoogle Scholar
- Rask-Andersen M, Alm’en MS, Schïoth HB. Trends in the exploitation of novel drug targets. Nat Rev Drug Discov. 2011;10(8):579–90.PubMedView ArticleGoogle Scholar
- Milligan JN, Jolly ER. Cercarial transformation and in vitro cultivation of Schistosoma mansoni schistosomules. J Vis Exp. 2011;54:4–7. doi:10.3791/3191.Google Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5. doi:10.1038/nmeth.2089.PubMedView ArticleGoogle Scholar
- Noldus LP, Spink AJ, Tegelenbosch RA. EthoVision: a versatile video tracking system for automation of behavioral experiments. Behav Res Methods Instruments Comp. 2001;33(3):398–414.View ArticleGoogle Scholar
- Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007;23(10):1289–91. doi:10.1093/bioinformatics/btm091.PubMedView ArticleGoogle Scholar
- Rouhana L, Weiss J, Forsthoefel DJ, Lee H, King RS, Inoue T, et al. RNA interference by feeding in vitro-synthesized double-stranded RNA to planarians: methodology and dynamics. Dev Dyn. 2013;242(6):718–30. doi:10.1002/dvdy.23950.PubMed CentralPubMedView ArticleGoogle Scholar
- Kao D, Felix D, Aboobaker A. The planarian regeneration transcriptome reveals a shared but temporally shifted regulatory program between opposing head and tail scenarios. BMC Genomics. 2013;14:797. doi:10.1186/1471-2164-14-797.PubMed CentralPubMedView ArticleGoogle Scholar
- Patocka N, Ribeiro P. The functional role of a serotonin transporter in Schistosoma mansoni elucidated through immunolocalization and RNA interference (RNAi). Mol Biochem Parasitol. 2013;187(1):32–42.PubMedView ArticleGoogle Scholar
- Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics (Oxford, England). 2014;30(9):1236–40. doi:10.1093/bioinformatics/btu031.View ArticleGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42(Database issue):199–205. doi:10.1093/nar/gkt1076.View ArticleGoogle Scholar
- Hetrick B, Han MS, Helgeson LA, Nolen BJ. Small molecules CK-666 and CK-869 inhibit actin-related protein 2/3 complex by blocking an activating conformational change. Chem Biol. 2013;20(5):701–12. doi:10.1016/j.chembiol.2013.03.019.PubMed CentralPubMedView ArticleGoogle Scholar
- Nogi T, Zhang D, Chan JD, Marchant JS. A novel biological activity of praziquantel requiring voltage-operated Ca2+ channel beta subunits: subversion of flatworm regenerative polarity. PLoS Negl Trop Dis. 2009;3(6):464. doi:10.1371/journal.pntd.0000464.View ArticleGoogle Scholar
- Crowther GJ, Shanmugam D, Carmona SJ, Doyle MA, Hertz-Fowler C, Berriman M, et al. Identification of attractive drug targets in neglected-disease pathogens using an in silico approach. PLoS Neglected Trop Dis. 2010;4(8):804. doi:10.1371/journal.pntd.0000804.View ArticleGoogle Scholar
- Wang III B, Collins JJ, Newmark PA. Functional genomic characterization of neoblast-like stem cells in larval Schistosoma mansoni. eLife. 2013;2:1–15. doi:10.7554/eLife.00768.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.