The primary cilium dampens proliferative signaling and represses a G2/M transcriptional network in quiescent myoblasts

Background Reversible cell cycle arrest (quiescence/G0) is characteristic of adult stem cells and is actively controlled at multiple levels. Quiescent cells also extend a primary cilium, which functions as a signaling hub. Primary cilia have been shown to be important in multiple developmental processes, and are implicated in numerous developmental disorders. Although the association of the cilium with G0 is established, the role of the cilium in the control of the quiescence program is still poorly understood. Results Primary cilia are dynamically regulated across different states of cell cycle exit in skeletal muscle myoblasts: quiescent myoblasts elaborate a primary cilium in vivo and in vitro, but terminally differentiated myofibers do not. Myoblasts where ciliogenesis is ablated using RNAi against a key ciliary assembly protein (IFT88) can exit the cell cycle but display an altered quiescence program and impaired self-renewal. Specifically, the G0 transcriptome in IFT88 knockdown cells is aberrantly enriched for G2/M regulators, suggesting a focused repression of this network by the cilium. Cilium-ablated cells also exhibit features of activation including enhanced activity of Wnt and mitogen signaling and elevated protein synthesis via inactivation of the translational repressor 4E-BP1. Conclusions Taken together, our results show that the primary cilium integrates and dampens proliferative signaling, represses translation and G2/M genes, and is integral to the establishment of the quiescence program.


Background
Primary cilia are microtubule-based membrane-encased organelles that function to receive and transduce signals from the extracellular milieu [1].These cellular antennae are intimately linked to the cell cycle, and barring rare instances [2,3], the extension of the cilium is restricted to cells in G0/G1, and suppressed in other cell cycle phases [4].The primary cilium is formed from and anchored by a basal body, which is essentially a sequestered centrosome that has been modified for this role.Upon cell cycle re-entry, cilia are actively dismantled prior to mitosis, releasing the basal body to function as a microtubule-organizing centre (MTOC) for spindle assembly.Multiple proteins that play important ciliary roles are also involved in cell cycle regulation: for example, Aurora Kinase A (AurA), Polo Like Kinase-1 (Plk1), and Anaphase Promoting Complex (APC) are all key regulators of centrosome function and ciliary disassembly, as well as mitotic progression (reviewed in [5]).
Thus, the current understanding is that the cilium exerts a check on cell cycle progression (reviewed in [6]) and that disassembly of the cilium is a prerequisite for cell cycle re-entry [7].The cilium is vital in embryonic patterning [8], and its role in terminally differentiated adult cells such as kidney epithelia is critical for their function [9].Mutations in ciliary genes give rise to a spectrum of complex diseases collectively termed ciliopathies [10], usually associated with increased cellular proliferation such as in polycystic kidneys, thus supporting the idea that cilia exert an inhibitory influence on proliferation.Additionally, while analysis thus far shows that ciliopathy patients do not display an increased risk of cancer, there is emerging evidence linking cilia with tumor suppressive mechanisms (Reviewed in [11]).However, the molecular mechanisms that couple ciliogenesis and cell cycle progression are still emerging and the functions of the cilium within the quiescent cell are largely unexplored.
Quiescence (G0) is a state of reversible cell cycle arrest in which cells have exited the cell cycle, but retain the ability to reenter in response to favourable signals and conditions.Previous studies [12][13][14][15] show that G0 can be distinguished from an extended G1 not only by the altered kinetics of entry into S phase, but by major alterations in the transcriptome.Indeed, analysis of G0 in multiple cell types points to the induction of a common cellular regenerative program in quiescence, suggesting a molecular basis for increased self-renewal [16].The quiescent state is critical for the function of adult stem cells, which contribute to regeneration and tissue homeostasis (reviewed in [17]).In skeletal muscle, homeostatic maintenance of differentiated tissue, as well as its regeneration following injury, are made possible by a small population of stem cells called muscle satellite cells (MuSC) [18].These muscle precursors persist in a quiescent undifferentiated state in adult muscle and are activated to re-enter the cell cycle when the tissue experiences damage.The majority of activated MuSC enter a program of proliferation and differentiation to regenerate damaged myofibers, while a minor subset enters a distinct program to self-renew and regenerate myofiberassociated quiescent MuSC, thus leading to functional tissue recovery [19].Although a number of signaling, metabolic, transcriptional and epigenetic regulatory circuits have been implicated in MuSC quiescence [15], the mechanisms that govern its establishment during early postnatal development and its maintenance during adult life are still being uncovered.Recently, quiescent MuSC were shown to possess primary cilia, which were implicated in the regulation of asymmetric cell division during cell cycle activation [20].However, the genetic and signaling networks potentially controlled by the cilium within quiescent MuSC are not known.
Here, we investigate the role of the primary cilium in quiescence using an established culture system to model alternate MuSC fates: asynchronously proliferating C2C12 myoblasts can be driven to enter either reversible [14,21] or terminal [22] cell cycle exit, by modulating culture conditions.We show that not only are primary cilia present on quiescent MuSC in vivo, and elaborated during quiescence in vitro, but are also transiently induced during early myogenic differentiation.When ciliogenesis is ablated using RNAi against IFT88, a key intraflagellar transport protein involved in ciliary extension, an altered quiescence-reactivation program is observed.Additionally, IFT88 knockdown in G0 leads to specific induction of G2/M transcriptional networks, and cells lacking cilia display enhanced activity of multiple cilium-associated signaling pathways, accompanied by an increase in protein synthesis, with a specific shift towards 4E-BP1-regulated translation.Intriguingly, knockdown cells do not show increased DNA synthesis or mitosis.These observations suggest that in quiescent cells, loss of this cellular antenna relieves the repression on mitogenic signaling but additional events are required to trigger the exit from quiescence.Our findings support the existence of complex interplay between the cilium cycle and the cell cycle at both transcriptional and translational levels.Overall, our study reveals that specific repression of G2/M genes by the primary cilium is integral to the establishment of the self-renewing quiescent state in culture, with implications for stem cell maintenance and function.

Primary cilia are dynamically regulated during the cell cycle in myogenic cells
Earlier, we showed that quiescent myoblasts in culture display a distinct transcriptional profile, with altered signaling modules, in particular, enhanced Wnt-TCF signaling [12,23].Interestingly, a number of ciliumassociated genes showed increased expression in G0.We compiled a putative "ciliome" consisting of 1896 ciliumrelated genes curated from available datasets from different cell types [24,25].Comparison of the "ciliome" with the "G0 altered" transcriptome consisting of 1747 annotated Differentially Expressed Genes derived from a microarray analysis of C2C12 myoblasts in proliferation vs. quiescence [12], yielded substantial overlap (345 genes [19%], p-value: 1.717854e-19) (Fig. S1A).Selected transcripts encoding cilium-associated proteins were validated as up-regulated in quiescent myoblasts (Fig. S1B).Together, these findings suggest a previously unreported induction of a ciliary transcriptional program in skeletal muscle myoblast quiescence.
To investigate whether the transcriptional induction of ciliary genes is accompanied by the elaboration of primary cilia in G0 myoblasts, we used immunofluorescence analysis (Fig. S2).C2C12 myoblasts grow asynchronously in adherent culture and can be triggered to enter reversible arrest (quiescence/G0) by culture in suspension, where deprivation of adhesion leads to cessation of cell division [21,26].G0 is rapidly reversed upon re-plating onto a culture surface, wherein cells synchronously re-enter the cell cycle.In actively proliferating cultures, 24% of cells were found to be ciliated (Fig. 1a, b).This population heterogeneity is consistent with the asynchronous distribution in different cell cycle phases and the extension of primary cilia only during a restricted window of the cell cycle, in G1/G0.Suspension-arrested myoblasts showed an increase in ciliation with ~60% of G0 cells marked by primary cilia, which is reversed upon reactivation into the cell cycle, wherein the frequency of ciliated cells returned to 23% within 24 h (Fig. 1a, b).In addition, G0 cells showed an increased average ciliary length: 70% of cilia were > 2 μm (Fig. 1c).By contrast, proliferating cells display shorter cilia: only ~30% of cilia were > 2 μm (Fig. 1c).
To establish the dynamics of cilium acquisition during cell cycle exit, we isolated cells at different time points after induction of quiescence and stained for acetylatedtubulin.Cells showed a rapid increase in ciliation during the early stages of quiescence entry to reach a maximum of ~60% ciliated cells at 24 h post suspension, after which there was little change (Fig. 1d) until the time of harvest (48 h).The increase in ciliated cells correlates with the reduction of markers of proliferation Ki67 and Cyclin A2, and induction of the growth arrest-specific (GAS) gene PDGFRα [27] (Fig. S3).We also analysed the dynamics of cilium loss during reactivation out of G0, and observed a rapid and complete loss of primary cilia within 30 min of re-plating.With progression past the G0-G1 transition (6 h) and entry into S phase (12 h), the frequency of ciliated cells gradually increased, reaching ~20% at 24 h post re-plating (Fig. 1e), by which time the culture resembles a pre-quiescent population of asynchronously proliferating myoblasts (Fig. 1b).These results establish that as in other cell types, quiescence in myoblasts is associated with a reversible extension of longer cilia on a larger proportion of cells than observed in proliferating cultures.The increased frequency and length of cilia in G0 myoblasts, and the rapid loss of ciliation within minutes of reactivation supports the notion of the cilium as an inhibitory influence on cell cycle re-entry and progression, as described in other cell types [7].

Myogenic differentiation involves a transient ciliation event
Myoblasts exit the cell cycle irreversibly during myogenic differentiation.We used immunofluorescence analysis to establish the status of ciliation during entry into terminal arrest.By contrast to reversibly arrested myoblasts, terminally differentiated myotubes lacked primary cilia (Fig. 1f).However, we observed a transient phase of ciliation in early differentiating cultures prior to myoblast fusion, with a steady increase in the frequency of ciliated cells from as early as 6 h in differentiation medium, to a maximum of 53% observed at 24 h.High levels of ciliation were still observed at 48 h when fusion of myoblasts into myotubes was evident (Fig. 1g).A recent report of this transient ciliation suggested that Hh signaling during the induction of myogenic differentiation is mediated by the primary cilium and that deregulation of this node is associated with tumorigenesis (rhabdomyosarcoma formation) [28].Interestingly, although the time window of transient ciliation corresponds to the increase in expression of Myogenin, a key regulator of differentiation (Fig. 1h), we found little overlap between cells expressing Myogenin and those bearing a primary cilium in differentiating cultures (Fig. 1i).Our observations suggest that establishment of the myogenic program is specifically excluded in cells bearing primary cilia.

Reversibly arrested reserve cells in differentiating cultures are ciliated
Cultures triggered for myogenic differentiation yield a heterogeneous population: at 5 days post serum withdrawal, while ~80% of nuclei may be found in multinucleated myotubes, the remaining ~20% constitute mononuclear "reserve cells" [29], a pool of quiescent undifferentiated cells that retain the ability to re-enter the cell cycle if proliferative conditions are restored [30].We separated out the mononucleated (undifferentiated) fraction from differentiating cultures and verified that these cells do not express Myogenin (Fig. S4).As seen in suspension-arrested G0 myoblasts, ~50% of unfused reserve cells are ciliated (Fig. 1j) and account for the number of ciliated cells observed in early triggered cultures.Taken together, these results show that the presence of primary cilia distinguishes reversibly arrested from irreversibly arrested cells.

Quiescent satellite cells in vivo are marked by primary cilia
Recently, primary cilia were shown to be preferentially present on quiescent MuSC and not on activated MuSC, where they play a role in asymmetric cell division during the early regenerative response [20].To investigate whether cilia are associated with reversibly versus irreversibly arrested muscle cells in vivo, we isolated single muscle fibers from the mouse Extensor Digitorum Longus (EDL) muscle and probed for the cilium by immunofluorescence analysis.We detected primary cilia on ~80% of Pax7 + quiescent MuSC but no cilia were found on the rest of the myofiber membrane, either associated with myonuclei (Pax7 − ) (Fig. 1k) or elsewhere.Thus, as seen in cultured myoblasts, most quiescent MuSC are ciliated, while the differentiated myofibers are not, indicating that the cilium distinguishes reversible from irreversible arrest in vivo as well.

Abrogation of ciliary extension leads to loss of canonical features of quiescence
Ciliary extension critically depends on intra-flagellar transport (IFT), the process by which material is transported on stable microtubule tracks in the cilium core.IFT is regulated by a number of genes, primarily encoding motor proteins and adaptors that link cargo to the cytoskeleton [31].Since ciliary extension was triggered in myoblasts early after receiving cues for cell cycle exit, we postulated that the primary cilium is involved in the regulation of cell cycle exit programs and cell fate determination.To explore this connection, we investigated the ability of cells to enter G0 when ciliogenesis was blocked by siRNAmediated suppression of a key anterograde transport protein IFT88.siRNA treatment resulted in a 60% reduction in the mRNA levels (Fig. 2a) and a 92% reduction in protein levels of IFT88 (Fig. 2b, d), compared to control.Knockdown of IFT88 led to significant suppression of ciliogenesis (Fig. 2c, e).Interestingly, we observed a stage-specific effect of cilium knockdown.Contrary to earlier reports where blocking ciliogenesis in HeLa cells results in hyperproliferation [32], cilium knockdown in myoblasts under growth conditions, did not lead to increased proliferation, as evidenced by unchanged frequencies of Ki67 + cells (Fig. 2f).However, when knockdown cells were triggered to enter quiescence, we observed an increased frequency of Ki67 + cells compared to control cells (Fig. 2g).Knockdown populations also showed a decreased frequency of cells expressing the quiescence marker, p27 (Fig. 2h, S10), as well as a reduction in p27 protein levels (Fig. 2i).Interestingly, cell cycle analysis revealed that myoblasts lacking cilia display an altered FACS profile with an increased proportion of G2/M cells (Fig. 2j, k).In order to verify that the cell cycle phenotype seen in quiescent myoblasts in IFT88KD was a true effect of ablated ciliogenesis rather than a non-ciliary effect of the protein as reported earlier [33], we performed RNAi experiments using two other candidate members of the IFT system that are involved in anterograde transport into primary cilia, IFT20 and KIF3A.We found that IFT20KD and KIF3AKD follow similar phenotypes as IFT88KD (Fig. 2l, m, Fig. S5).These results demonstrate that suppressing ciliary extension compromises entry into G0 suggesting that the cilium is required for establishment of myoblast quiescence.
(See figure on previous page.)Fig. 1 Primary cilia are dynamically regulated and are preferentially present on quiescent myoblasts in vitro and in vivo.a Immuno-detection of primary cilia on C2C12 myoblasts in different cellular states -proliferating (MB), quiescent (G0), and 24 h after reactivation into cell cycle from quiescence (R24).Antibodies against γ-tubulin (γ-tub) mark centrosome and Acetylated tubulin (Ac.Tub) mark cilia; DAPI was used to label DNA.While only a small proportion of MB is ciliated, G0 myoblasts show a higher frequency of ciliation, which is reversed at R24. Lower panel shows magnified views (boxed insets).(Scale bar, 10 μm.).b Quantification of data showed in (a).Values are mean ± s.e.m., N ≥ 3, * p < 0.05.c The length of cilia present on MB, G0, and R24 cultures was measured and the distribution plotted.Cilia elaborated by G0 myoblasts are longer (Average length 2.2 μm) than those in proliferating cultures (Average length 1.8 μm).The numbers of ciliated cells analyzed are MB (40), G0 (30), and R24 (54).d Increased frequency of ciliated cells observed as proliferating myoblasts withdraw into quiescence during a time course of suspension arrest (0-48 h).Values are mean ± s.e.m., N ≥ 3, * p < 0.05, ** p < 0.01, *** p < 0.005.e Rapid loss of cilia during reactivation of G0 cells by re-plating on to an adhesive substratum, and gradual increase in ciliation during cell cycle re-entry.Values are mean ± s.e.m., N = 3, * p < 0.05.f Immunofluorescence analysis of cilia in 5 day differentiated culture of C2C12 cells: cilia are marked with Acetylated tubulin (Ac.Tub), Myogenin marks differentiated cells, DNA is labeled with DAPI (Scale bar, 10 μm).Cilia are absent on multinucleated myotubes (open arrowheads), and preferentially present on mononuclear (unfused) cells that lack Myogenin (closed arrowheads).Zoomed insets are representative of (1) unfused myoblasts (Myogenin − Cilium + , and (2) multinucleated myotubes (Myogenin + Cilium − ).g Quantitative analysis of ciliation during a time course of differentiation (represented by image in f) reveals increased proportion of ciliated cells as myoblasts withdraw into terminal differentiation.Values are mean ± s.e.m., N = 3, * p < 0.05.h Increased proportion of Myogenin + cells during the period where ciliation was analyzed in (g) highlights the time course of increase in Myogenin expression during differentiation.Values are mean ± s.e.m., N = 3, * p < 0.05.i Quantitative analysis of cells that are both ciliated and Myogenin + shows that ciliated cells are largely negative for Myogenin, indicating exclusion from differentiation.Values are mean ± s.e.m., N = 2. j Primary cilia are present on reserve cells in differentiated cultures.5-day differentiated C2C12 cultures were mildly trypsinised to remove myotubes, enriching the adherent undifferentiated mononuclear reserve cells.Cilia were detected using γ-tubulin (γ tub) and Acetylated tubulin (Ac.Tub).(Scale bar, 10 μm).The graph shows quantification of the data represented in the image, values are mean ± s.e.m. 47% of reserve cells in differentiated cultures are ciliated.k Primary cilia are associated with Pax7 + satellite cells (SC) but not with Pax7 − myonuclei (MN).Representative immunofluorescence image of a single fiber isolated from mouse EDL muscle and stained for Pax7 to mark satellite cells, and Pericentrin (PCNT) and Acetylated tubulin (Ac.Tub) to mark primary cilia, DAPI marks nuclei (Scale bar, 10 μm).The graph shows the quantification of the data derived from 4 independent staining experiments on single fibers: each experiment used fibers isolated from one mouse and a total of 121 nuclei were evaluated for the presence/absence of an associated cilium.

Transcriptome profiling reveals a stage-specific cell cycle block in IFT88 knockdown cells
To determine the scale of changes that contribute to the altered quiescent state, we used global transcriptome profiling.Proliferating myoblasts (MB) were transfected with either IFT88 siRNA or a non-targeting control siRNA following which cells were transferred to suspension culture for 48 h to induce quiescence (G0) and subsequently reactivated for 2 h by re-plating (R2).Affymetrix microarrays were queried with RNA isolated from cells in these three different states and a comparative analysis was performed between control and knockdown samples for each cell state.
The transcriptional profile of IFT88KD myoblasts was consistent with the stage-specific cell cycle phenotype observed.Maximum number of differentially expressed genes were observed at conditions of G0 (253 genes upregulated, 157 genes down-regulated), while fewer genes were altered in proliferating conditions (51 up-regulated, 31 down-regulated) (Fig. 3a, b, Fig. S6).Thus, IFT88KD myoblasts cultured in suspension have a transcriptional profile distinct from control cells (Fig. 3a).Interestingly, while IFT88KD myoblasts showed substantial alteration of the transcriptome in reactivated cells (R2) (168 genes up-regulated and 85 genes down-regulated) (Fig. 3a, b), the G0 and R2 profiles were very similar in IFT88KD cells, suggesting a deviant reactivation program when ciliary signaling is dampened.Overall, these observations suggest that suppression of the cilium leads to a block beyond which cells neither entered nor exited quiescence normally.
To identify networks among the genes that were differentially regulated in cilium knockdown cells at G0 we used the STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database [36].The upregulated genes formed a tight network, corresponding to 2 major clusters representing cell cycle regulation and myogenic differentiation, and a minor cluster consisting of genes regulating ECM and related signaling (Fig. 3c).The down-regulated genes formed a more diffuse network of genes involved in the regulation of cell proliferation and response to stimulus (Fig. S7).
To gain further insight into the molecular pathways that may be altered in cells where ciliogenesis is blocked, we used Gene Set Enrichment Analysis (GSEA analysis) [37,38], which showed an enrichment of genes related to cell cycle in the transcriptome of IFT88 knockdown cells (Fig. 3d).Notably, we observed an enrichment of genes that define a proliferation signature (pSig) when compared with earlier established gene sets defining quiescence and proliferation [39] (Fig. 3d).Interestingly, the top ten gene sets with highest normalized enrichment (See figure on previous page.)Fig. 2 Abrogation of primary cilia has distinct effects in proliferating and quiescent myoblasts.C2C12 myoblasts were transfected with siRNAs targeting IFT88 to block ciliogenesis, and were analyzed for effects of knockdown on proliferation and quiescence.Cells transfected with nontargeting siRNA were used as control.a qRT-PCR analysis was used to detect levels of IFT88 mRNA in control and siRNA treated samples 48 h after transfection.Knockdown myoblasts show reduced levels of IFT88 mRNA.Values represented are mean ± s.e.m., N = 3, **p-value = 0.0062.b Western blotting analysis demonstrates efficient knockdown of IFT88 protein levels 48 h after transfection.c Primary cilia were visualized by immunofluorescence labeling of Acetylated tubulin (Ac.tubulin) in IFT88 knockdown myoblasts cultured in quiescence-inducing conditions.Representative immunofluorescence images are shown.A single cilium is present per cell in control cultures, whereas knockdown cells show greatly reduced cilia (accompanied by increased tubulin staining in the cell body).d Quantification of the reduction of IFT88 protein level shown in (b).Values indicated are mean ± s.e.m., N = 3, ***p-value = 0.0004.e Knockdown of IFT88 caused reduction in frequency of ciliated cells.Values represent mean ± s.e.m., N = 4, ** p-value = 0.034.f Proliferation of cells was analyzed by quantification of Ki67 + cells, detected by immunostaining.Targeting ciliogenesis had little effect in proliferating conditions (MB).Values represent mean ± s.e.m., N = 3, ns = not significant (pvalue = 0.959).g When induced to enter quiescence (G0), knockdown cells displayed increased frequency of Ki67 + compared with control siRNAtransfected cells.Values represent mean ± s.e.m., N = 5, ** p-value = 0.0015.h Ablation of ciliogenesis suppresses expression of quiescence marker p27.C2C12 myoblasts stably expressing the p27 sensor (mVenus-p27K − ) were transfected with siRNA targeting IFT88or control siRNA and then cultured in quiescence-inducing conditions.The frequency of p27 (mVenus) positive cells was measured by flow cytometry.Values represent mean ± s.e.m., N = 3, * p-value = 0.0234.i Western blotting revealed the reduction of the quiescence marker, p27 in IFT88KD cells under conditions of G0 arrest.The values are relative protein levels compared to Control, and represent mean of quantification from 2 experiments.j Flow cytometric analysis of IFT88 siRNA-transfected cells cultured in conditions that allow proliferation (MB) does not show a significant shift in cell cycle profile.When placed in conditions that induce quiescence (G0) in control cells, knockdown cells display an increased proportion of cells in G2/M.Representative profiles are shown where the X axis shows the DNA content estimated by florescence intensity of the DNA dye DRAQ5, and the Y axis represents cell numbers as a percentage of maximum (normalized to mode).k Quantification of experiment described in (j) showing shifts in proportion of cells in different cell cycle stages (G1, S, G2/M).Values represent mean ± s.e.m., N = 3, * p-value = 0.043.l Flow cytometric analysis of IFT20 and KIF3A siRNA-transfected cells cultured in conditions that induce quiescence (G0).Cells lacking cilia show an altered cell cycle profile with an increase in proportion of cells in G2/M.Representative profiles are shown where the X axis shows the DNA content estimated by florescence intensity of the DNA dye DRAQ5, and the Y axis represents cell numbers as a percentage of maximum (normalized to mode).m Ablation of ciliogenesis using other IFT targets suppresses expression of quiescence marker p27.C2C12 myoblasts stably expressing the p27 sensor (mVenus-p27K − ) were transfected with siRNA targeting IFT20 or KIF3A, or control siRNA and then cultured in quiescence-inducing conditions.The frequency of p27 (mVenus) positive cells was measured by flow cytometry.Values represent mean ± s.e.m., N = 3, ** p-value = 0.0001, n.s (p-value = 0.38).score (NES) that were enriched in IFT88KD at quiescence conditions, represent genes that are involved in G2/M phases (Fig. 3e), which is consistent with the increased population of cells with 4 N DNA content observed in the cell cycle profile (Fig. 2j, k).Accumulation of cells at G2/M phases might reflect non-ciliary roles of IFT88 [33], or a requirement for IFT88 in G2/M entry, with blocked cells yet to exit the cell cycle at the time of harvest.Therefore, we analyzed a time course of entry into quiescence to reveal any shifts in cell cycle exit profiles.Surprisingly, cells lacking cilia displayed similar exit kinetics to control cells and showed no significant change in the proportion of either S phase cells (EdU + ) or M phase cells (Histone3pS10 + ) (Fig. 3f) across the time course.Taken together, these results suggest that IFT88KD cells do not enter M phase of the cell cycle but are paused at G2.
To resolve the paradox of increased proliferative gene expression in cells displaying a non-proliferative phenotype, we analysed the expression levels of selected G2/M candidate genes that were up-regulated in knockdown cells.Control G0 cells showed 50-500 fold lower expression of transcripts encoding G2/M regulators AurA, Bub1, Cdc20, CenpF, FoxM1, Mcm3/5, and Plk4 compared to proliferating MB.Although there is increased expression of these G2/M marker genes in IFT88KD cells, the levels of these transcripts are still significantly lower than those observed in proliferating cells (10-50 fold lower) (Fig. 3g), suggesting that overall, the magnitude of altered expression may not achieve a threshold required for functional progression through G2 into M. Furthermore, GSEA analysis also revealed that the upregulated genes in IFT88KD cells were involved in cell cycle checkpoints, in particular G2/M checkpoints involving DNA damage and mitotic spindle assembly, functions associated with the centrosome.Collectively, these results show that ablation of ciliogenesis leads to deviation from a canonical quiescence program, resulting in cells that are blocked at G2 rather than at G0, as defined by their transcriptional profile.
The second signature of up-regulated genes in IFT88KD cells represented myogenic regulators and muscle structural genes.Since the cilium is transiently induced during early differentiation and appears to be retained only in the minor population of nondifferentiating reserve cells, the induction of myogenic genes in IFT88KD cells suggests that loss of the cilium de-repressed the differentiation program.The upregulation of two gene networks representing opposing cell fates (cell cycle [G2/M] vs. differentiation) is consistent with the repression of both these states by signaling from the cilium.

Myoblasts lacking cilia have reduced self-renewal potential
Induction of quiescence is associated with enhanced clonogenic self-renewal when cells are restored to proliferative conditions [16].Cells with impaired ciliogenesis showed a deviation from canonical features of G0 such as expression of p27 and 2 N DNA content.Therefore, we evaluated the self-renewal capability in IFT88KD myoblasts post-quiescence using colony formation (CFU assay) and found that knockdown myoblasts showed reduced CFU compared to control cells (Fig. 4a, S8A).The kinetics of cell cycle re-entry were also affected, with a significant decrease in EdU (See figure on previous page.)Fig. 3 IFT88 knockdown cells display an altered quiescence program enriched with G2/M signature pathways.Affymetrix array-based transcriptional profiling was performed on IFT88 siRNA and control siRNA transfected myoblasts at 3 different cell cycle states -proliferating (MB), quiescent (G0), and 2 h after reactivation into cell cycle from G0 (R2).a Two-dimensional Principal Component Analysis (PCA) plot [34] of microarray data showing that IFT88KD G0 cells have a transcriptional profile that is distinct from Control G0 cells.Samples from different cell states can be seen as distinct clusters while the replicates of each sample cluster together indicating that they behave similarly.Control samples are shown as green dots and IFT88KD samples are shown as red dots.Note that control and knockdown cells do not show differences in the MB state, but are well separated in both G0 and R2 states, indicating a stage-specific requirement for the cilium.b Size adjusted Venn diagram [35] depicting the number of differentially expressed genes in MB, G0 and R2 states.IFT88KD cells showed the greatest deviation from Control (maximum number of altered genes) in G0. c STRING analysis was used to visualize the networks formed by the altered genes at G0. Genes upregulated in IFT88KD cells displayed a strong interaction network, which comprised of three major clusters: (i) Cell cycle-related, (ii) Myogenic genes (iii) ECM-related.Disconnected nodes (genes) are not displayed.d GSEA was used to identify pathways of genes in the altered transcriptome of IFT88KD G0 myoblasts.The plots show enrichment of genes related to cell cycle and proliferation.The gene sets used as reference here are the 315 genes annotated in the GO term GO:0007049 (CELL CYCLE) and 338 genes identified and collated as a proliferation signature (pSig).e Gene sets with highest enrichment scores identified in GSEA are specific to G2/M phases.NES represents Normalized Enrichment Score.f IFT88KD cells show cell cycle exit kinetics similar to Control cells and do not participate in DNA synthesis (EdU incorporation, top panel) or mitosis (H3pS10, bottom panel), during a time course of suspension culture.Values represent mean ± s.e.m., N = 2. g Validation of microarray analysis: qRT PCR analysis was used to determine mRNA levels of selected G2/M genes identified as up-regulated in IFT88KD G0 cells at conditions of proliferation (MB) and quiescence (G0).The relative mRNA levels were calculated in comparison to MB for Control and IFT88KD cells.Control cells show a characteristic repression of expression of mitotic regulators in G0 (varying from 50 to 500-fold reduction).Although IFT88KD G0 cells display higher expression levels than Control G0 cells, they still show lower expression of these proliferative genes when compared to cycling cells (MB) (varying from 10 to 50-fold lower).Values represent mean ± s.Interestingly, similar to the cell cycle effect, this reduction in self-renewal is also stage-specific, with no significant impairment in CFU formation in proliferating cultures (Fig. S9).These results show that even though IFT88 is expressed in both proliferating and quiescent cells, compromising its expression specifically in conditions where the cilium is normally elaborated leads to decreased self-renewal ability, consistent with the failure to enter a canonical quiescence program due to suppressed ciliary extension.

Myoblasts lacking cilia exhibit enhanced signaling activity
The cilium is a known sensory hub that harbors receptors for multiple signaling pathways (reviewed in [1,40]).The enrichment of Wnt, Hh and mitogen receptors in the cilium is thought to enable growth factor induced reactivation out of G0.In our culture model, quiescence is triggered by the abrogation of adhesion-dependent signaling pathways [26,41].To elucidate the mechanism by which suppression of the primary cilium contributes to an altered quiescent program, we examined possible shifts in signaling cascades.Consistent with the notion of aberrant signaling, GSEA analysis of IFT knockdown cells showed an enrichment of genes related to ciliumdependent signaling pathways (Fig. 4c), including Notch, Hh, Wnt, and growth factor signaling.The primary cilium is known to show cell type-and conditionspecific influences in either promoting or dampening the activity of these pathways [42].Using a combination of reporter assays, qRT-PCR and western blot analysis, we assessed the activity of cilium-related pathways that have been previously implicated in G0 [12,13,21].We detected enhanced signaling through 3 specific pathways in IFT88 knockdown myoblasts.Wnt signaling was elevated as evidenced by increased Wnt-TCF reporter activity (TOPflash), and induction of the transcriptional effector, active (dephospho) β-catenin (Fig. 4d, e).Increased levels of IGFR protein and phosphorylation were also seen, as well as increased levels of a key mediator of the G0-G1 transition, PDGFRA (Fig. 4f).These results suggest that under quiescence conditions, the primary cilium functions to dampen multiple growth factor signaling pathways.We further examined whether the observed induction of upstream growth factor signaling events in knockdown cells led to enhanced activity at downstream signaling nodes.mTOR activity is an important integrator of growth factor signaling and functions by targeting protein synthesis [43].Suppression of ciliary extension and increased growth factor signaling resulted in increased mTOR phosphorylation (Fig. 4g), consistent with relay of an upstream signal such as the activation of PDGFR or IGFR.Two critical downstream targets of mTOR that directly affect translational activity are ribosomal protein S6 (rpS6), which is activated upon phosphorylation by S6 kinase (S6K), and 4E-BP1, a translational repressor which is inactivated upon phosphorylation (reviewed in [44]).Interestingly, while rpS6 did not show an increase in phosphorylation (Fig. 4g), IFT88KD myoblasts showed an increase in phosphorylation of 4E-BP1.
We next investigated whether the observed increase in the level of mTOR activity towards 4E-BP1 had phenotypic consequences i.e. increased translation.Indeed, we observed an appreciable increase in levels of protein (See figure on previous page.)Fig. 4 The altered quiescent state in IFT88KD is characterized by reduced self-renewal, deregulated reactivation kinetics and enhanced activity of proliferative signaling pathways.a IFT88KD G0 myoblasts were harvested from suspension and replated at clonal density for analysis of colony forming potential (% CFU).IFT88KD myoblasts showed reduced self-renewal when compared to control cells Values represent mean ± s.e.m., N = 4, *p = 0.0068.b EdU incorporation was estimated during a time course of reactivation (6, 12, 24 h) from quiescence (G0).Knockdown cells are able to return to the cell cycle, but at 24 h when control cultures are restored to pre-suspension levels of DNA synthesis, IFT88KD cells showed significantly fewer cells in S phase.Values represent mean ± s.e.m., N = 3, *p-value< 0.05.c-g Mitogenic signaling pathways are up-regulated by IFT88 knockdown.c Gene sets related to pathways that are enriched in the IFT88KD G0 myoblast transcriptome.IFT88KD G0 transcriptome shows enrichment for genes related to cell signaling.d Western blot analysis shows an increase in levels of the Wnt effector β-catenin (dephospho or active) upon ablation of ciliogenesis.The values are relative protein levels in IFT88KD compared to Control, and represent mean of densitometric quantification from 2 experiments.e Functional analysis: IFT88KD G0 cells show elevated signaling through the Wnt pathway.Quantification of increased Wnt signaling as indicated by higher TOPFlash activity (expressed as RLU/ μg protein).Values represent ± s.e.m., N = 2, p-value = 0.014.f IFT88KD show elevated growth factor signaling in conditions of suspension arrest.Western blot analysis shows increased growth factor receptor protein expression for IGFR and PDGFRA, as well as increased activation of IGFR (p-IGFR).The values are relative protein levels in IFT88KD compared to Control, and represent mean of densitometric quantification from 2 experiments.g Activation of selective arms of translation regulatory pathways: IFT88KD cells at G0 conditions show elevated levels of phosphorylated (active) mTOR.However, there is divergent response in the key mTOR targets; while ribosomal protein S6 shows decreased phosphorylation, the translational repressor 4E-BP1 shows increased phosphorylation consistent with loss of repressive activity and indicating preferential activity through this effector node.The values are relative protein levels in IFT88KD compared to Control, and represent mean of densitometric quantification from 2 experiments.h IFT88KD myoblasts show increased global levels of protein synthesis.Protein synthesis was detected in Control and IFT88KD myoblasts which were placed under suspension arrest, using the Click-iT® Plus OPP Alexa Fluor® 488 Protein Synthesis Assay Kit .(Representative images in top panel), and the fluorescence levels were quantified in over 50 cells per sample (bottom panel).p value < 0.0001 synthesis in IFT88KD cells, as evidenced by the increased incorporation of OPP into newly synthesized proteins (Fig. 4h).Taken together, this data suggests that the loss of the cilium channels mitogenic signals leading to translational control via one arm (4E-BP1) of the signaling pathway downstream of mTOR.4E-BP1 phosphorylation has been shown to be specifically elevated during myogenic differentiation [45], and during M phase in cycling cells [46], consistent with the enrichment of cells displaying these two expression profiles in IFT88 knockdown conditions.Thus, presence of the cilium in G0 maintains 4E-BP1 in its translationally repressive unphosphorylated state and dampens overall protein synthesis.
In summary, this work shows the preferential extension of the cilium in muscle cells during reversible cell cycle arrest, but not during irreversible arrest/differentiation.Disruption of ciliation leads to an altered cell cycle distribution in IFT88KD cells, characterized by a block in G2 rather than G0, an altered transcriptome showing enhanced expression of regulators of G2/M checkpoints as well as enhanced expression of myogenic differentiation markers, reduced levels of classical quiescence indicators, aberrant re-entry and reduced self-renewal ability.Although suppression of ciliary extension correlated with increased activity of cilium-mediated signaling pathways, some downstream signaling nodes were not activated, likely accounting for the absence of enhanced proliferation.We conclude that cilium extension in quiescent myoblasts contributes to the acquisition of key features of G0 including the characteristic quiescence transcriptional program which features suppression of both cell cycle and myogenic genes, as well as dampened signaling to the translational machinery (Fig. 5).
Fig. 5 A working model for primary cilium function in myoblast quiescence.Control myoblasts entering suspension induced quiescence elaborate primary cilia, which work to suppress proliferative signaling in G0, whereas myoblasts where ciliogenesis is blocked enter an altered quiescent state characterized by arrest in G2 and activation of G2/M checkpoint genes.

Discussion
In this study, we have explored the function of primary cilia in skeletal muscle myoblasts and report a novel role for the organelle in the establishment of quiescence.Existing reports linking the cilium to the cell cycle have focussed on transition from G1 to S phase [32], as well as re-entry into the cell cycle from a quiescent state [7].These studies have shown that resorption and modification of ciliary length and function contribute to cell cycle progression.Similarly, signaling pathways such as Hh and PDGF that act through the cilium also contribute to regulation of cell cycle progression (reviewed in [47]) and asymmetric cell division [20].Although the primary cilium is a quiescence-associated structure, there is surprisingly scanty evidence of the function of the cilium-centrosome axis in the establishment and maintenance of G0.Our study aims to bridge this gap, and provides evidence for repressive signaling from the cilium in establishment and/or maintenance of quiescence.
We show that as in other cells, primary ciliogenesis in myoblasts is induced during cell cycle exit.While this concurs with the well-established reciprocal relationship between the cilium cycle and cell cycle, we find that ciliation is sustained in quiescence, but lost in terminal differentiation, suggesting a role in distinguishing these non-dividing states.Increased ciliogenesis is observed within 24 h of receiving cell cycle exit cues, a time when expression of classical quiescence markers (such as p27) or differentiation markers (such as Myogenin) is still low.Thus, formation of the cilium occurs subsequent to receiving cell cycle exit cues (reduced adhesion or reduced mitogen availability), but prior to maturation of quiescent or differentiated states, implying a role for the signaling hub in this cell fate decision.
Our data also indicate that primary cilia mark only those cells that exit the cell cycle reversibly.Thus, in apparently homogeneous cultures induced to undertake myogenic differentiation, the primary cilium may help to select a sub-population (reserve cells) that resists differentiation cues and retains stem cell characteristics.We also observe heterogeneity in ciliation status within quiescent populations, suggesting that not all cells entering quiescence are equally capable of ciliogenesis.This may have implications for quiescent adult muscle stem cells, where some evidence exists indicating a correlation between heterogeneity and different propensities for activation [48].In conjunction with the finding that primary cilia are preferentially assembled on self-renewing MuSC in vivo ( [20]; this study), this observation provides interesting avenues for future exploration of ciliary roles in adult stem cell fate.The enhanced signaling activity in cilium-ablated cells may indicate that the cilium provides a zone of regulation for these pathways, possibly by physically localizing signaling intermediates or mechanical cues.These results further support earlier reports of cell type dependent effect of the primary cilium on signaling pathways such as Wnt and Hh (Reviewed in [42]), and that in some contexts, the cilium can dampen signaling.
Strikingly, when cells with impaired ciliogenesis were placed under conditions of suspension arrest, while they did not attain canonical quiescence, they also did not display continued cell proliferation as might be expected from increased mitogenic signaling.Instead, cells lacking cilia showed no increase in EdU incorporation (S phase), or phosphorylated H3 (M phase).This is also in contrast with an earlier report highlighting the role of IFT88 in mitotic spindle organization, which showed that IFT88 knockdown in HeLa cells resulted in mitotic catastrophe and delays [33].The observed accumulation of cells in the G2 phase in IFT88KD suggest either an incomplete reversal/residual retention of ciliary function, or the presence of an unreported role for the cilium in G2/M checkpoint control.The dual role of the centrosome/ basal body and the relationship of the centrosomecilium cycle with the cell cycle suggests possible mechanisms.When IFT is impaired, while ciliary extension is compromised, the centrosome-basal body transition may not be affected.Thus, the centrosome is likely still docked as a basal body and may not be free for its mitotic spindle functions, resulting in G2/M accumulation of cells unable to traverse that block.Alternatively, the centrosome and basal body may have differential affinity for localization for cell cycle regulatory proteins and disturbing this balance may result in aberrant regulation of cell cycle progression, since cell cycle regulators such as Cyclin B1, AurA and Plk1 are reported to localize to the centrosome or basal body and this localization is essential for their function [49,50].Since the centrosome alternates as a spindle organizer in M phase and the ciliary basal body in G0, it is possible that communication between the cilium and G2/M transcriptional networks is mediated by centrosomal components.Thus, the incomplete induction of G2/M transcriptional programs observed in the IFT88KD cells may reflect an additional requirement for transition to centrosomal function required for M phase entry.
The centrosome has also been linked with translational regulation of M phase regulators: Plk1was shown to co-localize with 4E-BP1 at centrosomes in mitotic cells, providing evidence for direct interaction and phosphorylation of 4E-BP1 by Plk1 [51].Interestingly, IFT88KD myoblasts induced to enter quiescence also show an increase in Plk1 expression as well as 4E-BP1 phosphorylation, accompanied by an increase in global protein synthesis.Considering emerging evidence of the functional association of translational machinery with the centrosome, our finding raises the possibility that the cilium-centrosome axis is involved in signaling between G2/M transcriptional programs and selective translation of 4E-BP1-dependent transcripts.
Our data shows that cells lacking cilia enter an alternate state of quiescence, characterized by a shift in transcriptional profile, as well as reduced ability to return to the cell cycle.These observations suggest that the primary cilium may also act as a transducer for quiescence signals and that the cilium-mediated establishment of the quiescence program is necessary for subsequent cell cycle re-entry and self-renewal.Disruption of primary ciliogenesis leads to a specific alteration of G2/M transcriptional networks in conditions of quiescence, suggesting a novel role for the control of mitotic progression by the cilium-centrosome axis.
The observed halt at G2 in cilium-ablated cells indicates the requirement of signals in addition to the canonical G2/M transcriptome for transition through mitosis, such as mechanical cues, which may in turn also be sensed, transduced and regulated by the cilium/ centrosome.While there is some evidence to support such non-ciliary roles for ciliary and centrosome proteins [33,46,[51][52][53], our findings may also suggest that retention of the centrosome as a basal body is monitored, and that complete dismantling of the cilium/basal body complex is required for the G2/M transition.

Conclusions
Taken together, our study shows that primary cilia are dynamically assembled when myoblasts receive cues for cell cycle exit, and strongly support the idea that the cilium is integral to the establishment of the quiescent state in skeletal muscle myoblasts, while also distinguishing reversible from irreversible arrest.The functions of the cilium in quiescence are mediated in part through the dampening of proliferative signaling, and suppression of cell cycle progression by influencing G2/M control mechanisms.

Methods
Cell culture C2C12 myoblasts were obtained originally from H. Blau, Stanford University and a sub-clone A2 was derived inhouse (Sachidanandan et al., 2002) and used for all experiments.Myoblasts were maintained in growth medium (GM; DMEM + 20% Fetal Bovine Serum (FBS) and antibiotics), ensuring that cells do not exceed 70-80% confluence before passaging.Differentiation was induced in low mitogen medium (DM: DMEM + 2% horse serum), for 5 days to form mature myotubes (MT).Synchronization in quiescence (G0) was induced by suspension culture of myoblasts for 48 h in 1.3% methylcellulose medium prepared with DMEM containing 20% FBS, 10 mM HEPES, and antibiotics, as described [21].G0 cells were reactivated into the cell cycle by re-plating at a sub-confluent density in GM and harvested at defined times (30 min -24 h) after activation.

Single muscle fiber analysis
Animal experiments were conducted with prior approval of Institute for Stem Cell Biology and Regenerative Medicine (InStem) and Center for Cellular and Molecular Biology (CCMB) Institutional Animal Ethics Committees and in accordance with Committee for the Purpose of Control and Supervision of Experiments on Animals (CPCSEA) guidelines.
Single skeletal muscle fibers were isolated using a technique modified from [55,56].Briefly, 6 to 8-week old male C57 BL/7 mice were euthanized by cervical dislocation, after which the Extensor Digitorum Longus (EDL) muscles were dissected from the both hindlimbs.The isolated muscle tissue was then treated with Collagenase Type 1 (Cat# LS4196 Worthington 400 U/ml final concentration) in DMEM for 1 h at 37 °C.The resultant dissociated muscle fibers were transferred into fresh DMEM and triturated to release individual muscle fibers, which were then washed through transfer to fresh medium, and then fixed with 4% paraformaldehyde for 15 min at RT.For immunostaining, single fibers were washed with PBS twice and mounted on charged slides (Cat#12-550-15, Fisher Scientific).The fibers were permeabilized with 0.5% Tween 20 in PBS for 1 h at RT followed by blocking with 2 mg/ml BSA in PBS, for 1 h at RT.Primary antibody incubations were performed overnight at 4 °C.Secondary antibody incubations were for 1 h at RT using fluorescently tagged antibodies from Invitrogen.All antibody dilutions were made in a solution of 1 mg/ml BSA in 0.25% Tween 20 in PBS.Antibody incubations were followed by three washes with blocking solution.4′, 6-Diamidino-2-Phenylindole (DAPI) (Cat# 32670 Sigma) was used to stain the DNA.Images were acquired using Zeiss LSM 510 Meta confocal microscope.Satellite cells were identified by their expression of nuclear Pax7 (Pax7+).Nuclei with no Pax7 expression (Pax7-) largely represent myonuclei, as practically no fiber-associated non-muscle cells co-purify in enzymatically disaggregated fiber preparations.

Cell cycle analysis
Adherent cells were trypsinized, washed in PBS and pelleted by centrifugation.Suspension-arrested cells were recovered from methylcellulose by dilution with PBS followed by centrifugation as described earlier.Cell pellets were dispersed in 0.75 ml of PBS, and fixed by drop wise addition into 80% ice-cold ethanol with gentle stirring, following which they were briefly washed with PBS and resuspended in PBS with 40 μM of the DNA dye DRAQ5™ (Cat.No DR50050, Biostatus) per 10 6 cells.Cell cycle analysis was performed on a FACS Caliber® Cytometer (Becton Dickenson) using CelQuest® software and analyzed using FlowJo® software.At least 10,000 cells were acquired for each sample.Forward scatter and side scatter were used to gate cell populations and doublets were removed from analysis.

Analysis of p27 levels
A stable C2C12 line expressing a p27 sensor (p27-mVenus, a fusion protein consisting of mVenus and a defective mutant of p27-CDKI, (p27K(−)) [57] was generated, transfected with either control or IFT88 siRNAs, and placed in suspension arrest.After recovery from methyl cellulose medium, these cells were assayed for proportion of mVenus positive cells by flow cytometry using a FACS Caliber cytometer (Becton Dickenson).CelQuest® software was used for acquisition and FlowJo® software was used to analyze the data.

EdU incorporation assay
Cells to be analyzed for DNA synthesis were pulsed with 10 μM EdU for 30 min, washed in PBS and fixed as described earlier.Samples were stored in PBS at 4 °C till further processing.For staining, cells were permeabilized and blocked in PBS having 10% FBS and 0.5% TX-100 and labeling was detected using Click-iT® imaging kit (Invitrogen) as per manufacturer's instructions.Samples were mounted in aqueous mounting agents with DAPI.Images were obtained using the Leica TCS SP8 confocal microscope.Minimum global changes in brightness or contrast were made, and composites were assembled using Fiji (ImageJ).

OPP incorporation assay
Cells to be analyzed for protein synthesis were pulsed with 20 μM OPP for 30 min, washed in PBS and fixed as described earlier.Samples were stored in PBS at 4 °C till further processing.For staining, cells were permeabilized and blocked in PBS having 10% FBS and 0.5% TX-100 and labeling was detected using Click-iT® imaging kit (Invitrogen) as per manufacturer's instructions.Samples were mounted in aqueous mounting agents with DAPI.Images were obtained using the Leica TCS SP8 confocal microscope.Image intensity was calculated using Fiji (ImageJ) software, and corrected mean intensity (CMI = total intensity of signal -(area of signal × mean background signal)) was determined.
Knockdown of gene expression using siRNA C2C12 cells were cultured in growth medium until 80% confluent.The cells were then trypsinized and plated on to tissue culture dishes at appropriate cell density according to size of the dish.Approximately 16 h post plating, the cells were transfected with siRNA using Lipofectamine RNAiMAX (Invitrogen) as per manufacturer's instructions.The cells were incubated with the RNA-Lipid complex for at least 18-24 h following which they were used for further experimental analysis.Both single siRNA as well as SMARTpools of 4 siRNAs were used for knockdown experiments.At least one replicate of all experiments, barring the microarray assay, was performed using SMARTpools.Typically, siRNA-mediated knockdown was in the range of 50-90% at RNA level.Details of siRNA used in this study are as follows.For IFT88 knockdown: 5′-GCUGUGAACUCGGAUAGAU-3′ (Eurogentec) or siGENOME Mouse Ift88 (21821) siRNA-SMARTpool (Dharmacon, M-050417-00-0010); For Control: Scrambled siRNA (Eurogentec SR-NP001-001) or siGENOME Non-Targeting Pool #1 (Dharmacon, D-00126-13-20).
Global transcription profiling using microarray 1 μg of total RNA isolated from Control or IFT88 KD cells was converted to cDNA using One-cycle labeling kit and amplified using IVT labeling kit following manufacturer's instructions (Affymetrix).The normalized cRNA was fragmented, hybridized to mouse Affymetrix Gene-chips (430A 2.0), washed, stained and scanned as per Affymetrix protocols.The experiment was repeated with three different biological replicates and data analyzed using Affymetrix Gene Chip operating software (GCOS).All the CEL (cell intensity) files generated by Expression Console were then loaded into the R Bioconductor "affy" package for microarray analysis [58].Briefly, the CEL files were normalized using Loess Normalization before proceeding to further analysis of differential expression between Control and IFT88KD samples for each specific cell state.Genes showing > 1.5-fold differential expression with P ≤ 0.05 were selected and a subset validated by real time Q-RT-PCR.The raw data is available on the GEO database (Series GSE110742).

STRING analysis
The lists of genes that were either upregulated or down-regulated in IFT88KD myoblasts under G0 conditions were used as input to search for networks in the STRING database using default settings with only connected nodes being represented in the network diagram.
GSEA Normalized expression values for Control (SCR) and IFT88KD were used as input dataset and the gene sets used were those from the C5 collection and those related to signaling, which are available on the Molecular Signatures Database (MSigDB), as well as a gene set identified from earlier published datasets (pSig).The enrichment analysis was carried out as described in the GSEA User guide (http://software.broadinstitute.org/gsea/doc/GSEAUserGuideFrame.html).
Luciferase activity was measured in lysates of TFC-1 and FFC-15 cells that were transfected with either IFT88 or negative control siRNAs.Assays were performed using the Luciferase Reporter Gene Assay Kit (Roche) to obtain data as relative light units (RLU).TCF activity was finally expressed as an Arbitrary Unit (AU), which was obtained after normalizing RLU to total protein estimated using BCA kit (Pierce).

Statistical analysis
Unless otherwise mentioned, all data represented are values derived from at least 3 biological replicates and is represented as mean ± s.e.m., analyzed using Student's ttest, where p < 0.05 was taken as significant.