Oncogenic RAS Generates Cancer Stem Cells in p53-Decient Fibroblasts Through SOX2 Induced By CDK1-Mediated Protein O-GlcNAcylation

Cancer stem cells (CSCs) have tumour initiation, self-renewal, and long-term tumour repopulation properties, and it is postulated that differentiated somatic cells can be reprogrammed to CSCs by oncogenic signals. We previously showed that oncogenic HRAS V12 conferred tumour initiation capacity in tumour suppressor p53-decient (p53 −/− ) primary mouse embryonic broblasts (MEFs) through transcription factor NF-κB-mediated enhancement of glucose uptake; however, the underlying mechanisms of RAS oncogene-induced CSC reprogramming have not been elucidated. Here, we found that the expression of the reprogramming factor SOX2 was induced by HRAS V12 in p53 −/− MEFs. Moreover, gene knockout studies revealed that SOX2 is an essential factor for the generation of CSCs by HRAS V12 in mouse and human broblasts. We demonstrated that HRAS V12 -induced cyclin-dependent kinase 1 (CDK1) activity and subsequent enhancement of protein O-GlcNAcylation were required for SOX2 induction and CSC generation in these broblasts and cancer cell lines containing RAS mutations. Moreover, the CDK inhibitor dinaciclib and O-GlcNAcylation inhibitor OSMI1 reduced the number of CSCs derived from these cells. Taken together, our results reveal a signalling pathway and mechanism for CSC generation by oncogenic RAS and suggest the possibility that this signalling pathway is a therapeutic target for CSCs. reprogramming into


Introduction
Cancers are derived from cancer stem cells (CSCs), also known as cancer-initiating cells, that form tumours in vivo 1 . CSCs comprise a subset of tumour cells with tumour formation capacity, self-renewal properties, and long-term tumour repopulation activity 2,3 . The molecular mechanisms that support the reprogramming events leading to the generation of CSCs are considered to be similar to those of induced pluripotent stem cells (iPSCs), including the same regulatory factors and control mechanisms 4 . Indeed, reprogramming factors, including SRY-box 2 (SOX2), octamer-binding transcription factor 4 (OCT4, also known as POU5F1), Krüppel-like factor 4 (KLF4), and MYC (collectively known as OKSM transcription factors), and chromatin modi ers are required for establishment of iPSCs are associated with cancer generation 4,5 . It has been shown that transient expression of OKSM factors induced cancer development in various tissues in vivo and was accompanied by global changes in DNA methylation status 6 . Alterations of genes involved in the regulation of cell proliferation, such as oncogenes and tumour suppressor genes, may lead to cancer development 7,8 . Therefore, it is possible that analyses of reprogramming factors and cell growth regulators will lead to the elucidation of the mechanism of CSC generation. However, cancer develops because of the acquisition of successive hallmark cancer capabilities in a multistep pathogenic process 7,8 , and it is di cult to analyse the role of reprogramming factors in CSC generation using established cancer cell lines.
P53 is the most frequently mutated gene in human cancer cells, and its gene product functions as a transcriptional activator of various target genes 9 . Accumulating evidence has shown that metabolic regulation by p53 is important for tumour suppression 10,11 . For example, p53 inhibits glycolysis, and the loss of p53 function changes the energy source from cellular respiration to glycolysis. This metabolic shift, known as the Warburg effect, favours tumour growth and is observed in many cancer cells 12 . In p53-de cient (p53 −/− ) mouse embryonic broblasts (MEFs), we previously found that oncogenic cell transformation and enhanced aerobic glycolysis were induced by an activating mutant of HRAS (HRAS V12 ) and were completely dependent on the transcription factor NF-κB 13 . Moreover, enhanced glycolysis-induced O-GlcNAcylation in HRAS V12 -expressing p53 −/− MEFs further activated NF-κB 14 , suggesting that the NF-κB/glycolysis activation loop in p53 −/− cells is involved in enhanced energy production and CSC generation. The loss of p53 markedly increases the e ciency of somatic cell reprogramming to iPSCs, suggesting the existence of a p53-mediated barrier system that inhibits this reprogramming 15−18 . In various stem cell populations, enhanced glucose uptake is commonly observed and is critical for the acquisition and maintenance of stemness 19 . Moreover, direct O-GlcNAcylation of OCT4 and SOX2 has been shown to regulate pluripotency and reprogramming of embryonic stem cells 20 .
In relation to these results, we found that IL-8 overexpression in colon and lung cancer cells enhanced O-GlcNAcylation, which was required for the production and maintenance of CSCs 21 . These results suggest that enhanced glycolysis and O-GlcNAcylation are involved in CSC generation and maintenance in p53deregulated cells.
Gain of function mutations in RAS genes are frequently found in human cancers and represent markers of poor prognosis in patients with certain cancers 22,23 . The RAS family in humans comprises three genes, which are HRAS, NRAS, and KRAS. The proteins function as small GTPases. The activation of RAS is tightly regulated by guanine nucleotide exchange factors, which promote GTP-bound active states, and GTP-bound RAS interacts with and activates effector enzymes, such as RAF kinases. In the RAS-RAF cascade, activated RAF kinase phosphorylates and activates mitogen-activated protein kinase kinases 1 and 2 (MEK1 and MEK2), which leads to phosphorylation and activation of extracellular signal-regulated kinase (ERK) 1 and 2, also known as mitogen-activated protein kinase (MAPK) 3 and 1, respectively. The RAS-RAF-MAPK pathway mediates various biological functions, such as cell proliferation, survival, and differentiation, through a large number of ERK1/2 substrates that include transcription factors and cell cycle regulators 24 . Oncogenic mutations in the RAS-RAF-MAPK pathway contribute to various malignant phenotypes that include invasion, metastasis, relapse, and angiogenesis 5 .
Studies have shown that RAS-MAPK signalling promotes stemness and maintenance of CSCs 25 , but the underlying mechanisms have not been completely elucidated. Based on these studies, we analysed the expression of reprogramming factors in oncogenic HRAS-expressing p53 −/− MEFs as a simple model system to understand the molecular mechanism of RAS-induced CSC generation. As a result, we found that SOX2 induction by CDK1-mediated protein O-GlcNAcylation was essential for generation and maintenance of CSCs in p53-de cient cells.

Results
Oncogenic RAS induces SOX2-initiated CSC generation. Previously, we reported that oncogenic HRAS (HRAS G12V mutant: HRAS V12 ) induced tumorigenic properties in p53 −/− MEFs 26 , suggesting that HRAS V12 promoted the generation of CSCs. To con rm this hypothesis, we generated HRAS V12 -expressing p53 −/− MEFs (Fig. 1A) and demonstrated that these cells formed tumours in nude mice (Fig. 1B). Previous studies have shown that CSCs display tumorigenic activity, and cancer cells with CSC-like properties will form spheres in low attachment culture conditions in medium containing growth factors 27,28 . As shown in Fig. 1C, approximately 0.4% of HRAS V12 -expressing p53 −/− MEFs developed spheres in low attachment plates. It has been reported that the reprogramming factors OCT4, KLF4, and SOX2 are important for generation and maintenance of CSCs and are used as CSC markers in various cancers [29][30][31] . Therefore, we analysed the mRNA expression levels of these factors and found that the sphere-forming cells exhibited enhanced expression of OCT4, KLF4, and SOX2 compared with that of parental adherent cells (Fig. 1D). These results strongly suggested that sphere-forming cells from HRAS V12 -expressing p53 −/ − MEFs were CSCs. We also compared expression levels of these genes between adherent HRAS V12expressing p53 −/− MEFs cells and p53 −/− MEFs control cells. Interestingly, although OCT4 expression was unchanged, the expression of SOX2 was markedly enhanced in HRAS V12 -expressing p53 −/− MEFs cells ( Fig. 1E and F). These results suggest that the expression of SOX2 as well as KLF4 was induced by HRAS V12 expression in p53 −/− MEFs, and OCT4 was induced in the process of CSC reprogramming.
It has been shown that tumorigenic properties are e ciently induced by co-expression of oncogenes in primary rodent cells but not in human cells, suggesting a mechanistic difference between human and rodent oncogenesis 32 . As previously shown, HRAS V12 and simian virus 40 large T antigen (SV40 LT), which functionally inactivates p53 and the retinoblastoma tumour suppressor protein, were unable to induce tumorigenic properties in human cells 32 ; SOX2 expression and sphere-forming activity were not induced by p53 knockout (KO) or HRAS V12 expression in the human lung broblast cell line TIG-3 ( Supplementary Fig. S1 online). However, TIG-3 cells were transformed by the combination of SV40 T-ag, c-MYC, and HRAS V12 (TIG-3-SMR cells) 33,34 , and sphere-forming cells were present in the TIG-3-SMR cell population ( Supplementary Fig. S1C online). Notably, SOX2 expression was also elevated in TIG-3-SMR cells when compared with control TIG-3 cells ( Supplementary Fig. S1D online). These results suggest that the induction of SOX2 by oncogenic signalling is regulated differently between mice and humans. However, the difference in tissue origin cannot be ruled out.
SOX2 expression is essential for HRAS V12 -induced development of CSC properties. Studies have shown that SOX2 functions as an oncogene, and expression is closely related to metastasis and relapse of many cancers 35,36 . Therefore, from these studies and the above results, it was possible that SOX2 functioned as an initiation factor for CSC reprogramming. To test this hypothesis, we deleted the SOX2 gene in HRAS V12 -expressing p53 −/− MEFs with three independent gRNAs using the CRISPR-Cas9 gene knockout system. Single cell-derived SOX2 KO clones from each gRNA were isolated, and immunoblotting and quantitative PCR (qPCR) showed a lack of SOX2 expression in each clone ( Supplementary Fig. S2A, B online). Inhibition of SOX2 expression was reported to attenuate cancer cell proliferation 37 . SOX2 KO cells showed suppressed cell growth properties; however, proliferation was not completely inhibited ( Fig. 2A).
Next, we analysed the effect of SOX2 KO on the CSC properties found in HRAS V12 -expressing p53 −/ − MEFs. The appearance of sphere-forming cells and an increase in anchorage-independent colonyforming cells correlate with tumourigenicity 38 ; however, these properties were not observed after SOX2 KO ( Fig. 2B and Supplementary Fig. S2C online). Furthermore, we analysed the effect of SOX2 KO on tumourinitiating activity in vivo. HRAS V12 -expressing p53 −/− MEFs with or without SOX2 KO were subcutaneously injected into nude mice and tumour growth was monitored. Consistent with data presented in Fig. 1B, we con rmed tumour progression in mice injected with the HRAS V12 -expressing p53 −/− MEFs (Fig. 2C). However, SOX2 KO in these cells completely inhibited tumour development (Fig. 2C, D). Because the cell growth rate was suppressed in vitro by SOX2 KO ( Fig. 2A), we measured tumorigenic activity for an extended period; however, we did not detect tumour formation up to 13 weeks after injection (Fig. 2C).
Moreover, forced expression of SOX2 in p53 −/− MEFs and p53 −/− TIG-3 cells promoted sphere formation and tumour development (Supplementary Fig. S2D-H online). These results suggest that SOX2 expression is required for CSC generation in p53-de cient cells. SOX2 expression is induced by the RAF-MAPK pathway in HRAS V12 -expressing p53 −/− MEFs. RAS proteins directly activate the downstream effectors RAF and PI3K followed by the downstream activation of MAPK and AKT pathways, respectively 24,39 . Therefore, we next examined whether these effectors were involved in promoting SOX2 expression. Constitutively active forms of RAF (BRAF V600E ) and PI3K (PI3K CAAX ) were stably expressed in p53 −/− MEFs. Although the increase in SOX2 mRNA expression was relatively weakly in PI3K CAAX -expressing cells, the levels of SOX2 and KLF4 mRNA in BRAF V600Eexpressing cells were similar to that in HRAS V12 -expressing cells (Fig. 3A). These results suggest that HRAS V12 induces the expression of SOX2 mRNA through RAF and its downstream effectors MEK/ERK. Furthermore, the expression of HRAS V12 enhanced the activating phosphorylation of ERK and both the mRNA and protein expression of SOX2. These enhancements were signi cantly reduced by the MEK inhibitor U0126 but not by the PI3K inhibitor LY294002 (Fig. 3B-D). Similar results were obtained after treatment with the MEK inhibitor PD184352 (Fig. 3E). Immuno uorescent staining con rmed that the number of SOX2-positive cells was reduced by the MEK inhibitor but not by the PI3K inhibitor (Fig. 3F).
Additionally, the number of sphere-forming cells was reduced after treatment of HRAS V12 -expressing p53 −/− MEFs with the MEK inhibitors (Fig. 3G). These results suggest that the RAF/MEK/ERK pathway is required for SOX2 induction and reprogramming of normal broblasts to CSCs. CDK1 activity is required for SOX2 induction and generation of CSCs in RAS-activated cells. Activation of ERK enhances expression of cyclin D1, which binds to and activates cyclin-dependent kinase 4 and 6 (CDK4/6) during the G1 phase of the cell cycle 40 . Aberrant regulation of the CDK4/6-cyclin D1 axis is associated with the development of metastatic melanoma and breast cancer 41,42 . From these ndings, it is possible that HRAS V12 induced SOX2 expression through ERK-mediated CDK4/6 activation. To test this possibility, we analysed the effect of the CDK4/6 inhibitor palbociclib 43 on SOX2 expression and con rmed that the expression of SOX2 mRNA and protein levels were slightly reduced in HRAS V12 -expressing p53 −/− MEFs after treatment with a high-dose of palbociclib ( Supplementary Fig. S3A, B online). Similar results were obtained with TIG-3-SMR cells ( Supplementary Fig. S3C online). These results suggest that CDK4/6 activity alone is insu cient for SOX2 induction. In addition to CDK4/6, the BRAF/MEK/MAPK pathway is also essential for some functions of CDK1 and CDK2 [44][45][46] . Recent studies have reported a role for CDK1 and CDK2 in the maintenance of stem cell pluripotency and CSC-like properties in breast cancer [47][48][49] . These reports suggested that CDK1 and/or CDK2 are required for SOX2 induction. Therefore, we examined the effects of the CDK1/2 inhibitors dinaciclib 50 and roscovitine 51 on SOX2 expression. We found that SOX2 mRNA and protein expression was strongly suppressed by dinaciclib (Fig. 4A, B) and roscovitine ( Supplementary Fig. S3D, E online) in HRAS V12 -expressing p53 −/− MEFs. Similar results for dinaciclib were observed in TIG-3-SMR cells (Fig. 4C). Moreover, the number of sphere-forming cells in HRAS V12 -expressing p53 −/− MEFs was largely suppressed by dinaciclib but only partially suppressed by palbociclib (Fig. 4D), which correlated with suppression of SOX2 expression levels ( Fig. 4A and Supplementary Fig. S3A online). To further analyse whether CDKs were essential for SOX2 expression, we analysed SOX2 levels in HRAS V12 -expressing p53 −/− MEFs after knockdown of CDK1 or CDK2 using siRNA. Immunoblot analysis showed that CDK1, but not CDK2, knockdown reduced the expression level of SOX2 protein (Fig. 4E). Treatment with dinaciclib, U0126, and PD184352, which suppressed the expression of SOX2 in this study, reduced the phosphorylation of CDK1 T161 , an indicator of CDK1 activation 52 ( Supplementary Fig. S4A-C online). These results indicate that CDK1 activity is required for SOX2 expression in HRAS V12 -expressing p53 −/− MEFs. Additionally, we analysed the effect of dinaciclib on human cancer cell lines, including KRAS-mutated colon cancer (HCT116, SW480, and DLD1) and lung cancer (H460 and A549) cells. Dinaciclib suppressed SOX2 mRNA expression in these cancer cells (Fig. 4F) and the number of sphere-forming cells from the HCT116, SW480, and H460 lines (Fig. 4G). These results suggest that CDK1-mediated SOX2-induction promotes the generation of CSCs in human cancer cells with RAS mutations. Therefore, we determined the O-GlcNAc levels in HRAS V12 -expressing p53 −/− MEFs and TIG-3-SMR cells and found elevated levels of protein O-GlcNAcylation compared with those in the respective control cells (Fig. 5A, B). Next, we analysed the cells after treatment with OSMI1, a cell-permeable, small molecule OGT inhibitor 37 . OSMI1 treatment suppressed total O-GlcNAcylation levels of proteins and SOX2 expression in these cells (Fig. 5A, B). In contrast, treatment of HRAS V12 -expressing p53 −/− MEFs with thiamet G, a speci c OGA inhibitor that increases O-GlcNAcylation 38 , enhanced total O-GlcNAcylation levels and SOX2 expression (Fig. 5C). Interestingly, the OSMI1-mediated reduction in SOX2 levels was attenuated by treatment with the proteasome inhibitor MG-132 (Fig. 5A, B). The mRNA levels of SOX2 in HRAS V12expressing p53 −/− MEFs were not signi cantly changed by OSMI1 and thiamet G (Fig. 5D), indicating that O-GlcNAcylation regulated SOX2 expression at the post-transcriptional level. Consistent with these results, the numbers of sphere-forming cells decreased and increased after treatment with OSMI1 and thiamet G, respectively (Fig. 5E). These results suggest that increased O-GlcNAcylation is required for SOX2 protein expression and sphere-forming activity in these cells.  Fig. S4C online). Furthermore, knockdown of CDK1, but not CDK2, with siRNA inhibited O-GlcNAcylation levels in these cells ( Supplementary Fig. S4D online). In contrast, treatment with palbociclib had no effect on O-GlcNAcylation levels ( Supplementary Fig. S3B, C online). In addition, SOX2 expression (as shown in Fig. 4F, G) and O-GlcNAcylation levels in KRASactivated cancer cells were suppressed by dinaciclib (Fig. 5F). These results suggest that RAS/RAF/MAPK pathway-induced CDK1 activation is important for induction of O-GlcNAcylation, and this activation pathway is required for SOX2 expression and subsequent CSC generation.

Discussion
Accumulating evidence has revealed that serial oncogenic mutations in stem cells and even in differentiated somatic cells may induce the generation of CSCs [54][55][56][57][58] . The dedifferentiation of somatic cells into CSCs is thought to be induced by a reprogramming mechanism similar to that observed in the production of iPSCs 59,60 . Viewed through the "hallmarks of cancer" 7,8 , CSCs may be generated by iPSC reprogramming factors that are induced by gain-of-function oncogenic and loss-of-function tumour suppressor gene mutations. Therefore, the HRAS V12 -expressing p53 −/− MEFs provided a simple and easyto-analyse model for evaluating the signalling system used by oncogenes to regulate reprogramming factors. In our previous studies, we found that the HRAS V12 mutation conferred tumour-initiating activity, a hallmark of CSCs, in p53 −/− MEFs 26 , and this phenomenon was completely dependent on NF-κBinduced aerobic glycolysis 13 . In this study, we found that HRAS V12 induced MAPK-CDK1 signal-dependent induction of SOX2 mRNA transcription and O-GlcNAcylation-mediated SOX2 protein accumulation. Although the transcriptional induction mechanism of SOX2 by CDK1 was not elucidated in this experimental system, CDK1-mediated induction of SOX2 has been analysed in other studies. For example, it was reported that CDK1-induced phosphorylation directly activated transcription factor CP2-like protein 1 (TFCP2L1), which is an activator of pluripotency-associated genes, including SOX2 61 . CDK1 phosphorylated and inhibited the histone lysine demethylase KDM5B, which is a transcriptional suppressor of the pluripotency genes SOX2 and NANOG 62 .
It has been shown that the hexosamine biosynthetic pathway shunts glycolysis toward the production of a key substrate for O-GlcNAcylation and is activated by enhanced glycolysis 63 . The RAS/MAPK pathway induces metabolic reprogramming, including enhanced glycolysis 64 . Therefore, in addition to NF-κBinduced enhancement of glycolysis in p53 −/− MEFs 13 , our present results suggest that the RAS/MAPK/CDK1 pathway further promoted glycolysis and resulted in enhanced protein O-GlcNAcylation. Moreover, we also found that O-GlcNAcylation inhibited proteasomal degradation of SOX2, suggesting that direct O-GlcNAcylation of SOX2 or other factor(s) is involved in the ubiquitinmediated degradation system. However, in our previous report and in the present study, we were unable to detect direct O-GlcNAcylation of SOX2 as demonstrated in pluripotent cells 65 . Although SOX2 expression was induced by mRNA levels in transformed cells (Fig. 1E), we found that O-GlcNAc modi cation regulated SOX2 at post-transcriptional levels ( Fig. 5A-D). For this issue, we indicated that deletion of SOX2 gene in HRAS V12 -expressing p53 −/− MEFs reduced O-GlcNAc levels ( Supplementary Fig. S4E online), and exogenous SOX2 expression induced protein O-GlcNAcylation in p53 −/− MEFs ( Supplementary Fig. S4F online). These data suggest that transcriptionally elevated SOX2 expression promotes O-GlcNAc levels, and enhanced O-GlcNAcylation could promote expression of SOX2 at posttranscriptional levels.
The tumour suppressor protein p53 is referred to as "the guardian of the genome" because it continuously surveys damaged genomic DNA and facilitates DNA repair 66 . Moreover, in response to oncogenic signalling, the p53 tumour surveillance system eliminates cells through p53-mediated apoptosis or induction of senescence 67 . In addition to these functions, accumulating evidence has indicated that regulation of cellular metabolism by p53 is involved in tumour suppression 68 . It was demonstrated that mice expressing acetylation-defective p53 mutants lost the ability to induce apoptosis and senescence but retained tumour suppressive function and the ability to modulate the expression of metabolic genes 69 . In p53 −/− MEFs, we previously found that oncogenic transformation by HRAS V12 was dependent on enhanced aerobic glycolysis through NF-κB-mediated induction of the glucose transporter GLUT3 13 . In this previous study, we found that GLUT3 (also known as SLC2A3) functioned as an oncogene because p53 −/− MEFs underwent oncogenic transformation after overexpression of GLUT3 without HRAS V12 . These results suggested that enforced glucose ux induced CSC generation in cells that lacked p53 function. Furthermore, we found that O-GlcNAcylation increased with accelerated glucose consumption in p53 −/− MEFs and was further enhanced by HRAS V12 expression 14 . In the present study, CSC generation in HRAS V12 -expressing p53 −/− MEFs was completely dependent on the induction of SOX2, and enhancement of O-GlcNAcylation by the HRAS V12 -MAPK-CDK1 pathway was important for SOX2 induction. Therefore, it is possible that glycolysis may be hampered by p53, which prevents excessive O-GlcNAcylation-mediated induction of SOX2. Indeed, it has been shown that oncogenes, such as HRAS V12 , activate p53 in untransformed cells 67 , and activation of p53 results in inhibition of glucose consumption and glycolysis 13,70 . Moreover, it has been reported that the OKSM reprogramming factors induce p53mediated senescence 71 , suggesting the possibility that these reprogramming factors limit their own function using p53-mediated negative feedback. Therefore, at least in MEFs, p53 may exert a tumour suppressor function through the regulation of O-GlcNAcylation by limiting glucose metabolism.
In MEFs, induction of SOX2 was essential for CSC reprogramming by HRAS V12 , but this pathway may not be the only signal for CSC generation. For example, in the context of in ammation-induced cancer, we found that the activated form of MYD88, a regulator of in ammatory signalling, promoted CSC generation in p53 −/− MEFs through activation of the NF-κB-HIF1-OCT4, but not the SOX2, induction pathway 72 . Furthermore, in the present study, expression of HRAS V12 alone was insu cient to induce SOX2 in human TIG-3 cells. Although it is possible that this difference was because of cell type speci city, differences between human and rodent cells cannot be ruled out. In contrast to human cells, primary rodent cells are e ciently converted to tumorigenic cells by the co-expression of oncogenes, suggesting a fundamental difference between humans and rodents 32 . Therefore, these results suggest the possibility that an unidenti ed tumour suppressive mechanism other than p53 exists in humans and prevents the induction of SOX2 by RAS. Further analyses are required to clarify this mechanism.

Immunoblot analysis
Cells were lysed with TNE buffer (10 mM Tris-HCl, pH 7.4; 1% NP-40; 150 mM NaCl; 1 mM EDTA; 1 mM DTT, and protease inhibitor cocktail) (Nacalai Tesque, Japan). Protein concentrations were measured using the Bradford assay. Cellular proteins (20 µg) were separated using sodium dodecyl sulfatepolyacrylamide gel electrophoresis and then transferred to polyvinylidene di uoride membranes (Merck). The membranes were probed with primary antibodies, followed by incubation with horseradish peroxidase-conjugated mouse or rabbit immunoglobulin G (GE Healthcare, England) and visualisation using Chemi-Lumi-One Super or Ultra assay kits (Nacalai Tesque). The protein bands were digitalised using the LAS-3000 mini image analyser (Fuji lm, Japan).

Quantitative real-time PCR
Total RNA was extracted using the NucleoSpin RNA kit (Macherey-Nagel, Germany) following the manufacturer's instructions. Double-stranded cDNA was prepared from total RNA using oligonucleotide (dT), random primers, and Superscript III reverse transcriptase (Invitrogen, Carlsbad, CA, USA).

Immuno uorescent analysis
The p53 −/− MEFs and HRAS V12 -expressing p53 −/− MEFs (4 × 10 4 ) were cultured on glass coverslips in 6well plates and treated with U0126 (10 µM) or LY294002 (20 µM) for 24 h. The cells were washed in phosphate-buffered saline (PBS), xed with 4% paraformaldehyde in PBS for 15 min, and permeabilised with 0.5% Triton X-100 in PBS for 10 min. Permeabilised cells were blocked with Blocking One Histo reagent (Nacalai Tesque) for 30 min and incubated with anti-SOX2 antibody for 1 h at room temperature. Alexa Flour 488-conjugated anti-rabbit antibody (Thermo Fisher Scienti c) was used as the secondary antibody for 1h at room temperature. VECTASHIELD Mounting Medium with DAPI (Vector Laboratories, Burlingame, CA, USA) was used to stain nuclei and for mounting the cells on slides. Images were acquired using a uorescence microscope (BioZero BZ-8100; Keyence, Japan) and analysed with BZ Analyzer software (Keyence).
Cell growth analysis HRAS V12 -expressing p53 −/− MEFs (1 × 10 5 ) expressing each SOX2 sgRNA separately were seeded in 6well plates. Cell numbers were determined using a Vi-CELL cell analyser (Beckman, Brea, CA, USA) on the indicated days after plating, and cell growth curves were created.
Sphere formation assay Cells (5 × 10 3 or 1 × 10 4 ) were plated in 6-well, ultra-low attachment plates and grown in serum-free DMEM/F12 medium containing epidermal growth factor (20 ng/mL) and basic broblast growth factor (10 ng/mL) for 7 days. The numbers of spheres were counted in each treatment group. Sphere images were captured using the BZ-8100 microscope.

Colony formation assay
A layer of 1.5% (weight/volume) agarose prepared in DMEM containing 10% FBS was added to the wells of 6-well plates. Agarose (0.6%) containing 3 × 10 4 HRAS V12 -expressing p53 −/− MEFs containing each SOX2 sgRNA was added to the top of the rst layer. After 30 days, each well was stained with 0.005% crystal violet (Sigma-Aldrich), and the colonies were counted.

Animal experiments and cell line xenografts
The animal experiment protocol was approved by the Ethics Committee on Animal Experiments of Nippon Medical School (ethics approval number 26-020, . Animal experiments were carried out in accordance with the guidelines for Animal Experiments of Nippon Medical School and the guidelines of The Law and Noti cation of the Government of Japan as well as the ARRIVE guidelines. Mice were maintained at 20-24 °C in a facility with a 12 h light/ dark cycle and 40%-70% humidity. The mice were allowed free access to water and standard MF laboratory mouse chow (Oriental Yeast Co., ltd. Tokyo, Japan) and housed at a maximum number of ve per cage. All mice were checked for stress each day. For the xenograft experiments, 5-week-old male BALB/cAJcl-nu/nu mice were purchased from CLEA Japan, Inc. (Tokyo, Japan) and assigned at random to the experiments. These mice were subcutaneously injected with 1 × 10 5 cells of p53 −/− MEFs, HRAS V12 -expressing p53 −/− MEFs or HRAS V12 -expressing p53 −/ − MEFs containing SOX2 gRNA. The number of mice used are indicated for each experiment. Tumour growth was monitored every 3 days for 3 weeks and each week thereafter by caliper measurements, and tumour size was determined using the following formula: (Length × Width 2 )/2. At the end of the experiments (3 or 13 weeks post-injection), mice were euthanized by cervical dislocation, then each tumour were removed and weighed, and also collected for further experiments.

Statistical analysis
All experiments were repeated at least three times independently. Data are presented as means ± standard deviation (SD). Statistical analysis was performed with the Student's t-test using Microsoft Excel (Microsoft, DC, USA). P < 0.05 was considered statistically signi cant. Figure 1 Oncogenic RAS induces SOX2-initiated CSC generation. A, Vector only or the HRASV12 mutant was transduced into p53−/− mouse embryonic broblasts (MEFs) by retroviral infection for 2 days, and cells were selected with hygromycin for 3 days. Expression of the HRASV12 mutant was con rmed by immunoblotting analysis. B, p53−/− MEFs (CTL) or HRASV12-expressing p53−/− MEFs (1 × 105) were subcutaneously injected into immunode cient mice (n = 5 per group). Tumour weights were measured 3 weeks after injection. **P < 0.01.