BACKGROUND
the freshwater planarian schmidtea mediterranea has the remarkable ability to regenerate any part of its body following amputation, and uninjured animals constantly replace cells in the differentiated tissues through a homeostatic process .
planarians owe their remarkable regenerative capacity to a pool of adult stem cells called neoblasts  <cit> .
recent studies have detected heterogeneity within the neoblast population  <cit> , demonstrating that some of these cells are pluripotent and capable of producing differentiated cells of any type  <cit> , while others are specialized progenitors of particular cell types .
following injury, the neoblasts divide and migrate toward the wound site  <cit> , where the regeneration blastema forms in the event of tissue loss.
neoblast progeny then differentiate into the appropriate cell types and integrate with the pre-existing tissue, and regeneration culminates with full restoration of function.
many of the genes underlying these processes are evolutionarily conserved ; thus, planarians are an excellent model for the study of tissue regeneration and adult stem cell maintenance and function, and studying planarian homologs of genes involved in human diseases can provide insight into those genes’ functions that could be applied toward finding new treatments.
benefits to using planarians as a model system are the ability to do high throughput screens and the availability of protocols to analyze gene expression and function.
available tools include microarrays and mrna sequencing to study changes in gene expression, optimized in situ hybridization protocols to identify patterns of expression in specific tissues, and rna interference  for assessing function.
using these techniques, planarian researchers have begun characterizing genes necessary for stem cell function and tissue regeneration.
for example, high-throughput sequencing of transcriptomes of sorted stem cells and post-mitotic progeny has identified many genes expressed in neoblasts that may contribute to maintaining their pluripotency .
other screens identified genes required for the initial response to wounding  <cit>  or for regeneration of specific tissue types .
finally, smaller scale studies have shown that patterning of the regenerating tissues is guided in large part by conserved developmental pathways, including wnt/β-catenin, bmp, and slit-robo signaling .
questions still remain, however, regarding the full complement of genes required for planarian regeneration and how various cell types pattern relative to one another and the non-regenerating tissue.
we performed a functional genomics screen aimed at identifying previously uncharacterized genes involved in planarian regeneration.
we amputated one side of the head and used microarrays to analyze changes in gene expression in the regeneration blastema as well as in the non-regenerating side of the head.
in contrast to other screens that have examined differential gene expression during regeneration of the entire head  <cit> , our strategy allowed us to additionally look for genes required for repair  of tissues, including the central nervous system  and for signaling between the new and old tissue.
we next performed large scale in situ hybridization experiments to examine the expression patterns of the differentially expressed genes, and based on these patterns and homology, we chose  <dig> genes to knock down by rnai.
we identified  <dig> genes required for various aspects of planarian regeneration, five of which have also been functionally characterized by other labs .
the  <dig> newly characterized genes included those encoding homologs of the rna binding protein sart <dig> and chromatin remodeler smarcc- <dig>  which were necessary for neoblast maintenance, smed-ddc, which was needed for production of photoreceptor pigment, and smed-pgdb <dig>  which resulted in loss of the pharynx upon knockdown.
we also identified a homolog of nuclear transport factor importin-α and found it to be necessary for stem cell function, differentiation, and patterning of regenerating tissues relative to the midline.
this is notable because it suggests that regulated nuclear import of proteins may play a role in regeneration.
the results of this screen further our understanding of the molecular control of planarian regeneration by pinpointing additional genes involved in the animal’s remarkable ability to replace its tissues.
we have created a publicly accessible database to house our in situ and rnai results, making these data available to the planarian research community.
methods
planarian care
asexual schmidtea mediterranea from the ciw <dig> strain were maintained in 1x montjuïc salts and fed homogenized calf liver.
animals were starved for 1 week prior to use in experiments.
we used worms 1–2 mm in length for in situ hybridization and 3–5 mm for rnai experiments.
clones and accession numbers
constructs available in a collection of cdna clones  <cit>  were pulled from glycerol stocks, and inserts were subcloned into pjc <dig>   <cit>  for use in rnai experiments.
other sequences were directionally cloned into pjc <dig>  using gene-specific primers to amplify from cdna.
accession numbers or primers for each clone used in the screen are listed in additional file  <dig> and additional file  <dig>  clones used for in situ markers or follow-up experiments include smed-pc <dig> , smed-cintillo , smed-inx , smed-slit , smed-laminin  and smed-ima- <dig>   and  .
microarrays
samples were collected from intact and regenerating animals as described in the results section .
four separate biological replicates were processed for control and each of the blastema samples across the time course, and two were processed for the “opposite side” samples.
total rna was extracted with trizol reagent , treated with turbo dna-free dnase , then further purified using the qiagen rneasy minelute cleanup kit.
cdna was synthesized with the superscript double-stranded cdna synthesis kit .
cdna labeling and array hybridization were carried out by the genechip core at the university of california san diego.
each sample was hybridized to a separate nimblegen array  containing probes for  <dig>  ests representing approximately  <dig>  unique loci, with probes spotted in duplicate to provide technical replication.
microarray data analysis
nimblescan software was used to perform robust multi-array analysis  for quantile normalization and probe-level summarization of the expression data, with normalization occurring first between arrays of the same sample type and then across all samples.
background correction was applied during normalization.
probe sets with values below  <dig>  in all samples were omitted from further analysis, and expression of ests representing the same genomic locus was averaged to give one value per locus for each sample.
the limma bioconductor program  <cit>  was used to fit the expression data to a linear model and make pairwise contrasts between the control and each of the other six sample types .
the benjamini and hochberg correction  <cit>  was applied to adjust p-values to control the false discovery rate for multiple comparisons .
the microarray data have been deposited in the gene expression omnibus .
heat map generation
gene cluster  <dig>   <cit>  was used to produce the heat map of microarray data for the set of differentially expressed genes.
mean expression values for each gene in each sample type were log transformed, and the data were adjusted to center the genes relative to the overall mean.
a self-organizing map  was generated using the default settings to organize both genes and arrays and followed by hierarchical clustering using the average linkage method.
whole-mount in situ hybridization
riboprobes containing dig-11-utp were synthesized with t <dig> or t <dig> polymerase  from pcr product templates.
animals were fixed with formaldehyde and processed for whole-mount in situ hybridization as previously described  <cit>  using an intavis insitupro vs liquid handling robot.
worms for γ-irradiation experiments were fixed 3 days after a 60 gy treatment in a jl shepherd mark i cesium- <dig> irradiator; control and irradiated samples were processed side-by-side through fixation and hybridization and were developed for the same length of time.
antibody staining
worms stained with anti-arrestin  were fixed with formaldehyde using the same protocol as for in situ hybridization.
those stained with anti-phospho-histone h <dig>   were fixed with carnoy’s solution.
following fixation, staining was carried out as previously described  <cit> .
rna interference
planarians were fed bacteria induced to express double-stranded rna  against the gene of interest as previously described  <cit> .
bacterially expressed gfp dsrna was used as a negative control.
for the rnai screen, feedings were administered twice a week for three or six feedings, and animals were amputated pre-pharyngeally on the day after the last feeding, allowed to regenerate for 14 days, and then fixed for staining.
the number of feedings varied in other experiments as described in the results.
identification and phylogenetic analysis of importin-α proteins
to identify planarian homologs of importin-α proteins, importin-α sequences from species listed in fig.
<dig> were obtained from the uniprot database  <cit>  and used for tblastn searches against the planarian genome  <cit>  .
conserved protein domains, including the importin-β binding domain and armadillo repeats were identified using prosite  <cit> .
alignment of the amino acid sequences was performed using clustalw within mega <dig> software  <cit>  using the default settings.
the phylogenetic tree was constructed from this alignment by the neighbor-joining method in mega <dig> using the jones-taylor-thornton substitution model and settings for uniform mutation rates and pairwise gap deletion.
the reliability of the tree was tested by the bootstrap method with  <dig> replicates.
imaging
a leica m205c microscope outfitted with leica dfc <dig> camera was used to image live animals and those stained by colorimetric in situ hybridization.
fluorescent samples were imaged with a zeiss axiocam mrm camera mounted on a zeiss stereo lumar.v <dig> stereomicroscope or a zeiss axio observer.z <dig> inverted microscope outfitted with an apotome for optical sectioning.
cell counting
the number of ph <dig> positive cells was counted by hand from a single focal plane image of the whole worm and normalized to the area inside a hand-drawn outline of the worm calculated by axiozoom.
counts were performed using imagej  <dig> software  <cit> .
database creation
a relational database was developed using mysql to house multi-format data, and a web interface was created to intuitively interact with the database.
the database includes a customizable blast search , export and search functions, and image-viewing capability.
each record is uniquely identified by the est contig name.
annotation , expression, and phenotype data are searchable via the web interface, and search results can be downloaded into a comma-separated value file   <cit> .
RESULTS
changes in gene expression during head repair and regeneration
as a first step in our screen, we measured changes in gene expression on days two, three, and four of regeneration following amputation of one side of the head.
we chose to focus on the head and these time points because this is a region and period where many different cell and tissue types differentiate and pattern relative to one another; we were also interested in studying the regeneration and patterning of the cns.
the nascent cns can be visualized with the neural marker anti-synapsin by day three of regeneration and has connected to the cephalic ganglion on the non-amputated side of the head by day four .
removing only one side of the head allowed us to examine mechanisms of repair in addition to those required for de novo regeneration of the entire head; we also aimed to identify molecular interactions between the blastema and the non-regenerating tissue that was left behind on the opposite side.
figure 1b outlines our sample collection procedure.
we amputated half of the head, saved the tissue that was removed as the non-regenerating control sample, and then collected samples from the blastema and “opposite side” on days two through four of regeneration.
we extracted rna from each of these samples, reverse transcribed the rna to cdna and submitted these samples to the genechip microarray core at university of california san diego for labeling and hybridization to custom nimblegen microarrays  <cit> .
analysis of the resulting data revealed  <dig> genes with significant  changes in gene expression in one or more of the blastema or opposite side samples compared to control .
of these,  <dig> genes were upregulated in the blastema and  <dig> were downregulated.
for the majority of the upregulated genes, expression peaked on day two and then declined over the rest of the time course, but a cluster of  <dig> genes remained highly upregulated throughout .
on the opposite side of the head, we found  <dig> genes with increased expression and  <dig> with decreased expression; for  <dig> of these genes, the change in expression was unique to the opposite side tissue and not mirrored in the blastema.fig.
<dig> identification of genes differentially expressed during head repair and regeneration.
a half-head regeneration.
anti-synapsin antibody was used to label the central nervous system of intact animals and animals that had regenerated for two to 4 days following half-head amputation.
samples were imaged from the ventral side with anterior to the left, and the amputated portion of the head in regenerates is toward the bottom.
cg = cephalic ganglia; vnc = ventral nerve cords.
asterisks mark the pharynx, and triangles point to where the nervous system is regenerating in the blastema.
scale bar =  <dig>  mm.
b sample collection strategy.
one side of the head was removed from intact animals and saved as the non-regenerating control sample.
after two to four days of regeneration, samples were collected from the blastema and non-regenerating opposite side of the head.
c heat map summary of changes in gene expression in the blastema and other side of the head during a time course of half-head regeneration.
low expression is represented in blue and high in red.
control = non-regenerating samples.
day  <dig> b, day  <dig> b, and day  <dig> b = blastema samples collected on day two, three or four of regeneration.
day  <dig> os, day  <dig> os, and day  <dig> os = samples collected from the opposite side of the head on day two, three, or four of regeneration.
d categorization of differentially expressed genes based on the function of their homologs in other species using clusters of orthologous groups.
red bars show the number of upregulated genes in each functional group in either the blastema or opposite side tissue relative to control.
blue bars represent downregulated genes
to begin characterizing the differentially expressed genes, we grouped them into functional categories based on homology using eukaryotic clusters of orthologous groups  .
approximately half of the genes  could be definitively assigned to a kog group, and the others were unclassifiable due to low sequence conservation.
the genes upregulated in the blastema were divided across many functional groups.
the most highly represented category was “posttranslational modification, protein turnover, chaperones”, which included homologs of meprin a metalloproteinase and other proteases, proteasome subunits, and heat shock proteins.
the category of “replication, recombination and repair” was the next most common and contained several dna replication licensing factors , exo- and endonucleases involved in dna repair, and cell cycle protein cdc <dig>  among others; representation of genes in this category is likely due to the inclusion of some postblastema tissue in our samples.
notably, the “transcription” category, which was of interest to us because it contains proteins such as transcription factors and chromatin regulators likely to play important roles in controlling changes in gene expression during regeneration, was also highly represented among the differentially expressed genes.
nearly half  of the annotated downregulated genes and almost as many upregulated genes were identified as functioning in “signal transduction mechanisms”.
several genes potentially involved in neural function were among those downregulated in the signal transduction class, including homologs of three human gaba receptors  and a sodium-dependent noradrenaline transporter, slc6a <dig>
in samples from the side of the head opposite to the blastema of regenerates, differentially-expressed genes again fell into diverse functional categories.
genes upregulated in the opposite side tissue but not in the blastema included a homolog of kcnn <dig>  a small conductance calcium-activated potassium channel with neuronal function in humans.
“signal transduction” was again the most highly represented functional category among the downregulated genes and included wnt signaling protein frizzled and homologs of an acetylcholine receptor and a serine/threonine kinase.
given that most cell types in the animal are contained within the samples we collected, we expected to find a broad range of functional categories represented in our expression data.
identification of genes expressed in the blastema, neoblasts, and cns
based on our expectation that genes involved in head regeneration and repair would be expressed in the blastema, neoblasts, or cns, we further analyzed the genes identified from our microarray study by whole-mount in situ hybridization .
we performed these experiments to determine which cell and tissue types they are expressed in and to validate our microarray results.
additional file  <dig> provides a list of the genes tested and a summary of the results.
among  <dig> genes that were identified as upregulated in the blastema by our microarray analysis, we found  <dig> that were expressed at a higher level in the blastema than in the non-regenerating tissue by wish.
for some, the increased expression was tightly confined to the blastema itself , whereas others were strongly expressed both in the blastema and in the area beneath it  .
sixty genes appeared strongly upregulated in the blastema based on in situ hybridization, but surprisingly only one of these fell within the group of genes showing the largest fold-change by microarray .
this can be explained by noting that in non-regenerating animals, many of the genes with the largest fold-change were detected almost exclusively in the secretory cells surrounding the pharynx - any expression in the blastema represented a large change compared to the control head region, even if the resulting expression was not particularly strong.
taken together, the in situ staining validated the microarray results with regard to blastema expression.fig.
<dig> identification of genes expressed in tissues of interest and validation of microarray results.
representative examples from the wish screen.
a blastema expression.
the lower half of the head of the animal on the right in each pair was removed two to four days before fixation, on the day of greatest expression change for each gene.
arrows indicate the blastema.
b volcano plot of p-value vs. fold change for blastema to control comparison colored to indicate genes with no increased expression in the blastema , moderate blastema expression  or strong blastema expression .
c neoblast expression.
the animals on the right were irradiated  3 days before staining to destroy the neoblasts.
the animals on the left are non-irradiated.
d volcano plot colored to indicate neoblast expression; green marks genes detected in neoblasts and red marks non-neoblast genes.
genes were scored as expressed in neoblasts if staining was reduced following irradiation, as for histone h2a and jerky homolog-like.
e cns expression.
arrows and arrowheads indicate expression in the cephalic ganglia and ventral nerve cords, respectively.
f volcano plot colored to indicate genes with no cns expression , weak cns expression  or strong cns expression .
g venn diagram showing overlap of expression between tissues of interest.
h examples of genes not expressed in the blastema, neoblasts, or cns, but in the epidermis , intestine , secretory cells , or other discrete cells .
i examples of wish to genes with proposed opposite side of the head expression.
the lower half of the head of the animal on the right in each pair was amputated  <dig> to 4 days before fixation.
the arrow indicates cephalic ganglia expression.
animals were imaged ventrally, anterior to the left.
scale =  <dig>  mm
to identify genes expressed in the neoblasts, we performed wish in control animals and planarians treated with γ-irradiation to destroy the neoblasts  <cit> .
we found  <dig> genes with reduced staining in the irradiated samples.
some of these genes were expressed only in the neoblasts , whereas others were also strongly expressed in other tissues , including the cns and intestine .
the vast majority of neoblast-expressed genes were also expressed in the blastema , and many were found among the genes identified by microarray to be upregulated in the blastema .
we identified  <dig> genes expressed in the cns exclusively or in combination with other tissues.
some genes were expressed in the cns at levels similar to their expression in other tissues and were classified as simply having cns expression, whereas others were more highly expressed in the cns and were classified as “strong cns” genes .
genes with strong cns expression were enriched among those that appeared downregulated in blastema samples compared to control on the microarrays .
this is likely not due to actual downregulation but more representative of genes that are highly expressed in the cephalic ganglia of intact animals and that are not strongly upregulated in the blastema.
among the genes with strong cns expression, some appeared pan-neuronal  whereas others were expressed in a subset of neuronal cells or tissues  .
besides genes expressed in the blastema, neoblasts and/or the cns,  <dig> genes were found primarily in other tissues, including the epidermis, intestine, secretory cells, or other discrete cells .
we also performed in situ hybridization to  <dig> genes that were identified by microarray to have changed in expression only on the opposite side of the head  during regeneration.
of these, seven were identified as potentially upregulated in the opposite side tissue.
we did not observe an obvious increase in expression in the opposite side of the head over that of the surrounding tissue, although there may have been a more widespread increase in expression that included this region.
the staining pattern for several of these genes  was similar to that of the ca2+activated k+channel shown in fig.
2i, with expression mainly in the mesenchyme.
none of these seven genes were expressed in the cns.
in contrast, the majority of genes  identified as downregulated in the opposite side tissue showed strong cns expression, as exemplified by contig3034_se <dig> .
in total, we characterized the expression patterns of  <dig> genes, most of which had not been previously reported.
we have created a publicly accessible database to house our in situ data.
images of the expression patterns of all genes analyzed in our screen are available at http://planaria.sdsu.edu.
a targeted rnai screen to identify genes required for regeneration
we chose  <dig> genes from those differentially expressed during head regeneration to knock down by rnai .
we gave higher priority to genes that were expressed in tissues of interest  from our wish experiments or that we hypothesized would function in head/cns regeneration based on their homology.
planarians were fed double-stranded rna  twice per week for either three or six feedings and then amputated anterior to the pharynx and observed through 14 days of regeneration .
knockdown of  <dig> genes resulted in phenotypes , which included stem cell loss, impaired regeneration, effects on patterning or differentiation, and defects affecting specific tissues such as the photoreceptors and pharynx.
the results for genes in each of the phenotype categories are described below.fig.
<dig> representative images of phenotypes observed following rnai knockdown.
a rnai feeding schedules.
phenotypes in panels , , and  appeared after three feedings, panel  after  <dig> feedings, and the remainder after six feedings.
b-c loss of neoblasts in smed-slbp.
arrow indicates reduced-size blastema.
anti-ph <dig> labels the mitotic neoblasts.
live animals were imaged on day seven of regeneration, and stained animals were fixed on day  <dig>  d-e phenotypes of smed-dkc <dig>  arrows indicate missing blastemas.
d = day seven following pre-pharyngeal amputation, e = ventral curling and lysis on day  <dig>  f impaired regeneration following six smed-rbap46/48- <dig> dsrna feedings.
animals imaged on day seven.
g abnormally elongated photoreceptor pigment  in day seven regenerates following three feedings against smed-rbap46/48- <dig>  h asymmetric photoreceptors of smed-mcm <dig>  the red triangle indicates underdeveloped photoreceptor in day six regenerates.
i reduced blastema  and underdeveloped photoreceptors  in day six smed-fen <dig> regenerates.
j reduced blastema  and forked tail  in day seven smed-ptbp <dig> regenerates.
k
smed-tph animals lacking photoreceptor pigment on day  <dig> of regeneration.
red triangles indicate missing pigment.
l reduction of photoreceptor pigment in smed-ddc.
red triangles indicate photoreceptor pigment cups.
m elongated body shape following extended smed-tph rnai.
n edema following smed-gas <dig> rnai.
the animal shown is affected in the posterior.
o-q
smed-pgbd <dig> phenotypes.
white arrows indicate dorsal humps in uninjured animals  and day  <dig> regenerates .
the blue triangle indicates a post-pharyngeal lesion, and the green arrowhead indicates lateral bulging.
q wish to pharynx marker laminin in smed-pgbd <dig>  upper animals are shown ventrally, and lower animals were imaged from the side, with dorsal toward the top.
anterior is to the left.
yellow arrows and arrowheads indicate the base of the pharynx and its ventral opening, respectively.
brackets indicate the pharynx in controls.
scale =  <dig>  mm for all panels
smed–slbp
smed–sart3
smed–smarcc–1
smed–cycd
smed–espl1
smed–rrm2b–1
smed–rrm2b–2
smed–pdss2/dlp1
smed–dkc1
smed–emg1/nep1
smed–lig1
smed–prim2
smed–mcm7
pl08006b2e08
smed–rbap46/48–2
smed–mcm2
smed–ptbp1
smed–fen–1
smed–morf4l1/mrg–1
smed–ddc
smed–tph
smed–gas8
smed–pgbd4
smed–b9d2
smed–ima–1
only genes producing a phenotype are shown; see additional file  <dig> for a full list of genes tested.
gene names were assigned based on homology from blastx searches against the ncbi database.
e–values are the lower value between blastx of the cloned est or a longer sequence from published transcriptomes  <cit>  against the corresponding human protein.
kog group and functional category assignments were made using the eukaryotic clusters of orthologous groups database  <cit> .
an asterisk  indicates that loss of stem cells was verified by staining with anti–phospho–histone h3
differentially expressed genes required for neoblast survival/maintenance
knockdown of  <dig> genes resulted in phenotypes associated with loss of the stem cells, which include impaired blastema formation, head regression, lesions, ventral curling and lysis .
staining with anti-phospho-histone h <dig>  revealed a striking reduction in the number of mitotic cells relative to controls following rnai against four of the genes that had the “loss of stem cells” phenotype ; staining was not performed for the other genes in this group, in most cases due to death of the animals by lysis prior to the end of the two week observation period.
all  <dig> of the genes in this phenotypic class were upregulated in the blastema based on both microarray and wish data, and  <dig> of them  were expressed in the stem cells.
these data suggest that these genes play important roles in the survival or maintenance of planarian stem cells.
genes required for blastema formation and regenerative capacity
knockdown of five other genes resulted in a general reduction in regeneration ability without the dramatic loss of stem cells observed in the previous group.
for each of these genes, blastemas developed but were smaller than those in controls, and photoreceptors were delayed in formation or underdeveloped relative to controls at the same time point .
there was a three- to four-fold decrease in the number of mitotic cells following knockdown of smed-rbap46/48- <dig> , n =  <dig> animals/gene) and smed-mcm <dig> , n =  <dig> animals/gene).
the impaired regeneration phenotype was accompanied by additional defects in smed-rbap46/48- <dig>  smed-mcm <dig> and smed-ptbp <dig> animals.
in contrast to the underdeveloped photoreceptors observed following six dsrna feedings targeting smed-rbap46/48- <dig>  shortening the treatment to three feedings uncovered a defect in photoreceptor morphology in which the pigment cups appeared abnormally elongated at early stages of regeneration .
in some smed-mcm <dig> animals , the photoreceptors were asymmetric in size or only one photoreceptor fully developed .
finally, following knockdown of smed-ptbp <dig>  head fragments regenerating a new tail sometimes developed a forked blastema with two separate points rather than the typical tapered tail , and both head and trunk fragments displayed abnormal inching movements.
the knockdown results for this group of genes as a whole suggest they may act in cell differentiation and/or tissue patterning during regeneration.
other defects and effects on specific tissues or cell types
knockdown of five other genes resulted in phenotypes distinct from those described above in that they did not affect general stem cell function or blastema formation but instead had more specific effects.
two of these genes, dopa decarboxylase  and tryptophan 5-hydroxylase , were required for production of photoreceptor pigment; in trunks regenerating their heads, smed-tph animals lacked photoreceptor pigment entirely , and the pigment cups in smed-ddc animals were visible as faint brown spots rather than the usual dark black color .
both smed-tph and smed-ddc are expressed in the pigment cup   <cit> .
in addition to the photoreceptor phenotype described above, extended knockdown of smed-tph caused the worms to become longer and thinner .
the length-to-width ratio of non-regenerating smed-tph animals after  <dig> dsrna feedings was  <dig>  ±  <dig>  compared to  <dig>  ±  <dig>  for gfp controls .
these longer smed-tph animals also displayed abnormal inchworm-like movement rather than the usual gliding.
rnai against growth arrest-specific protein  <dig>  led to edema , a phenotype typically associated with dysfunction of the protonephridia.
finally, we observed movement defects  following knockdown of a homolog of b <dig> domain-containing protein  <dig> .
this phenotype was most noticeable shortly after the final dsrna feeding in both uninjured and regenerating animals and seemed to wear off over time after the rnai feedings stopped.
knockdown of a protein with weak homology to piggybac transposable element-derived protein  <dig>  led to the formation of a raised hump over the pharynx and a single lesion at the posterior end of the pharynx on the dorsal side of the animal .
in some  regenerating animals, we also observed lateral bulges near the pharynx .
staining with the pharynx marker laminin in smed-pgbd <dig> trunk fragments after 14 days regeneration or in uninjured animals following a period of 10 days starvation after the final dsrna feeding revealed that the pharynx was lost after this treatment .
the observation that knockdown of genes in this group leads to specific defects rather than a general loss of regenerative capacity suggests that their products may act in pathways that direct regeneration or maintenance of particular tissues or cell types, including the photoreceptors, ciliated cells and the pharynx.
planarians have two homologs of importin-α with distinct expression patterns and function
the final gene that produced a phenotype upon knockdown in our screen was a homolog of the nuclear transport factor importin-α that we named smed-ima- <dig>  knockdown of this gene led to defects in differentiation and patterning of regenerating tissues along with a decrease in the number of mitotic cells .
importin-α homologs in other species aid in the transport of nls-containing proteins into the nucleus by acting as an adaptor between the target proteins and importin-β, which in turn interacts with the nuclear pore complex  <cit> .
intrigued by the possibility that regulated nuclear import of proteins could be important for planarian regeneration, we decided to characterize the importin-α gene family in s. mediterranea.
in addition to smed-ima- <dig>  we identified another member of this family in the s. mediterranea genome, which we named smed-ima- <dig> .
each of these proteins contains two armadillo repeats, which may mediate protein-protein interactions, and an importin-β binding domain , which is required for the interaction between importin-α and importin-β .
a handful of other transcripts  shared some sequence similarity with this family but did not contain the importin-β binding domain  and could not be detected by in situ hybridization, suggesting that they are not functional homologs.
in contrast to smed-ima- <dig>  which was predominantly expressed in the stem cell population , smed-ima- <dig> was ubiquitously expressed but still present in the stem cells based on its reduction following γ-irradiation .
smed-ima- <dig> was strongly upregulated in the region beneath the blastema during regeneration.fig.
<dig> analysis of importin-α homologs in s. mediterranea.
a phylogenetic analysis of importin-α  homologs.
each protein is labeled with uniprot accession number, gene name, and species.
species abbreviations are as follows: human = homo sapiens, ponab = pongo abelii, bovin = bos taurus, mouse = mus musculus, danre = danio rerio, rat = rattus norvegicus, chick = gallus, yeast = saccharomyces cerevisiae, schpo = schizosaccharomyces pombe, dicdi = dictyostelium discoideum, arath = arabidopsis thaliana, orysj = oryza sativa subspecies japonica, sollc = solanum lycopersicum, smed = schmidtea mediterranea, caeel = caenorhabditis elegans, drome = drosophila melanogaster, xenla = xenopus laevis.
arrows indicate smed-ima- <dig> and smed-ima- <dig>  b domain structure of schmidtea mediterranea importin-α proteins.
ibb = importin-β binding domain.
arm = armadillo repeats.
c-d whole-mount in situ hybridization to smed-ima- <dig> and smed-ima- <dig>  the irradiated worms were treated with 60 gy γ-irradiation 3 days prior to fixation to destroy the stem cells.
the regenerating animals were amputated to remove the half of the head oriented toward the bottom of the picture two days prior to fixation.
arrows indicate the blastema.
animals were imaged from the ventral side with anterior to the left.
scale =  <dig>  mm
the two importin-α homologs had distinct functions based on rnai knockdown experiments.
knockdown of smed-ima- <dig> caused a rapid loss of stem cells accompanied by ventral curing and lysis , suggesting that smed-ima- <dig> may serve an essential function in bulk nuclear import of nls-containing proteins.
the smed-ima- <dig> rnai phenotype was more specific; when animals were fed smed-ima- <dig> dsrna four times over two weeks then amputated pre-pharangeally, the photoreceptors often appeared closer together than in gfp controls  .
in other animals, one photoreceptor was very small or absent , or a single cycloptic photoreceptor appeared near the center of the head  .
the abnormality of the photoreceptors extended to their neuronal connections , which we visualized by staining rnai-treated animals with an antibody against arrestin  <cit>  .
the cephalic ganglia in smed-ima- <dig> planarians were reduced in size, with an average area of  <dig>  ±  <dig>  mm2/mm animal length compared to  <dig>  ±  <dig>  mm2/mm animal length in controls  .
smed-ima- <dig> cephalic ganglia also appeared less developed than in controls and were also collapsed toward the midline  .
the smaller size of the cephalic ganglia may have been due to a defect in differentiation; this defect is further illustrated by a reduction in the number of regenerated sensory neurons, from  <dig>  ±  <dig>  cintillo positive cells/mm length in control to  <dig>  ±  <dig>  in smed-ima- <dig>  .
the regenerating tails of trunk fragments also developed abnormally in many cases, with the blastema growing asymmetrically slanted or forked  .
interestingly, we did not observe similar defects in the tail blastema morphology of regenerating head fragments.fig.
<dig>
smed-ima- <dig> is required for normal stem cell function and regeneration.
a-f images of control and smed-ima- <dig> animals fed bacterially expressed dsrna targeting each gene four times over 2 weeks then amputated transversely both pre- and post-pharyngeally.
animals were imaged or fixed on day  <dig> of regeneration.
anterior is toward the top.
live animals were imaged from the dorsal side and all others were imaged ventrally.
scale bars =  <dig>  mm for  and ,  <dig>  mm for  and ,  <dig>  mm for .
a live animals following rnai treatment.
dashed lines indicate amputation sites.
red triangles indicate photoreceptors forming abnormally close to the midline, and the black arrow indicates a mis-positioned and underdeveloped photoreceptor.
the white arrow indicates forking of the tail blastema.
b in situ hybridization to the neuronal marker smed-pc <dig>  white arrow indicates small cephalic ganglia collapsed toward the midline.
c staining with anti-arrestin antibody  to mark photoreceptor neurons.
arrows indicate aberrant neuronal projections.
d in situ hybridization to smed-cintillo, which labels sensory neurons.
e in situ hybridization to smed-inx to label the intestine.
arrows mark space between the two posterior intestinal branches.
f in situ hybridization to midline marker smed-slit.
g uninjured animals stained with anti-phospho-histone h <dig> antibody  to mark mitotic cells following six dsrna feedings.
anterior is to the left
the abnormal positioning of the photoreceptors and cephalic ganglia with respect to the midline led us to investigate whether smed-ima- <dig> caused other defects in midline patterning.
knockdown of smed-slit, a known regulator of midline patterning, leads to fusion of the two posterior branches of the intestine during regeneration  <cit> , however knockdown of smed-ima- <dig> did not lead to the same phenotype .
we also assayed smed-slit mrna expression by wish in smed-ima- <dig> knockdown animals and did not find any overt difference in its pattern or levels .
therefore, the midline collapse phenotype of smed-ima- <dig> rnai does not appear to be caused by a defect in smed-slit regulation.
we also observed a reduction in the number of mitotic cells in smed-ima- <dig> animals compared to controls.
after six feedings of dsrna over three weeks, the number of ph3+ cells in intact animals was reduced from 258 ±  <dig> cells/mm <dig> in gfp controls to 162 ±  <dig> cells/mm <dig> in smed-ima- <dig>  .
worms treated with dsrna against smed-ima- <dig> for extended periods of time  began to show other signs of stem cell dysfunction, including head regression and lysis.
the rnai phenotypes of the two smed-importin-α genes suggest that regulated nuclear import of proteins is a key factor in stem cell function.
our functional genomics screen identified several previously uncharacterized genes required for neoblast maintenance and blastema formation and others affecting regeneration or homeostatic maintenance of specific cell and tissue types.
this, along with our implication of a nuclear import factor in patterning of regenerating tissues provides new insights into the molecular basis of planarian regeneration.
discussion
we performed a screen to identify genes required for stem cell-based tissue replacement, approaching the problem from multiple angles.
first, we identified genes differentially expressed during regeneration of one half of the head.
we removed only one side of the head in order to look for genes required not only for de novo head regeneration, but also for those involved in repair and reconnection to tissue that remained following injury.
we examined expression over a time course from two to four days after amputation to go beyond the initial wound response to include the period when patterning and reconnection of tissues such as the cns occur.
working from the list of differentially expressed genes, we identified genes expressed in tissues of interest to us  and then examined the function of selected genes by rnai.
in the course of the screen, we identified genes involved in all aspects of regeneration, from neoblast maintenance to differentiation and patterning.
many of the genes required for neoblast maintenance are genes with essential cellular functions, such as ribosome biogenesis and dna replication.
others, however, may play a more specific role in cell cycle regulation or maintaining pluripotency.
for example, sart <dig> , an rna binding nuclear protein involved in pre-mrna splicing  <cit> , tumor rejection  <cit> , and hiv- <dig> replication  <cit> , has also been implicated in maintaining pluripotency of human embryonic stem cells  <cit> .
in addition, smed-smarcc- <dig>  which wenemoser et al.
<cit>  identified as upregulated as part of the initial wound response, is predicted based on homology to human smarcc <dig> to play a role in chromatin-mediated regulation of gene expression.
smarcc <dig> is part of the npbaf and nbaf chromatin remodeling complexes, which are required for proliferation of neural progenitors and dendrite growth in differentiated neurons, respectively  <cit> .
the loss of neoblasts phenotype in smed-smarcc- <dig> animals suggests it may act in a broader capacity in maintaining pluripotency and mitotic activity of planarian neoblasts.
the phenotypes we observed for other genes in the screen were also consistent with their putative function based on homology.
for example, ddc is required in other species for synthesis of melanin  <cit> , the primary pigment of the planarian photoreceptor  <cit> .
smed-ddc is expressed in the photoreceptors , and knockdown leads to a reduction in photoreceptor pigment, consistent with a function in melanin production.
rnai against smed-tph led to loss of photoreceptor pigment and changes in body morphology and movement, with animals becoming longer and thinner over an extended feeding scheme and losing their ability to glide smoothly.
our observation that smed-tph is involved in pigmentation is consistent with the reported role of smed-tph in eye melanogenesis  <cit> .
additionally, because smed-tph is required for serotonin biosynthesis , the smed-tph rnai movement defect is in line with previous observations that the coordinated wave-like motion of the ventral cilia of planarians is controlled by serotonin  <cit> .
along these same lines, the movement defects we observed following knockdown of smed-b9d <dig> in planarians are likely directly linked to a conserved function in the production of the ventral cilia the animals use for gliding.
b9d <dig> homologs play an important role in ciliogenesis in other species  <cit> , and loss-of-function mutations in this gene lead to meckel syndrome in humans, which is characterized by renal cysts and other abnormalities related to cilia dysfunction  <cit> .
further investigation of the potential role of smed-b9d <dig> in planarian ciliogenesis could yield information useful for understanding and treating this disease.
communication between the blastema and non-regenerating tissue
we identified  <dig> genes that displayed differential expression in comparisons between samples collected from the non-regenerating side of the head during regeneration and control samples from non-regenerating animals.
these genes are potentially involved in communication between the blastema and surrounding tissue to instruct the development of missing tissues in injured animals.
we knocked down expression of  <dig> of these genes, seven of which were predicted to be upregulated and seven downregulated during regeneration based on the microarrays.
among these experiments, we only observed a phenotype following knockdown of smed-tph, and this phenotype  was not consistent with smed-tph acting as a coordinator in cellular communication between the blastema and non-regenerating tissues.
our lack of success finding any genes specifically involved in this process may have been due to our focus on only genes with a significant change in expression during regeneration; if the genes we aimed to find were already expressed in the non-regenerating animals to coordinate normal cell turnover and did not react to the onset of regeneration with a change in expression, then we would have missed them in our screen.
expanding the screen to additional genes or using other criteria  may eventually uncover genes involved in the communication process.
nuclear import as a regulatory mechanism during regeneration
nuclear import of proteins that contain a nuclear localization signal  is facilitated by members of the karyopherin protein family, importin-α and importin-β.
in the classical nuclear import pathway, the nlss of target proteins are recognized and bound by importin-α, which serves as an adapter between the nls and importin-β  <cit> .
transport of the cargo into the nucleus is then accomplished through interactions of importin-β with the nuclear pore complex.
many species express several distinct importin-α proteins with variations in binding affinity for different nls sequences  <cit> .
this specificity allows individual importin-α family members to govern localization of unique sets of target proteins.
the ability to grant or limit access to the nucleus to other proteins such as transcription factors allows importin-α proteins to serve as regulatory molecules that affect cell fate  <cit> .
for example, control of the nuclear localization of transcription factor oct3/ <dig> by importin-α <dig> affects maintenance of mouse es cells in a pluripotent state, and control of brn <dig> and sox <dig> localization by importin-α <dig> regulates their differentiation toward a neural fate  <cit> .
additionally, the observation that different importin-α subtypes with distinct target specificities are expressed in various temporal and spatial patterns during mouse spermatogenesis has led to the proposal that regulated nuclear import plays a role in this process  <cit> .
rnai knockdown of the two importin-α homologs in schmidtea mediterranea led to distinct phenotypes.
we found that smed-ima- <dig> is required for general stem cell survival, while smed-ima- <dig> is important for stem cell maintenance, differentiation, and patterning of regenerating tissues.
we hypothesize that these proteins mediate nuclear localization of distinct groups of target proteins to achieve stem cell maintenance and regulate cellular differentiation and tissue patterning during regeneration.
interestingly, labbé et al.
previously demonstrated that knockdown of nuclear pore complex components nup <dig> and nup <dig> in planarians leads to defects in tissue homeostasis  <cit> .
together these results establish planarians as a model system for further investigation of the role of nuclear import in development, and identifying specific target proteins will be the next step toward a greater understanding of this process in vivo.
planarian in situ database
we have created an online database to make our in situ hybridization data available to the public.
it contains the sequence, annotation, expression, and rnai phenotype data, as well as links to images of in situ hybridization experiments where available.
a custom blast page was created to allow for searches against the est contig sequences.
the database is available at http://planaria.sdsu.edu.
CONCLUSIONS
we have identified new genes involved in various aspects of planarian regeneration and stem cell function.
most notably, characterization of smed-ima- <dig> as a key player in these processes opens the door to the study of controlled nuclear import of macromolecules as a regulatory mechanism in regeneration.
additionally, the online database we created to house our in situ data should serve as a valuable resource for the planarian research community.
availability of supporting data
the data sets supporting the results of this article are available in the gene expression omnibus  repository, accession gse <dig>  and at http://planaria.sdsu.edu/.
additional files
additional file 1:
identifies the ests matched to each gene id and provides their sequence.
additional file 2:
lists 1) all differentially expressed genes by gene id with blast2go and kog annotation for each gene; 2) normalized average expression levels for each gene tested in each sample type; 3) adjusted
p
-values
for each comparison performed with differentially expressed genes highlighted in red.
where a significant change was observed, the “expression change” column provides a description of the direction of the change  and in which tissue type  the change was detected.
additional file 3:
provides identifiers for 5′ and 3′ est sequence reads or primer sequences for clones used to make riboprobes.
“blastema”, “neoblasts” and “cns” columns indicate whether or not each gene was expressed in each of those tissues based on wish.
additional file 4:
lists all genes for which rnai experiments were performed, with identifiers for the 5′- and 3′-end sequence reads or primer sequences for each clone used.
the genes for which knockdown produced a phenotype are marked with the names they are referred to by in table  <dig>
additional file 5:
shows in situ hybridization to the genes  for which knockdown resulted in a phenotype.
animals were imaged from the ventral side and are shown with anterior to the left.
inset images for ddc and tph show photoreceptor expression imaged dorsally with anterior toward the top.
additional file 6:
shows univariate scatter plots of measurements/counts taken on regenerating
smed-ima-1
animals.
open circles represent values for individual animals and horizontal lines indicate the mean for each group.
6a shows differences in the size of the cephalic ganglia on day ten of regeneration in smed-ima- <dig> compared to gfp controls, and 6b shows differences in the number of sensory neurons  on the same day.
abbreviations
armarmadillo repeat
cgcephalic ganglia
cnscentral nervous system
dsrnadouble-stranded rna
ibbimportin-β binding domain
nlsnuclear localization signal
ph3phospho-histone h3s10
rnairna interference
smedgdschmidtea mediterranea genome database
vncventral nerve cords
wishwhole-mount in situ hybridization
competing interests
the authors declare that they have no competing interests.
authors’ contributions
ah designed, performed and supervised all the experimental work, analyzed the data and wrote the manuscript.
jmh, mwc, kgr, ca and cjs performed experiments.
mh created the database.
rmz conceived and supervised the project, and edited the manuscript.
all authors have read and approved the final manuscript.

 the freshwater planarian schmidtea mediterranea has the remarkable ability to regenerate any part of its body following amputation, and uninjured animals constantly replace cells in the differentiated tissues through a homeostatic process. planarians owe their remarkable regenerative capacity to a pool of adult stem cells called neoblasts. recent studies have detected heterogeneity within the neoblast population, demonstrating that some of these cells are pluripotent and capable of producing differentiated cells of any type, while others are specialized progenitors of particular cell types. following injury, the neoblasts divide and migrate toward the wound site, where the regeneration blastema forms in the event of tissue loss. neoblast progeny then differentiate into the appropriate cell types and integrate with the pre-existing tissue, and regeneration culminates with full restoration of function. many of the genes underlying these processes are evolutionarily conserved; thus, planarians are an excellent model for the study of tissue regeneration and adult stem cell maintenance and function, and studying planarian homologs of genes involved in human diseases can provide insight into those genes’ functions that could be applied toward finding new treatments. benefits to using planarians as a model system are the ability to do high throughput screens and the availability of protocols to analyze gene expression and function. available tools include microarrays and mrna sequencing to study changes in gene expression, optimized in situ hybridization protocols to identify patterns of expression in specific tissues, and rna interference for assessing function. using these techniques, planarian researchers have begun characterizing genes necessary for stem cell function and tissue regeneration. for example, high-throughput sequencing of transcriptomes of sorted stem cells and post-mitotic progeny has identified many genes expressed in neoblasts that may contribute to maintaining their pluripotency. other screens identified genes required for the initial response to wounding or for regeneration of specific tissue types. finally, smaller scale studies have shown that patterning of the regenerating tissues is guided in large part by conserved developmental pathways, including wnt/β-catenin, bmp, and slit-robo signaling. questions still remain, however, regarding the full complement of genes required for planarian regeneration and how various cell types pattern relative to one another and the non-regenerating tissue. we performed a functional genomics screen aimed at identifying previously uncharacterized genes involved in planarian regeneration. we amputated one side of the head and used microarrays to analyze changes in gene expression in the regeneration blastema as well as in the non-regenerating side of the head. in contrast to other screens that have examined differential gene expression during regeneration of the entire head, our strategy allowed us to additionally look for genes required for repair of tissues, including the central nervous system and for signaling between the new and old tissue. we next performed large scale in situ hybridization experiments to examine the expression patterns of the differentially expressed genes, and based on these patterns and homology, we chose genes to knock down by rnai. we identified genes required for various aspects of planarian regeneration, five of which have also been functionally characterized by other labs. the newly characterized genes included those encoding homologs of the rna binding protein sart and chromatin remodeler smarcc- which were necessary for neoblast maintenance, smed-ddc, which was needed for production of photoreceptor pigment, and smed-pgdb which resulted in loss of the pharynx upon knockdown. we also identified a homolog of nuclear transport factor importin-α and found it to be necessary for stem cell function, differentiation, and patterning of regenerating tissues relative to the midline. this is notable because it suggests that regulated nuclear import of proteins may play a role in regeneration. the results of this screen further our understanding of the molecular control of planarian regeneration by pinpointing additional genes involved in the animal’s remarkable ability to replace its tissues. we have created a publicly accessible database to house our in situ and rnai results, making these data available to the planarian research community. methods planarian care asexual schmidtea mediterranea from the ciw strain were maintained in 1x montjuïc salts and fed homogenized calf liver. animals were starved for 1 week prior to use in experiments. we used worms 1–2 mm in length for in situ hybridization and 3–5 mm for rnai experiments. clones and accession numbers constructs available in a collection of cdna clones were pulled from glycerol stocks, and inserts were subcloned into pjc for use in rnai experiments. other sequences were directionally cloned into pjc using gene-specific primers to amplify from cdna. accession numbers or primers for each clone used in the screen are listed in additional file and additional file clones used for in situ markers or follow-up experiments include smed-pc, smed-cintillo, smed-inx, smed-slit, smed-laminin and smed-ima- and. microarrays samples were collected from intact and regenerating animals as described in the results section. four separate biological replicates were processed for control and each of the blastema samples across the time course, and two were processed for the “opposite side” samples. total rna was extracted with trizol reagent, treated with turbo dna-free dnase, then further purified using the qiagen rneasy minelute cleanup kit. cdna was synthesized with the superscript double-stranded cdna synthesis kit. cdna labeling and array hybridization were carried out by the genechip core at the university of california san diego. each sample was hybridized to a separate nimblegen array containing probes for ests representing approximately unique loci, with probes spotted in duplicate to provide technical replication. microarray data analysis nimblescan software was used to perform robust multi-array analysis for quantile normalization and probe-level summarization of the expression data, with normalization occurring first between arrays of the same sample type and then across all samples. background correction was applied during normalization. probe sets with values below in all samples were omitted from further analysis, and expression of ests representing the same genomic locus was averaged to give one value per locus for each sample. the limma bioconductor program was used to fit the expression data to a linear model and make pairwise contrasts between the control and each of the other six sample types. the benjamini and hochberg correction was applied to adjust p-values to control the false discovery rate for multiple comparisons. the microarray data have been deposited in the gene expression omnibus. heat map generation gene cluster was used to produce the heat map of microarray data for the set of differentially expressed genes. mean expression values for each gene in each sample type were log transformed, and the data were adjusted to center the genes relative to the overall mean. a self-organizing map was generated using the default settings to organize both genes and arrays and followed by hierarchical clustering using the average linkage method. whole-mount in situ hybridization riboprobes containing dig-11-utp were synthesized with t or t polymerase from pcr product templates. animals were fixed with formaldehyde and processed for whole-mount in situ hybridization as previously described using an intavis insitupro vs liquid handling robot. worms for γ-irradiation experiments were fixed 3 days after a 60 gy treatment in a jl shepherd mark i cesium- irradiator; control and irradiated samples were processed side-by-side through fixation and hybridization and were developed for the same length of time. antibody staining worms stained with anti-arrestin were fixed with formaldehyde using the same protocol as for in situ hybridization. those stained with anti-phospho-histone h were fixed with carnoy’s solution. following fixation, staining was carried out as previously described. rna interference planarians were fed bacteria induced to express double-stranded rna against the gene of interest as previously described. bacterially expressed gfp dsrna was used as a negative control. for the rnai screen, feedings were administered twice a week for three or six feedings, and animals were amputated pre-pharyngeally on the day after the last feeding, allowed to regenerate for 14 days, and then fixed for staining. the number of feedings varied in other experiments as described in the results. identification and phylogenetic analysis of importin-α proteins to identify planarian homologs of importin-α proteins, importin-α sequences from species listed in fig.  were obtained from the uniprot database and used for tblastn searches against the planarian genome. conserved protein domains, including the importin-β binding domain and armadillo repeats were identified using prosite. alignment of the amino acid sequences was performed using clustalw within mega software using the default settings. the phylogenetic tree was constructed from this alignment by the neighbor-joining method in mega using the jones-taylor-thornton substitution model and settings for uniform mutation rates and pairwise gap deletion. the reliability of the tree was tested by the bootstrap method with replicates. imaging a leica m205c microscope outfitted with leica dfc camera was used to image live animals and those stained by colorimetric in situ hybridization. fluorescent samples were imaged with a zeiss axiocam mrm camera mounted on a zeiss stereo lumar.v stereomicroscope or a zeiss axio observer.z inverted microscope outfitted with an apotome for optical sectioning. cell counting the number of ph positive cells was counted by hand from a single focal plane image of the whole worm and normalized to the area inside a hand-drawn outline of the worm calculated by axiozoom. counts were performed using imagej software. database creation a relational database was developed using mysql to house multi-format data, and a web interface was created to intuitively interact with the database. the database includes a customizable blast search, export and search functions, and image-viewing capability. each record is uniquely identified by the est contig name. annotation, expression, and phenotype data are searchable via the web interface, and search results can be downloaded into a comma-separated value file.  changes in gene expression during head repair and regeneration as a first step in our screen, we measured changes in gene expression on days two, three, and four of regeneration following amputation of one side of the head. we chose to focus on the head and these time points because this is a region and period where many different cell and tissue types differentiate and pattern relative to one another; we were also interested in studying the regeneration and patterning of the cns. the nascent cns can be visualized with the neural marker anti-synapsin by day three of regeneration and has connected to the cephalic ganglion on the non-amputated side of the head by day four. removing only one side of the head allowed us to examine mechanisms of repair in addition to those required for de novo regeneration of the entire head; we also aimed to identify molecular interactions between the blastema and the non-regenerating tissue that was left behind on the opposite side. figure 1b outlines our sample collection procedure. we amputated half of the head, saved the tissue that was removed as the non-regenerating control sample, and then collected samples from the blastema and “opposite side” on days two through four of regeneration. we extracted rna from each of these samples, reverse transcribed the rna to cdna and submitted these samples to the genechip microarray core at university of california san diego for labeling and hybridization to custom nimblegen microarrays. analysis of the resulting data revealed genes with significant changes in gene expression in one or more of the blastema or opposite side samples compared to control. of these, genes were upregulated in the blastema and were downregulated. for the majority of the upregulated genes, expression peaked on day two and then declined over the rest of the time course, but a cluster of genes remained highly upregulated throughout. on the opposite side of the head, we found genes with increased expression and with decreased expression; for of these genes, the change in expression was unique to the opposite side tissue and not mirrored in the blastema.fig. identification of genes differentially expressed during head repair and regeneration. a half-head regeneration. anti-synapsin antibody was used to label the central nervous system of intact animals and animals that had regenerated for two to 4 days following half-head amputation. samples were imaged from the ventral side with anterior to the left, and the amputated portion of the head in regenerates is toward the bottom. cg = cephalic ganglia; vnc = ventral nerve cords. asterisks mark the pharynx, and triangles point to where the nervous system is regenerating in the blastema. scale bar =   mm. b sample collection strategy. one side of the head was removed from intact animals and saved as the non-regenerating control sample. after two to four days of regeneration, samples were collected from the blastema and non-regenerating opposite side of the head. c heat map summary of changes in gene expression in the blastema and other side of the head during a time course of half-head regeneration. low expression is represented in blue and high in red. control = non-regenerating samples. day b, day b, and day b = blastema samples collected on day two, three or four of regeneration. day os, day os, and day os = samples collected from the opposite side of the head on day two, three, or four of regeneration. d categorization of differentially expressed genes based on the function of their homologs in other species using clusters of orthologous groups. red bars show the number of upregulated genes in each functional group in either the blastema or opposite side tissue relative to control. blue bars represent downregulated genes to begin characterizing the differentially expressed genes, we grouped them into functional categories based on homology using eukaryotic clusters of orthologous groups. approximately half of the genes could be definitively assigned to a kog group, and the others were unclassifiable due to low sequence conservation. the genes upregulated in the blastema were divided across many functional groups. the most highly represented category was “posttranslational modification, protein turnover, chaperones”, which included homologs of meprin a metalloproteinase and other proteases, proteasome subunits, and heat shock proteins. the category of “replication, recombination and repair” was the next most common and contained several dna replication licensing factors, exo- and endonucleases involved in dna repair, and cell cycle protein cdc among others; representation of genes in this category is likely due to the inclusion of some postblastema tissue in our samples. notably, the “transcription” category, which was of interest to us because it contains proteins such as transcription factors and chromatin regulators likely to play important roles in controlling changes in gene expression during regeneration, was also highly represented among the differentially expressed genes. nearly half of the annotated downregulated genes and almost as many upregulated genes were identified as functioning in “signal transduction mechanisms”. several genes potentially involved in neural function were among those downregulated in the signal transduction class, including homologs of three human gaba receptors and a sodium-dependent noradrenaline transporter, slc6a in samples from the side of the head opposite to the blastema of regenerates, differentially-expressed genes again fell into diverse functional categories. genes upregulated in the opposite side tissue but not in the blastema included a homolog of kcnn a small conductance calcium-activated potassium channel with neuronal function in humans. “signal transduction” was again the most highly represented functional category among the downregulated genes and included wnt signaling protein frizzled and homologs of an acetylcholine receptor and a serine/threonine kinase. given that most cell types in the animal are contained within the samples we collected, we expected to find a broad range of functional categories represented in our expression data. identification of genes expressed in the blastema, neoblasts, and cns based on our expectation that genes involved in head regeneration and repair would be expressed in the blastema, neoblasts, or cns, we further analyzed the genes identified from our microarray study by whole-mount in situ hybridization. we performed these experiments to determine which cell and tissue types they are expressed in and to validate our microarray results. additional file provides a list of the genes tested and a summary of the results. among genes that were identified as upregulated in the blastema by our microarray analysis, we found that were expressed at a higher level in the blastema than in the non-regenerating tissue by wish. for some, the increased expression was tightly confined to the blastema itself, whereas others were strongly expressed both in the blastema and in the area beneath it. sixty genes appeared strongly upregulated in the blastema based on in situ hybridization, but surprisingly only one of these fell within the group of genes showing the largest fold-change by microarray. this can be explained by noting that in non-regenerating animals, many of the genes with the largest fold-change were detected almost exclusively in the secretory cells surrounding the pharynx - any expression in the blastema represented a large change compared to the control head region, even if the resulting expression was not particularly strong. taken together, the in situ staining validated the microarray results with regard to blastema expression.fig. identification of genes expressed in tissues of interest and validation of microarray results. representative examples from the wish screen. a blastema expression. the lower half of the head of the animal on the right in each pair was removed two to four days before fixation, on the day of greatest expression change for each gene. arrows indicate the blastema. b volcano plot of p-value vs. fold change for blastema to control comparison colored to indicate genes with no increased expression in the blastema, moderate blastema expression or strong blastema expression. c neoblast expression. the animals on the right were irradiated 3 days before staining to destroy the neoblasts. the animals on the left are non-irradiated. d volcano plot colored to indicate neoblast expression; green marks genes detected in neoblasts and red marks non-neoblast genes. genes were scored as expressed in neoblasts if staining was reduced following irradiation, as for histone h2a and jerky homolog-like. e cns expression. arrows and arrowheads indicate expression in the cephalic ganglia and ventral nerve cords, respectively. f volcano plot colored to indicate genes with no cns expression, weak cns expression or strong cns expression. g venn diagram showing overlap of expression between tissues of interest. h examples of genes not expressed in the blastema, neoblasts, or cns, but in the epidermis, intestine, secretory cells, or other discrete cells. i examples of wish to genes with proposed opposite side of the head expression. the lower half of the head of the animal on the right in each pair was amputated to 4 days before fixation. the arrow indicates cephalic ganglia expression. animals were imaged ventrally, anterior to the left. scale =   mm to identify genes expressed in the neoblasts, we performed wish in control animals and planarians treated with γ-irradiation to destroy the neoblasts. we found genes with reduced staining in the irradiated samples. some of these genes were expressed only in the neoblasts, whereas others were also strongly expressed in other tissues, including the cns and intestine. the vast majority of neoblast-expressed genes were also expressed in the blastema, and many were found among the genes identified by microarray to be upregulated in the blastema. we identified genes expressed in the cns exclusively or in combination with other tissues. some genes were expressed in the cns at levels similar to their expression in other tissues and were classified as simply having cns expression, whereas others were more highly expressed in the cns and were classified as “strong cns” genes. genes with strong cns expression were enriched among those that appeared downregulated in blastema samples compared to control on the microarrays. this is likely not due to actual downregulation but more representative of genes that are highly expressed in the cephalic ganglia of intact animals and that are not strongly upregulated in the blastema. among the genes with strong cns expression, some appeared pan-neuronal whereas others were expressed in a subset of neuronal cells or tissues. besides genes expressed in the blastema, neoblasts and/or the cns, genes were found primarily in other tissues, including the epidermis, intestine, secretory cells, or other discrete cells. we also performed in situ hybridization to genes that were identified by microarray to have changed in expression only on the opposite side of the head during regeneration. of these, seven were identified as potentially upregulated in the opposite side tissue. we did not observe an obvious increase in expression in the opposite side of the head over that of the surrounding tissue, although there may have been a more widespread increase in expression that included this region. the staining pattern for several of these genes was similar to that of the ca2+activated k+channel shown in fig. 2i, with expression mainly in the mesenchyme. none of these seven genes were expressed in the cns. in contrast, the majority of genes identified as downregulated in the opposite side tissue showed strong cns expression, as exemplified by contig3034_se. in total, we characterized the expression patterns of genes, most of which had not been previously reported. we have created a publicly accessible database to house our in situ data. images of the expression patterns of all genes analyzed in our screen are available at http://planaria.sdsu.edu. a targeted rnai screen to identify genes required for regeneration we chose genes from those differentially expressed during head regeneration to knock down by rnai. we gave higher priority to genes that were expressed in tissues of interest from our wish experiments or that we hypothesized would function in head/cns regeneration based on their homology. planarians were fed double-stranded rna twice per week for either three or six feedings and then amputated anterior to the pharynx and observed through 14 days of regeneration. knockdown of genes resulted in phenotypes, which included stem cell loss, impaired regeneration, effects on patterning or differentiation, and defects affecting specific tissues such as the photoreceptors and pharynx. the results for genes in each of the phenotype categories are described below.fig. representative images of phenotypes observed following rnai knockdown. a rnai feeding schedules. phenotypes in panels,, and appeared after three feedings, panel after feedings, and the remainder after six feedings. b-c loss of neoblasts in smed-slbp. arrow indicates reduced-size blastema. anti-ph labels the mitotic neoblasts. live animals were imaged on day seven of regeneration, and stained animals were fixed on day d-e phenotypes of smed-dkc arrows indicate missing blastemas. d = day seven following pre-pharyngeal amputation, e = ventral curling and lysis on day f impaired regeneration following six smed-rbap46/48- dsrna feedings. animals imaged on day seven. g abnormally elongated photoreceptor pigment in day seven regenerates following three feedings against smed-rbap46/48- h asymmetric photoreceptors of smed-mcm the red triangle indicates underdeveloped photoreceptor in day six regenerates. i reduced blastema and underdeveloped photoreceptors in day six smed-fen regenerates. j reduced blastema and forked tail in day seven smed-ptbp regenerates. k smed-tph animals lacking photoreceptor pigment on day of regeneration. red triangles indicate missing pigment. l reduction of photoreceptor pigment in smed-ddc. red triangles indicate photoreceptor pigment cups. m elongated body shape following extended smed-tph rnai. n edema following smed-gas rnai. the animal shown is affected in the posterior. o-q smed-pgbd phenotypes. white arrows indicate dorsal humps in uninjured animals and day regenerates. the blue triangle indicates a post-pharyngeal lesion, and the green arrowhead indicates lateral bulging. q wish to pharynx marker laminin in smed-pgbd upper animals are shown ventrally, and lower animals were imaged from the side, with dorsal toward the top. anterior is to the left. yellow arrows and arrowheads indicate the base of the pharynx and its ventral opening, respectively. brackets indicate the pharynx in controls. scale =   mm for all panels smed–slbp smed–sart3 smed–smarcc–1 smed–cycd smed–espl1 smed–rrm2b–1 smed–rrm2b–2 smed–pdss2/dlp1 smed–dkc1 smed–emg1/nep1 smed–lig1 smed–prim2 smed–mcm7 pl08006b2e08 smed–rbap46/48–2 smed–mcm2 smed–ptbp1 smed–fen–1 smed–morf4l1/mrg–1 smed–ddc smed–tph smed–gas8 smed–pgbd4 smed–b9d2 smed–ima–1 only genes producing a phenotype are shown; see additional file for a full list of genes tested. gene names were assigned based on homology from blastx searches against the ncbi database. e–values are the lower value between blastx of the cloned est or a longer sequence from published transcriptomes against the corresponding human protein. kog group and functional category assignments were made using the eukaryotic clusters of orthologous groups database. an asterisk indicates that loss of stem cells was verified by staining with anti–phospho–histone h3 differentially expressed genes required for neoblast survival/maintenance knockdown of genes resulted in phenotypes associated with loss of the stem cells, which include impaired blastema formation, head regression, lesions, ventral curling and lysis. staining with anti-phospho-histone h revealed a striking reduction in the number of mitotic cells relative to controls following rnai against four of the genes that had the “loss of stem cells” phenotype; staining was not performed for the other genes in this group, in most cases due to death of the animals by lysis prior to the end of the two week observation period. all of the genes in this phenotypic class were upregulated in the blastema based on both microarray and wish data, and of them were expressed in the stem cells. these data suggest that these genes play important roles in the survival or maintenance of planarian stem cells. genes required for blastema formation and regenerative capacity knockdown of five other genes resulted in a general reduction in regeneration ability without the dramatic loss of stem cells observed in the previous group. for each of these genes, blastemas developed but were smaller than those in controls, and photoreceptors were delayed in formation or underdeveloped relative to controls at the same time point. there was a three- to four-fold decrease in the number of mitotic cells following knockdown of smed-rbap46/48-, n =  animals/gene) and smed-mcm, n =  animals/gene). the impaired regeneration phenotype was accompanied by additional defects in smed-rbap46/48- smed-mcm and smed-ptbp animals. in contrast to the underdeveloped photoreceptors observed following six dsrna feedings targeting smed-rbap46/48- shortening the treatment to three feedings uncovered a defect in photoreceptor morphology in which the pigment cups appeared abnormally elongated at early stages of regeneration. in some smed-mcm animals, the photoreceptors were asymmetric in size or only one photoreceptor fully developed. finally, following knockdown of smed-ptbp head fragments regenerating a new tail sometimes developed a forked blastema with two separate points rather than the typical tapered tail, and both head and trunk fragments displayed abnormal inching movements. the knockdown results for this group of genes as a whole suggest they may act in cell differentiation and/or tissue patterning during regeneration. other defects and effects on specific tissues or cell types knockdown of five other genes resulted in phenotypes distinct from those described above in that they did not affect general stem cell function or blastema formation but instead had more specific effects. two of these genes, dopa decarboxylase and tryptophan 5-hydroxylase, were required for production of photoreceptor pigment; in trunks regenerating their heads, smed-tph animals lacked photoreceptor pigment entirely, and the pigment cups in smed-ddc animals were visible as faint brown spots rather than the usual dark black color. both smed-tph and smed-ddc are expressed in the pigment cup. in addition to the photoreceptor phenotype described above, extended knockdown of smed-tph caused the worms to become longer and thinner. the length-to-width ratio of non-regenerating smed-tph animals after dsrna feedings was  ±  compared to  ±  for gfp controls. these longer smed-tph animals also displayed abnormal inchworm-like movement rather than the usual gliding. rnai against growth arrest-specific protein led to edema, a phenotype typically associated with dysfunction of the protonephridia. finally, we observed movement defects following knockdown of a homolog of b domain-containing protein. this phenotype was most noticeable shortly after the final dsrna feeding in both uninjured and regenerating animals and seemed to wear off over time after the rnai feedings stopped. knockdown of a protein with weak homology to piggybac transposable element-derived protein led to the formation of a raised hump over the pharynx and a single lesion at the posterior end of the pharynx on the dorsal side of the animal. in some regenerating animals, we also observed lateral bulges near the pharynx. staining with the pharynx marker laminin in smed-pgbd trunk fragments after 14 days regeneration or in uninjured animals following a period of 10 days starvation after the final dsrna feeding revealed that the pharynx was lost after this treatment. the observation that knockdown of genes in this group leads to specific defects rather than a general loss of regenerative capacity suggests that their products may act in pathways that direct regeneration or maintenance of particular tissues or cell types, including the photoreceptors, ciliated cells and the pharynx. planarians have two homologs of importin-α with distinct expression patterns and function the final gene that produced a phenotype upon knockdown in our screen was a homolog of the nuclear transport factor importin-α that we named smed-ima- knockdown of this gene led to defects in differentiation and patterning of regenerating tissues along with a decrease in the number of mitotic cells. importin-α homologs in other species aid in the transport of nls-containing proteins into the nucleus by acting as an adaptor between the target proteins and importin-β, which in turn interacts with the nuclear pore complex. intrigued by the possibility that regulated nuclear import of proteins could be important for planarian regeneration, we decided to characterize the importin-α gene family in s. mediterranea. in addition to smed-ima- we identified another member of this family in the s. mediterranea genome, which we named smed-ima-. each of these proteins contains two armadillo repeats, which may mediate protein-protein interactions, and an importin-β binding domain, which is required for the interaction between importin-α and importin-β. a handful of other transcripts shared some sequence similarity with this family but did not contain the importin-β binding domain and could not be detected by in situ hybridization, suggesting that they are not functional homologs. in contrast to smed-ima- which was predominantly expressed in the stem cell population, smed-ima- was ubiquitously expressed but still present in the stem cells based on its reduction following γ-irradiation. smed-ima- was strongly upregulated in the region beneath the blastema during regeneration.fig. analysis of importin-α homologs in s. mediterranea. a phylogenetic analysis of importin-α homologs. each protein is labeled with uniprot accession number, gene name, and species. species abbreviations are as follows: human = homo sapiens, ponab = pongo abelii, bovin = bos taurus, mouse = mus musculus, danre = danio rerio, rat = rattus norvegicus, chick = gallus, yeast = saccharomyces cerevisiae, schpo = schizosaccharomyces pombe, dicdi = dictyostelium discoideum, arath = arabidopsis thaliana, orysj = oryza sativa subspecies japonica, sollc = solanum lycopersicum, smed = schmidtea mediterranea, caeel = caenorhabditis elegans, drome = drosophila melanogaster, xenla = xenopus laevis. arrows indicate smed-ima- and smed-ima- b domain structure of schmidtea mediterranea importin-α proteins. ibb = importin-β binding domain. arm = armadillo repeats. c-d whole-mount in situ hybridization to smed-ima- and smed-ima- the irradiated worms were treated with 60 gy γ-irradiation 3 days prior to fixation to destroy the stem cells. the regenerating animals were amputated to remove the half of the head oriented toward the bottom of the picture two days prior to fixation. arrows indicate the blastema. animals were imaged from the ventral side with anterior to the left. scale =   mm the two importin-α homologs had distinct functions based on rnai knockdown experiments. knockdown of smed-ima- caused a rapid loss of stem cells accompanied by ventral curing and lysis, suggesting that smed-ima- may serve an essential function in bulk nuclear import of nls-containing proteins. the smed-ima- rnai phenotype was more specific; when animals were fed smed-ima- dsrna four times over two weeks then amputated pre-pharangeally, the photoreceptors often appeared closer together than in gfp controls. in other animals, one photoreceptor was very small or absent, or a single cycloptic photoreceptor appeared near the center of the head. the abnormality of the photoreceptors extended to their neuronal connections, which we visualized by staining rnai-treated animals with an antibody against arrestin. the cephalic ganglia in smed-ima- planarians were reduced in size, with an average area of  ±   mm2/mm animal length compared to  ±   mm2/mm animal length in controls. smed-ima- cephalic ganglia also appeared less developed than in controls and were also collapsed toward the midline. the smaller size of the cephalic ganglia may have been due to a defect in differentiation; this defect is further illustrated by a reduction in the number of regenerated sensory neurons, from  ±  cintillo positive cells/mm length in control to  ±  in smed-ima-. the regenerating tails of trunk fragments also developed abnormally in many cases, with the blastema growing asymmetrically slanted or forked. interestingly, we did not observe similar defects in the tail blastema morphology of regenerating head fragments.fig. smed-ima- is required for normal stem cell function and regeneration. a-f images of control and smed-ima- animals fed bacterially expressed dsrna targeting each gene four times over 2 weeks then amputated transversely both pre- and post-pharyngeally. animals were imaged or fixed on day of regeneration. anterior is toward the top. live animals were imaged from the dorsal side and all others were imaged ventrally. scale bars =   mm for and,  mm for and,  mm for. a live animals following rnai treatment. dashed lines indicate amputation sites. red triangles indicate photoreceptors forming abnormally close to the midline, and the black arrow indicates a mis-positioned and underdeveloped photoreceptor. the white arrow indicates forking of the tail blastema. b in situ hybridization to the neuronal marker smed-pc white arrow indicates small cephalic ganglia collapsed toward the midline. c staining with anti-arrestin antibody to mark photoreceptor neurons. arrows indicate aberrant neuronal projections. d in situ hybridization to smed-cintillo, which labels sensory neurons. e in situ hybridization to smed-inx to label the intestine. arrows mark space between the two posterior intestinal branches. f in situ hybridization to midline marker smed-slit. g uninjured animals stained with anti-phospho-histone h antibody to mark mitotic cells following six dsrna feedings. anterior is to the left the abnormal positioning of the photoreceptors and cephalic ganglia with respect to the midline led us to investigate whether smed-ima- caused other defects in midline patterning. knockdown of smed-slit, a known regulator of midline patterning, leads to fusion of the two posterior branches of the intestine during regeneration, however knockdown of smed-ima- did not lead to the same phenotype. we also assayed smed-slit mrna expression by wish in smed-ima- knockdown animals and did not find any overt difference in its pattern or levels. therefore, the midline collapse phenotype of smed-ima- rnai does not appear to be caused by a defect in smed-slit regulation. we also observed a reduction in the number of mitotic cells in smed-ima- animals compared to controls. after six feedings of dsrna over three weeks, the number of ph3+ cells in intact animals was reduced from 258 ±  cells/mm in gfp controls to 162 ±  cells/mm in smed-ima-. worms treated with dsrna against smed-ima- for extended periods of time began to show other signs of stem cell dysfunction, including head regression and lysis. the rnai phenotypes of the two smed-importin-α genes suggest that regulated nuclear import of proteins is a key factor in stem cell function. our functional genomics screen identified several previously uncharacterized genes required for neoblast maintenance and blastema formation and others affecting regeneration or homeostatic maintenance of specific cell and tissue types. this, along with our implication of a nuclear import factor in patterning of regenerating tissues provides new insights into the molecular basis of planarian regeneration. discussion we performed a screen to identify genes required for stem cell-based tissue replacement, approaching the problem from multiple angles. first, we identified genes differentially expressed during regeneration of one half of the head. we removed only one side of the head in order to look for genes required not only for de novo head regeneration, but also for those involved in repair and reconnection to tissue that remained following injury. we examined expression over a time course from two to four days after amputation to go beyond the initial wound response to include the period when patterning and reconnection of tissues such as the cns occur. working from the list of differentially expressed genes, we identified genes expressed in tissues of interest to us and then examined the function of selected genes by rnai. in the course of the screen, we identified genes involved in all aspects of regeneration, from neoblast maintenance to differentiation and patterning. many of the genes required for neoblast maintenance are genes with essential cellular functions, such as ribosome biogenesis and dna replication. others, however, may play a more specific role in cell cycle regulation or maintaining pluripotency. for example, sart, an rna binding nuclear protein involved in pre-mrna splicing, tumor rejection, and hiv- replication, has also been implicated in maintaining pluripotency of human embryonic stem cells. in addition, smed-smarcc- which wenemoser et al. identified as upregulated as part of the initial wound response, is predicted based on homology to human smarcc to play a role in chromatin-mediated regulation of gene expression. smarcc is part of the npbaf and nbaf chromatin remodeling complexes, which are required for proliferation of neural progenitors and dendrite growth in differentiated neurons, respectively. the loss of neoblasts phenotype in smed-smarcc- animals suggests it may act in a broader capacity in maintaining pluripotency and mitotic activity of planarian neoblasts. the phenotypes we observed for other genes in the screen were also consistent with their putative function based on homology. for example, ddc is required in other species for synthesis of melanin, the primary pigment of the planarian photoreceptor. smed-ddc is expressed in the photoreceptors, and knockdown leads to a reduction in photoreceptor pigment, consistent with a function in melanin production. rnai against smed-tph led to loss of photoreceptor pigment and changes in body morphology and movement, with animals becoming longer and thinner over an extended feeding scheme and losing their ability to glide smoothly. our observation that smed-tph is involved in pigmentation is consistent with the reported role of smed-tph in eye melanogenesis. additionally, because smed-tph is required for serotonin biosynthesis, the smed-tph rnai movement defect is in line with previous observations that the coordinated wave-like motion of the ventral cilia of planarians is controlled by serotonin. along these same lines, the movement defects we observed following knockdown of smed-b9d in planarians are likely directly linked to a conserved function in the production of the ventral cilia the animals use for gliding. b9d homologs play an important role in ciliogenesis in other species, and loss-of-function mutations in this gene lead to meckel syndrome in humans, which is characterized by renal cysts and other abnormalities related to cilia dysfunction. further investigation of the potential role of smed-b9d in planarian ciliogenesis could yield information useful for understanding and treating this disease. communication between the blastema and non-regenerating tissue we identified genes that displayed differential expression in comparisons between samples collected from the non-regenerating side of the head during regeneration and control samples from non-regenerating animals. these genes are potentially involved in communication between the blastema and surrounding tissue to instruct the development of missing tissues in injured animals. we knocked down expression of of these genes, seven of which were predicted to be upregulated and seven downregulated during regeneration based on the microarrays. among these experiments, we only observed a phenotype following knockdown of smed-tph, and this phenotype was not consistent with smed-tph acting as a coordinator in cellular communication between the blastema and non-regenerating tissues. our lack of success finding any genes specifically involved in this process may have been due to our focus on only genes with a significant change in expression during regeneration; if the genes we aimed to find were already expressed in the non-regenerating animals to coordinate normal cell turnover and did not react to the onset of regeneration with a change in expression, then we would have missed them in our screen. expanding the screen to additional genes or using other criteria may eventually uncover genes involved in the communication process. nuclear import as a regulatory mechanism during regeneration nuclear import of proteins that contain a nuclear localization signal is facilitated by members of the karyopherin protein family, importin-α and importin-β. in the classical nuclear import pathway, the nlss of target proteins are recognized and bound by importin-α, which serves as an adapter between the nls and importin-β. transport of the cargo into the nucleus is then accomplished through interactions of importin-β with the nuclear pore complex. many species express several distinct importin-α proteins with variations in binding affinity for different nls sequences. this specificity allows individual importin-α family members to govern localization of unique sets of target proteins. the ability to grant or limit access to the nucleus to other proteins such as transcription factors allows importin-α proteins to serve as regulatory molecules that affect cell fate. for example, control of the nuclear localization of transcription factor oct3/ by importin-α affects maintenance of mouse es cells in a pluripotent state, and control of brn and sox localization by importin-α regulates their differentiation toward a neural fate. additionally, the observation that different importin-α subtypes with distinct target specificities are expressed in various temporal and spatial patterns during mouse spermatogenesis has led to the proposal that regulated nuclear import plays a role in this process. rnai knockdown of the two importin-α homologs in schmidtea mediterranea led to distinct phenotypes. we found that smed-ima- is required for general stem cell survival, while smed-ima- is important for stem cell maintenance, differentiation, and patterning of regenerating tissues. we hypothesize that these proteins mediate nuclear localization of distinct groups of target proteins to achieve stem cell maintenance and regulate cellular differentiation and tissue patterning during regeneration. interestingly, labbé et al. previously demonstrated that knockdown of nuclear pore complex components nup and nup in planarians leads to defects in tissue homeostasis. together these results establish planarians as a model system for further investigation of the role of nuclear import in development, and identifying specific target proteins will be the next step toward a greater understanding of this process in vivo. planarian in situ database we have created an online database to make our in situ hybridization data available to the public. it contains the sequence, annotation, expression, and rnai phenotype data, as well as links to images of in situ hybridization experiments where available. a custom blast page was created to allow for searches against the est contig sequences. the database is available at http://planaria.sdsu.edu.  we have identified new genes involved in various aspects of planarian regeneration and stem cell function. most notably, characterization of smed-ima- as a key player in these processes opens the door to the study of controlled nuclear import of macromolecules as a regulatory mechanism in regeneration. additionally, the online database we created to house our in situ data should serve as a valuable resource for the planarian research community. availability of supporting data the data sets supporting the results of this article are available in the gene expression omnibus repository, accession gse and at http://planaria.sdsu.edu/. additional files additional file 1: identifies the ests matched to each gene id and provides their sequence. additional file 2: lists 1) all differentially expressed genes by gene id with blast2go and kog annotation for each gene; 2) normalized average expression levels for each gene tested in each sample type; 3) adjusted p -values for each comparison performed with differentially expressed genes highlighted in red. where a significant change was observed, the “expression change” column provides a description of the direction of the change and in which tissue type the change was detected. additional file 3: provides identifiers for 5′ and 3′ est sequence reads or primer sequences for clones used to make riboprobes. “blastema”, “neoblasts” and “cns” columns indicate whether or not each gene was expressed in each of those tissues based on wish. additional file 4: lists all genes for which rnai experiments were performed, with identifiers for the 5′- and 3′-end sequence reads or primer sequences for each clone used. the genes for which knockdown produced a phenotype are marked with the names they are referred to by in table  additional file 5: shows in situ hybridization to the genes for which knockdown resulted in a phenotype. animals were imaged from the ventral side and are shown with anterior to the left. inset images for ddc and tph show photoreceptor expression imaged dorsally with anterior toward the top. additional file 6: shows univariate scatter plots of measurements/counts taken on regenerating smed-ima-1 animals. open circles represent values for individual animals and horizontal lines indicate the mean for each group. 6a shows differences in the size of the cephalic ganglia on day ten of regeneration in smed-ima- compared to gfp controls, and 6b shows differences in the number of sensory neurons on the same day. abbreviations armarmadillo repeat cgcephalic ganglia cnscentral nervous system dsrnadouble-stranded rna ibbimportin-β binding domain nlsnuclear localization signal ph3phospho-histone h3s10 rnairna interference smedgdschmidtea mediterranea genome database vncventral nerve cords wishwhole-mount in situ hybridization competing interests the authors declare that they have no competing interests. authors’ contributions ah designed, performed and supervised all the experimental work, analyzed the data and wrote the manuscript. jmh, mwc, kgr, ca and cjs performed experiments. mh created the database. rmz conceived and supervised the project, and edited the manuscript. all authors have read and approved the final manuscript. 
in an effort to better understand the molecular mechanisms underlying planarian regeneration, we performed a functional genomics screen aimed at identifying genes involved in this process in schmidtea mediterranea. neoblastspatterningimportinfunctional genomicsschmidtea mediterraneaissue-copyright-statement© the author 2015 we used microarrays to detect changes in gene expression in regenerating and non-regenerating tissues in planarians regenerating one side of the head and followed this with high-throughput screening by in situ hybridization and rnai to characterize the expression patterns and function of the differentially expressed genes. we also found that a homolog of the nuclear transport factor importin-α plays a role in stem cell function and tissue patterning, suggesting that controlled nuclear import of proteins is important for regeneration. we have additionally created an online database to house our in situ and rnai data to make it accessible to the planarian research community. planarians are renowned for their regenerative capacity and are an attractive model for the study of adult stem cells and tissue regeneration. along with five previously characterized genes, we identified additional genes necessary for stem cell maintenance or general regenerative capability or that caused tissue-specific defects upon knockdown. through this work, we described the roles of several previously uncharacterized genes in planarian regeneration and implicated nuclear import in this process.
 planarians are renowned for their regenerative capacity and are an attractive model for the study of adult stem cells and tissue regeneration. in an effort to better understand the molecular mechanisms underlying planarian regeneration, we performed a functional genomics screen aimed at identifying genes involved in this process in schmidtea mediterranea. methods we used microarrays to detect changes in gene expression in regenerating and non-regenerating tissues in planarians regenerating one side of the head and followed this with high-throughput screening by in situ hybridization and rnai to characterize the expression patterns and function of the differentially expressed genes.  along with five previously characterized genes, we identified additional genes necessary for stem cell maintenance or general regenerative capability or that caused tissue-specific defects upon knockdown. we also found that a homolog of the nuclear transport factor importin-α plays a role in stem cell function and tissue patterning, suggesting that controlled nuclear import of proteins is important for regeneration.  through this work, we described the roles of several previously uncharacterized genes in planarian regeneration and implicated nuclear import in this process. we have additionally created an online database to house our in situ and rnai data to make it accessible to the planarian research community. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords neoblastspatterningimportinfunctional genomicsschmidtea mediterraneaissue-copyright-statement© the author 2015 