Abstract
Oncogenic KRAS is the key driver oncogene for several of the most aggressive human cancers. One key feature of oncogenic KRAS expression is an early increase in cellular reactive oxygen species (ROS) which promotes cellular transformation if cells manage to escape cell death, mechanisms of which remain incompletely understood. Here, we identify that expression of oncogenic as compared to WT KRAS in isogenic cellular systems renders cells more resistant to ferroptosis, a recently described type of regulated necrosis. Mechanistically, we find that cells with mutant KRAS show a specific lack of ferroptosis-induced lipid peroxidation. Interestingly, KRAS-mutant cells upregulate expression of ferroptosis suppressor protein 1 (FSP1). Indeed, elevated levels of FSP1 in KRAS-mutant cells are responsible for mediating ferroptosis resistance and FSP1 is upregulated as a consequence of MAPK and NRF2 pathway activation downstream of KRAS. Strikingly, FSP1 activity promotes cellular transformation in soft agar and its overexpression is sufficient to promote spheroid growth in 3D in KRAS WT cells. Moreover, FSP1 expression and its activity in ferroptosis inhibition accelerates tumor onset of KRAS WT cells in the absence of oncogenic KRAS in vivo. Consequently, we find that pharmacological induction of ferroptosis in pancreatic organoids derived from the LsL-KRASG12D expressing mouse model is only effective in combination with FSP1 inhibition. Lastly, FSP1 is upregulated in non-small cell lung cancer (NSCLC), colorectal cancer (CRC) and pancreatic ductal adenocarcinoma (PDAC) as compared to the respective normal tissue of origin and correlates with NRF2 expression in PDAC patient datasets. Based on these data, we propose that KRAS-mutant cells must navigate a ferroptosis checkpoint by upregulating FSP1 during tumor establishment. Consequently, ferroptosis-inducing therapy should be combined with FSP1 inhibitors for efficient therapy of KRAS-mutant cancers.
Similar content being viewed by others
Introduction
The Ras proto-oncogenes (HRAS, NRAS and KRAS) are amongst the most frequently mutated genes across human cancers [1,2,3]. KRAS in particular is mutated in lung and pancreatic ductal adenocarcinoma (PDAC) as well as colorectal cancer. Point mutations within KRAS favor its active, GTP-bound state [2, 3]. Thereby, oncogenic forms of KRAS constitutively signal through the mitogen-activated protein (MAPK) pathway, PI3K and Rac1 signaling pathways endowing them with a variety of advantages including evasion of extrinsic apoptosis [4, 5]. It is established that KRAS-mediated cellular transformation requires the generation of reactive oxygen species (ROS) through elevated expression of NADPH oxidase 1 (Nox1) [6]. Yet, it is poorly understood how cells expressing mutated KRAS can mitigate the problem of ROS-induced cell death. We recently showed that cells expressing oncogenic KRAS upregulate the cystine/glutamate antiporter xCT (SLC7A11) upon hydrogen peroxide stimulation to promote cellular transformation [7]. Interestingly, xCT has been shown to protect cells from ferroptosis, an iron-dependent type of regulated necrosis executed by the accumulation of lipid ROS [8]. Cells are protected from ferroptosis by glutathione peroxidase 4 (GPX4) [9] which depends on glutathione (GSH) as an electron donor to reduce lipid hydroperoxides. GSH synthesis is coupled to the availability of intracellular cysteine which can be generated from cystine imported via xCT [10]. In addition, recent studies indicate that the Coenzyme Q10 (COQ10) oxidoreductase ferroptosis suppressor protein 1 (FSP1, formerly AIFM2) protects cells from ferroptosis through the generation of the lipid radical-trapping agent ubiquinol [11, 12]. Although ferroptosis has been suggested to represent a vulnerability in HRAS-mutant cells [13], GPX4 deletion in pancreatic intraepithelial neoplasia (PanINs) in genetically engineered mouse models driven by KRASG12D did not effectively kill PanINs [14] strongly suggesting KRASG12D-driven PanINs to be protected from ferroptosis through unknown mechanisms.
Here, through the use of various independent isogenic cellular models expressing near-endogenous levels of oncogenic or wild type (WT) KRAS, we demonstrate that oncogenic forms of KRAS render cells more resistant to ferroptosis through NRF2-mediated FSP1 upregulation in vitro and during tumor initiation in vivo.
Materials and methods
Cell lines
The panel of âRaslessâ mouse embryonic fibroblasts (MEFs) reconstituted with various oncogenic KRAS mutations (RPZ26216, RPZ25854, RPZ26198, RPZ26186, RPZ26425, RPZ26299, RPZ26295) was generated and kindly provided by the RAS Initiative at the Frederick National Laboratory for Cancer Research (FNLCR), US. Independently, Rasless MEFs were also obtained from M. Barbacid to generate KRAS point mutants from bulk sorting without deletion of the endogenous floxed KRAS allele. All MEFs were grown in Dulbeccoâs modified Eagleâs (DMEM)â+âGlutaMAX⢠medium (Gibco) with 4âµg/ml of blasticidin. NIH-3T3 cells expressing KRASG12V were generated and described previously [7]. NIH-3T3 stably expressing 4OHT-inducible HRASG12V and freshly isolated KRASG12D-inducible mouse embryonic fibroblasts (MEFs) were cultured in DMEM (Gibco) supplied with 1% L-Glutamine (Sigma) and 1% Sodium Pyruvate (Sigma). Human pancreatic duct epithelial cells (HPDE) were kindly provided by A. Trauzold (University of Kiel) and cultured in 75% RPMI 1640/ medium in presence of 25% keratinocyte growth medium 2 (PromoCell). The human non-small cell lung cancer (NSCLC) cell line A549 and mouse Lewis lung carcinoma cell line (3LL) were kindly provided by Prof. Julian Downward and cultured in RPMI 1640 medium (Gibco). HEK-293T cells were cultured in DMEM medium (Gibco). All media were supplemented with 10% fetal calf serum (FCS) (Sigma Aldrich) and 1000âU/mL of both penicillin and streptomycin (Pen/Strep) (Sigma Aldrich). All cells were kept at 37â°C with 5% CO2 and tested for mycoplasma at regular intervals (mycoplasma barcodes, Eurofins Genomics).
Reagents
Blasticidin (AppliChem GmbH), RSL3 (Selleckchem), ML210 (Tocris), ML162 (Caymann), Erastin (Biomol), Sulfasalazine (SAS) (MedChemExpress), Imidazole Ketone Erastin (IKE) (Sellekchem), Ferrostatin-1 (Sigma Aldrich), Liproxstatin-1 (Biozol), Necrostatin-1s (Abcam), zVAD (ENZO), iFSP1 (Cayman Chemicals), Tert-butylhydroquinone (TBHQ) (Sellekchem), AMG510 (MedChemExpress), ARS1620 (Chemgood), PD184352 (Sigma Aldrich), MK2206 (Sellekchem), DRAQ7 (Biolegend), BODIPY C11 (Invitrogen), H2DCFDA (Invitrogen), Dharmafect I (Dharmacon), Puromycin (Sigma), Doxycycline hydrochloride (Alfa Aesar), 4-hydroxy-tamoxifen (4OHT) (Sigma), Polybrene (Merck), CaCl2 (Sigma Aldrich), HBS (Sigma Aldrich), propidium iodide (Sigma).
Antibodies
Ras (clone RAS10, #05-516; Millipore, 1:1000), GPX4 (Abcam, ab41787, 1:2,000), xCT (Abcam, ab37185, 1:2000), Ã-Actin (Sigma, A1978, 1:10,000), GAPDH (Cell Signaling, #97166, 1:2000), FSP1 (previously described [11], kindly provided by M. Conrad, undiluted hybridoma supernatant), p44/42 MAPK (Erk1/2) (Cell Signaling, #9102 1:1000), phospho-p44/42 MAPK (Erk1/2) (Thr202/Tyr204) (Cell Signaling, #4370, 1:1000), HRP-conjugated secondary antibodies: goat-anti-mouse-HRP (Linaris GmBH, 20400-1âmg, 1:10,000), goat-anti-rabbit-HRP (Linaris GmBH, 20402-1âmg, 1:10,000), goat-anti-rat-HRP (Sigma, A9037-1âml, 1:10,000).
Plasmids
The packaging plasmids pCMV-VSV-G (#8454), pCMV-VSV-G (#8454), pMDLg/pRRE (#12251) were obtained from Addgene, P442-empty vector and P442-PLI-AIFM2-WT was kindly provided by J. P. Friedmann-Angeli, pLKO.1-empty vector and pLKO.1-shFSP1 were purchased from Merck (NM_153779/TRCN0000112139/pLKO.1). pCW-Puro-KRASG12D to generate doxycycline KRASG12D-inducible HPDE cells was cloned from pCW (addgene #50661 [15]) by replacing the existing Cas9 gene by human KRASG12D cDNA.
siRNA transfections
Two hundred microliters Opti-MEM (Gibco) and 1.5âµL Dharmafect Reagent I (Dharmacon) were mixed and incubated for 5â10âmin at room temperature. 2.2âµL of siRNA (Stock 20âmM) (Dharmacon) were added to the mixture and incubated for another 30âmin at room temperature. After incubation, 200âµL of the mixture was added to each well (6-well) plate and cells were plated on top. Knockdowns were incubated for 48â72âh, as indicated.
Cell viability assays
For this assay, 5000 or 10,000 (for iFSP1â±âRSL3 viability assays) cells were plated per 96-well plate 24âh before treatment. Cell viability was determined by Cell Titer Blue assay according to the manufacturerâs instructions (Promega).
Cell death assays (flow cytometry)
One day before treatment 45,000 cells (clonal and bulk-sorted MEFs), 55,000 cells (cells expressing either P442-empty vector or P442-PLI-AIFM2-WT) or 50,000 cells (3LLs) were plated in each well of a 24-well plate. For FSP1 siRNA knockdown, 40,000 cells were seeded 48âh before treatment. To determine cell death, adherent and detached cells were harvested and stained with propidium iodide (PI) (1âµg/ml) (Sigma Aldrich) in PBS (Thermo Fisher) supplemented with 2% FBS. PI-positive cells were quantified by flow cytometry using an LSR-FACS Fortessa (BD Bioscience) and FlowJo software (BD Bioscience).
Live cell imaging (IncuCyte)
Five thousand, 7500 or 10,000 cells per 96-well plate, 55,000 cells per 24-well plate or 300,000 cells per 6-well plate were seeded 24âh in advance, respectively. For KEAP1 siRNA knockdown, 20,000 cells were seeded in a 24-well plate on top of the transfection mix and incubated for 48âh followed by treatments for 24âh. Upon treatment, (Ferrostatin-1 [1 or 5âµM], RSL3 [0.1âµM or 1âµM], iFSP1 [10âµM], Erastin [0.37âµM], Sulfasalazine (SAS) [0.17âmM], Imidazole Ketone Erastin (IKE) [1.11âµM], ML210 [0.37âµM], ML162 [1.11âµM], TBHQ [25ânM]) cells were imaged using the 10à objective within the IncuCyte live cell imager (Sartorius). For dead cell quantification, 100ânM DRAQ7 (Thermofisher) were added to each well. For lipid ROS determination, cells were stained with 5âµM BODIPY C11. Cells were imaged for indicated timepoints every 2âh. Analysis for confluence, DRAQ7-positive (dead) or BODIPY C11-positive cells was performed using the Software IncuCyte 2021A (Sartorius).
Quantitative PCR
For KRASWT and KRASG12D comparison, 300,000 cells were seeded per well in a 6-well plate and RNA was extracted 24âh later. For MEK and AKT inhibition treatment experiments 200,000 cells were seeded in a 6-well plate one day in advance followed by treatments for 48âh. For TBHQ treatment experiments, 150,000 of KRASWT cells were seeded in a 6-well plate a day in advance followed by treatment for 24âh. For KEAP1 knockdowns, 200,000 of KRASWT cells were seeded in a 6-well plate on top of the transfection mix and incubated for 72âh. For LsL-KRASG12D-inducible MEFs, 200,000 cells were seeded in 1âµg/ml 4OHT in a 6-well plate and incubated for 4 or 5 days. For doxycycline-inducible KRASG12D HPDE cells, 450,000 cells were seeded in 0.5âµg/mL doxycycline in a 6-well plate and incubated for 72âh. For 4OHT-inducible HRASG12V NIH-3T3 cell, 35,000 cells were seeded for 72âh in a 6-well plate and 100ânM 4OHT was added and incubated for another 48âh.
For total RNA isolation, the NucleoSpin RNA kit (740955.5, Macherey-Nagel) was used according to the manufacturerâs instructions. Next, isolated RNA was reverse transcribed into cDNA using the LunaScript RT SuperMix Kit (E3010L, NEB). For quantitative PCR, 5âµl of Power SYBR GREEN PCR Master Mix (4368702, Thermo Fisher) was mixed with 2âµl of nuclease-free water (NEB), 1âµl (10âµM) of primer mix (forward and reverse primers) (see Supplementary Table 1 for primers used) and 2âµl of cDNA (5âµg/µl). Real-time qPCR was performed in triplicate or in quadruplicate on the Quant Studio5 qRT PCR cycler and results were normalized to the expression of the house-keeping gene indicated. Actin, Rplp0, Rpl13a or 18S were used as house-keeping gene controls as indicated.
Lipid ROS quantification (flow cytometry)
Thirty-five thousand cells per well were plated in a 24-well plate 24âh before treatment. Lipid ROS levels were quantified by BODIPY C11 (Invitrogen) staining. To this end, cells were stained using 5âµM BODIPY C11 during the last 30âmin of treatment incubation. Mean fluorescence intensity (MFI) was determined by flow cytometry using an LSR-FACS Fortessa (BD Bioscience) and FlowJo software (BD Bioscience). Flow cytometry data were collected from at least 5000 cells.
General ROS quantification (flow cytometry)
Fifty-five thousand cells per well were plated in a 24-well plate 24âh before treatment. Cells were incubated with 20âµM H2DCFDF (Invitrogen) per well to stain cellular ROS. Mean fluorescence intensity (MFI) was determined by flow cytometry using an LSR-FACS Fortessa (BD Bioscience) and FlowJo software (BD Bioscience). Flow cytometry data were collected from at least 5000 cells.
NADPH Assay
Twenty thousand cells per 96-well plate were seeded in advance. NADP/NADPH was determined by NADP/NADPH-Glo⢠Detection Reagent assay according to the manufacturerâs instructions (Promega).
Lipidomics to determine oxidized lipids and levels of total phospholipids
Mass spectrometry experiments to determine total phospholipids and oxidized lipids were performed as described previously [16]. In brief, levels of oxidized phosphatidylcholine (PC) and phosphatidylethanolamine (PE) species were determined by Liquid Chromatography coupled to Electrospray Ionization Tandem Mass Spectrometry (LC-ESI-MS/MS). Oxidized PC and PE species were quantified by normalizing their peak areas to those of the internal standards. Glycerophospholipids (PC and PE, including ether-linked species) in cells were analyzed by Nano-Electrospray Ionization Tandem Mass Spectrometry (Nano-ESI-MS/MS) with direct infusion of the lipid extract (Shotgun Lipidomics). The protein content of the homogenate was routinely determined using bicinchoninic acid. Endogenous glycerophospholipids were quantified by referring their peak areas to those of the internal standards. The calculated glycerophospholipid amounts were normalized to the protein content of the cell homogenate.
Western blotting
For KRASWT and KRASG12D 300,000 cells were seeded one day in advance in a 6-well plate before cells were treated for another 5âh with or without RSL3. For FSP1 knockdowns, 75,000 cells were seeded in a 6-well plate on top of the transfection mix and incubated for 72âh. For KRASWT and KRASG12D cells expressing either empty Vector or FSP1-WT 300,000 cells were seeded in a 6-well plate for 24âh. For LsL-KRASG12D-inducible MEFs, cells were incubated for 72âh with or without 1âµg/ml 4OHT before 250,000 cells were seeded into a 6-well plate for indicted timepoints. For MEKi and AKTi treatment experiments 150,000 cells were seeded one day in advance in a 6-well plate before cells were treated for 72âh. Cell lysates were prepared in lysis buffer (30âmM Tris-HCl [pH 7.4], 120âmM NaCl, 2âmM EDTA, 2âmM KCl, 1% Triton X-100, 1à COMPLETE protease-inhibitor and phosphatase- inhibitor cocktail). Lysate concentrations were adjusted to equal protein concentrations using the bicinchoninic acid (BCA) protein assay (Biorad). Equal amounts of protein were mixed with a final concentration of 1à reducing sample buffer (Invitrogen) and 200âmM DTT (VWR). Samples were heated to 95â°C for 10âmin, separated via gel electrophoresis and transferred to nitrocellulose membranes (Biorad). Membranes were blocked in PBS with 0.1% Tween 20 (PBST) (VWR) with 5% (w/v) dried milk powder (AppliChem) for at least 30âmin and incubated with primary antibodies over night. After washing with PBST, membranes were incubated with horse radish peroxidase (HRP)-coupled secondary antibodies (Biotium) diluted 1:10,000 for at least 1âh at room temperature. After another washing step, bound antibodies were detected using chemiluminescent Amersham ECL Prime Western Blotting Detection Reagent (Cytiva) or SuperSignal⢠West Femto Maximum Sensitivity Substrate (Thermo Fisher). X-ray films CL-XPosure⢠(Thermo Fisher) or the FUSION Solo S system and software (Vilber) were used to develop the membranes.
RNA sequencing
For RNA sequencing, 70,000 cells per well (6-well plate) of either KRASWT or KRASG12D were seeded 24âh in advance. The next day, cells were washed with PBS and RNA was isolated using the NucleoSpin RNA kit (740955.5, Macherey-Nagel) according to the manufacturerâs instructions. cDNA libraries amplified from the 3â² UTR were generated from total RNA using the Lexogen QuantSeq kit (Lexogen, Austria) according to the standard protocol and sequenced with a 50-bp single-end protocol on Illumina HiSeq4000 sequencer (Illumina, USA). The raw-sequencing data was aligned to the respective mouse reference genomes and quantified prior to differential expression analyses. Raw FPKM values of each transcript were transformed by log2 (FPKMâ+â0.01). Data processing and statistical analyses were performed using Microsoft Excel (Microsoft, USA) and Instant Clue software [17] which performed a hierarchical clustering to classify the experiments and generate a heatmap for the visualization of different RNA expression.
Generation of FSP1-overexpressing cells
To generate stable cells overexpressing FSP1, viral particles were produced in HEK-293T cells. HEK-293T cells were plated in a 6-well plate the day before transfection. For 6 à 6 wells of a 6-well plate 5âµg of each packaging plasmid and 10âµg of transfer vector plasmid were mixed together. Fifty microliters of 250âmM CaCl2 and 444âµl H2O were added to the plasmid mixture and mixed well by pipetting. For the formation of calcium-phosphate-DNA co-precipitate, the plasmid transfection mix (~500âµl) was carefully dropped into 500âµl of 2à HBS buffer under constant vortexing. The precipitate was incubated for 20âmin at room temperature and added dropwise into freshly replaced media without Pen/Strep. After 6â8âh of transfection medium was aspirated, and fresh normal medium was added. The following two days, virus-containing supernatant was harvested and filtered with 0.45âμm sterile syringe filter. Fresh medium was always added again on the cells. Virus harvest was centrifuged, and supernatant was collected and stored at â80â°C. For transduction of the KRASWT and KRASG12D cells viral supernatant was added to wells containing cells with 6âµg/ml polybrene and centrifuged for 45âmin at 2500ârpm at 30â°C. Cells were incubated afterwards at 37â°C and 5% CO2 until they were confluent and selected for positively-transfected cells using Puromycin (KRASWT 1âμg/ml and KRASG12D 2.5âμg/ml) for 7 days.
Generation of stable FSP1 knockdown cells
For stable transduction of FSP1 knockdown cells, lentiviral supernatant with an shFSP1 transfer plasmid was produced as described above and KRASG12D cells were transduced with the virus. After selection with Puromycin (2.5âμg/ml) for 7 days, knockdown was validated using qPCR.
Generation of LsL-KRASG12D-inducible MEFs
Mouse embryonic fibroblasts (MEFs) were generated from E13.5 mouse embryos by standard Jacks lab procedure (https://jacks-lab.mit.edu/protocols/making_mefs). MEFs with positive genotypes for 4-hydroxytamoxifen (4OHT)-inducible Cre (primer: #1 GCG GTC TGG CAG TAA AAA CTA TC, #2 GTG AAA CAG CAT TGC TGT CAC TT) and LsL-KRasG12D (primer: #1 gtc ttt ccc cag cac agt gc, #2 ctc ttg cct acg cca cca gct c; #3 agc tag cca cca tgg ctt gag taa gtc tgc a) were used for the experiments. To obtain Cre expression cells were treated with 1âµg/ml 4OHT.
Generation of HRASG12V-inducible NIH-3T3
NIH-3T3 stably expressing 4OHT-inducible HRASG12Vgenerated as previously described [18]. Activation of Ras was induced by exposure of cells to 100ânM 4OHT, with 4OHT-supplemented medium being refreshed every two days.
Generation of KRASG12D- inducible human pancreatic duct epithelial cells (HPDE)
Lentiviral particles were produced in HEK-293T cells transfected with packaging plasmids pCMV-VSV-6, pMDLg and pRSV-Rev as well as 10âµg pCW-Puro-KRASG12D (a newly cloned pCW (addgene #50661 [15]) backbone replacing the existing Cas9 by human KRASG12D cDNA) using a 1:1 mixture of 2à HBS and 250âmM CaCl2 in DMEM F12 medium with 10% FCS, 1% Pen/Strep and 1% Glutamine. After overnight incubation, medium was replaced for HPDE medium (see above), and medium containing viral particles was subsequently harvested after 48âh. HPDE cells were infected in a 6-well plate at 30% confluency by replacing the medium by viral supernatant after adding polybrene at a final concentration of 8âµg/mL. After 3 days, medium was replaced by selection medium containing 0.5âµg/mL puromycin. For KRASG12D induction, HPDE cells were treated with 0.5âµg/ml doxycycline for 72âh.
Spheroid assays
To generate spheroid cultures, 5000 cells per 96-well were plated in a 96 ultra-low attachment multi-well plate (Corning) in 100âµl media containing 4 % Matrigel (Corning). After 24âh 2à iFSP1 [10âµM or 20âµM], 2à Fer1 [2.5âµM] or 2à Liproxstatin-1 [0.3âµM] treatment in 100âµl media was added to the cells and incubated for 9 (A549) or 14 days (MEFs). Pictures of the spheroids were taken with the BZ-X800E microscope (Keyence). Spheroid assay and organoid assay colony area and brightness were analyzed using the BZ-H4M/Measurement Application Software (Keyence).
Soft agar colony formation assays
Cells were plated in 6-well plates at 8000 cells per well. Equal volumes of culturing medium and agarose were used such that the final concentrations were DMEM, 10% bovine serum (for 3T3 KRAS) or RPMI, 10% fetal bovine serum (for A549) and 0.25% agarose for the top layer or 0.4% agarose for the bottom layer, respectively. Where indicated, DMSO, iFSP1 (10 or 20âµM) or Fer1 [5âµM] was added to the top agar layer. Cells were fed twice a week with 1âmL of corresponding DMSO or iFSP1 treated medium onto the top layer. Colonies were allowed to form over the course of 18 to 30 days, following which they were imaged and quantified using ImageJ.
Tumor xenograft studies
Mice were maintained on a 12-h light/dark cycle with water and food ad libitum throughout the duration of the project. Mouse embryonic fibroblast (MEF) cell lines (5âÃâ105 cells either KRASWT e.V., KRASWT FSP1-WT, KRASG12D e.V. or KRASG12D shFSP1) were injected in 200âµl PBS into both flanks of 8â10 weeks old male NMRI-Foxn1 nu/nu mice (Janvier). Mice were not randomized. A group size of at least 10 tumors per condition was assumed to achieve significantly different results (pâ=â0.05) with a power of 80%. For that, cells were harvested from plates using trypsin and washed five times with PBS to remove residual FCS. Mice injected with KRASWT e.V. or KRASWT FSP1-WT were assigned to either vehicle or Liproxstatin-1 treatment groups once tumors reached a minimum size of 2.5 à 2.5âmm. For two consecutive weeks, mice were injected 5 à per week either with vehicle (PBS with 1% DMSO) or Liproxstatin-1 (10âmg/kg). Tumor size was tracked by caliper measurements and volume was calculated as (length à width à width)/2. People performing tumor measurements and calculating tumor volume were blinded to the group allocation. Mice were sacrificed at the end of the treatment and fresh-frozen tumor tissue was used for further analysis.
Protein extraction from fresh-frozen tumors
For protein isolation, 20â30âmg of fresh-frozen tumor tissue were mixed in a peqlab vial with the adequate number of ceramic beads and 500âµl IP-lysis buffer (30âmM Tris-HCl, 120âmM NaCl, 2âmM EDTA, 2âmM KCl, 1% Triton-X-100, pH 7.4, Protease and Phosphatase inhibitor (Roche)). For lysis, samples were homogenized for 2 à 30âs using the Precellys 24-dual homogenisator (Peqlab). Samples were centrifuged at 14,000 RPM for 20âmin at 4â°C and then further used for western blotting.
Isolation and treatment of murine pancreatic organoids
The pancreas was isolated from PDX1-Cre KRAS G12D mice, washed with cold mouse wash medium (DMEM high glucose + Pen/Strep + 1% FCS) and cut into 1â2âmm pieces using scalpels. Pancreatic pieces were transferred into 50âml falcons containing 10âml of mouse digestion medium (200âml Mouse wash medium + 25âmg Collagenase P; Sigma-Aldrich #C9407â+â25âmg Dispase II; Thermo Fisher #17105041) and shaken at 130ârpm at 37â°C for 20âmin. The supernatant was transferred to a Petri dish containing 10âml of mouse washing medium to obtain the first wash fraction. Ten milliliters of mouse digestion medium were added to the remaining pancreatic pieces. Cycles of shaking at 130ârpm at 37â°C for 10âmin were repeated until wash fractions with mainly pancreatic ducts and almost no acinar cells were observed. Wash fractions enriched with ducts were combined and spun down at 1200ârpm for 5âmin. The obtained duct pellet was resuspended in 200âµl of ice-cold Matrigel (Growth Factor Reduced, Phenol Red Free; Corning #356231) and a 30âµl dome was seeded into the middle of a well in a prewarmed 24-well plate. Ducts were consecutively diluted in Matrigel to obtain a cellular density with the most favorable conditions for organoids growth. The plate containing domes was placed into the cell culture incubator for 10â15âmin for the Matrigel to set before 500âµl of PancreaCult⢠Organoid Growth Medium (Mouse) (Stemcell #06040) was added. The organoids were maintained in cell culture with splitting once a week and twice-weekly medium change. For organoids treatments, single cells from 5 days old PDX1-Cre KRAS G12D organoids were isolated according to Boj et al. [19] and Huch et al. [20] and seeded at 200 cells in 100âµl medium/well in a 96-well plate covered in Matrigel:DPBS (1:1). Organoids from single cells were left to grow for 6 days, then treated with 50âµl of PancreaCult medium containing the indicated treatments for 2 days before microscopic pictures were taken using the BZ-X800E microscope (Keyence).
Analysis software and bioinformatic analysis
Heatmaps visualizing ferroptosis and KEGG pathway component expression were generated using Instant Clue software [17]. FACS data were analyzed and quantified using the FlowJo 10.4.2 software. Cell Titer Blue viability assays and qPCR results were analyzed using Excel. Lipidomics measurements were analyzed by MultiQuant 3.0.2 software (SCIEX). IncuCyte experiments were analyzed by using the Software IncuCyte 2021A. Soft Agar colonies were imaged and quantified using ImageJ. Spheroid assay and organoid assay colony area and brightness were analyzed using the BZ-H4M/Measurement Application Software (Keyence). Figures were assembled and data plotted and analyzed using GraphPad Prism 7 for Mac OS X.
Quantification, statistical analysis and reproducibility
GraphPad Prism 7 software (GraphPad Software Inc.) was used for Mac OS X to execute statistical analysis. For comparison between two conditions two-tailed t-tests were performed and for comparison between multiple samples two-way ANOVA and Tukeyâs post test for multiple comparisons were used. All data are presented as meanâ±âSEM of at least three independent biological replicates. From at least three independent experiments all means are calculated and plotted. Biological replicates gave comparable results, and no technical or biological replicates were excluded. In the respective figure legends statistical tests are declared. The following p value cut-offs were used for all tests: ****pâ<â0.0001, ***pâ<â0.001, **pâ<â0.01, *pâ<â0.05, nspâ>â0.05. Representative western blots are shown.
Results
Endogenous-level expression of oncogenic KRAS protects from ferroptosis
In order to interrogate the influence of oncogenic KRAS expression in genetically defined isogenic cellular systems, we made use of N- and HRAS-deficient mouse embryonic fibroblasts (MEFs) containing a LoxP-flanked KRAS gene as well as tamoxifen-inducible Cre recombinase. After induction of Cre recombinase- which renders these cells âRaslessâ [21]- cells were reconstituted with comparable expression levels of either wild type KRAS 4B or commonly mutated forms of KRAS 4B (hereafter referred to as KRAS; cell line panel available from the Ras initiative at the NIH national cancer institute, US). Strikingly, when treating this cell line panel with the GPX4 small molecule inhibitor RSL3, both WT KRAS-expressing clones were killed within 24âh of treatment whilst all cells expressing oncogenic variants of KRAS were more resistant (Fig. 1a). Moreover, co-treatment with the ferroptosis-selective antioxidant ferrostatin-1 (Fer-1) [8], blocked cell death induced in KRAS WT cells confirming the induction of ferroptosis (Fig. 1a). Importantly, this phenotype was not caused by varying levels in KRAS or expression of GPX4, as all cell lines expressed comparable levels of both proteins (Fig. 1b). Moreover, treatment with other class II ferroptosis-inducing compounds (FINs) which directly inhibit GPX4 equally led to a more drastic loss of viability in KRAS WT cells as compared to oncogenic KRAS while for class I FINs no significant difference could be observed (Supplementary Fig. 1a). In addition, loss of cell confluence induced by RSL3 as well as cell death was consistently less pronounced over time also in live cell imaging kinetic experiments in a representative KRASG12D-mutated MEF line (Supplementary Fig. 1b, c). As these MEFs are generated from single cellular clones, we independently generated Rasless MEFs from bulk-sorted populations to exclude clonal effects. These express near endogenous levels of FLAG-tagged KRAS behind GFP and an internal ribosomal entry site (IRES) and were enriched via bulk sorting of GFP+ cells (Fig. 1c). Importantly, expression of oncogenic KRASG12D but not KRAS WT equally rendered cells more resistant to ferroptotic cell death in bulk-sorted MEFs (Fig. 1d). Of note, ferrostatin-1 but not the caspase inhibitor zVAD or the RIPK1 inhibitor nec-1s could block cell death induced by RSL3 indicating ferroptotic cell death. Moreover, WT KRAS-expressing cells also died more rapidly upon RSL3 treatment than KRAS-mutant cells (Fig. 1e, f). In order to determine whether direct signaling from oncogenic KRAS is responsible for increased ferroptosis resistance we made use of the fact that effective small molecule inhibitors against KRASG12C have recently been developed [22,23,24,25]. Indeed, treating either KRASG12C-expressing MEFs or the established Lewis lung carcinoma cell line (3LL) - in which NRAS was knocked out (deltaNRAS86) leaving them only with an activating point mutation in KRASG12C [26] - with the KRASG12C inhibitor AMG510 sensitized these cells to ferroptosis (Supplementary Fig. 1d, e). Taken together, we find that expression of near-endogenous levels of oncogenic KRAS renders cells more resistant to ferroptotic cell death in various isogenic cellular and experimental setting.
KRAS-mutated cells are protected from ferroptosis-induced lipid peroxidation
A major hallmark of ferroptosis is a lipid ROS-dependent lipid peroxidation chain reaction [27] which oxidizes phosphatidylethanolamine (PE) and phosphatidylcholine (PC) species containing arachidonic (AA) and adrenic acid (AdA) [28, 29]. Since we observed that cells expressing various forms of mutated KRAS were more resistant to ferroptosis than KRAS WT cells, we next determined the extent of lipid ROS accumulation. Indeed, both WT KRAS clones readily accumulated lipid ROS 5âh after stimulation with RSL3 while all KRAS mutants tested did not show lipid ROS accumulation at this time (Fig. 2a, b). Bulk-sorted MEFs expressing flag-tagged mutant KRASG12D but not WT KRAS, similarly presented with decreased lipid ROS accumulation upon inhibition of GPX4 (Supplementary Fig. 2a). Moreover, KRAS-mutated cells also showed a decrease in the accumulation of oxidized BODIPY C11 indicative of the presence of lipid ROS in time-lapse imaging (Fig. 2c). These data suggested that oncogenic KRAS limits the propagation of lipid ROS and thereby acute lipid peroxidation upon induction of ferroptosis. As the extent of lipid peroxidation during ferroptosis is coupled to cellular amounts of AA-containing polyunsaturated fatty acid (PUFA) PE and PC species including ether-linked PUFAs [30], we first measured basal levels of diacylglycerol (DAG) and ether-linked PE and PC PUFAs in KRAS WT as compared to mutant cells using mass spectrometry. While levels of most PUFA species were comparable, a few PC species were elevated in KRAS WT cells (Supplementary Fig. 2bâe). Next, we determined the levels of phospholipid oxidation upon ferroptosis induction in Rasless MEFs expressing either KRAS WT or KRASG12D by mass spectrometry. Although we observed elevated basal lipid ROS in KRAS-mutated cells (Supplementary Fig. 2f) -likely as a result of elevated basal levels of total ROS due to oncogene expression (Supplementary Fig. 2g) - specific induction of lipid peroxidation upon GPX4 inhibition was absent in KRAS-mutated cells (Fig. 2d). Interestingly, total ROS was nevertheless readily induced in KRAS-mutated cells by GPX4 inhibition (Supplementary Fig. 2g) suggesting the protection from ROS to be specific towards lipid ROS. Together, these data establish that mutant KRAS endows cells with superior capacity to protect cells from a ferroptosis-specific increase in lipid peroxidation.
Elevated levels of FSP1 protect KRAS-mutated cells from ferroptosis
In order to determine the mechanism by which KRAS-mutated cells might buffer acute lipid peroxidation, we performed comparative 3ⲠRNA sequencing of KRAS WT and KRASG12D-mutated MEFs. Interestingly, when analyzing the top 1000 upregulated genes in KRAS-mutated cells for functional association networks using STRING, we obtained a significant enrichment of the ferroptosis pathway along with several other metabolic pathways (Fig. 3a). As the list of genes annotated in the KEGG ferroptosis pathway does not contain more recently discovered regulators of ferroptosis, we manually extended this list (KEGG+) and analyzed expression of these genes in our comparative dataset. Strikingly, AIFM2 mRNA, recently renamed as ferroptosis suppressor protein 1 (FSP1) due to its ferroptosis protective activity [11, 12], was upregulated in KRAS-mutated cells within a cluster of genes (Fig. 3b). Within several ferroptosis regulatory genes of the dataset, FSP1 was upregulated significantly (Supplementary Fig. 3a). Moreover, FSP1 mRNA upregulation in KRAS-mutated cells could also be confirmed by quantitative real-time PCR (qPCR) (Fig. 3c). Importantly, FSP1 was basally upregulated also on protein level in KRAS-mutated cells and, unlike xCT, a recently identified target gene further upregulated upon H2O2 [7], was not further increased upon stimulation with RSL3 (Fig. 3d). As FSP1 has been shown to render cells more resistant to ferroptosis, we next tested whether its elevated expression in KRAS-mutated cells was responsible for mediating increased ferroptosis resistance of KRAS-mutated cells. Indeed, FSP1 suppression was sufficient to sensitize KRAS-mutated cells to ferroptosis (Fig. 3e, f). Vice versa, overexpression of FSP1 was sufficient to render KRAS WT cells as resistant to ferroptosis as KRAS-mutated control cells (Fig. 3g, h). As FSP1 is an NADH ubiquinone oxidoreductase, FSP1 activity requires NADPH as an electron source [11, 12]. Therefore, we also tested whether basal levels of the FSP1 cofactor NADPH would differ in KRAS-mutated as compared to WT cells and thereby contribute to differential activity, yet this was not the case (Supplementary Fig. 3b). To next determine whether FSP1 activity may also protect KRAS-mutated cells from ferroptosis, we employed a recently developed small molecule inhibitor against FSP1 [11]. Strikingly, co-incubation with this inhibitor (iFSP1) reverted ferroptosis resistance endowed by oncogenic KRAS expression (Supplementary Fig. 3c, d). Of note, iFSP1 also slightly sensitized KRAS WT cells, yet due to the fact WT cells were already very sensitive, the relative sensitization observed was much stronger for KRAS-mutated cells. Thus, KRAS-mutated cells display increased ferroptosis resistance due to elevated levels of FSP1.
Oncogenic KRAS upregulates FSP1 via NRF2 and the MAPK pathway
To thoroughly test a direct mechanistic link of oncogenic KRAS and FSP1 induction, we generated primary MEFs with inducible expression of KRASG12D from its endogenous locus similar to an established approach [31] (LsL-KRASG12D-inducible MEFs). As expected, induction of KRASG12D by 4-hydroxytamoxifen (4OHT) treatment led to enhanced basal phosphorylation of ERK which was further increased upon refeeding with FCS (Fig. 4a). Moreover, KRASG12D induction readily elevated expression of the established MAPK target gene dual specificity phosphatase 6 (DUSP6) along with the antioxidant transcription factor nuclear factor erythroid 2-related factor 2 (NRF2) and its bona fide target genes glutamate-cysteine ligase catalytic subunit (GCLC) and heme oxigenase-1 (HO-1) (Fig. 4b). Strikingly, induction of KRASG12D was indeed sufficient to also directly induce FSP1 expression (Fig. 4b). In order to validate direct induction of FSP1 by oncogenic RAS, we validated this finding in human pancreatic duct epithelial cells (HPDE) with doxycycline-inducible expression of KRASG12D and HRASG12V-inducible NIH-3T3 cells (Supplementary Fig. 4a, b). Given that oncogenic KRAS is known to induce NRF2 [31] and FSP1 was very recently shown to be a direct transcriptional target of NRF2 [32] in NSCLC mutated in the NRF2 inhibitor kelch-like ECH-associated protein 1 (KEAP1), we next tested whether activating endogenous NRF2 would be sufficient in our cells to induce FSP1. Strikingly, silencing of KEAP1 readily induced FSP1 along with NRF2 target genes (Fig. 4c). Moreover, KEAP1 knockdown rendered KRAS WT cells resistant to ferroptosis, which could partially be reverted by iFSP1 treatment (Fig. 4d).
Furthermore, treatment with the chemical NRF2 activator tert-butylhydrochinone (TBHQ) equally induced NRF2 target genes along with FSP1 (Fig. 4e). While silencing of NRF2 was toxic to MEFs (a toxicity blockable by Fer-1), we mined publicly available datasets from the KRAS-mutated NSCLC cell line A549 in which NRF2 was knocked down (Fig. 4f) [33]. Indeed, upon NRF2 silencing, FSP1 expression was also significantly decreased (Fig. 4g) confirming NRF2-mediated FSP1 regulation to be present in KRAS-mutated cells. These data suggest that oncogenic KRAS-induced NRF2 directly leads to elevated transcription of FSP1 thereby protecting them from ferroptosis. To also determine which other major KRAS effector pathway may upregulate FSP1, KRAS-mutated cells were treated with the MEK inhibitor PD184352 to block the MAPK arm downstream of KRAS as well as the AKT inhibitor MK2206. Interestingly, FSP1 mRNA was reduced only under MEK but not AKT inhibition indicating that MAPK pathway activation is responsible for elevated FSP1 levels in KRAS-mutated cells (Supplementary Fig. 4c). Of note, MEK inhibition in WT cells did not significantly regulate FSP1 mRNA despite effectively blunting expression of the established MAPK pathway target gene DUSP6 indicating preferential regulation in KRAS-mutated cells (Supplementary Fig. 4d). In support of this, MEK inhibition resulted in dose-dependent reduction in FSP1 protein levels along with decreased phosphorylation of ERK (Supplementary Fig. 4e). Furthermore, gene-set enrichment analysis of genes co-expressed with FSP1 in the lung adenocarcinoma (LUAD) TCGA dataset were significantly enriched for the MAPK pathway as indicated by enrichment of the RAF and MEK pathways (Supplementary Fig. 4f). Of note, NRF2 has been reported to be phosphorylated and activated by MAPK signaling [34]. Therefore, we also measured expression of the NRF2 target gene GCLC under MEK as compared to AKT inhibition. Indeed, GCLC was reduced by inhibition of MEK but not AKT (Supplementary Fig. 4g). Taken together, our data propose that FSP1 is upregulated in KRAS-mutated cells as a direct result of MAPK and NRF2 pathway activation.
FSP1 aids KRAS-mediated cellular transformation and promotes tumor onset in vivo
One hallmark of oncogenic KRAS is its capacity to mediate cellular transformation. Therefore, we hypothesized that FSP1 activity may aid cellular transformation capacity of KRAS-mutated cells. Indeed, 3T3 cells transformed by KRASG12V expression presented with decreased colony formation in soft agar in the presence of iFSP1. Yet, this activity of FSP1 was not due to ferroptosis protection as co-treatment with Fer-1 did not rescue decreased colony formation (Fig. 5a). Moreover, human KRAS-mutated A549 cells equally showed decreased soft agar colony formation in the presence of iFSP1 (Fig. 5b) suggesting that elevated FSP1 expression additionally may promote cellular transformation of KRAS-mutated cells independently of ferroptosis protection. Recently, in 3D Matrigel-based spheroid assays, ferroptosis was shown to occur in spheroid centers thereby limiting their growth [35]. Therefore, we next tested the extent of spheroid formation in KRAS-mutated as compared to WT cells in this experimental system. As expected, KRAS-mutated cells formed spheroids much more efficiently than KRAS WT cells although a few colonies could be detected likely due to some extent of spontaneous transformation enabling continuous proliferation of KRAS WT MEFs (Fig. 5c). Here, overexpression of FSP1 was sufficient to allow for spheroid growth in KRAS WT cells to a similar extent as KRAS-mutant cells. While FSP1 overexpression in KRAS-mutant cells led to a slight decrease in colony formation, importantly, iFSP1 incubation readily reverted spheroid formation enabled by FSP1 overexpression in both cases (Fig. 5c). Moreover, iFSP1 also significantly impacted spheroid formation of human A549 cells yet again, this activity was not caused by protecting from ferroptosis (Fig. 5d, Supplementary Fig. 5a). Thereby, our data support the concept that cellular transformation endowed by KRAS may at least in part depend upon FSP1 expression and activity, yet this activity is independent of its role in ferroptosis protection. Next, we aimed to test whether FSP1 expression was sufficient to impact tumor initiation. To this end, we transplanted mice with KRAS WT cells with either control or FSP1 overexpression in comparison to KRAS-mutant control and cells with short hairpin RNA (shRNA)-mediated FSP1 silencing (Supplementary Fig. 5b). As expected, KRAS-mutated cells presented with earlier tumor onset than KRAS WT cells. Yet strikingly, FSP1 expression was sufficient in KRAS WT cells to significantly accelerate tumor onset and increase tumor incidence closer to the rates of KRAS-mutant tumors (Fig. 5e). Moreover, treatment of mice bearing KRAS WT tumors with the ferroptosis-selective inhibitor Liproxstatin-1 accelerated tumor onset to similar levels as FSP1 overexpression suggesting ferroptosis to be indeed responsible for suppression of tumor initiation capacity of WT cells in vivo. Importantly, Liproxstatin-1 treatment of FSP1-overexpressing WT tumors did not result in additional promotion of tumor onset indicating FSP1 to promote KRAS WT tumor initiation by protecting from ferroptosis. Interestingly, and, similar to our findings in spheroid assays, tumor volumes of tumors arising also increased when expressing exogenous FSP1, yet this increase was independent of ferroptosis protection as Liproxstatin-1 treatment did not affect tumor volumes (Supplementary Fig. 5c). Vice versa, FSP1 silencing in KRAS-mutated tumors delayed tumor onset albeit not to the levels of KRAS WT tumors. Therefore, oncogenic KRAS partially promotes its early tumor onset through FSP1-mediated ferroptosis protection. Moreover, FSP1 expression alone is sufficient to promote tumor initiation in the absence of oncogenic KRAS by suppressing ferroptosis in vivo.
Oncogenic mutations in KRAS are most frequently observed in patients with pancreatic ductal adenocarcinoma (PDAC). Yet, most small molecules inducing ferroptosis do not display pharmacokinetics and solubilities suitable for in vivo use yet. Therefore, we next generated organoids from mice developing pancreatic intraepithelial neoplasia (PanINs) as a result of KRASG12D expression from the endogenous promotor [36] and treated them either with RSL3 alone or in combination with iFSP1. Strikingly, RSL3 treatment alone was insufficient to induce ferroptosis in pancreatic organoids but the combination with iFSP1 led to effective killing of pancreatic organoids expressing KRASG12D (Fig. 5f). Based on these data, we propose that breaking ferroptosis resistance through the use of FSP1 inhibitors might be a particularly potent treatment strategy against KRAS-driven cancers.
FSP1 expression is upregulated in KRAS-mutated cancers and correlates with poor outcome in PDAC patients
In order to test whether FSP1 expression may be upregulated in KRAS-driven cancer, we probed publicly available tumor (TCGA) and normal (GTEX) datasets for FSP1 expression in colorectal cancer, lung adenocarcinoma and pancreatic adenocarcinoma in comparison to their respective normal tissue of origin using Gene Expression Profiling Interactive Analysis (GEPIA). Interestingly, FSP1 was significantly overexpressed in all three tumor types as compared to the respective normal tissues (Fig. 6a). While in pancreatic cancer the vast majority of patients present with activating KRAS mutations and a KRAS WT group is therefore difficult to obtain, in non-small cell lung cancer (NSCLC) about half of the patients are usually WT. Strikingly, when dividing an NSCLC dataset (GSE31852) by KRAS-mutation status, FSP1 expression was significantly upregulated in KRAS-mutated patient material (Fig. 6b). Moreover, in two independent PDAC datasets [37, 38] with included adjacent normal tissue, FSP1 expression was significantly upregulated in PDAC over normal pancreas and FSP1 mRNA correlated with NRF2 mRNA within these two datasets (Fig. 6câf). Moreover, we performed gene-set enrichment analysis (GSEA) on transcription factor motifs of genes co-expressed with FSP1 within the TCGA PDAC dataset. Strikingly, NRF2 was amongst the top 10 enriched motifs (Fig. 6g) suggesting FSP1 upregulation in PDAC patients to be a result of NRF2-mediated transcription. Lastly, high FSP1 expression in PDAC patients showed a strong trend towards drastically shortened relapse-free survival (Fig. 6h).
Taken together, our data establish that endogenous levels of oncogenic KRAS expression render cells more resistant to ferroptosis by upregulating FSP1 through the NRF2 and MAPK pathway allowing for a superior capacity to buffer acute lipid peroxidation during tumor initiation. Hence, only combined targeting of GPX4 and FSP1 is effective at killing KRAS-driven pancreatic organoids and FSP1 is upregulated in human KRAS-driven cancers. Based on these data, we propose that pro-ferroptotic therapy for KRAS-driven cancers should include inhibition of FSP1 in order to achieve efficient tumor cell killing.
Discussion
In the present study, we identify that expression of near-endogenous levels of oncogenic KRAS renders cells more resistant to ferroptosis through elevated expression of FSP1. Of note, in earlier studies overexpression of oncogenic HRAS in fibroblasts as compared to cells expressing empty vector led to a sensitization to erastin-induced cell death [13, 39], later found to be ferroptotic due to erastin-mediated targeting of xCT [10]. What might seem as a discrepancy might in fact represent two distinct stages of cellular transformation. Upon acute overexpression of an oncogenic RAS variant, cellular levels of ROS are known to be upregulated due to NOX1 induction [6]. Moreover, in line with an earlier study by Yang et al. [40] we observed that a cluster of genes involved in iron uptake (TFRC, STEAP3) was in fact upregulated in KRAS-mutated cells (Fig. 3b), possibly to feed an increased requirement for iron in the mitochondrial respiratory chain. This may suggest that RAS expression would fuel lipid peroxidation due to elevated basal ROS which is also what was observed by Yang et al. 2008. Yet, chronic elevated levels of ROS are known to activate NRF2. In fact, endogenous expression levels as opposed to overexpression of oncogenic KRAS were shown to effectively induce NRF2 activation [31]. Along these lines, a recent study found that NRF1 and NRF2 protect cells from ferroptosis through distinct and independent mechanisms [41]. Moreover, NRF2 was recently shown to induce FSP1 transcription via NRF2 in KEAP1-mutant NSCLC [32]. In keeping with these results, we now find that cells expressing oncogenic KRAS directly induce FSP1 expression through MAPK-NRF2 pathway activation. Of note, class I FINs, which lead to GSH depletion, did not show a significant difference between KRAS-mutated and -WT cells. This observation may suggest that GSH depletion triggers general ROS accumulation which in turn activates NRF2 and via this route may upregulate FSP1 also in WT cells thereby neutralizing the difference in killing. Hence, we propose a model wherein effective activation of NRF2-mediated transcription might be decisive for whether oncogenic KRAS expression renders cells more or less sensitive to ferroptosis.
An interesting possibility is that the Nonsense-Mediated Decay (NMD) pathway might contribute to the regulation of FSP1 mRNA. NMD can regulate a number of perfectly functional transcripts many of which usually have an abnormally long 3â²UTR. Indeed, FSP1 murine transcript variants contain several NMD-inducing features including an abnormally long 3â²UTR in one of them. Interestingly, FSP1 mRNA was found to be upregulated in murine embryonic stem cells upon NMD inhibition [42], and oncogenic KRAS might inhibit NMD via several non-exclusive mechanisms. For instance, KRASG12D is known to activate p38 MAPK [43,44,45] and p38 MAPK has been shown to inhibit NMD [46]. Additionally, ER stress is well-known to inhibit NMD [47] and has been observed in KRASG12D-expressing cells [48]. It is therefore tempting to speculate whether cells expressing oncogenic KRAS may inhibit NMD which in turn could result in FSP1 mRNA stabilization.
Expression of oncogenic KRAS creates selective metabolic addiction to nucleotide synthesis via the pentose phosphate pathway (PPP) [49]. Interestingly, high cellular levels of NADPH, the product of the PPP, were identified as markers of resistance to ferroptosis [50]. NADPH in turn is an important electron donor for a variety of cellular enzymes including FSP1. Consequently, high-level exogenous overexpression of FSP1, as obtained in our experiments, may impact cellular NADPH levels and, with that, be problematic for metabolic NADPH addiction of KRAS-mutant cells [49]. In line with this, FSP1 overexpression in KRAS-mutant cells decreased their capacity to form spheroids. Moreover, inhibition of the PPP was shown to reduce soft agar colony formation of transformed 3T3 MEFs [51]. In favor of another ferroptosis-independent function for FSP1 expression in cancer, we found that tumor volume growth in WT cells in vivo was promoted by FSP1 expression but not Liproxstatin-1 treatment. This growth promotion may be facilitated by FSP1- generated NAD+ which was shown to promote glycolysis [52], a mode of energy generation advantageous for hypoxic tumors.
Interestingly, two recent studies identified GTP cyclohydrolase-1 (GCH1) and its products tetrahydrobiopterin /dihydrobiopterin (BH4/BH2) to act as potent cellular antioxidants protecting from ferroptosis in the absence of GPX4 [53, 54]. Yet, cells not expressing GCH1 seem to solely depend upon FSP1 for the generation of endogenous radical-trapping agents, a fact used for screening for novel FSP1 inhibitors [55]. Of note, GCH1 mRNA expression was barely detectable in our cellular systems and also not influenced by KRAS-mutation status in contrast to FSP1 (data not shown).
In a genetically engineered mouse model of KRAS-driven PDAC, inducible whole-body deletion of xCT led to significant tumor regression [56]. However interestingly, cancer-associated fibroblast (CAF)-restricted deletion of xCT in the very same mouse model was sufficient to achieve a strong anti-tumor effect [57]. These data together with the fact that GPX4 deletion within PanINs was insufficient to trigger ferroptosis in pancreatic cancer [14], support the idea that KRAS-mutated cells have evolved an additional layer of protection against ferroptotic cell death. Our data propose that elevated FSP1 expression in KRAS-mutated cells is, at least in part, responsible for this protection. Based on these considerations, we propose that combined induction of ferroptosis and FSP1 inhibition should be considered for therapeutic strategies developed against KRAS-mutated cancers.
Data availability
The data and material that support the findings of this study are available from the corresponding author upon reasonable request.
References
Bos JL. ras oncogenes in human cancer: a review. 1989;49:4682â9.
Downward J. Targeting RAS signalling pathways in cancer therapy. Nat Rev Cancer. 2003;3:11â22.
Downward J. Targeting RAS and PI3K in lung cancer. Nat Med. 2008;14:1315â6.
Hoogwater FJH, Nijkamp MW, Smakman N, Steller EJA, Emmink BL, Westendorp BF, et al. Oncogenic K-Ras turns death receptors into metastasis-promoting receptors in human and mouse colorectal cancer cells. Gastroenterology. 2010;138:2357â67.
Karstedt S, von Conti A, Nobis M, Montinaro A, Hartwig T, Lemke J, et al. Cancer cell-autonomous TRAIL-R signaling promotes KRAS-driven cancer progression, invasion, and metastasis. Cancer Cell. 2015;27:561â73.
Mitsushita J, Lambeth JD, Kamata T. The superoxide-generating oxidase Nox1 is functionally required for Ras oncogene transformation. Cancer Res. 2004;64:3580â5.
Lim JKM, Delaidelli A, Minaker SW, Zhang H-F, Colovic M, Yang H, et al. Cystine/glutamate antiporter xCT (SLC7A11) facilitates oncogenic RAS transformation by preserving intracellular redox balance. Proc Natl Acad Sci USA 2019;1092:201821323.
Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149:1060â72.
Yang WS, SriRamaratnam R, Welsch ME, Shimada K, Skouta R, Viswanathan VS, et al. Regulation of ferroptotic cancer cell death by GPX4. Cell. 2014;156:317â31.
Dixon SJ, Patel DN, Welsch M, Skouta R, Lee ED, Hayano M, et al. Pharmacological inhibition of cystine-glutamate exchange induces endoplasmic reticulum stress and ferroptosis. eLife. 2014;3:e02523.
Doll S, Freitas FP, Shah R, Aldrovandi M, Silva MC da, Ingold I, et al. FSP1 is a glutathione-independent ferroptosis suppressor. Nature. 2019;149:1060â20.
Bersuker K, Hendricks J, Li Z, Magtanong L, Ford B, Tang PH, et al. The CoQ oxidoreductase FSP1 acts parallel to GPX4 to inhibit ferroptosis. Nature. 2019;575:688â92.
Yagoda N, Rechenberg M, von Zaganjor E, Bauer AJ, Yang WS, Fridman DJ, et al. RAS-RAF-MEK-dependent oxidative cell death involving voltage-dependent anion channels. Nature. 2007;447:864â8.
Dai E, Han L, Liu J, Xie Y, Zeh HJ, Kang R, et al. Ferroptotic damage promotes pancreatic tumorigenesis through a TMEM173/STING-dependent DNA sensor pathway. Nat Commun. 2020;11:6339â11.
Wang T, Wei JJ, Sabatini DM, Lander ES. Genetic screens in human cells using the CRISPR-Cas9 system. Science. 2014;343:80â84.
Bebber CM, Thomas ES, Stroh J, Chen Z, Androulidaki A, Schmitt A, et al. Ferroptosis response segregates small cell lung cancer (SCLC) neuroendocrine subtypes. Nat Commun. 2021;12:2048â19.
Nolte H, MacVicar TD, Tellkamp F, Krüger M. Instant clue: a software suite for interactive data visualization and analysis. Sci Rep. 2018;8:12648.
Hosogane M, Funayama R, Nishida Y, Nagashima T, Nakayama K. Ras-induced changes in H3K27me3 occur after those in transcriptional activity. PLoS Genet. 2013;9:e1003698.
Boj SF, Hwang C-I, Baker LA, Chio IIC, Engle DD, Corbo V, et al. Organoid models of human and mouse ductal pancreatic cancer. Cell. 2015;160:324â38.
Huch M, Bonfanti P, Boj SF, Sato T, Loomans CJM, Wetering Mvande, et al. Unlimited in vitro expansion of adult biâpotent pancreas progenitors through the Lgr5/Râspondin axis. Embo J. 2013;32:2708â21.
Drosten M, Dhawahir A, Sum EYM, Urosevic J, Lechuga CG, Esteban LM, et al. Genetic analysis of Ras signalling pathways in cell proliferation, migration and survival. EMBO J. 2010;29:1091â104.
Ostrem JM, Peters U, Sos ML, Wells JA, Shokat KM. K-Ras(G12C) inhibitors allosterically control GTP affinity and effector interactions. Nature. 2013;503:548â51.
Vasta JD, Peacock DM, Zheng Q, Walker JA, Zhang Z, Zimprich CA, et al. KRAS is vulnerable to reversible switch-II pocket engagement in cells. Nat Chem Biol. 2022;18:596â604.
Janes MR, Zhang J, Li L-S, Hansen R, Peters U, Guo X, et al. Targeting KRAS Mutant Cancers with a Covalent G12C-Specific Inhibitor. Cell. 2018;172:578â581.e17.
Canon J, Rex K, Saiki AY, Mohr C, Cooke K, Bagal D, et al. The clinical KRAS(G12C) inhibitor AMG 510 drives anti-tumour immunity. Nature. 2019;56:779â21.
Molina-Arcas M, Moore C, Rana S, Maldegem F van, Mugarza E, Romero-Clavijo P, et al. Development of combination therapies to maximize the impact of KRAS-G12C inhibitors in lung cancer. Sci. Transl. Med. 2019;11.
Angeli JPF, Krysko DV, Conrad M. Ferroptosis at the crossroads of cancer-acquired drug resistance and immune evasion. Nat Rev Cancer. 2019;15:348.
Kagan VE, Mao G, Qu F, Angeli JPF, Doll S, Croix CS, et al. Oxidized arachidonic and adrenic PEs navigate cells to ferroptosis. Nat Chem Biol. 2017;13:81â90.
Wiernicki B, Dubois H, Tyurina YY, Hassannia B, Bayir H, Kagan VE, et al. Excessive phospholipid peroxidation distinguishes ferroptosis from other cell death modes including pyroptosis. Cell Death Dis. 2020;11:922â11.
Zou Y, Henry WS, Ricq EL, Graham ET, Phadnis VV, Maretich P, et al. Plasticity of ether lipids promotes ferroptosis susceptibility and evasion. Nature. 2020;585:603â8.
DeNicola GM, Karreth FA, Humpton TJ, Gopinathan A, Wei C, Frese K, et al. Oncogene-induced Nrf2 transcription promotes ROS detoxification and tumorigenesis. Nature. 2011;475:106â9.
Koppula P, Lei G, Zhang Y, Yan Y, Mao C, Kondiparthi L, et al. A targetable CoQ-FSP1 axis drives ferroptosis- and radiation-resistance in KEAP1 inactive lung cancers. Nat Commun. 2022;13:2206.
Namani A, Cui QQ, Wu Y, Wang H, Wang XJ, Tang X. NRF2-regulated metabolic gene signature as a prognostic biomarker in non-small cell lung cancer. Oncotarget. 2017;8:69847â62.
Sun Z, Huang Z, Zhang DD. Phosphorylation of Nrf2 at multiple sites by MAP kinases has a limited contribution in modulating the Nrf2-dependent antioxidant response. PLoS ONE. 2009;4:e6588.
Takahashi N, Cho P, Selfors LM, Kuiken HJ, Kaul R, Fujiwara T, et al. 3D culture models with CRISPR screens reveal hyperactive NRF2 as a prerequisite for spheroid formation via regulation of proliferation and ferroptosis. Mol Cell. 2020;80:828â844.e6.
Hingorani SR, Petricoin EF, Maitra A, Rajapakse V, King C, Jacobetz MA, et al. Preinvasive and invasive ductal pancreatic cancer and its early detection in the mouse. Cancer Cell. 2003;4:437â50.
Pei H, Li L, Fridley BL, Jenkins GD, Kalari KR, Lingle W, et al. FKBP51 affects cancer cell response to chemotherapy by negatively regulating Akt. Cancer Cell. 2009;16:259â66.
Janky R, Binda MM, Allemeersch J, broeck AV den, Govaere O, Swinnen JV, et al. Prognostic relevance of molecular subtypes and master regulators in pancreatic ductal adenocarcinoma. Bmc Cancer. 2016;16:632.
Dolma S, Lessnick SL, Hahn WC, Stockwell BR. Identification of genotype-selective antitumor agents using synthetic lethal chemical screening in engineered human tumor cells. Cancer Cell. 2003;3:285â96.
Yang WS, Stockwell BR. Synthetic lethal screening identifies compounds activating iron-dependent, nonapoptotic cell death in oncogenic-RAS-harboring cancer cells. Chem Biol. 2008;15:234â45.
Forcina GC, Pope L, Murray M, Dong W, Abu-Remaileh M, Bertozzi CR, et al. Ferroptosis regulation by the NGLY1/NFE2L1 pathway. Proc Natl Acad Sci USA 2022;119:e2118646119.
Hurt JA, Robertson AD, Burge CB. Global analyses of UPF1 binding and function reveal expanded scope of nonsense-mediated mRNA decay. Genome Res. 2013;23:1636â50.
Munoz L, Yeung YT, Grewal T. Oncogenic Ras modulates p38 MAPK-mediated inflammatory cytokine production in glioblastoma cells. Cancer Biol Ther. 2016;17:355â63.
Houdt WJ van, Bruijn MT de, Emmink BL, Raats D, Hoogwater FJH, Rinkes IHMB, et al. Oncogenic K-ras activates p38 to maintain colorectal cancer cell proliferation during MEK inhibition. Cell Oncol. 2010;32:245â57.
Coelho MA, TÅecesson S, de C, Rana S, Zecchin D, Moore C, et al. Oncogenic RAS signaling promotes tumor immunoresistance by stabilizing PD-L1 mRNA. Immunity. 2017;47:1083â1099.e6.
Nickless A, Cheruiyot A, Flanagan KC, Piwnica-Worms D, Stewart SA, You Z. p38 MAPK inhibits nonsense-mediated RNA decay in response to persistent DNA damage in noncycling cells. J Biol Chem. 2017;292:15266â76.
Li Z, Vuong JK, Zhang M, Stork C, Zheng S. Inhibition of nonsense-mediated RNA decay by ER stress. Rna. 2017;23:378â94.
Ramadori G, Konstantinidou G, Venkateswaran N, Biscotti T, Morlock L, Galié M, et al. Diet-induced unresolved ER stress hinders KRAS-driven lung tumorigenesis. Cell Metab. 2015;21:117â25.
Santana-Codina N, Roeth AA, Zhang Y, Yang A, Mashadova O, Asara JM, et al. Oncogenic KRAS supports pancreatic cancer through regulation of nucleotide synthesis. Nat Commun. 2018;9:4945â13.
Shimada K, Hayano M, Pagano NC, Stockwell BR. Cell-line selectivity improves the predictive power of pharmacogenomic analyses and helps identify NADPH as biomarker for ferroptosis sensitivity. Cell Chem Biol. 2016;23:225â35.
Kuo W, Lin J, Tang TK. Human glucoseâ6âphosphate dehydrogenase (G6PD) gene transforms NIH 3T3 cells and induces tumors in nude mice. Int J Cancer. 2000;85:857â64.
Nguyen HP, Yi D, Lin F, Viscarra JA, Tabuchi C, Ngo K, et al. Aifm2, a NADH oxidase, supports robust glycolysis and is required for cold- and diet-induced thermogenesis. Mol Cell. 2020;77:600â617.e4.
Kraft VAN, Bezjian CT, Pfeiffer S, Ringelstetter L, Müller C, Zandkarimi F, et al. GTP cyclohydrolase 1/Tetrahydrobiopterin counteract ferroptosis through lipid remodeling. Acs Cent Sci. 2020;6:41â53.
Soula M, Weber RA, Zilka O, Alwaseem H, La K, Yen F, et al. Metabolic determinants of cancer cell sensitivity to canonical ferroptosis inducers. Nat Chem Biol. 2020;16:1351â60.
Yoshioka H, Kawamura T, Muroi M, Kondoh Y, Honda K, Kawatani M, et al. Identification of a small molecule that enhances ferroptosis via inhibition of ferroptosis suppressor protein 1 (FSP1). Acs Chem Biol. 2022;17:483â91.
Badgley MA, Kremer DM, Maurer HC, DelGiorno KE, Lee H-J, Purohit V, et al. Cysteine depletion induces pancreatic tumor ferroptosis in mice. Science. 2020;368:85â89.
Sharbeen G, McCarroll JA, Akerman A, Kopecky C, Youkhana J, Kokkinos J, et al. Cancer-associated fibroblasts in pancreatic ductal adenocarcinoma determine response to SLC7A11 inhibition. Cancer Res. 2021;81:3461â79.
Lánczky A, GyÅrffy B. Web-based survival analysis tool tailored for medical research (KMplot): development and implementation. J Med Internet Res. 2021;23:e27633.
Acknowledgements
Susanne Brodesser for advice on lipidomics, M. Conrad for providing the FSP1 antibody, J.P. Friedmann-Angeli for providing FSP1 plasmids and advice on Liproxstatin-1 in vivo usage, Dennis Plencker for advice on pancreatic organoid cultures, A. Trauzold and J. Downward for providing cells (HPDE, A549, respectively), Nina Wobst for additional validation in human cell lines, and Karin Schlegelmilch for solving cloning problems.
Funding
Work in the von Karstedt lab was funded through a Max Eder junior research group grant (701125509) by the German cancer aid, a project grant (A06) by the center for molecular medicine cologne (CMMC), an eMed consortium grant by the BMBF (InCa-01ZX1901A) and a collaborative research center grant on cell death (CRC1403, project ID 414786233) and predictability in evolution (CRC1310, project ID 325931972) funded by the Deutsche Forschungsgemeinschaft (DFG). JB was funded through a Mildred Scheel Junior Research Group Grant (ID 70113307) by the German cancer aid, JKML was funded by the Volkswagen foundation (Experiment!; A131649-97391) and by the German research foundation (Walter Benjamin Program; LI3844/1-1). Work in the Leprivier lab was funded by the German research foundation (LE 3751/2-1), the German Cancer Aid (70112624), and the Dr. Rolf M. Schwiete Stiftung (2020-018). Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
Conceptualization: SvK; visualization, FM; methodology: FM, JKML, CMB, ES, ST, AD, JS, JuBe, FIY, KN, JB, SvK; investigation: FM, JKML, CMB, ES, ST, AD, JS, JuBe, FIY, JB, SvK; formal analysis: FM, JKML, GL and SvK; writing-original draft: SvK; writing-review and editing: FM, JKML, ES, ST, FIY, SvK, LTF, GL, JB; funding acquisition; SvK; project administration: SvK; supervision: SvK, GL.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethics approval
Our study makes use of publicly available RNA-sequencing datasets as cited in the results section. As such, prior ethics approval has been obtained for these studies and no additional approval is required. All mouse experiments were conducted in accordance with an Institutional Animal Care and Use Committee (IACUC). All people involved in animal experiments received prior training and have passed the additionally required personal licensing course (FELASA B). All animal experiments were approved by the local authorities (LANUV, North-Rhine-Westphalia, Germany) and performed under license number 81-02.04.2017.A477.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Edited by G. Melino
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Müller, F., Lim, J.K.M., Bebber, C.M. et al. Elevated FSP1 protects KRAS-mutated cells from ferroptosis during tumor initiation. Cell Death Differ 30, 442â456 (2023). https://doi.org/10.1038/s41418-022-01096-8
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41418-022-01096-8
This article is cited by
-
Ferroptotic therapy in cancer: benefits, side effects, and risks
Molecular Cancer (2024)
-
Ferroptosis in cancer: from molecular mechanisms to therapeutic strategies
Signal Transduction and Targeted Therapy (2024)
-
Ferroptosis as a new tool for tumor suppression through lipid peroxidation
Communications Biology (2024)
-
New insights into crosstalk between Nrf2 pathway and ferroptosis in lung disease
Cell Death & Disease (2024)
-
FSP1 inhibition enhances olaparib sensitivity in BRCA-proficient ovarian cancer patients via a nonferroptosis mechanism
Cell Death & Differentiation (2024)