|Home | About | Journals | Submit | Contact Us | Français|
The Hedgehog (Hh) and Wnt signal transduction pathways are master regulators of embryogenesis and tissue renewal and represent anticancer therapeutic targets. Using genome-wide RNA interference screening in murine cultured cells, we established previously unknown associations between these signaling pathways and genes linked to developmental malformations, diseases of premature tissue degeneration, and cancer. We identified functions in both pathways for the multitasking kinase Stk11 (also known as Lkb1), a tumor suppressor implicated in lung and cervical cancers. We found that Stk11 loss resulted in disassembly of the primary cilium, a cellular organizing center for Hh pathway components, thus dampening Hh signaling. Loss of Stk11 also induced aberrant signaling through the Wnt pathway. Chemicals that targeted the Wnt acyltransferase Porcupine or that restored primary cilia length by inhibiting the tubulin deacetylase HDAC6 (histone deacetylase 6) countered deviant pathway activities driven by Stk11 loss. Our study demonstrates that Stk11 is a critical mediator in both the Hh and the Wnt pathways, and our approach provides a platform to support the development of targeted therapeutic strategies.
Studies spanning nearly a half-century aimed at elucidating the underlying mechanisms of embryonic development in metazoans have yielded a handful of secreted intercellular signaling molecules that determine cell fate. These include members of the Hedgehog (Hh) and Wnt families, which additionally maintain homeostatic tissue renewal in adult animals and frequently contribute to degenerative disease and cancer (1, 2). Despite some successes in pharmacologically influencing the Hh and Wnt pathways (3, 4), the effective deployment of small molecules that target these pathways is limited by our understanding of pathway complicity in disease etiology.
The coordinated action of several well-defined pathway components transduces signals from Hh or Wnt into cellular responses specific to each ligand (5, 6). Essential to the fidelity of these signaling processes is the dynamic redistribution of pathway components to subcellular compartments that function as signaling relay points upon pathway stimulation. In the Hh pathway, the antenna-like primary cilium functions as an assembly area for pathway regulators (7) and contains the molecular machinery that mediates constitutive proteolytic processing of Gli transcriptional activators into repressor molecules (8–10). Activation of Hh-mediated response first entails engagement of Hh with the Dally-like protein (Dlp), Cell adhesion molecule–related/down-regulated by oncogenes (Cdo), and Patched (Ptch) receptors (5). In response to Ptchmediated reception of Hh, the Smoothened (Smo) transmembrane effector accumulates in the primary cilium, disengages Gli protein processing, and induces Gli-dependent transcription of target genes (11). Wnt is a ligand that stimulates multiple downstream pathways, those involving activation of the transcriptional coactivator β-catenin and those independent of the β-catenin effector (so-called canonical or noncanonical Wnt pathway responses, respectively) (2). Common to the known Wnt pathways is the activation of the Dishevelled (Dvl) signaling molecule, which typically entails its phosphorylation and redistribution from the cytoplasm to membrane-bound receptors in response to Wnt stimulation (2).
We performed a genome-wide RNA interference (RNAi) screen in cultured mouse cells to identify genes involved in both Hh and Wnt signaling. From this gene set, we discovered unanticipated associations between these pathways and genes previously linked to disease. We demonstrated that the tumor suppressor serine-threonine kinase and cell polarity regulator Stk11 (also known as Lkb1) functions as a component of both Hh and Wnt pathways. We found that Stk11 influenced the response to Hh by controlling the length of the primary cilium and Wnt responses by gating the abundance of phosphorylated Dvl protein. We formulated chemically based approaches to impede the altered cell signaling in cancerous cells lacking Stk11, which suggests that our screening strategy may be useful in defining molecular targets for therapeutic approaches for other diseases associated with multifunctional genes.
To identify previously unknown Hh pathway components, we performed two parallel genome-wide RNAi screens in cultured cells by transiently introducing gene-specific small interfering RNA (siRNA) pools (four non-redundant siRNAs) from two mouse siRNA libraries into 3T3-ShhFL cells (Fig. 1A; fig. S1, A to F; and tables S1 and S2). The 3T3-ShhFL cells enable monitoring of Sonic hedgehog (Shh)–induced transcriptional response in a cell-autonomous manner with a well-established Hh pathway reporter (the GliBS reporter) (12). Because these cells produce Shh, targeting of genes encoding proteins involved in either Shh production or downstream signaling can reduce reporter activity. With this system, we screened a combined total of 39,000 siRNA pools, with ~15,700 of these mutually targeting the same genes with mostly nonoverlapping siRNA sequences (13) (Fig. 1A). Focusing on positive pathway regulators (siRNAs that reduced reporter activity greater than two SDs from the mean signal), we identified 535 genes of interest (tables S3 and S4). Recovery of known Hh pathway components such as Smo, the seven-transmembrane pathway effector, and components of the primary cilium was improved with the use of both siRNA libraries, suggesting an advantage in using multiple siRNA libraries in RNAi-based screening strategies (Fig. 1B). We were also intrigued by the sparseness of overlapping genes identified from the initial parallel screening of two siRNA libraries and further interrogated this phenomenon with additional RNAi-based studies using conditions that minimize inherent noise associated with high-throughput screening platforms (fig. S1, G and H). These additional studies revealed that the integration of RNAi-derived data sets is limited by the effectiveness of the siRNA design and current high-throughput screening formats and showed the need to screen multiple siRNA libraries to reduce false-negative discovery rates.
We retested the siRNAs against the 535 genes of interest in ShhLightII cells stimulated with exogenous Shh to enrich for genes that act downstream of Shh protein production (Fig. 1C). ShhLightII cells are NIH 3T3 cells that do not endogenously produce Shh and are stably transfected with the GliBS reporter. We counterscreened the resulting hits against the Wnt/β-catenin pathway in NIH 3T3 cells with a Wnt-responsive luciferase reporter [SuperTopFlash (STF)] (Fig. 1, C and D). Genes when targeted with siRNAs that decreased activities of both pathways may influence cell viability and were not considered for further testing. We cross-referenced the remaining genes with the OMIM (Online Mendelian Inheritance in Man) database to identify previously unknown associations between diseases and the Hh and Wnt pathways (Fig. 1, C and D).
We confirmed the specificity of siRNA pools targeting the disease-associated genes by evaluating the ability of individual siRNAs from the pool to recapitulate the original observation (fig. S2 and table S5). Genes already validated by identification from two siRNAs pools were not retested. Using a minimum threshold of three distinct siRNAs inducing a shared phenotypic outcome as the criteria for high-confidence genes with minimal influence from RNAi off-targeting (14), we retained a total of 10 genes associated with either malformation, degenerative disease, or cancer (Table 1). The heparan sulfate proteoglycan glypican 1 (Gpc1) is a homolog of Dlp that functions as part of the receptor complex for Hh in Drosophila (15). Prkar1a, the regulatory subunit of protein kinase A (PKA), was previously identified as critical to Hh pathway response in a kinome-wide RNAi screen (14). We also identified two gene products linked to ciliopathies: retinaldehyde-binding protein 1 (Rlbp1) and intraflagellar transport 57 (Ift57) (16, 17). With the exception of Smo and Gli2, the Hh-related functions of the remaining genes [the serine-threonine kinase Stk22b; diacyl-glycerol kinase θ (Dgkq); and the achalasia, adrenocortical insufficiency, alacrima syndrome (Aaas) gene] are unknown. The Stk11 tumor suppressor controls microtubule dynamics, metabolism, and cell polarity but has not been previously implicated in Hh-mediated signaling (18).
Among the genes identified as potential Hh pathway regulators, Prkar1a and Stk11 are both recognized as tumor suppressors (Table 1). By analyzing the transcriptional response to Shh in mouse embryonic fibroblasts (MEFs) null for PRKAR1a or STK11, we confirmed that loss of either of these two molecules decreased Hh pathway response (Fig. 2, A to C). Whereas the requirement for Prkar1a in Hh-induced expression of Ptch1 has been reported, its mechanism of action within the pathway has not been well interrogated. In accordance with the role of PKA in Gli processing, loss of Prkar1a results in an increase in the amount of Gli3 repressor (Gli3R), presumably as a result of unregulated PKA activity (19) (Fig. 2D). STK11−/− MEFs also exhibited increased Gli3R abundance relative to wild-type MEFs that cannot be accounted for by changes in the Gli3 expression (Fig. 2E and fig. S3A). Unlike in the PRKAR1A-null cells, the production of Gli3R in STK11-null cells was reduced by Shh protein or a Smo agonist (SAG), suggesting that Smo function was intact (Fig. 2E). Pulse-chase analysis of the rate of Gli3R formation after addition of either SAG or Hh-conditioned medium confirmed that the rate of Gli3R formation was increased by nearly threefold in STK11-null MEFs relative to that in wild-type MEFs (Fig. 2F and fig. S3B). Furthermore, in the presence of a proteasome inhibitor, which blocks Gli3R production, the stability of Gli3R does not appear to be increased in STK11-null MEFs when compared to that in wild-type MEFs. Thus, alteration of the Gli3R turnover rate does not appear to contribute to the accumulation of Gli3R observed in STK11-null MEFs (Fig. 2G).
To further understand the mechanism underlying the Stk11-induced increase in the abundance of Gli3R, we examined the primary cilium, an organelle essential for proteolytic processing of Gli molecules into transcriptional repressors (7). The cilia in STK11−/− cells were on average one-half the length of their counterparts in wild-type cells, suggesting a defect in achieving the normal length of this organelle in STK11−/− cells (Fig. 3A and fig. S3C). Examination of cilia in STK11−/− cells by transmission electron microscopy confirmed the establishment of a basal body, the cellular anchor of cilia, and the ability of STK11−/− cells to elaborate rudimentary cilia (Fig. 3B). Because the dynamics of cilia formation and disassembly are influenced by the cell cycle, and STK11−/− MEFs double more rapidly than do wild-type MEFs (20) (fig. S3D), we analyzed the propensity of the wild-type and STK11−/− MEFs to exhibit growth arrest under conditions typically used to examine primary cilia (low serum, high cell density). Both cell lines achieved cell cycle arrest at confluency, suggesting that the differences in primary cilium length are not attributable to changes in the duration of cell cycle phases (Fig. 3C).
The tubulin deacetylase HDAC6 (histone deacetylase 6) promotes ciliary disassembly and thus influences cilia length (21). Addition of suberoylanilide hydroxamic acid (SAHA), an inhibitor of deacetylases including HDAC6 (22), increased cilia length in STK11−/− cells but not in wild-type cells, suggesting that STK11−/− cells undergo excessive cilia disassembly (Fig. 3, D and E). Consistent with this model, the HDAC6 inhibitor tubacin, but not the HDAC1 inhibitor MS-275, restored cilia length in STK11−/− cells (Fig. 3E and fig. S3E). The extension of primary cilia in SAHA-treated cells also correlated with restoration of Gli3R to an amount similar to that in wild-type cells (Fig. 3F). Concomitantly, SAHA increased the abundance of acetylated tubulin, which is decreased in STK11−/− cells, compared to wild-type cells (Fig. 3F). Considering that Gli proteins must cycle through the primary cilium to acquire PKA-dependent phosphorylation before proteolytic processing in the cytoplasm (23–25), the increased abundance of the GliR molecules may be due to decreased transit time of full-length Gli molecules in the shortened cilia. However, we cannot rule out the contribution of other Stk11-dependent mechanisms that may directly influence Gli processing.
We also examined in STK11−/− cells the ability of Hh pathway components to accumulate in the primary cilium in response to Shh (26, 27). Consistent with the ability of Hh and SAG to block Gli3 processing (Fig. 2E), the translocation of Smo into the cilium and accumulation of Gli2 protein at the tip of the cilium in response to Shh appeared normal in the STK11−/−cells (fig. S4, A and B). Thus, the excessive accumulation of Gli3R in STK11−/− cells is not likely due to altered subcellular localization of Smo or Gli proteins. We also confirmed, using a membrane-associated protein marker (Smo), that cilia length is decreased in STK11−/− cells (fig. S4C).
We also identified Stk11 as a suppressor of Wnt-mediated responses (Fig. 1D). Examination of biochemical changes associated with Wnt/β-catenin pathway activity in STK11−/− MEFs revealed that the amount of phosphorylated Dvl2 was increased compared to that in wild-type MEFs (Fig. 4A and fig. S5A). This phosphorylation was not sensitive to HDAC6 inhibitors, suggesting that in MEFs, acetylated tubulin abundance does not affect Dvl protein phosphorylation (fig. S5B). Phosphorylation of Dvl2 by casein kinase 1 (CK1) family members at consensus CK1 phosphorylation sites contributes to both β-catenin–dependent and noncanonical Wnt signal transduction (28). We interrogated altered pathway activity with a synthetic Wnt pathway inhibitor that targets the pathway at the level of Wnt production [inhibitor of Wnt production (IWP), an inhibitor of the acyltransferase Porcupine] (29) (Fig. 4B). The IWP compound restored the ratio of phosphorylated to unphosphorylated Dvl proteins in STK11−/− MEFs to that observed in wild-type MEFs, suggesting that the increased Dvl activity is dependent on production of endogenous Wnt protein (Fig. 4C and fig. S5C). Whereas specificity studies associated with the IWP compounds have been previously described (29), we provide here additional support for the mode of action attributed to these compounds (fig. S5, D and E). Consistent with the ligand-dependent activity of Dvl in STK11−/−MEFs, inhibition of CK1, a kinase essential to ligand-mediated Dvl activity (28), reduced Dvl phosphorylation (Fig. 4D), whereas an inhibitor of PKA had little effect on Dvl phosphorylation (fig. S5F).
Dvl proteins function as cytoplasmic routers that can transduce both β-catenin–dependent and –independent Wnt-mediated responses (30). RNAi targeting of Dvl2 reduced the rate of growth in STK11−/− MEFs, suggesting that Dvl2 plays a key role in promoting aberrant cell proliferation (Fig. 4E). In response to Wnt, Dvl adopts a diffuse cytoplasmic distribution that contrasts with its localization to punctate structures in the absence of Wnt (31) (fig. S6). We found that there was abundant diffusely localized Dvl2 [fused to green fluorescent protein (GFP)] in STK11−/− MEFs even in the absence of added Wnt, and this contrasted with the distribution of a GFP-tagged Axin1, a binding partner of Dvl (31), which retained its typical punctate cellular distribution (Fig. 4F). Although IWP restored ratios of phosphorylated to unphosphorylated Dvl2 in the STK11−/− cells to that observed in wild-type MEFs (Fig. 4C), it failed to restore punctate Dvl distribution (Fig. 4F). Together, these findings suggest that loss of Stk11 results in altered Dvl protein distribution that facilitates Wnt-dependent Dvl phosphorylation.
We investigated the role of Stk11 in vivo by analyzing zebrafish embryos treated with morpholino oligonucleotides (MOs) targeting Stk11. Consistent with reduced Hh signaling, the embryos had U-shaped somites when analyzed at 24 hours after fertilization (Fig. 5A). We also noted that some treated embryos exhibited shortened tails, a phenotype that may arise from altered noncanonical Wnt-mediated responses (32) (Fig. 5A and fig. S7). Engrailed 1a (Eng1a) expression, which is associated with Hh signaling (33), was also reduced in the Stk11 morphants in the midbrain-hindbrain boundary and muscle pioneers (Fig. 5B). We also noted little change in β-catenin–dependent Wnt signaling (axin2 expression) in zebrafish from the same developmental stage, suggesting potential differences in Stk11-mediated regulation of Wnt/β-catenin signaling across species or between developing and adult tissues (Fig. 5C). Together, these data are consistent with a role for Stk11 in both Hh and Wnt signaling in vivo but do not rule out the possibility that other functions of Stk11, such as its activity in cell polarity, also contribute to the observed developmental defects.
Whereas the Stk11-associated Peutz-Jeghers hereditary tumor syndrome is a relatively rare disease, sporadic cancers, such as those from the lung and cervix, frequently harbor mutations in STK11 (34, 35). We investigated the relevance of Stk11 loss in Wnt signaling in a panel of cells derived from cervical carcinoma. Consistent with our findings in MEFs, cervical carcinoma cell lines, previously described as lacking Stk11 (35), all exhibited enhanced Dvl2 phosphorylation relative to normal endometrium-derived cells (Endo cells) (Fig. 6A). A similar link between Stk11 status and Dvl phosphorylation was also present in cell lines derived from another Stk11-associated cancer, non–small cell lung cancer (NSCLC) (fig. S8A). Notably, cells from this cancer type lacking STK11 generally exhibit poor ciliogenesis, a phenotype that can be partially rescued by reintroduction of Stk11 DNA (fig. S8, B and C).
Despite the observation that the STK11-null cervical carcinoma cell lines uniformly exhibited an increase in Dvl phosphorylation, only a subset is positive for c-Myc protein, which is encoded by a target gene of the Wnt/β-catenin pathway (Fig. 6A). For example, HeLa and SiHa cells both lack Stk11, but only in HeLa cells was c-Myc detected. This suggests that the outcome of STK11 loss is context specific with respect to aberrant engagement of Wnt signaling.
We confirmed the role of Stk11 in controlling Wnt/β-catenin pathway activity in HeLa cells. Consistent with the presence of c-Myc in these cells, we observed increased activity of the STF reporter in HeLa cells relative to that in Endo and SiHa cells (Fig. 6B). This reporter activity was sensitive to the IWP compound that disrupts Wnt ligand production, as well as a chemical that inhibits Wnt/β-catenin pathway activity by stabilizing Axin, a negative regulator of β-catenin [the inhibitor of Wnt response (IWR)] (Fig. 6B). We also noted in HeLa, but not SiHa cells, the presence of Axin2 and phosphorylated LRP6, both hallmarks of Wnt/β-catenin pathway activity (fig. S9).
Introduction of Stk11 into HeLa cells reduced the Wnt-mediated transcriptional response (Fig. 6C), partially restored the punctate distribution of Dvl2 (Fig. 6D) (Axin1 distribution shown in fig. S10), and promoted the appearance of nonphosphorylated Dvl2 (Fig. 6E). Expression of a kinase-dead (KD) form of Stk11 did not result in an increase in the abundance of unphosphorylated Dvl2, suggesting that kinase activity of Stk11 is essential to its ability to restrain Dvl activity (Fig. 6E). These observations, together with the ability of Stk11 to reduce cell growth (either by stalling proliferation or by inducing cell death) in HeLa cells (Fig. 6F), support a role for Stk11 in controlling cell growth or viability by influencing Wnt/β-catenin pathway activity. This cell growth regulation is not likely to involve Hh signaling, given that HeLa cells rarely form primary cilia (36). Engagement of Dvl as a consequence of Stk11 loss in cancerous cells may thus give rise to β-catenin–dependent or –independent Wnt-mediated responses, or both. Furthermore, the type of response is likely dictated by the expression pattern of Wnt family members and their receptors (37).
We used an RNAi screening approach in murine cultured cells to establish previously unknown relationships between two therapeutically relevant signal transduction pathways and genes of importance in diseases of cell growth, tissue degeneration, and malformation. We investigated the mechanism of action for Stk11 in both the Hh and the Wnt pathways and biochemically established a role for Prkar1a, the regulatory subunit of PKA, in governing GliR formation. The use of two readouts (Hh and Wnt pathway reporters) to elucidate cellular mechanisms of Stk11 function suggests that expanded data sets encompassing analysis of additional pathways would improve our capability to assign gene function to cellular processes.
Proteolytic processing of Gli proteins entails their cycling through the primary cilium to acquire PKA-dependent phosphorylation (23–25). However, this process has been difficult to study, given that genetic approaches to disrupt intraciliary trafficking through targeted deletion of intraflagellar trafficking (IFT) complex components typically result in disruption of ciliogenesis as well as Gli processing (38–40). Our observations in STK11-null cells suggest that the length of the primary cilium controls the rate of Gli processing. The ability of the tubulin deacetylase inhibitors SAHA and tubacin to influence cilia length in STK11−/−, but not wild-type, MEFs argues that Stk11 plays an important regulatory role in limiting ciliary disassembly and that chemicals targeting HDAC6 may be useful in restoring some aspects of Stk11-induced pathology (Fig. 6G).
Our study suggests that sequestration of Dvl away from Wnt receptors represents an important checkpoint in Wnt pathway responses and that this checkpoint may be lost in Stk11-associated tumors. About two-thirds of cervical carcinomas exhibit increased abundance of Dvl protein (41), suggesting that increased Dvl activity may be a common defect in this cancer type. Given the role of Dvl in mediating a diverse array of Wnt-dependent responses, our findings reveal a mechanistic basis underlying deviant Wnt-dependent transcriptional response and possibly other Wnt-associated cellular derangements in cancerous cells lacking Stk11 (42). Thus, inhibitors of Wnt production, such as the Porcn inhibitors, may be useful in some Stk11-associated cancers, whereas inhibitors of canonical Wnt responses, such as those that induce Axin stability, may be additionally relevant in a subset of such cancers (Fig. 6G).
Given the tremendous resources devoted to achieving chemical control of the Hh and Wnt pathways for therapeutic goals, our screening strategy potentially broadens the number of diseases that may be managed with modulators of these pathways. We focused here on the role of Hh and Wnt signaling in Stk11-associated tumors; however, further studies are required to understand the nature of other disease-pathway relationships established in our screen. Our identification of chemical means to counter some of the effects resulting from loss of a multitasking gene product, such as Stk11, suggests that screening additional signal transduction pathways with this approach would further facilitate the deconstruction and, ultimately, management of complex diseases.
3T3-ShhFL cells used in the primary screen were generated by the stable transfection of NIH 3T3 cells expressing Shh protein (provided by P. A. Beachy) with a Hh-responsive reporter [8x-Gli-BS-Luc (43)], and a Renilla luciferase (RL) reporter, pRL-SV40 (Promega). ShhLightII cells (12), HeLa cells, NIH 3T3 cells, and DLD-1 cells were purchased from the American Type Culture Collection. Wild-type MEFs and STK11−/− MEFs were provided by N. Bardeesy, IFT88−/− MEFs by A. Liu, and PRKAR1A−/− MEFs by L. Kirschner. HBEC cells were provided by J. Shay. Lung cancer cell lines H1819, H460, A549, and H23 were provided by J. Minna. Endo and SiHa cells were provided by D. Castrillon.
The human STK11 (NM_000455) expression construct was purchased from Origene. KD STK11 (STK11 K78M) was engineered with polymerase chain reaction (PCR)–based cloning and mutagenesis strategies. The STF reporter plasmid was provided by R. T. Moon. The pRK5-ShhN and Wnt3A plasmids were provided by P. A. Beachy. Axin1-GFP and Dvl2-GFP plasmids were provided by M. Bienz.
SAG was purchased from Alexis Biochemicals, D4476 from Tocris, and Sant1 from Sigma. IWP and IWR compounds targeting the Wnt pathway components Porcupine and Tankyrase, respectively, were synthesized de novo as previously described (29). The HDAC inhibitor SAHAwas provided by C. Chen. Tubacin (a specific HDAC6 inhibitor) and MS-275 (HDAC1-3 inhibitor) were provided by R. Mazitschek and S. Schreiber. 5-(Tetradecyloxy)-2-furoic acid (TOFA) was purchased from Cayman Chemicals.
The antibody that recognizes Gli3 was obtained from Genentech (6F5), B. Wang, or R&D Systems. The antibody against Prkar1a was purchased from BD Biosciences. The antibody that recognizes Smo was provided by J. Kim and P. A. Beachy. Antibodies against acetylated tubulin, β-catenin, and β-actin were purchased from Sigma. Antibodies against Stk11, Dvl2, Axin1, Axin2, p-ERK, p-LRP6, LRP6, Dvl3, and α-tubulin were purchased from Cell Signaling Technologies. Antibodies against Dvl1 and Tnks 1/2 were purchased from Santa Cruz Biotechnology. The antibody that recognizes GFP was purchased from MBL.
Nonsilencing siRNAs (AllStars Negative Control siRNA) were purchased from Qiagen. Sequences of siRNAs used in secondary tests are in table S5.
Mouse genome-wide siRNA libraries, produced by Dharmacon and Qiagen, were screened by transfecting siRNA pools targeting a single gene into 3T3-ShhFL cells in 96-well format with the X-tremeGENE transfection reagent (Roche), according to the manufacturer’s protocol. The final siRNA concentration was 57 nM. Twenty-four hours after transfection, cells were placed under low-serum conditions (0.5% calf serum). Seventy-two hours after transfection, firefly luciferase (FL) and RL activities were assessed with the Dual-Luciferase kit (Promega). To account for edge effects and variation across plates, we normalized luciferase measurements to the average luciferase value for each well position across all plates and to the average luciferase measurement of each plate. Samples for which the triplicate data were inconsistent (SD >0.1) were removed from consideration. z scores were calculated by determining the number of SDs that the luciferase measurements for each siRNA pool fell away from the mean values for all the siRNA pools in each primary screen. siRNA pools that affected the RL control reporter (|z| >2.5) were removed from consideration. siRNA pools for which FL z ≤ −2.0 (97.73% confidence) were further evaluated by secondary screens.
The exogenous Shh secondary screen was performed essentially as above by the transfection of ShhLightII cells with siRNAs and, 24 hours later, the addition of ShhN-conditioned media, collected from human embryonic kidney (HEK) 293–ShhN cells. Average FL measurements ≤75% of FL for nonsilencing siRNA controls and SD <0.15 were considered for further testing. The Wnt counterscreen was performed by transfection of NIH 3T3 cells with the STF reporter, the RL reporter, the Wnt3A plasmid, and siRNAs using the Effectene transfection reagent (Qiagen), according to the manufacturer’s protocol. Final siRNA concentration was 25 nM. The low-throughput counterscreen (cross-library test) was performed by hand in 96-well format with only the center wells of the plate. The RNAi specificity test was performed in the same manner as the exogenous Shh secondary test with individual siRNAs (table S5).
Reverse transcription PCR (RT-PCR) analysis was performed as previously described (19). Primer design for RT-PCR was as follows: GapDH forward (F), ATCCTGCACCACCAACT; GapDH reverse (R), TGCCTGCTTCA-CCACCTT; Ptch1 F, ACTGTCCAGCTACCCCAATG; and Ptch1 R, CAT-CATGCCAAAGAGCTCAA. Quantitative PCR (qPCR) analysis primers were as follows: GapDH F, ACATCGCTCAGACACCATG; GapDH R, TGTAGTTGAGGTCAATGAAGGG; Gli3 F, GGGATTCCGACAGTT-CTGAACC; Gli3 R, CTGGGGAGGTCTTCATCAGGC; Ptch F, CGTCAGAAGATAGGAGAAGAGGC; and Ptch R, GTAGCACAAA-TGTTCCAACTTCC.
Gli3R formation was inhibited by addition of either Shh-conditioned medium or SAG (0.2 μM) for 24 hours. Conditioned medium or SAG was then removed and replaced with a medium containing Smo antagonist (SANT, 2.5 μM) to inhibit Shh pathway response and to allow reinitiation of Gli3R formation. Cellular lysate was then collected at indicated time points. To examine the rate of Gli3R turnover, we treated cells with MG132 to inhibit proteasome-dependent Gli processing. Cellular lysate was then collected at various time points. Samples were analyzed by Western blot, and Gli3R protein was quantified with ImageJ or Gene Tools software.
Wild-type zebrafish embryos were injected with 3 pmol of MO targeting STK11 (sequence: GAGATCCGCGCCCACGCTCATCTTT) or standard control (Gene Tools) at the one- to two-cell stage. Whole-mount in situ hybridization was performed at 24 hours after fertilization with digoxigenin-labeled antisense RNA probes generated against Eng1A (BC080209; base pairs 3 to 1433), Axin2 (BC045281; base pairs 341 to 1144), or MyoD (NM_131262; base pairs 316 to 927).
Cells were fixed in 2.5% glutaraldehyde in 0.1 M sodium cacodylate buffer, postfixed in buffered 1% osmium tetroxide, en bloc–stained in 2% uranyl acetate, dehydrated with a graded series of ethanol, and embedded in EMbed-812 resin. Thin sections were cut on a Leica Ultracut UCT ultramicrotome and then stained with 2% uranyl acetate and lead citrate. Images were acquired on an FEI Tecnai G2 Spirit electron microscope equipped with a LaB6 source and operating at 120 kV.
For studies involving transfection of siRNAs, DharmaFECT 3 (Dharmacon), Effectene (Qiagen), or X-tremeGENE (Roche) reagent was used per the manufacturer’s protocol. For studies involving transfection of complementary DNAs (cDNAs), FuGENE HD (Roche) or Effectene (Qiagen) was used according to the manufacturer’s instructions.
Cell lysates for Western blot analysis were generated with either phosphate-buffered saline (PBS)/1% NP-40 or RIPA (radioimmunoprecipitation assay) buffer (for detection of Gli proteins). Western blot results were quantified with ImageJ or Gene Tools software. For imaging of cilia, cells were plated at high density and low serum on polylysinecoated coverslips (BD Biosciences). In experiments with SAHA or IWP, cells were treated for 24 hours. For Smo localization experiments, cells were treated with ShhN-conditioned media for 4 hours before fixation. Cells were fixed in either 4% formaldehyde or ice-cold methanol and incubated with antibody against acetylated tubulin and then secondary antibodies (mouse fluorescein isothiocyanate–conjugated or Alexa Fluor 488–conjugated antibodies) diluted in PBS, 0.2% Triton X-100, and 5% goat serum. Cilia were imaged and measurements were taken with Adobe Photoshop.
Cancer cell growth studies were performed as previously described (29). Briefly, cells were seeded at clonal density (1500 cells in 24-well format) and cellular adenosine 5′-triphosphate (ATP) levels were measured 6 days later with CellTiter-Glo reagent (Promega).
Fig. S1. Additional material relating to the RNAi screen.
Fig. S2. Confirmatory siRNA assays for gene function in Hh pathway response.
Fig. S3. Additional characterization of STK11−/− cells.
Fig. S4. Hh-induced Smo or Gli protein accumulation in the primary cilium is not affected in STK11−/− cells.
Fig. S5. Dvl phosphorylation is sensitive to IWP but not SAHA.
Fig. S6. Dvl2 undergoes a change in distribution in response to Wnt protein stimulation in MEFs.
Fig. S7. The defect in posterior morphogenesis and somite development is apparent in STK11 morphants stained for MyoD expression, which labels somitic tissue.
Fig. S8. Loss of STK11 in cancerous cells induces Dvl protein phosphorylation and compromises primary cilia formation.
Fig. S9. Cell line–dependent responses to Wnt ligand in cervical carcinoma cells lacking Stk11.
Fig. S10. Axin1 forms punctate structures in HeLa cells.
Table S1. Dharmacon library primary screen results (Excel file).
Table S2. Qiagen library primary screening results (Excel file).
Table S3. Genes identified by two siRNA pools from the primary screen (Excel file).
Table S4. Counterscreen results (Excel file).
Table S5. siRNA sequences targeting disease genes (Excel file).
We thank N. Bardeesy, J. Kim, P. Beachy, A. Liu, L. Kirschner, D. Castrillon, R. Mazitschek, S. Schreiber, J. Shay, J. Minna, C. Chen, and S. Scales for reagents. We also thank B. Perkins, C. Gilpin, K. Luby-Phelps, T. Carroll, L. Huang, and U. Eskiocak for technical assistance. We thank J. K. Chen and colleagues for useful discussions. Finally, we thank members of the Lum lab and W. Snell for critical reading of this manuscript.
Funding: This work was supported by the NIH (1R01GM076398), the American Cancer Society (RSG GMC-07-062-01), the Welch Foundation (I-1665), and an endowment from Virginia Murchinson Linthicum to L.L.
Author contributions: L.S.J., X.W., M.E.D., C.-W.F., O.K., B.C., and W.T. performed all the experiments. L.S.J., X.W., M.E.D., C.-W.F., O.K., J.F.A., and L.L. designed the experiments and analyzed the data. B.W. provided reagents. L.S.J. and L.L. wrote the paper.
Competing interests: L.L., B.C., W.T., and M.E.D. are listed as inventors in a pending patent application associated with the IWR/IWP compounds described in this study.