The CCR4-NOT Complex Mediates Deadenylation and Degradation of Stem Cell mRNAs and Promotes Planarian Stem Cell Differentiation

Post-transcriptional regulatory mechanisms are of fundamental importance to form robust genetic networks, but their roles in stem cell pluripotency remain poorly understood. Here, we use freshwater planarians as a model system to investigate this and uncover a role for CCR4-NOT mediated deadenylation of mRNAs in stem cell differentiation. Planarian adult stem cells, the so-called neoblasts, drive the almost unlimited regenerative capabilities of planarians and allow their ongoing homeostatic tissue turnover. While many genes have been demonstrated to be required for these processes, currently almost no mechanistic insight is available into their regulation. We show that knockdown of planarian Not1, the CCR4-NOT deadenylating complex scaffolding subunit, abrogates regeneration and normal homeostasis. This abrogation is primarily due to severe impairment of their differentiation potential. We describe a stem cell specific increase in the mRNA levels of key neoblast genes after Smed-not1 knock down, consistent with a role of the CCR4-NOT complex in degradation of neoblast mRNAs upon the onset of differentiation. We also observe a stem cell specific increase in the frequency of longer poly(A) tails in these same mRNAs, showing that stem cells after Smed-not1 knock down fail to differentiate as they accumulate populations of transcripts with longer poly(A) tails. As other transcripts are unaffected our data hint at a targeted regulation of these key stem cell mRNAs by post-transcriptional regulators such as RNA-binding proteins or microRNAs. Together, our results show that the CCR4-NOT complex is crucial for stem cell differentiation and controls stem cell-specific degradation of mRNAs, thus providing clear mechanistic insight into this aspect of neoblast biology.


Introduction
Post-transcriptional control is central for the regulation of gene expression in stem cells [1]. A key post-transcriptional process is mRNA degradation [2][3][4] the control of which is believed to be as important as transcriptional regulation [5,6]. Although transcriptional regulation has been extensively studied, less is known about the developmental and physiological roles of mRNA degradation in stem cells, which are thought to involve the same RNA binding proteins [7] that act together to coordinate many complex aspects of mRNA biology, one of which is degradation. mRNA degradation starts with deadenylation (i.e. shortening of the poly(A) tail) [2,8]. This affects gene expression both by decreasing translational activity and mRNA stability [9,10]. The major deadenylase in eukaryotes is the CCR4-NOT complex [11][12][13], which is also involved in regulating several other aspects of mRNA metabolism, such as mRNA export [14,15], translation [15] and transcription itself [11,[16][17][18].
In yeast, the CCR4-NOT complex is composed of nine different subunits [11]: two deadenylases (Ccr4p and Pop2p/Caf1p), five Not proteins (Not1p-Not5p), Caf40p and Caf130p. Among them, Not1p, a 240 kDa protein, is thought to act as a scaffold and is the only subunit required for yeast viability [11,19]. Most of the subunits of the yeast complex are conserved across metazoans [20][21][22], with the exception of Not5p and Caf130p. In mammals two paralogous genes with mutually exclusive expression patterns encode each deadenylase of the complex [23]. Furthermore, the two deadenylase subunits Caf1 and Ccr4 regulate distinct sets of mRNAs [24,25].
Little is known, however, about the role of the CCR4-NOT complex in stem cells. It was found, for instance, that different components are important in maintaining mouse and human ESC identity [38], but the mechanisms remain largely unexplored. The freshwater planarian Schmidtea mediterranea is an emerging model for stem cell biology [39][40][41]. Its striking regeneration capacities are sustained by the presence of the neoblasts, a population of pluripotent stem cells that not only drive regeneration but sustain constant homeostatic cell turnover as well [42]. Planarian neoblasts can be eliminated by irradiation and are amenable to RNAi-mediated gene knock down. Furthermore, their abundance in the organism allows the quantitative evaluation of phenotypes.
Here, we use S. mediterranea to investigate the role of the CCR4-NOT complex in stem cell regulation by characterizing the function of the Smed-not1 gene, which encodes for the homologue of Not1p in yeast and CNot1 in mammals. We report that its knock down specifically affects the differentiation and selfmaintenance capabilities of neoblasts. Smed-not1 knock down results in a progressive increase in the levels of several neoblast transcripts, and we demonstrate that this increase is stem cell specific. Finally, we observe that these same mRNAs have stem cell specific increases in the frequency of long poly(A) tails after Smed-not1 knock down, showing that the observed increases in mRNA levels in stem cells are likely a consequence of decreased targeted deadenylation by the CCR4-NOT complex. Our findings highlight a likely central role for poly(A) tail length regulation in orchestrating pluripotent stem cell differentiation.

Results
In silico characterization of the CCR4-NOT complex in S. mediterranea We identified the different subunits of the CCR4-NOT complex in the planarian species S. mediterranea by TBLASTN searches in the S. mediterranea genome and in our reference transcriptome assembly [50,53] and other genomic and transcriptomic resources [58][59][60] (Table S1). We identified homologues of all metazoan genes known to encode subunits of the CCR4-NOT complex. Two components of the CCR4-NOT complex, the deadenylases Ccr4 and Caf1 were previously described in the planarian species Dugesia japonica [45]. However, no phenotype was reported for these enzymatic components after RNAi-mediated knock down and we also observed no strong phenotype ( Figure S1) for Smed-not6 (the Ccr4 orthologue, Table S1), Smed-not7A and Smed-not7B (the two paralogues of Caf1 in S. mediterranea, Table S1). We instead chose to focus on Smed-not1, as Not1 is believed to act as the central scaffolding protein in the complex and it is the only component of the complex essential for viability in yeast [11,19].

Smed-not1 is expressed in CNS and throughout the parenchyma in an irradiation-sensitive manner
We investigated the expression pattern of Smed-not1 by whole mount in situ hybridization (WMISH). We observed broad expression of this key component of the core deadenylation complex, including expression throughout the parenchyma and the central nervous system (CNS) ( Figure 1A, top panel). This pattern suggested to us that Smed-not1 may be expressed in neoblasts, since they are distributed throughout the parenchyma. To test this we monitored parenchymal expression after irradiation to remove neoblasts and observed that the parenchymal component of expression disappeared over a period of 5 days after irradiation ( Figure 1A, middle and bottom panels). All neoblasts disappear 24-48 hours after lethal irradiation, and consequently the expression of neoblast specific genes disappears over a similar period. As controls we analyzed the expression of Smedtud-1, the orthologue of the previously described Schmidtea polychroa Tudor gene Spoltud-1 [48] (Figure 1B), Smed-vasa-1, a Vasa orthologue of S. mediterranea [54,61] (Figure 1C), and Smed-pcna, the orthologue of the PCNA gene described in Dugesia japonica [62] ( Figure 1D). As expected, the neoblast-specific signals of all three disappeared by day 3 post-irradiation, while irradiation insensitive expression in differentiated cells of the CNS remained for Smedtud-1 and Smedvasa-1 ( Figure 1B-D). Smed-eye53, a marker expressed in differentiated cells [63], was used as a control to demonstrate that gene expression in post-mitotic cells is not ablated by irradiation ( Figure 1E). Since Smed-not1 hybridization signals present in the parenchyma were reduced but did not completely disappear by day 3 post-irradiation, it is likely that Smed-not1 is expressed in neoblasts and their recent progeny, but also other post-mitotic cells. These data indicates that Smed-not1 is broadly expressed throughout the planarian body, consistent with a housekeeping function, in a pattern that includes neoblasts, their progeny as well as differentiated cells. Investigation of Not1 expression in previous transcriptome based studies of mRNAs expressed in neoblasts is consistent with this expression pattern ( Figure S2).

Author Summary
Although transcriptional regulation in stem cells is a very active subject, much less is known about how posttranscriptional mechanisms of gene expression affect stem cells. Here, we use freshwater planarians in order to address this question. Planarians have a striking regenerative capacity driven by a population of pluripotent stem cells, the neoblasts. Control of both proliferation and differentiation is thought to rely heavily on post-transcriptional mechanisms, but their precise role is unknown. Poly(A) tail length regulation is an important mechanism of post-transcriptional control of gene expression as changes can be very rapid, and longer poly(A) tails are linked to increased mRNA stability and translational activity. We investigated the role of the CCR4-NOT complex, the major deadenylating complex in eukaryotes, by knocking down its main scaffolding subunit called Not1. Neoblasts in knock down animals are unable to differentiate and accumulate mRNAs with longer poly(A) tails. Our results show that the CCR4-NOT complex is needed for the targeted degradation of specific mRNAs expressed in stem cells, and the failure of this process likely prevents neoblasts from differentiating. These results reveal a new functional aspect of the CCR4-NOT complex and offer a mechanistic insight into the regulation of planarian stem cells.
Smed-not1 is required for planarian regeneration and homeostatic cell turnover We then analyzed the function of Smed-not1 by RNAi experiments. All Smed-not1(RNAi) animals displayed abrogated regeneration capacities and eventually died. They were able to produce both anterior and posterior regeneration blastemas, but never completed the regenerative process ( Figure 2B, vs. Figure 2A). In order to test if the formation of a regeneration blastema depended on the time of transection, we cut animals at 1, 3, 5, 10 and 15 days after Smed-not1(RNAi) treatment ( Figure S3). We found that animals were able to produce a regeneration blastema at all-time points, however animals cut earlier produced larger blastemas. All blastemas of Smed-not1(RNAi) worms eventually regressed. The ability of Smed-not1(RNAi) animals to produce a large regeneration blastema at early time points after RNAi suggests that mitotic neoblasts, the source of blastema cells, are still present and proliferating.
We then analyzed the phenotype of intact Smed-not1(RNAi) animals. We found that Smed-not1(RNAi) worms presented head regression ( Figure 2D, vs. Figure 2C) [43,44,64,65], a symptom of interrupted homeostatic cell turnover. A survival curve during which the onset of tissue homeostasis defects was recorded (N = 40) demonstrated temporal phenotypic variability ( Figure 2E). In Smed-not1(RNAi) animals these defects were seen first at 15 days after RNAi, and in the majority of animals after 20 days ( Figure 2C). Variable degrees of head regression were also observed after 20 days of RNAi ( Figure 2C, also see Figure S3). By day 22 all animals had defects, showing complete penetrance of Smed-not1 RNAi. All animals died by day 36 after dsRNA delivery, with the majority of deaths occurring between day 26 and day 34 ( Figure 2E). All control(RNAi) animals survived without any defect for .35 days, however. Together, these results demonstrate that Smed-not1 is needed for regeneration and homeostatic cell turnover in S. mediterranea.

Smed-not1(RNAi) animals maintain proliferative neoblasts
Next, we analyzed the mitotic marker phospho-histone-H3 (h3p) in Smed-not1(RNAi) animals at 5, 10, 15 and 20 days after RNAi. Up to 15 days all (N = 7 per time point) had mitotic neoblasts comparable in numbers to those of control(RNAi) animals ( Figure 3A). After 20 days all animals still had mitotic cells, but at variable levels, in agreement with our earlier phenotypic characterisation ( Figure 3B-E). Most had normal levels ( Figure 3C, vs. Figure 3B), although animals with more severe head regression defects showed a visible reduction in the mitotic levels ( Figure 3D-E, vs. Figure 3B-C), but overall the reduction in mitoses was not statistically significant. These experiments show that Smed-not1(RNAi) worms have mitotic cells, even as head regression defects progress. Significantly, similar defects are seen in irradiated worms only weeks after complete loss of mitotic activity. Therefore, we interpret our data as showing that effects on stem cell proliferation are not the primary cause of the Smed-not1(RNAi) phenotype, instead implicating neoblast differentiation impairment as responsible for regenerative failure, head regression and other defects. By performing FACS experiments we found that Smednot(RNAi) led to a moderate reduction of the sorted irradiation sensitive X1 cells, which primarily contains neoblasts, at 15 days but not at 10 days ( Figure 3F. Figure S4A-B). Even though no significant decrease of h3p cells was detected at this time point we interpret our FACS data as more sensitive and conclude that both methods consistently detect large numbers of proliferating neoblasts 15 days after Smed-not1 dsRNA administration, further implicating neoblast differentiation defects instead. In order to monitor the behaviour of neoblasts and their postmitotic progeny during progression of the Smed-not1 knock down phenotype we analysed the expression of neoblast and progeny markers [66]. Smedwi-1, a marker of neoblasts, Smed-nb.21.11e, a marker of early neoblast progeny, and Smed-agat-1, a marker of late neoblast progeny, were analyzed in control(RNAi) worms ( Figure 4A-C) and Smed-not1(RNAi) worms after 10 ( Figure 4D-F), 15 ( Figure 4G-I), and 20 ( Figure 4J-L, Figure S5A-F) days of RNAi. Only one time point (10 days) is shown for control(RNAi) worms, since no differences were observed between time points. Smedwi-1 expression was qualitatively the same after 10 and 15 days of Smed-not1(RNAi), ( Figure 4D, G, vs. Figure. 4A), but clearly reduced to a variable extent after 20 days ( Figure 4J vs. Figure 4A, D and G; Figure S5A vs. Figure S5D); some animals had nearly normal expression while Smedwi-1 expression was severely reduced in those with the most severe head regression defects. These results, like those above, suggest that prominent stem cell loss occurs only when Smed-not1(RNAi) animals begin to regress the head and to die, again implicating differentiation impairment instead of proliferation or self-renewal as a primary cause for regenerative failure.
Expression of Smed-nb.21.11e, a marker of early neoblast progeny, looked broadly equivalent to that of control worms in animals fixed 10 days after Smed-not1 dsRNA delivery ( Figure 4E, vs. Figure 4B) but a clearly visible progressive decrease in the number of Smed-nb.21.11e-positive cells was detected after 15 ( Figure 4H, vs. Figure 4B) or 20 ( Figure 4K, vs. Figure 4B) days after RNAi. After 20 days of Smed-not1 knock down animals had only a few remaining Smed-nb.21.11e-positive cells ( Figure S5E). These results show that clearly visible decreases in early progeny cell number precede the prominent decrease in neoblasts themselves. When we checked the expression of Smed-agat-1, a marker of late neoblast progeny, we observed a consistent qualitative increase in Smed-agat-1-positive cells in Smed-not1(RNAi) worms after 10 days ( Figure 4F vs. Figure 4C). At later time points, however, the number of Smed-agat-1-positive cells also progressively declined ( Figure 4I, L, vs. Figure 4C). Again, we observed a considerable variability in Smed-not1(RNAi) worms after 20 days ( Figure S5F). However, in all animals with a clear decrease in the number of Smed-agat-1-positive cells this was particularly apparent in the anterior region, a characteristic feature of Smed-agat-1positive cell depletion upon neoblast elimination by irradiation or perturbation by RNAi [50,64,65,67].
Therefore, neoblasts are abundant 15 days after RNAi, and only clearly start to be depleted later, coinciding with the onset of head regression defects. Similar defects are observed after irradiation, however, these take .10 days to manifest after elimination of mitotic activity. In contrast Smed-not1(RNAi) animals display these defects when neoblasts are still present. These results suggest that a primary defect in neoblast differentiation, rather than neoblast maintenance, causes failure in tissue homeostasis. In support of this idea, altered stem cell progeny numbers precede the disappearance of Smedwi-1 signals and mitotic activity. This alteration can be observed as early as 10 days after RNAi in the case of Smed-agat-1-positive cells and 15 days for Smed-nb.21.11epositive cells, which are clearly depleted at this time point.
To further test this, we compared the dynamics of the neoblast and progeny cell markers in Smed-not1(RNAi) animals to those of Smedwi-2(RNAi) animals, in which neoblast differentiation is abrogated [44]. Smedwi-1, Smed-nb.21.11e and Smed-agat-1 were expressed in Smedwi-2(RNAi) animals with very similar dynamics to Smed-not1(RNAi) animals ( Figure S5G-R). Taken together, these Smed-not1(RNAi) animals have abnormal numbers of Smed-agat-1 transcripts and Smed-agat-1-positive cells In order to achieve a quantitative measure of stem cell progeny mRNA levels in Smed-not1(RNAi) animals we performed quantitative real time PCR experiments (qRT-PCR) of Smed-nb.21.11e and Smed-agat-1 transcripts in Smed-not1(RNAi) worms. We focused on earlier time points of 10 and 15 days after RNAi. qRT-PCR showed that Smed-nb.21.11e levels in whole animals were similar to those of control(RNAi) after 10 and 15 days of dsRNA administration ( Figure 5A). However, Smed-agat-1 mRNA levels increased by day 10 and were almost two-fold higher after 15 days. Since this result did not correlate with what we observed by colorimetric WMISH, we quantified Smed-agat-1-positive cells by fluorescent WMISH (FWMISH). A significant increase in the numbers of Smed-agat-1-positive cells was found in Smed-not1(RNAi) animals 10 days after dsRNA administration ( Figure 5B), but their numbers declined to control levels after 15 days. These data confirmed the qualitative data from colorimetric WMISH ( Figure 4). The distribution of Smed-agat-1-positive cells in Smed-not1(RNAi) animals 15 days after dsRNA delivery was different from controls, with less Smed-agat-1-positive cells in the anterior part of the worm ( Figure 5C). This revealed a stark discordance between the number of Smed-agat-1-positive cells and the level of Smed-agat-1 mRNA. We conclude that Smed-agat-1 transcripts are accumulating in decreasing numbers of Smed-agat-1-positive cells after 15 days of Smed-not1 knock down, and that each Smed-agat-1-positive cell contains an increased average number of Smed-agat-1 transcripts. A similar, but less pronounced, process could also explain reduced numbers of Smed-nb.21.11e positive cells and discordant stable levels of this transcript, which do not drop by day 15.
Given that the CCR4-NOT complex is known to regulate gene expression through its deadenylating activity we wished to ascertain if it could be directly responsible of the discordance between Smed-agat-1 positive cell number and mRNA levels. If an increase in mRNA levels is caused by impaired deadenylation and subsequent degradation we would expect to observe increased frequency of long poly(A) tail lengths. Using a poly(A) tail length (PAT) assay [68][69][70][71], in whole worms we observed that this was the case for both Smed-agat-1 and Smed-nb.21.11.e ( Figure 5D), with the first giving starker differences, while control mRNAs Smed-eif-3, Smed-mhc, or Smed-ef-2 were only mildly affected or not affected at all. In addition a spike-in control of exogenous mRNA showed equal poly(A) tail length distribution across samples ( Figure 5D), showing that the differences observed are present in our different planarian mRNA samples and not introduced by the PAT assay technique.

Smed-not1(RNAi) animals have increased levels of transcripts expressed in stem cells with increased frequency of long poly(A) tails
Next, we found that WMISH analysis of neoblast markers, Smedtud-1, Smed-vasa-1 and Smed-pcna, was suggestive of increased levels of these transcripts in Smed-not1(RNAi) worms, with qualitatively visible differences after both 10 and 15 days after RNAi treatment ( Figure 6A-C). These data further demonstrate that neoblast maintenance is not affected by Smed-not1 knock down at these time points.
Since WMISH does not provide a quantitative measure of mRNA levels we quantified these differences by qRT-PCR experiments on RNA from Smed-not1(RNAi) animals. Smedtud-1, Smed-vasa-1 and Smed-pcna all progressively increased to levels approximately 50% and 100% higher than those of control(RNAi) after 10 and 15 days, respectively ( Figure 6D). Smedwi-1 transcript levels were also significantly increased in whole animals after 15 days ( Figure 6D). Collectively, these results demonstrate that Smed-not1 RNAi knock down leads to an increase in mRNA levels in several genes expressed in neoblasts and their progeny.
Given the known conserved function of the CCR4-NOT complex in regulating mRNA levels through targeted deadenylation we performed PAT assays on the set of neoblast markers and on the samples above ( Figure 5D). In all cases Smed-not1 knock down resulted in increased average poly(A) tail length ( Figure 6E), demonstrating that increased transcript levels correlate with a failure in deadenylation. These data confirm that knock down of the CCR4-NOT complex subunit Smed-not1 leads to increased transcript levels of genes known to be key to Upregulation and increased polyadenylation of Smedtud-1, Smed-vasa-1 and Smed-pcna occurs specifically in stem cells Since many neoblast mRNAs are also expressed in differentiated cells (e. g. Smedtud-1 and Smed-vasa-1 are also prominently expressed in the CNS, Figure 1B-C), the increased mRNA levels we detected could arise from a response in differentiated cells alone, stem cells alone or from both differentiated and stem cells. To distinguish these possibilities, we planned to use an irradiation approach to remove all neoblasts and then measure transcript levels in Smed-not1(RNAi) and control(RNAi) worms by qRT-PCR. We reasoned that if transcript accumulation was indeed limited to stem cells then mRNA levels of these transcripts after irradiation would be equal in both irradiated control(RNAi) and irradiated Smed-not1(RNAi) samples.
Smedtud-1 and Smed-vasa-1 are expressed in the CNS to levels that amount respectively to roughly 70% and 40% of their total expression, according to our previous neoblast profiling by a combinatorial RNA-seq, RNA interference and irradiation approach [50]. We confirmed this by qRT-PCRs ( Figure S5A) and observed that Smed-pcna expression amounts to roughly only 10% of its normal expression 24 hours after irradiation (Figure S6A), indicating that 24 hours of irradiation suffice to eliminate around 90% of neoblasts. Conversely, most neoblast progeny survive beyond 1 day post-irradiation, as measured by our qRT-PCR experiments with the markers Smed-nb.21.11e and Smed-agat1 ( Figure S6B) and consistent with previously published data [50,66], making then 24 hours after irradiation the ideal time point to perform our experiment.
Therefore, in order to find out if the increased levels of neoblast mRNAs originate in neoblasts or instead in the CNS or elsewhere, we then used this irradiation approach in Smednot-1(RNAi) animals and compared them to control(RNAi) animals. We irradiated both Smednot-1(RNAi) and control(RNAi) animals at 9 and 14 days after dsRNA administration, 24 hours before the data collection time points of 10 and 15 days respectively. We then examined the expression pattern of Smedtud-1 by WMISH. This confirmed our qRT-PCR experiments, and was consistent with our previous WMISH results ( Figure 7A). Non irradiated Smed-not1(RNAi) animals showed qualitatively more intense expression of Smedtud-1, However, all animals that were irradiated 24 hours prior to fixation showed an identical expression pattern of Smedtud-1, with similar levels of signal detected only in the CNS to their control(RNAi) irradiated counterparts ( Figure 7A). This experiment shows that Smedtud-1 is not ectopically expressed in other tissues or organs after Smed-not1 knock down, since this ectopic expression should be visible either in non irradiated or irradiated samples, and suggests instead that the increased levels of Smedtud-1 come from an accumulation of this transcript in neoblasts.
We performed qRT-PCR and confirmed that Smedtud-1, Smedvasa-1 and Smed-pcna mRNA levels increased progressively in Smed-not1(RNAi) animals ( Figure 7B, left) but both control(RNAi) and Smed-not1(RNAi) animals irradiated 24 hours previously contained similar levels of all three transcripts ( Figure 7B, right). Similar to results for wild type worms ( Figure S6A), the levels of the Smedtud-1, Smed-vasa-1 and Smed-pcna transcripts were respectively around 70%, 40% and 10% of the expression in non-irradiated control(RNAi) animals. Therefore we conclude that the overexpression observed in Smed-not1(RNAi) animals for these three transcripts disappears 24 hours after irradiation, and is therefore located in irradiation-sensitive neoblasts, rather than elsewhere in the body.
To implicate the CCR4-NOT complex-mediated targeted mRNA degradation directly in this effect we performed PAT assays on non irradiated and irradiated samples and observed that the distribution of poly(A) tail lengths of all these neoblast mRNAs was increased in Smed-not1(RNAi) animals, but this effect was removed by irradiation ( Figure 7C and D). In the case of neoblast specific Smedwi-1 and Smed-pcna no signal was detected after irradiation ( Figure 7C). In the case of Smed-vasa-1 and Smedtud-1 signal was detectable after irradiation, reflecting expression in the CNS, but these transcripts did not show any increase in poly(A) tail length distribution after Smed-not1 knock down ( Figure 7D). These data confirm that effects on mRNA levels for these important neoblast genes are confined to stem cells. As an additional control we also measured Smed-nb.21.11e and Smed-agat-1 poly(A) lengths ( Figure 7E). 24 hours after irradiation cells expressing these transcripts are still present and the increase in poly(A) tail length caused by Smed-not1 knock down is still evident. This confirms that irradiation itself does not cause the absence of detected differences in poly(A) tail length per se, but by removing the cycling neoblasts. The poly(A) tail length distribution of control mRNAs were not affected by irradiation ( Figure 7F).
As a further proof that the effects we observe are confined to stem cells we looked at gene expression, poly(A) tail length and the efficacy of Smed-not1 knock down across FACS cellular compart-ments [51,72]. Progressive increases in neoblast gene mRNA levels were confined to X1 and X2 populations of sorted cells from Smed-not1(RNAi) animals ( Figure 8A, X1 and X2). Both of these sorted fractions contain neoblasts to different extents. Levels of these mRNAs were not increased in irradiation insensitive (Xins) sorted fractions from Smed-not1(RNAi) animals compared to controls ( Figure 8A, Xins). This fraction contains primarily differentiated cells including CNS cells which also express Smedtud-1 and Smedvasa-1.
One possibility for the specificity we observe is that Smed-not1 knock down has a high efficacy in stem cells but not in post-mitotic cells, as it has been suggested recently for Smed-bruno-like knock down [38,73]. We performed qRT-PCR measurement of Smed-not1 transcript levels and found that Smed-not1 mRNA is consistently expressed in all cellular fractions (consistent with Figure S2 and Figure 1A) and that Smed-not1 knock down significantly depletes Smed-not1 in all sorted compartments compared to control(RNAi) animals ( Figure 8B). Consequently the observed neoblast specificity is unlikely due to an absence of knock down in differentiated cells.
To finally link the mechanism of CCR4-NOT complexmediated deadenylation to increased mRNA levels we also checked poly(A) tail lengths in FACS sorted cells. For Smedtud-1, Smed-vasa-1, Smed-pcna and Smedwi-1 we observed a progressive increase in long poly(A) tails in the X1 compartment of stem cells ( Figure 8C and D). This trend was also observed for these transcripts in X2 cells ( Figure 8C and D). No poly(A) tail signal was detectable in Xins fractions for Smed-pcna and Smedwi-1, consistent with their low abundance in this fraction, and poly(A) tail length was not affected for Smedtud-1 and Smed-vasa-1 in Xins cells ( Figure 8C and D). These data are in agreement with our irradiation based experiments (Figure 7). The poly(A) tail lengths of control mRNAs Smed-ef-2 and Smed-mhc were only mildly or not affected by Smed-not1 knock down ( Figure 8E). These mild differences are likely due to CCR4-NOT mediated deadenylation, but again seem to be restricted to X1 and X2 cells, whereas Xins cells remain unaffected. To which extent all, most or only a subset of neoblast transcripts are affected after Smed-not1 knock down remains an open question.
Taken together, these results clearly demonstrate that Smed-not1 knock down induces a prominent increase of key transcripts expressed in neoblasts, that this accumulation occurs specifically in neoblasts, and that it is associated with an increased frequency of long poly(A) tails of these transcripts specifically in neoblasts. Furthermore, these results strongly suggest that the neoblastspecific increase of mRNA levels of genes such as Smedtud-1, Smedvasa-1 and Smed-pcna may be responsible for the impaired differentiation capacities of neoblasts observed in Smed-not1(RNAi) animals. It is likely that other genes expressed in neoblasts are similarly upregulated in Smed-not1(RNAi) animals and contribute to differentiation defects. Our results suggest a mechanism by which the differentiation capacities of neoblasts are dependent on CCR4-NOT mediated degradation of specific neoblast mRNAs.

Discussion
Planarians are an emerging in vivo model for stem cell biology because of their unique stem cell population. In this study we used the planarian S. mediterranea as a model system to establish a function for the CCR4-NOT complex in stem cell regulation. We found that Smed-not1 knock down abrogated regeneration and impaired homeostatic cell turnover. Interestingly, Smed-not1 knock down primarily affects the stem cell compartment of S. mediterranea rather than inducing more widespread effects, even though the CCR4-NOT complex is the major deadenylating complex in eukaryotes and regulates at least 85% of mRNAs in yeast [74].
While we observed a stark and specific effect of Smed-not1 knock down on deadenylation, it is still possible that other functions of the CCR4-NOT complex might also be impaired. The CCR4-NOT complex is involved in several steps of RNA metabolism [12,13,16] and further work is therefore needed to elucidate which ones are also at work in stem cells.
We observed effective gene knock down of Smed-not1 even in differentiated cell fractions, but specific effects on neoblast transcripts were limited to stem cells. This fact suggests that targeted deadenylation by either RNA-binding proteins or miRNAs is providing specificity and is therefore central to stem cell differentiation and self-renewal properties. Consistently, several studies have implicated the CCR4-NOT complex in mRNA-specific deadenylation via targeted recruitment of the CCR4-NOT complex by RNA-binding proteins [26,27,30,31], which are in turn known to be highly enriched and functionally important in neoblasts [43][44][45][46][48][49][50]52,75]. It is possible to hypothesize that disruption of the CCR4-NOT complex via knock down of its scaffolding subunit might impair the protein-protein interactions needed to tether deadenylating activity to transcripts targeted by RNA-binding proteins or miRNAs for degradation, while general non-targeted deadenylation would remain relatively unaffected, therefore not causing a broader metabolic failure at earlier time points. In fact, the miRNA silencing complex protein GW182 specifically interacts with the Not1 subunit of the CCR4-NOT complex through specific and conserved domains [36,37].
Despite the specificity of the effects seen in neoblast transcripts, mRNAs expressed elsewhere were also found to be affected. Both Smed-agat-1 and, to a lesser extent, Smed-nb.21.11e were also affected. These results show that CCR4-NOT deadenylating activity is present in cell types other than neoblasts and that specificity is not due to restriction of activity to stem cells. Furthermore, the effects on well described neoblast progeny markers suggest that Smed-not1 knock down likely influences several steps of cellular differentiation that may all contribute to the observed effects on homeostasis and regeneration.
The Smed-not1 phenotype is progressive with respect to both the decreasing capacity of the animals to produce blastema cells and by the accumulation of mRNAs in stem cells. The phenotype results from a drop in neoblast progeny numbers, followed by stem cell loss. Similar neoblast and progeny dynamics have been shown by us in Smedwi-2(RNAi) organisms here and by another group in Smed-CHD4(RNAi) organisms [64]. Both Smedwi-2 and Smed-CHD4 knock downs initially affect neoblast differentiation rather than their self-renewal and proliferative capabilities [44,64] which are only affected at later time points, similar to Smed-not1 knock down. The ultimate cause for stem cell loss after the impairment of neoblast differentiation is unknown, and likely to be a broad failure in homeostasis as organs and tissues fail. However, for more than 15 days after Smed-not1 dsRNA administration proliferating neoblasts are detectable in large numbers while regeneration is abrogated, suggesting neoblast loss is not a primary cause for the regeneration defect, instead implicating impaired neoblast differentiation capacities. Smed-not1 knock down induces an increased frequency of long poly(A) tails of Smedtud-1, Smed-vasa-1, Smed-pcna and Smedwi-1 mRNAs as these same mRNAs accumulate. This regulation occurs specifically in neoblasts, rather than in the CNS, where two of these transcripts are also expressed. However, it is likely that Smed-not1 knock down affects many more transcripts that contribute to failure in stem cell differentiation. Due to the lack of specific antibodies it is difficult to evaluate if the differences observed in transcript abundance and polyadenylation state affect the abundance of the proteins that these transcripts encode. However, taking into account that all four transcripts expressed in neoblasts analyzed accumulate in these stem cells, it is reasonable to conclude that Smed-not1 knock down induces transcriptome-wide changes in stem cell expression patterns, and that these changes will likely affect protein levels.
Our results offer a new mechanistic insight into post-transcriptional regulation in neoblasts and its targets. After depletion of a post-transcriptional regulator many transcripts accumulate without being degraded, and this likely prevents neoblast differentiation, which needs the effective removal of these transcripts. The deadenylating activity of the CCR4-NOT complex is clearly central to this process. For genes like Smedtud-1 and Smed-vasa-1, which are expressed in neoblasts and the CNS, we observe accumulation and increased frequency of long poly(A) tails of the transcript only in neoblasts. This suggests that the deadenylation of these mRNAs is regulated specifically during the onset of differentiation and requires the targeted recruitment of the CCR4-NOT complex by RNA-binding proteins, as has been described in other organisms [26,27,30,31,76]. Interestingly, several RNA binding proteins and post-transcriptional regulators have already been described as crucial for neoblast function [43,44,46,48], and some of them have been already been functionally linked with the CCR4-NOT complex in other model organisms. Future research will help elucidate the mechanisms by which these proteins orchestrate planarian stem cell processes.
The CCR4-NOT complex has been shown to mediate the repressive function of both miRNAs and piRNAs [32][33][34][35]. Small RNAs are believed to be very important regulators of mammalian stem cells [77] and neoblasts [78,79]. The Smed-not1 RNAi phenotype is very similar to those of Smedwi-2 and Smedwi-3, two Piwi proteins involved in piRNA regulation [44,80]. Furthermore, several studies have highlighted the presence of miRNAs highly enriched in stem cells [78,79,81]. Future research will help in understanding if these phenotypic similarities reflect a functional link between Piwi proteins, piRNAs, miRNAs and the CCR4-NOT complex in planarian stem cells. Our results highlight the importance of the CCR4-NOT complex in the regulation of stem cells, the fact that post-transcriptional regulation of gene expression is a key element in the regulation of pluripotency, and that planarians will provide an excellent platform for these studies.

Organisms
Planarians of the asexual strain of S. mediterranea were kept and used as previously described [82].

Sequences
The putative members of the S. mediterranea CCR4-NOT complex and other transcripts were identified in published S. mediterranea transcriptomic and genomic sequences [53,[58][59][60] and the longest transcripts for Smed-not-1 confirmed by PCR and sequencing. The putative members of the S. mediterranea CCR4-NOT complex were identified by TBLASTN searches in the current assembly of the S. mediterranea genome and in the available transcriptomic data. In order to determine the number of loci for each of the components, the different transcripts identified were mapped to the S. mediterranea genome. The genomic sequence encoding Smed-not1 was found split in two contigs (v31.001778 and v31.002774), the existence of one single transcript for these two genomic contigs was confirmed by PCR using the primers 59-CATCGCAACAATGGAGAGAA-39 and 59-ATTTGAGCTG-TATGGGCGAT-39, each mapping to one of the two contigs. These PCR experiments revealed as well the existence of a 3 kb region not present in the S. mediterranea genomic data.
The full sequence of Smed-not1 was obtained by de novo assembling the transcript from the raw transcriptomic data, using the known transcriptomic and genomic data as a guide. The Smed-not1 sequence has been deposited in Genbank (accession KF781122).

RNAi
RNAi experiments were performed as previously described [82]. Control(RNAi) worms were injected with dsRNA encoding for GFP, a gene not present in the S. mediterranea genome. dsRNA encoding for Smed-not1 was prepared by in vitro transcription of a region of the Smed-not1 gene. Briefly, an amplicon was generated from S. mediterranea reverse transcribed RNA with the primers 59-GGCCGCGGTGTCCAAGAAAAAGCAAGTCAG-39 and 59-GCCCCGGCCAGCTGGCGTCAGTTTAGTGAA-39, containing a 59 adaptor sequence for the subsequent addition of T7 promoter. The product of this PCR was gel purified and subjected to another step of amplification with the primers 59-GAGAATTC-TAATACGACTCACTATAGGGCCGCGG-39 and 59-AGG-GATCCTAATACGACTCACTATAGGCCCCGGC-39 with the purpose of attaching T7 promoter sequences to both ends of the amplicon. The product of this PCR was further purified with Purelink DNA purification columns (Invitrogen) and used as a template for in vitro transcription using T7 RNA polymerase (Roche). The product of the in vitro transcription was treated with Turbo DNAse (Ambion), phenol extracted, precipitated in ethanol in presence of sodium acetate, glycogen and EDTA and resuspended in water. dsRNA encoding GFP for use as a negative control was similarly prepared from a vector encoding the GFP gene. The final concentration of the injected solution was 1 mg/ ml. Animals were injected with a Nanoject II (Drummond) for three consecutive days and monitored or used for experiments in the subsequent days. Day 1 after RNAi is considered to be in all experiments the first day after the third dsRNA injection. Alive animals were imaged in a Zeiss Discovery V8 with a Zeiss AxioCam MRC camera.

Irradiation
Irradiation was performed as previously described [79]. Animals were placed in a sealed c-ray source and administered an irradiation dose of 100 Gy.
In situ hybridization, immunohistochemistry and imaging WMISH, FWMISH and IHC were performed and imaged as previously described [82]. When qualitative differences are shown, animals were processed and monitored in parallel. Riboprobes were generated by in vitro transcription of PCR products generated as described above, with only one T7 promoter linked to the 39 end of the amplicon, and in the presence of digoxigeninlabelled UTP (Roche). The products of in vitro transcription reactions were then treated with Turbo DNAse (Ambion). Riboprobes were then precipitated in ethanol in the presence of LiCl and glycogen and resuspended in 50% formamide in TE buffer, 0.01% Tween.
The following primers were used: Briefly, animals were killed in a 2%HCl solution, fixed in Carnoy's solution, bleached in a 8% H 2 O 2 /methanol, rehydrated, permeabilized with Proteinase K (Sigma), treated with 0.25% and 0.5% acetic anhydride in 0.1M triethanolamine pH 7.6, prehybridized and hybridized with digoxigenin labelled riboprobes (0.2 ng/ml, O/N at 56uC) They were then washed in buffers of increasing stringency, immunolabelled with anti-digoxigeninalkaline-phosphatase antibody (Roche) and developed in the presence of NBT and BCIP (Roche). For FWMISH, an antidigoxigenin-peroxidase antibody (Roche) was used and the signal was developed with the Tyramide Signal Amplification kit (Perkin Elmer). WMISH specimens were imaged on a Zeiss Discovery V8 equipped with a Zeiss AxioCam MRC camera. FWMISH specimens were imaged in a Leica SP3 confocal.
Quantification of Smed-agat-1-positive cells Animals (N = 8 per time point and treatment) were stained by FWMISH of Smed-agat-1 and monitored by confocal microscopy. Cell counts were performed in z-projections of both the dorsal and ventral sides of the animals. 130 squares corresponding to 250 mm 2 in both the dorsal and ventral parts of the animals were selected for counting along the length of the animals. Counts were performed using ImageJ software. Error bars represent standard error of the mean. Statistical significance was analyzed by Student's T test by comparing values from each sample to its respective control sample.

Immunohistochemistry
Whole mount immunohistochemistry of phospho-histone-3 was performed as previously described [83]. Essentially, animals were fixed as above, blocked in a 1% BSA/PBS 0.3% triton X-100 solution, incubated overnight with anti-phospho-histone-3 (Millipore, 1/500 dilution) and anti-phospho-tyrosine (Cell Signaling, 1/200 dilution), washed and incubated in Alexa Fluor 488 and 568 secondary antibodies (Molecular Probes, 1/400 dilution). Animals were then washed, mounted in 70% glycerol/PBS and imaged with a Leica SP3 confocal and a Leica MZ16F fluorescence stereomicroscope and a Leica DFC 300Fx camera. Countings of phospho-histone-3 were performed with ImageJ software and normalized to the total area of the sample.

qRT-PCR
qRT-PCR experiments were performed as previously described [48] with modifications. Essentially, total RNA from samples of 5 animals was extracted with Trizol reagent (Invitrogen) according to the manufacturer's instructions, and cDNAs were synthesized with SuperScriptIII Reverse Transcriptase (Invitrogen). qRT-PCR experiments were then performed using the Absolute qPCR SYBR Green Master Mix (Thermo Scientific). Experiments were performed on three biological replicates per time point and treatment. Each biological replicate was technically replicated three times in each reaction, each reaction was repeated three times. The gene Smed-ef-2 was used for normalization.
The following primers were used:

FACS
Fluorescence activated cell sorting (FACS) of planarian samples was performed as previously described [51,72]. Analysis was performed with FlowJo.

PAT assays
PAT assays were performed as previously described [69,71] with minor modifications. A total of 400-1000 ng of total RNA extracted from five animals were used per each time point, replicate and treatment, except for FACS samples, in which 40 animals were used per time point and treatment and PAT reactions were performed with 100 ng of total RNA. Three biological replicates were analyzed per time point and treatment, and technically replicated at least twice, except for FACS samples, which were only technically replicated. C. elegans total RNA was spiked-in as a control and assayed with a primer for cpg-2.
RNA samples were incubated with 0.3 mg of 59-phosphorylated oligo d(T) in a total volume of 8 ml, and heat denatured for 5 min at 65uC. Then, the following mixture was added: 4 ml of Super Script II First Strand Buffer, 0.5 ml of 0.1 M DTT, 2.25 ml of 10 mM ATP, 0.125 ml of RNAsin (Promega), 1.25 ml of 1 mM (each) dNTPs (Promega), and 1 ml of 2000 units/ml T4 DNA Ligase (New England Biolabs). The volume was then brought to a total of 20 ml with water, and the samples were incubated at 42uC for 30 min. in order to allow the saturation of polyA tails with oligo d(T). Then, 0.5 ml of an oligod(T)-anchor 100 mM was added , and the mixture was incubated for 2 hours at 12uC to allow ligation of oligo d(T) molecules. Then, samples were pre-warmed at 42uC and 1 ml of SuperScript II (Invitrogen) was added. Finally the samples were incubated for 1 hour at 42uC and the SuperScript enzyme was heat inactivated at 70uC for 30 min.
These cDNAs were then used in PCR reactions, with fresh aliquots of the anchor primer and gene specific forward primers designed close to the 39 region of the mRNA. The sequences of the mRNAs tested were obtained in our reference transcriptome (4). The following primers were used: Smedtud  21.11e and Smed-agat-1 (B) in animals 1, 3 and 5 days after irradiation, normalized expression and relative to non irradiated samples. Error bars represent standard deviation. Animals 1 day after irradiation have around 10% of Smed-pcna transcripts of non-irradiated controls, and this number further decreases 3 and 5 days after irradiation. However, the expression of Smedtud-1 and Smedvas-1 mRNAs only decreases to around 70% and 40% respectively of the level of non irradiated controls, reflecting expression that does not localize to neoblasts and is therefore not eliminated by irradiation. Similar to Smed-pcna, the level of Smedwi-1 transcripts decreases to around 10% of the expression in non irradiated controls and becomes almost undetectable later. The levels of the progeny specific mRNAs Smed-nb.21.11e and Smed-agat-1 decrease progressively at later time points of irradiation. Therefore, around 90% of the neoblast specific transcripts are eliminated only 1 day after irradiation while most of the expression of progeny specific transcripts is still detected and the non-neoblast expression of Smedtud-1 and Smedvas-1 localized in the CNS is not eliminated by irradiation. (TIF)

Table S1
In silico search of CCR4-NOT complex components in S. mediterranea. Summary of the CCR4-NOT complex components found in silico in S. mediterranea. Each of the described components of the yeast, Drosophila melanogaster and human CCR4-NOT complexes is indicated. The column ''S. mediterranea'' indicates the given name for each of the genes, the column ''Contig Smed genome'' indicates the genomic contigs in which each locus was found, followed by the transcriptomic datasets from Blythe et al. and Adamidi et al. The column ''D. japonica'' indicates the accession numbers of the CCR4-NOT components previously described in this planarian species. Nine different components of the CCR4-NOT complex were found both in genomic and transcriptomic sequences, corresponding to the orthologues of all metazoan CCR4-NOT complex components. Similar to humans, two paralogues of the not7/caf1 gene were found (Smed-not7A and Smed-not7B) in both genomic and transcriptomic sequences and two additional genomic loci encoding two similar versions of an additional not7/caf1 gene were found (Smed-not7C.1 and Smed-not7C.2). However, the transcripts encoded by these two genomic loci were not found in transcriptomic datasets, and therefore they are possible pseudogenes. No orthologue of the yeast specific not5 was found, but one orthologue of the metazoan specific not10 was found (Smed-not10). We found several transcripts mapping to the same genomic locus for most of the genes, encoding different regions of the gene or different splicing variants. The Smed-not1 gene was split in two different contigs (v31.001778 and v31.002774) encoding respectively the 59 and 39 regions of the same gene. PCR experiments confirmed that they correspond to the same transcript. (PDF)